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

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

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

revision 1637 by ksteube, Mon Jul 14 05:34:59 2008 UTC revision 1661 by gross, Mon Jul 21 22:08:27 2008 UTC
# Line 63  double Paso_FCTransportProblem_getSafeTi Line 63  double Paso_FCTransportProblem_getSafeTi
63                 for (i=0;i<n_rows;++i) {                 for (i=0;i<n_rows;++i) {
64                      rtmp1=fctp->main_diagonal_low_order_transport_matrix[i];                      rtmp1=fctp->main_diagonal_low_order_transport_matrix[i];
65                      rtmp2=fctp->lumped_mass_matrix[i];                      rtmp2=fctp->lumped_mass_matrix[i];
66                      if ( (rtmp1<0 && rtmp2>0.) || (rtmp1>0 && rtmp2<0.) ) {                      if ( (rtmp1<0 && rtmp2>0.) || (rtmp1>0 && rtmp2<0) ) {
67                          dt_max_loc=MIN(dt_max_loc,-rtmp2/rtmp1);                          dt_max_loc=MIN(dt_max_loc,-rtmp2/rtmp1);
68                      }                      }
69                  }                  }
# Line 339  void Paso_SolverFCT_solve(Paso_FCTranspo Line 339  void Paso_SolverFCT_solve(Paso_FCTranspo
339                         }                         }
340                         Paso_Coupler_startCollect(u_m_coupler,u_m);                         Paso_Coupler_startCollect(u_m_coupler,u_m);
341                         #ifdef PASO_MPI                         #ifdef PASO_MPI
342                  double local_norm[2], norm[2];
343                            local_norm[0]=norm_u;                            local_norm[0]=norm_u;
344                            local_norm[1]=norm_du;                            local_norm[1]=norm_du;
345                        MPI_Allreduce(local_norm,norm, 2, MPI_DOUBLE, MPI_MAX, fctp->mpi_info->comm);                        MPI_Allreduce(local_norm,norm, 2, MPI_DOUBLE, MPI_MAX, fctp->mpi_info->comm);

Legend:
Removed from v.1637  
changed lines
  Added in v.1661

  ViewVC Help
Powered by ViewVC 1.1.26