Changeset 12106

Timestamp:

04/23/12 11:03:24 (13 years ago)

Author:

habbalf

Message:

Create dense matrix if sparsity=1

File:

• issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp

@@ -68 +68 @@
 
         #ifdef _HAVE_PETSCDEV_
-        MatCreateAIJ(MPI_COMM_WORLD,m,n,M,N,d_nz,NULL,o_nz,NULL,&outmatrix); 
+        if(sparsity==1){
+                MatCreateDense(MPI_COMM_WORLD,m,n,M,N,NULL,&outmatrix); 
+        }
+        else{
+                MatCreateAIJ(MPI_COMM_WORLD,m,n,M,N,d_nz,NULL,o_nz,NULL,&outmatrix); 
+        }
         #else
         MatCreateMPIAIJ(MPI_COMM_WORLD,m,n,M,N,d_nz,NULL,o_nz,NULL,&outmatrix);