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

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

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

revision 700 by gross, Thu Apr 6 00:13:40 2006 UTC revision 1028 by gross, Wed Mar 14 00:15:24 2007 UTC
# Line 49  void Paso_Solver_ILU_free(Paso_Solver_IL Line 49  void Paso_Solver_ILU_free(Paso_Solver_IL
49  Paso_Solver_ILU* Paso_Solver_getILU(Paso_SystemMatrix * A,bool_t verbose) {  Paso_Solver_ILU* Paso_Solver_getILU(Paso_SystemMatrix * A,bool_t verbose) {
50    dim_t n=A->num_rows;    dim_t n=A->num_rows;
51    dim_t n_block=A->row_block_size;    dim_t n_block=A->row_block_size;
52    index_t num_colors=0;    index_t num_colors=0, *mis_marker=NULL;
53    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;
54    register double mainA11,mainA12,mainA13,mainA21,mainA22,mainA23,mainA31,mainA32,mainA33;    register double mainA11,mainA12,mainA13,mainA21,mainA22,mainA23,mainA31,mainA32,mainA33;
55    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;
# Line 58  Paso_Solver_ILU* Paso_Solver_getILU(Paso Line 58  Paso_Solver_ILU* Paso_Solver_getILU(Paso
58    /* allocations: */      /* allocations: */  
59    Paso_Solver_ILU* out=MEMALLOC(1,Paso_Solver_ILU);    Paso_Solver_ILU* out=MEMALLOC(1,Paso_Solver_ILU);
60    if (Paso_checkPtr(out)) return NULL;    if (Paso_checkPtr(out)) return NULL;
61    index_t* mis_marker=TMPMEMALLOC(n,index_t);    mis_marker=TMPMEMALLOC(n,index_t);
62    out->colorOf=MEMALLOC(n,index_t);    out->colorOf=MEMALLOC(n,index_t);
63    out->factors=MEMALLOC(A->len,double);    out->factors=MEMALLOC(A->len,double);
64    out->main_iptr=MEMALLOC(n,index_t);    out->main_iptr=MEMALLOC(n,index_t);

Legend:
Removed from v.700  
changed lines
  Added in v.1028

  ViewVC Help
Powered by ViewVC 1.1.26