/[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 1363 by gross, Mon Dec 17 03:42:03 2007 UTC revision 1364 by gross, Mon Dec 17 07:22:45 2007 UTC
# Line 39  void Paso_FCTransportProblem_free(Paso_F Line 39  void Paso_FCTransportProblem_free(Paso_F
39             Paso_SystemMatrix_free(in->flux_matrix);             Paso_SystemMatrix_free(in->flux_matrix);
40             Paso_MPIInfo_free(in->mpi_info);             Paso_MPIInfo_free(in->mpi_info);
41    
42               MEMFREE(in->u);
43             MEMFREE(in->lumped_mass_matrix);             MEMFREE(in->lumped_mass_matrix);
44             MEMFREE(in->row_sum_flux_matrix);             MEMFREE(in->row_sum_flux_matrix);
45             MEMFREE(in->colorOf);             MEMFREE(in->colorOf);
# Line 106  Paso_FCTransportProblem* Paso_FCTranspor Line 107  Paso_FCTransportProblem* Paso_FCTranspor
107           out->main_iptr=MEMALLOC(n,index_t);           out->main_iptr=MEMALLOC(n,index_t);
108           out->lumped_mass_matrix=MEMALLOC(n,double);           out->lumped_mass_matrix=MEMALLOC(n,double);
109           out->row_sum_flux_matrix=MEMALLOC(n,double);           out->row_sum_flux_matrix=MEMALLOC(n,double);
110             out->u=MEMALLOC(n,double);
111    
112           if ( ! (Paso_checkPtr(out->colorOf) || Paso_checkPtr(out->main_iptr) ||           if ( ! (Paso_checkPtr(out->colorOf) || Paso_checkPtr(out->main_iptr) ||
113                   Paso_checkPtr(out->lumped_mass_matrix) || Paso_checkPtr(out->row_sum_flux_matrix)) ) {                   Paso_checkPtr(out->lumped_mass_matrix) || Paso_checkPtr(out->row_sum_flux_matrix) || Paso_checkPtr(out->u)) ) {
114                            
115               printf("Paso_SolverFCT_getFCTransportProblem: Revise coloring!!\n");               printf("Paso_SolverFCT_getFCTransportProblem: Revise coloring!!\n");
116               Paso_Pattern_color(pattern->mainPattern,&(out->num_colors),out->colorOf);               Paso_Pattern_color(pattern->mainPattern,&(out->num_colors),out->colorOf);
# Line 141  Paso_FCTransportProblem* Paso_FCTranspor Line 143  Paso_FCTransportProblem* Paso_FCTranspor
143       return NULL;       return NULL;
144    }    }
145  }  }
146    
147    void Paso_FCTransportProblem_checkinSolution(Paso_FCTransportProblem* in, double* u) {
148        dim_t i, n;
149      
150        n=Paso_FCTransportProblem_getTotalNumRows(in);
151        #pragma omp parallel for schedule(static) private(i)
152        for (i = 0; i < n; ++i) {
153             in->u[i]=u[i];
154        }
155    }

Legend:
Removed from v.1363  
changed lines
  Added in v.1364

  ViewVC Help
Powered by ViewVC 1.1.26