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

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

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

revision 4796 by jfenwick, Thu Feb 6 06:12:20 2014 UTC revision 4797 by caltinay, Mon Mar 24 05:42:42 2014 UTC
# Line 45  void Paso_Solver_ILU_free(Paso_Solver_IL Line 45  void Paso_Solver_ILU_free(Paso_Solver_IL
45    
46  /*   constructs the incomplete block factorization            */  /*   constructs the incomplete block factorization            */
47    
48  Paso_Solver_ILU* Paso_Solver_getILU(Paso_SparseMatrix * A,bool verbose) {  Paso_Solver_ILU* Paso_Solver_getILU(paso::SparseMatrix * A,bool verbose) {
49    const dim_t n=A->numRows;    const dim_t n=A->numRows;
50    const dim_t n_block=A->row_block_size;    const dim_t n_block=A->row_block_size;
51    const index_t* colorOf = Paso_Pattern_borrowColoringPointer(A->pattern);    const index_t* colorOf = Paso_Pattern_borrowColoringPointer(A->pattern);
52    const dim_t num_colors = Paso_Pattern_getNumColors(A->pattern);    const dim_t num_colors = Paso_Pattern_getNumColors(A->pattern);
53    const index_t *ptr_main = Paso_SparseMatrix_borrowMainDiagonalPointer(A);    const index_t *ptr_main = paso::SparseMatrix_borrowMainDiagonalPointer(A);
54    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;
55    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;
56    register index_t i,iptr_main,iptr_ik,k,iptr_kj,j,iptr_ij,color,color2, iptr;    register index_t i,iptr_main,iptr_ik,k,iptr_kj,j,iptr_ij,color,color2, iptr;
# Line 322  Paso_Solver_ILU* Paso_Solver_getILU(Paso Line 322  Paso_Solver_ILU* Paso_Solver_getILU(Paso
322     vector is available.     vector is available.
323  */  */
324    
325  void Paso_Solver_solveILU(Paso_SparseMatrix * A, Paso_Solver_ILU * ilu, double * x, const double * b) {  void Paso_Solver_solveILU(paso::SparseMatrix * A, Paso_Solver_ILU * ilu, double * x, const double * b) {
326       register dim_t i,k;       register dim_t i,k;
327       register index_t color,iptr_ik,iptr_main;       register index_t color,iptr_ik,iptr_main;
328       register double S1,S2,S3,R1,R2,R3;       register double S1,S2,S3,R1,R2,R3;
# Line 330  void Paso_Solver_solveILU(Paso_SparseMat Line 330  void Paso_Solver_solveILU(Paso_SparseMat
330       const dim_t n_block=A->row_block_size;       const dim_t n_block=A->row_block_size;
331       const index_t* colorOf = Paso_Pattern_borrowColoringPointer(A->pattern);       const index_t* colorOf = Paso_Pattern_borrowColoringPointer(A->pattern);
332       const dim_t num_colors = Paso_Pattern_getNumColors(A->pattern);       const dim_t num_colors = Paso_Pattern_getNumColors(A->pattern);
333       const index_t *ptr_main = Paso_SparseMatrix_borrowMainDiagonalPointer(A);       const index_t *ptr_main = paso::SparseMatrix_borrowMainDiagonalPointer(A);
334            
335            
336       /* copy x into b */       /* copy x into b */

Legend:
Removed from v.4796  
changed lines
  Added in v.4797

  ViewVC Help
Powered by ViewVC 1.1.26