/[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 1362 by gross, Mon Dec 17 02:28:16 2007 UTC revision 1363 by gross, Mon Dec 17 03:42:03 2007 UTC
# Line 40  void Paso_FCTransportProblem_free(Paso_F Line 40  void Paso_FCTransportProblem_free(Paso_F
40             Paso_MPIInfo_free(in->mpi_info);             Paso_MPIInfo_free(in->mpi_info);
41    
42             MEMFREE(in->lumped_mass_matrix);             MEMFREE(in->lumped_mass_matrix);
43               MEMFREE(in->row_sum_flux_matrix);
44             MEMFREE(in->colorOf);             MEMFREE(in->colorOf);
45             MEMFREE(in->main_iptr);             MEMFREE(in->main_iptr);
46             MEMFREE(in);             MEMFREE(in);
# Line 87  Paso_FCTransportProblem* Paso_FCTranspor Line 88  Paso_FCTransportProblem* Paso_FCTranspor
88       if (Paso_checkPtr(out)) return NULL;       if (Paso_checkPtr(out)) return NULL;
89    
90       out->theta=theta;       out->theta=theta;
91         out->valid_matrices=FALSE;
92       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);
93       Paso_SystemMatrix_allocBuffer(out->transport_matrix);       Paso_SystemMatrix_allocBuffer(out->transport_matrix);
94       out->flux_matrix=Paso_SystemMatrix_alloc(matrix_type,pattern,block_size,block_size);       out->flux_matrix=Paso_SystemMatrix_alloc(matrix_type,pattern,block_size,block_size);
# Line 95  Paso_FCTransportProblem* Paso_FCTranspor Line 97  Paso_FCTransportProblem* Paso_FCTranspor
97       out->colorOf=NULL;       out->colorOf=NULL;
98       out->main_iptr=NULL;       out->main_iptr=NULL;
99       out->lumped_mass_matrix=NULL;       out->lumped_mass_matrix=NULL;
100         out->row_sum_flux_matrix=NULL;
101    
102       if (Paso_noError()) {       if (Paso_noError()) {
103           n=Paso_SystemMatrix_getTotalNumRows(out->transport_matrix);           n=Paso_SystemMatrix_getTotalNumRows(out->transport_matrix);
# Line 102  Paso_FCTransportProblem* Paso_FCTranspor Line 105  Paso_FCTransportProblem* Paso_FCTranspor
105           out->colorOf=MEMALLOC(n,index_t);           out->colorOf=MEMALLOC(n,index_t);
106           out->main_iptr=MEMALLOC(n,index_t);           out->main_iptr=MEMALLOC(n,index_t);
107           out->lumped_mass_matrix=MEMALLOC(n,double);           out->lumped_mass_matrix=MEMALLOC(n,double);
108             out->row_sum_flux_matrix=MEMALLOC(n,double);
109    
110           if ( ! (Paso_checkPtr(out->colorOf) || Paso_checkPtr(out->main_iptr) || Paso_checkPtr(out->lumped_mass_matrix)) ) {           if ( ! (Paso_checkPtr(out->colorOf) || Paso_checkPtr(out->main_iptr) ||
111                     Paso_checkPtr(out->lumped_mass_matrix) || Paso_checkPtr(out->row_sum_flux_matrix)) ) {
112                            
113               printf("Paso_SolverFCT_getFCTransportProblem: Revise coloring!!\n");               printf("Paso_SolverFCT_getFCTransportProblem: Revise coloring!!\n");
114               Paso_Pattern_color(pattern->mainPattern,&(out->num_colors),out->colorOf);               Paso_Pattern_color(pattern->mainPattern,&(out->num_colors),out->colorOf);

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

  ViewVC Help
Powered by ViewVC 1.1.26