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

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

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

revision 2659 by artak, Thu Sep 10 03:54:50 2009 UTC revision 2661 by artak, Fri Sep 11 00:59:59 2009 UTC
# Line 43  void Paso_Solver_AMG_System_free(Paso_So Line 43  void Paso_Solver_AMG_System_free(Paso_So
43          Paso_SparseMatrix_free(in->block1);          Paso_SparseMatrix_free(in->block1);
44          Paso_SparseMatrix_free(in->block2);          Paso_SparseMatrix_free(in->block2);
45          Paso_SparseMatrix_free(in->block3);          Paso_SparseMatrix_free(in->block3);
   
46          MEMFREE(in);          MEMFREE(in);
47       }       }
48  }  }
# Line 147  Paso_Solver_AMG* Paso_Solver_getAMG(Paso Line 146  Paso_Solver_AMG* Paso_Solver_getAMG(Paso
146       out->solver=NULL;       out->solver=NULL;
147       /*out->GS=Paso_Solver_getGS(A_p,verbose);*/       /*out->GS=Paso_Solver_getGS(A_p,verbose);*/
148       out->level=level;       out->level=level;
149         out->n=n;
150         out->n_block=n_block;
151            
152       if (level==0 || n<=options->min_coarse_matrix_size) {       if (level==0 || n<=options->min_coarse_matrix_size) {
153           out->coarsest_level=TRUE;           out->coarsest_level=TRUE;
# Line 185  Paso_Solver_AMG* Paso_Solver_getAMG(Paso Line 186  Paso_Solver_AMG* Paso_Solver_getAMG(Paso
186             #pragma omp parallel for private(i) schedule(static)             #pragma omp parallel for private(i) schedule(static)
187             for (i = 0; i < n; ++i) counter[i]=mis_marker[i];             for (i = 0; i < n; ++i) counter[i]=mis_marker[i];
188    
            out->n=n;  
            out->n_block=n_block;  
189             out->n_F=Paso_Util_cumsum(n,counter);             out->n_F=Paso_Util_cumsum(n,counter);
190             out->mask_F=MEMALLOC(n,index_t);             out->mask_F=MEMALLOC(n,index_t);
191             out->rows_in_F=MEMALLOC(out->n_F,index_t);             out->rows_in_F=MEMALLOC(out->n_F,index_t);

Legend:
Removed from v.2659  
changed lines
  Added in v.2661

  ViewVC Help
Powered by ViewVC 1.1.26