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

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

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

revision 3120 by gross, Mon Aug 30 10:48:00 2010 UTC revision 3158 by gross, Mon Sep 6 06:09:11 2010 UTC
# Line 27  Line 27 
27  #include "SparseMatrix.h"  #include "SparseMatrix.h"
28  #include "Solver.h"  #include "Solver.h"
29  #include "BlockOps.h"  #include "BlockOps.h"
30    #include "PasoUtil.h"
31    
32  void Paso_SparseMatrix_invMain(Paso_SparseMatrix * A_p, double* inv_diag, int* pivot) {  void Paso_SparseMatrix_invMain(Paso_SparseMatrix * A_p, double* inv_diag, int* pivot) {
33     /*   inv_diag=MEMALLOC( A->numRows * A_p-> block_size,double);     /*   inv_diag=MEMALLOC( A->numRows * A_p-> block_size,double);
# Line 129  void Paso_SparseMatrix_applyBlockMatrix( Line 130  void Paso_SparseMatrix_applyBlockMatrix(
130     pivot=MEMALLOC( A->numRows * A->row_block_size */     pivot=MEMALLOC( A->numRows * A->row_block_size */
131     dim_t n=A_p->numRows;     dim_t n=A_p->numRows;
132     dim_t n_block=A_p->row_block_size;     dim_t n_block=A_p->row_block_size;
133         Paso_Copy(n_block*n, x,b);
134     Paso_BlockOps_allMV(n_block,n,block_diag,pivot,x,b);     Paso_BlockOps_allMV(n_block,n,block_diag,pivot,x);
135  }  }
136    

Legend:
Removed from v.3120  
changed lines
  Added in v.3158

  ViewVC Help
Powered by ViewVC 1.1.26