/[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 1410 by gross, Thu Feb 7 04:24:00 2008 UTC revision 1552 by gross, Thu May 8 08:52:41 2008 UTC
# Line 34  void Paso_FCTransportProblem_free(Paso_F Line 34  void Paso_FCTransportProblem_free(Paso_F
34       if (in!=NULL) {       if (in!=NULL) {
35          in->reference_counter--;          in->reference_counter--;
36          if (in->reference_counter<=0) {          if (in->reference_counter<=0) {
   
37             Paso_SystemMatrix_free(in->transport_matrix);             Paso_SystemMatrix_free(in->transport_matrix);
38             Paso_SystemMatrix_free(in->mass_matrix);             Paso_SystemMatrix_free(in->mass_matrix);
39             Paso_SystemMatrix_free(in->iteration_matrix);             Paso_SystemMatrix_free(in->iteration_matrix);
# Line 92  Paso_FCTransportProblem* Paso_FCTranspor Line 91  Paso_FCTransportProblem* Paso_FCTranspor
91       out->dt_max=LARGE_POSITIVE_FLOAT;       out->dt_max=LARGE_POSITIVE_FLOAT;
92       out->valid_matrices=FALSE;       out->valid_matrices=FALSE;
93       out->transport_matrix=Paso_SystemMatrix_alloc(matrix_type,pattern,block_size,block_size);       out->transport_matrix=Paso_SystemMatrix_alloc(matrix_type,pattern,block_size,block_size);
      Paso_SystemMatrix_allocBuffer(out->transport_matrix);  
94       out->mass_matrix=Paso_SystemMatrix_alloc(matrix_type,pattern,block_size,block_size);       out->mass_matrix=Paso_SystemMatrix_alloc(matrix_type,pattern,block_size,block_size);
95       out->iteration_matrix=NULL;       out->iteration_matrix=NULL;
96    
# Line 164  void Paso_FCTransportProblem_checkinSolu Line 162  void Paso_FCTransportProblem_checkinSolu
162      }      }
163      #ifdef PASO_MPI      #ifdef PASO_MPI
164           local_u_min=u_min;           local_u_min=u_min;
165           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, in->mpi_info->comm);
166      #endif      #endif
167      in->u_min=u_min;      in->u_min=u_min;
168      #pragma omp parallel for schedule(static) private(i)      #pragma omp parallel for schedule(static) private(i)

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

  ViewVC Help
Powered by ViewVC 1.1.26