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

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

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

revision 1841 by gross, Fri Oct 3 03:57:52 2008 UTC revision 1975 by artak, Thu Nov 6 03:07:02 2008 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_SparseMatrix * A,bool_t verbose) {  Paso_Solver_ILU* Paso_Solver_getILU(Paso_SparseMatrix * A,bool_t verbose) {
50    dim_t n=A->numRows;    dim_t n=A->numRows;
51    dim_t n_block=A->row_block_size;    dim_t n_block=A->row_block_size;
   index_t num_colors=0, *mis_marker=NULL;  
52    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;
53    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;
54    register index_t i,iptr_main,iptr,iptr_ik,k,iptr_kj,j,iptr_ij,color,color2;    register index_t i,iptr_main,iptr,iptr_ik,k,iptr_kj,j,iptr_ij,color,color2;
55    double time0,time_color,time_fac;    double time0=0,time_color=0,time_fac=0;
56    /* allocations: */      /* allocations: */  
57    Paso_Solver_ILU* out=MEMALLOC(1,Paso_Solver_ILU);    Paso_Solver_ILU* out=MEMALLOC(1,Paso_Solver_ILU);
58    if (Paso_checkPtr(out)) return NULL;    if (Paso_checkPtr(out)) return NULL;

Legend:
Removed from v.1841  
changed lines
  Added in v.1975

  ViewVC Help
Powered by ViewVC 1.1.26