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

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

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

revision 1973 by ksteube, Thu Sep 25 23:11:13 2008 UTC revision 1974 by jfenwick, Thu Nov 6 02:40:10 2008 UTC
# Line 247  err_t Paso_Solver_PCG( Line 247  err_t Paso_Solver_PCG(
247             /* v=A*p */             /* v=A*p */
248             Performance_stopMonitor(pp,PERFORMANCE_SOLVER);             Performance_stopMonitor(pp,PERFORMANCE_SOLVER);
249             Performance_startMonitor(pp,PERFORMANCE_MVM);             Performance_startMonitor(pp,PERFORMANCE_MVM);
250         Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(ONE, A, p,ZERO,v);         Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(PASO_ONE, A, p,PASO_ZERO,v);
251             Performance_stopMonitor(pp,PERFORMANCE_MVM);             Performance_stopMonitor(pp,PERFORMANCE_MVM);
252             Performance_startMonitor(pp,PERFORMANCE_SOLVER);             Performance_startMonitor(pp,PERFORMANCE_SOLVER);
253    
# Line 333  err_t Paso_Solver_PCG( Line 333  err_t Paso_Solver_PCG(
333              sum_5 = 0;              sum_5 = 0;
334                  #pragma omp parallel private(i0, istart, iend, ipp, ss, gamma_1,gamma_2)                  #pragma omp parallel private(i0, istart, iend, ipp, ss, gamma_1,gamma_2)
335                  {                  {
336                    gamma_1= ( (ABS(sum_3)<= ZERO) ? 0 : -sum_4/sum_3) ;                    gamma_1= ( (ABS(sum_3)<= PASO_ZERO) ? 0 : -sum_4/sum_3) ;
337                    gamma_2= ONE-gamma_1;                    gamma_2= PASO_ONE-gamma_1;
338                    ss=0;                    ss=0;
339                    #ifdef USE_DYNAMIC_SCHEDULING                    #ifdef USE_DYNAMIC_SCHEDULING
340                        #pragma omp for schedule(dynamic, 1)                        #pragma omp for schedule(dynamic, 1)

Legend:
Removed from v.1973  
changed lines
  Added in v.1974

  ViewVC Help
Powered by ViewVC 1.1.26