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

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

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

revision 1407 by gross, Mon Feb 4 06:45:48 2008 UTC revision 1410 by gross, Thu Feb 7 04:24:00 2008 UTC
# Line 156  void Paso_FCTransportProblem_checkinSolu Line 156  void Paso_FCTransportProblem_checkinSolu
156      {      {
157           local_u_min=0.;           local_u_min=0.;
158           #pragma omp for schedule(static) private(i)           #pragma omp for schedule(static) private(i)
159           for (i=0;i<n;++i) {           for (i=0;i<n;++i) local_u_min=MIN(local_u_min,u[i]);
             in->u[i]=u[i];  
             local_u_min=MIN(local_u_min,u[i]);  
          }  
160           #pragma omp critical           #pragma omp critical
161           {           {
162              u_min=MIN(u_min,local_u_min);              u_min=MIN(u_min,local_u_min);
# Line 170  void Paso_FCTransportProblem_checkinSolu Line 167  void Paso_FCTransportProblem_checkinSolu
167           MPI_Allreduce(&local_u_min,&u_min, 1, MPI_DOUBLE, MPI_MIN, fctp->mpi_info->comm);           MPI_Allreduce(&local_u_min,&u_min, 1, MPI_DOUBLE, MPI_MIN, fctp->mpi_info->comm);
168      #endif      #endif
169      in->u_min=u_min;      in->u_min=u_min;
170        #pragma omp parallel for schedule(static) private(i)
171        for (i=0;i<n;++i) {
172            in->u[i]=u[i]-u_min;
173        }
174  }  }

Legend:
Removed from v.1407  
changed lines
  Added in v.1410

  ViewVC Help
Powered by ViewVC 1.1.26