Changeset 13869

Timestamp:

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

Author:

Mathieu Morlighem

Message:

NEW: added AllocationInfo method for matrices

Location:

issm/trunk-jpl/src/c

Files:

• classes/matrix/Matrix.h (modified) (1 diff)
• toolkits/petsc/objects/PetscMat.cpp (modified) (2 diffs)
• toolkits/petsc/objects/PetscMat.h (modified) (1 diff)

issm/trunk-jpl/src/c/classes/matrix/Matrix.h

@@ -189 +189 @@
                 }
                 /*}}}*/
+                /*FUNCTION AllocationInfo{{{*/
+                void AllocationInfo(void){
+                        _assert_(this);
+                        if(type==PetscMatType){
+                                #ifdef _HAVE_PETSC_
+                                this->pmatrix->AllocationInfo();
+                                #endif
+                        }
+                        else if(type==SeqMatType){
+                                //this->smatrix->AllocationInfo();
+                        }
+                        else _error_("Matrix type: " << type << " not supported yet!");
+                }/*}}}*/
                 /*FUNCTION Assemble{{{*/
                 void Assemble(void){

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){

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

@@ -40 +40 @@
                 /*}}}*/
                 /*PetscMat specific routines {{{*/
+                void AllocationInfo(void);
                 void Echo(void);
                 void Assemble(void);