Changeset 13036 for issm/trunk-jpl/src/c/classes/matrix/Matrix.cpp
- Timestamp:
- 08/14/12 16:58:21 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/matrix/Matrix.cpp
r12987 r13036 51 51 this->pmatrix=new PetscMat(M,N); 52 52 #else 53 _error 2_("Petsc matrix format not usable, as Petsc has not been compiled!");53 _error_("Petsc matrix format not usable, as Petsc has not been compiled!"); 54 54 #endif 55 55 } … … 57 57 this->smatrix=new SeqMat(M,N); 58 58 } 59 else _error 2_("Matrix type: " << type << " not supported yet!");59 else _error_("Matrix type: " << type << " not supported yet!"); 60 60 61 61 } … … 74 74 this->pmatrix=new PetscMat(M,N,sparsity); 75 75 #else 76 _error 2_("Petsc matrix format not usable, as Petsc has not been compiled!");76 _error_("Petsc matrix format not usable, as Petsc has not been compiled!"); 77 77 #endif 78 78 } … … 80 80 this->smatrix=new SeqMat(M,N,sparsity); 81 81 } 82 else _error 2_("Matrix type: " << type << " not supported yet!");82 else _error_("Matrix type: " << type << " not supported yet!"); 83 83 } 84 84 /*}}}*/ … … 96 96 this->pmatrix=new PetscMat(serial_mat,M,N,sparsity); 97 97 #else 98 _error 2_("Petsc matrix format not usable, as Petsc has not been compiled!");98 _error_("Petsc matrix format not usable, as Petsc has not been compiled!"); 99 99 #endif 100 100 } … … 102 102 this->smatrix=new SeqMat(serial_mat,M,N,sparsity); 103 103 } 104 else _error 2_("Matrix type: " << type << " not supported yet!");104 else _error_("Matrix type: " << type << " not supported yet!"); 105 105 106 106 } … … 119 119 this->pmatrix=new PetscMat(M,N,connectivity,numberofdofspernode); 120 120 #else 121 _error 2_("Petsc matrix format not usable, as Petsc has not been compiled!");121 _error_("Petsc matrix format not usable, as Petsc has not been compiled!"); 122 122 #endif 123 123 } … … 125 125 this->smatrix=new SeqMat(M,N,connectivity,numberofdofspernode); 126 126 } 127 else _error 2_("Matrix type: " << type << " not supported yet!");127 else _error_("Matrix type: " << type << " not supported yet!"); 128 128 129 129 } … … 136 136 delete this->pmatrix; 137 137 #else 138 _error 2_("Petsc matrix format not usable, as Petsc has not been compiled!");138 _error_("Petsc matrix format not usable, as Petsc has not been compiled!"); 139 139 #endif 140 140 } … … 142 142 delete this->smatrix; 143 143 } 144 else _error 2_("Matrix type: " << type << " not supported yet!");144 else _error_("Matrix type: " << type << " not supported yet!"); 145 145 146 146 } … … 159 159 this->smatrix->Echo(); 160 160 } 161 else _error 2_("Matrix type: " << type << " not supported yet!");161 else _error_("Matrix type: " << type << " not supported yet!"); 162 162 163 163 } … … 175 175 } 176 176 else{ 177 _error 2_("Matrix type: " << type << " not supported yet!");177 _error_("Matrix type: " << type << " not supported yet!"); 178 178 } 179 179 } … … 192 192 norm=this->smatrix->Norm(norm_type); 193 193 } 194 else _error 2_("Matrix type: " << type << " not supported yet!");194 else _error_("Matrix type: " << type << " not supported yet!"); 195 195 196 196 return norm; … … 208 208 this->smatrix->GetSize(pM,pN); 209 209 } 210 else _error 2_("Matrix type: " << type << " not supported yet!");210 else _error_("Matrix type: " << type << " not supported yet!"); 211 211 212 212 } … … 223 223 this->smatrix->GetLocalSize(pM,pN); 224 224 } 225 else _error 2_("Matrix type: " << type << " not supported yet!");225 else _error_("Matrix type: " << type << " not supported yet!"); 226 226 227 227 } … … 238 238 this->smatrix->MatMult(X->svector,AX->svector); 239 239 } 240 else _error 2_("Matrix type: " << type << " not supported yet!");240 else _error_("Matrix type: " << type << " not supported yet!"); 241 241 242 242 } … … 257 257 output->smatrix=this->smatrix->Duplicate(); 258 258 } 259 else _error 2_("Matrix type: " << type << " not supported yet!");259 else _error_("Matrix type: " << type << " not supported yet!"); 260 260 261 261 return output; … … 275 275 output=this->smatrix->ToSerial(); 276 276 } 277 else _error 2_("Matrix type: " << type << " not supported yet!");277 else _error_("Matrix type: " << type << " not supported yet!"); 278 278 279 279 … … 292 292 this->smatrix->SetValues(m,idxm,n,idxn,values,mode); 293 293 } 294 else _error 2_("Matrix type: " << type << " not supported yet!");294 else _error_("Matrix type: " << type << " not supported yet!"); 295 295 } 296 296 /*}}}*/ … … 307 307 } 308 308 else{ 309 _error 2_("Matrix type: " << type << " not supported yet!");310 } 311 312 } 313 /*}}}*/ 309 _error_("Matrix type: " << type << " not supported yet!"); 310 } 311 312 } 313 /*}}}*/
Note:
See TracChangeset
for help on using the changeset viewer.