PGI User Forum
 SearchSearch   MemberlistMemberlist     RegisterRegister   ProfileProfile    Log inLog in 

CUDA-x86.

How can I change the non-rectangular loops to rectangular

 
Post new topic   Reply to topic    PGI User Forum Forum Index -> Accelerator Programming
View previous topic :: View next topic  
Author Message
Teslalady



Joined: 16 Mar 2012
Posts: 75

PostPosted: Mon Jun 04, 2012 8:17 am    Post subject: How can I change the non-rectangular loops to rectangular Reply with quote

As I know that all loops must be rectangular, but in my code, the loop is non-rectangular,how can I change it to rectangular loop, can you advise me?

My code is as below:

Code:
for(i=0; i<n-1; i++)
        {
       
                for(j=i+1; j<n; j++)
                {
                        for(k=0,sum=0; k<n; k++)
                        {
                                if(k != i)
                                {
                                   sum += l[j][k]*u[k][i];
                                }
                        }
                        l[j][i] = (float)((a[j][i]-sum)/u[i][i]);
                }

                for(j=i+1; j<n; j++)
                {
                        for(k=0,sum=0; k<n; k++)
                        {
                                if(k != i+1)
                                {
                                  sum += l[i+1][k]*u[k][j];
                                }
                        }
                        u[i+1][j] = (float)((a[i+1][j]-sum));
                }
        }
Back to top
View user's profile
mkcolg



Joined: 30 Jun 2004
Posts: 6129
Location: The Portland Group Inc.

PostPosted: Mon Jun 04, 2012 8:53 am    Post subject: Reply with quote

Make it rectangular then add an if statement to ignore part some of the iterations.

Code:
for(i=0; i<n-1; i++)
        {
                for(j=0; j<n; j++)
                {
                    if (j > i) {
                        for(k=0,sum=0; k<n; k++)
                        {
                                if(k != i)
                                {
                                   sum += l[j][k]*u[k][i];
                                }
                        }
                        l[j][i] = (float)((a[j][i]-sum)/u[i][i]);
                     }
                }

                for(j=0; j<n; j++)
                {
                    if (j > i) {
                        for(k=0,sum=0; k<n; k++)
                        {
                                if(k != i+1)
                                {
                                  sum += l[i+1][k]*u[k][j];
                                }
                        }
                        u[i+1][j] = (float)((a[i+1][j]-sum));
                    }
                }
        }


Hope this helps,
Mat
Back to top
View user's profile
Display posts from previous:   
Post new topic   Reply to topic    PGI User Forum Forum Index -> Accelerator Programming All times are GMT - 7 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © phpBB Group