Changeset 12903
- Timestamp:
- 08/04/12 09:17:13 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/Solverx/Solverx.cpp
r12900 r12903 30 30 31 31 /*According to matrix type, use specific solvers: */ 32 if(Kff->type==PetscMatType){32 switch(Kff->type){ 33 33 #ifdef _HAVE_PETSC_ 34 PetscVec* uf0_vector = NULL; 35 PetscVec* df_vector = NULL; 36 if(uf0) uf0_vector = uf0->pvector; 37 if(df) df_vector = df->pvector; 38 39 SolverxPetsc(&uf->pvector,Kff->pmatrix,pf->pvector,uf0_vector,df_vector,parameters); 40 #else 41 _error_("ISSM should be compiled with Petsc support for this operation!"); 34 case PetscMatType: 35 PetscVec* uf0_vector = NULL; 36 PetscVec* df_vector = NULL; 37 if(uf0) uf0_vector = uf0->pvector; 38 if(df) df_vector = df->pvector; 39 SolverxPetsc(&uf->pvector,Kff->pmatrix,pf->pvector,uf0_vector,df_vector,parameters); 40 break; 42 41 #endif 42 case SeqMatType: 43 SolverxSeq(&uf->svector,Kff->smatrix,pf->svector); 44 break 45 default: 46 _error2_("Matrix type: " << Kff->type << " not supported yet!"); 43 47 } 44 else if(Kff->type==SeqMatType){45 SolverxSeq(&uf->svector,Kff->smatrix,pf->svector);46 }47 else _error2_("Matrix type: " << Kff->type << " not supported yet!");48 48 49 49 /*Assign output pointers: */
Note:
See TracChangeset
for help on using the changeset viewer.