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

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

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

revision 414 by jgs, Wed Nov 9 02:02:19 2005 UTC revision 415 by gross, Wed Jan 4 05:37:33 2006 UTC
# Line 326  void Paso_Solver_solveILU(Paso_Solver_IL Line 326  void Paso_Solver_solveILU(Paso_Solver_IL
326          /* x_F=invA_FF*b_F  */          /* x_F=invA_FF*b_F  */
327          Paso_Solver_applyBlockDiagonalMatrix(n_block,ilu->n_F,ilu->inv_A_FF,ilu->A_FF_pivot,ilu->x_F,ilu->b_F);          Paso_Solver_applyBlockDiagonalMatrix(n_block,ilu->n_F,ilu->inv_A_FF,ilu->A_FF_pivot,ilu->x_F,ilu->b_F);
328          /* b_C=b_C-A_CF*x_F */          /* b_C=b_C-A_CF*x_F */
329          Paso_SystemMatrix_MatrixVector(-1.,ilu->A_CF,ilu->x_F,1.,ilu->b_C);          Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(-1.,ilu->A_CF,ilu->x_F,1.,ilu->b_C);
330          /* x_C=ILU(b_C)     */          /* x_C=ILU(b_C)     */
331          Paso_Solver_solveILU(ilu->ILU_of_Schur,ilu->x_C,ilu->b_C);          Paso_Solver_solveILU(ilu->ILU_of_Schur,ilu->x_C,ilu->b_C);
332          /* b_F=b_F-A_FC*x_C */          /* b_F=b_F-A_FC*x_C */
333          Paso_SystemMatrix_MatrixVector(-1.,ilu->A_FC,ilu->x_C,1.,ilu->b_F);          Paso_SystemMatrix_MatrixVector_CSR_OFFSET0(-1.,ilu->A_FC,ilu->x_C,1.,ilu->b_F);
334          /* x_F=invA_FF*b_F  */          /* x_F=invA_FF*b_F  */
335          Paso_Solver_applyBlockDiagonalMatrix(n_block,ilu->n_F,ilu->inv_A_FF,ilu->A_FF_pivot,ilu->x_F,ilu->b_F);          Paso_Solver_applyBlockDiagonalMatrix(n_block,ilu->n_F,ilu->inv_A_FF,ilu->A_FF_pivot,ilu->x_F,ilu->b_F);
336          /* x<-[x_F,x_C]     */          /* x<-[x_F,x_C]     */

Legend:
Removed from v.414  
changed lines
  Added in v.415

  ViewVC Help
Powered by ViewVC 1.1.26