Changeset 13869 for issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp

Timestamp:

11/02/12 14:42:38 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added AllocationInfo method for matrices

File:

• issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp ¶

@@ -40 +40 @@
 }
 /*}}}*/
-/*FUNCTION PetscMat(IssmDouble* serial_mat,int M,int N,IssmDouble sparsity){{{*/
+/*FUNCTION PetscMat::PetscMat(IssmDouble* serial_mat,int M,int N,IssmDouble sparsity){{{*/
 PetscMat::PetscMat(IssmDouble* serial_mat,int M,int N,IssmDouble sparsity){
 
@@ -77 +77 @@
 
 /*PetscMat specific routines: */
+/*FUNCTION PetscMat::AllocationInfo{{{*/
+void PetscMat::AllocationInfo(void){
+
+        MatInfo info;
+        MatGetInfo(this->matrix,MAT_GLOBAL_SUM,&info);
+        _pprintLine_("=========================== Stiffness matrix allocation info ===========================");
+        _pprintLine_("");
+        _pprintLine_(" nz_allocated: "<<info.nz_allocated);
+        _pprintLine_(" nz_used     : "<<info.nz_used);
+        _pprintLine_(" nz_unneeded : "<<info.nz_unneeded<<" ("<<double(info.nz_unneeded)/double(info.nz_allocated)*100.<<"%)");
+        _pprintLine_("");
+        _pprintLine_("========================================================================================");
+}
+/*}}}*/
 /*FUNCTION PetscMat::Echo{{{*/
 void PetscMat::Echo(void){