- Timestamp:
- 10/11/12 11:23:47 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp
r13596 r13622 46 46 int* idxm=NULL; 47 47 int* idxn=NULL; 48 48 49 49 if(M)idxm=xNew<int>(M); 50 50 if(N)idxn=xNew<int>(N); … … 52 52 for(i=0;i<M;i++)idxm[i]=i; 53 53 for(i=0;i<N;i++)idxn[i]=i; 54 55 54 56 55 this->matrix=NewMat(M,N,sparsity,IssmComm::GetComm()); … … 66 65 /*FUNCTION PetscMat::PetscMat(int M,int N, int connectivity, int numberofdofspernode){{{*/ 67 66 PetscMat::PetscMat(int M,int N, int connectivity,int numberofdofspernode){ 68 67 69 68 this->matrix=NewMat(M,N,connectivity,numberofdofspernode,IssmComm::GetComm()); 70 69 … … 99 98 IssmDouble PetscMat::Norm(NormMode mode){ 100 99 101 102 100 IssmDouble norm=0; 103 101 _assert_(this->matrix); 104 102 MatNorm(this->matrix,ISSMToPetscNormMode(mode),&norm); 105 103 106 104 return norm; 107 105
Note:
See TracChangeset
for help on using the changeset viewer.