/[escript]/trunk/paso/src/GMRES2.c
ViewVC logotype

Diff of /trunk/paso/src/GMRES2.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3013 by gross, Thu Apr 22 05:59:31 2010 UTC revision 3014 by gross, Wed Apr 28 04:05:21 2010 UTC
# Line 111  err_t Paso_Solver_GMRES2( Line 111  err_t Paso_Solver_GMRES2(
111                          normv2=Paso_l2(n,v[k],F->mpi_info);                          normv2=Paso_l2(n,v[k],F->mpi_info);
112                          h[INDEX2(k,k-1,l)]=normv2;                          h[INDEX2(k,k-1,l)]=normv2;
113                      }                      }
114    /*
115  {int p,q;  {int p,q;
116     for (p=0;p<k+1;p++){     for (p=0;p<k+1;p++){
117     for (q=0;q<k+1;q++)printf("%e ",h[INDEX2(p,q,l)]);     for (q=0;q<k+1;q++)printf("%e ",h[INDEX2(p,q,l)]);
# Line 118  err_t Paso_Solver_GMRES2( Line 119  err_t Paso_Solver_GMRES2(
119        
120     }     }
121  }  }
122    */
123                     /*                     /*
124                      *   watch out for happy breakdown                      *   watch out for happy breakdown
125                      */                      */
# Line 151  err_t Paso_Solver_GMRES2( Line 153  err_t Paso_Solver_GMRES2(
153        /*        /*
154         * all done and ready for the forward substitution:         * all done and ready for the forward substitution:
155        */        */
 printf("k =  %d\n", k);  
156    
157        for (i=k-1;i>=0;--i) {        for (i=k-1;i>=0;--i) {
158             for (j=i+1;j<k;j++) {             for (j=i+1;j<k;j++) {
 printf("%d %d : %e\n", i,j,h[INDEX2(i,j,l)]);  
159                g[i]-=h[INDEX2(i,j,l)]*g[j];                g[i]-=h[INDEX2(i,j,l)]*g[j];
160             }             }
161             g[i]/=h[INDEX2(i,i,l)];             g[i]/=h[INDEX2(i,i,l)];

Legend:
Removed from v.3013  
changed lines
  Added in v.3014

  ViewVC Help
Powered by ViewVC 1.1.26