/[escript]/branches/windows_from_1456_trunk_1580_merged_in/paso/src/PCG.c
ViewVC logotype

Diff of /branches/windows_from_1456_trunk_1580_merged_in/paso/src/PCG.c

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

revision 1592 by phornby, Wed Jun 4 07:20:59 2008 UTC revision 1593 by phornby, Thu Jun 5 11:01:24 2008 UTC
# Line 154  err_t Paso_Solver_PCG( Line 154  err_t Paso_Solver_PCG(
154         #endif         #endif
155      }      }
156      num_iter=0;      num_iter=0;
157    
158        /* PGH */
159        /* without this we get a use of an unititialised var below */
160        tau = 0;
161    
162      /* start of iteration */      /* start of iteration */
163      while (!(convergeFlag || maxIterFlag || breakFlag)) {      while (!(convergeFlag || maxIterFlag || breakFlag)) {
164             ++(num_iter);             ++(num_iter);
165    
166               /* PGH */
167               /* The next lines were commented out before I got here */
168             /* v=prec(r)  */             /* v=prec(r)  */
169               /* tau=v*r; */
170               /* leading to the use of an unititialised var below */
171    
172             Performance_stopMonitor(pp,PERFORMANCE_SOLVER);             Performance_stopMonitor(pp,PERFORMANCE_SOLVER);
173             Performance_startMonitor(pp,PERFORMANCE_PRECONDITIONER);             Performance_startMonitor(pp,PERFORMANCE_PRECONDITIONER);
174             Paso_Solver_solvePreconditioner(A,v,r);             Paso_Solver_solvePreconditioner(A,v,r);
175             Performance_stopMonitor(pp,PERFORMANCE_PRECONDITIONER);             Performance_stopMonitor(pp,PERFORMANCE_PRECONDITIONER);
176             Performance_startMonitor(pp,PERFORMANCE_SOLVER);             Performance_startMonitor(pp,PERFORMANCE_SOLVER);
177             /* tau=v*r    */  
178         sum_1 = 0;         sum_1 = 0;
179             #pragma omp parallel private(i0, istart, iend, ipp, ss)             #pragma omp parallel private(i0, istart, iend, ipp, ss)
180             {             {
Legend:
Removed from v.1592  
changed lines
  Added in v.1593
  ViewVC Help
Powered by ViewVC 1.1.26