Changeset 11739

Timestamp:

03/19/12 20:58:58 (13 years ago)

Author:

Eric.Larour

Message:

Some legitimate demands if the vector is a NULL pointer

File:

• issm/trunk-jpl/src/c/objects/Numerics/Vector.cpp (modified) (5 diffs)

issm/trunk-jpl/src/c/objects/Numerics/Vector.cpp

@@ -112 +112 @@
 
         #ifdef _HAVE_PETSC_
-        VecView(this->vector,PETSC_VIEWER_STDOUT_WORLD);
+        if(!this->vector){
+                printf("Vector size: 0\n");
+        }
+        else VecView(this->vector,PETSC_VIEWER_STDOUT_WORLD);
         #else
         this->vector->Echo();
@@ -209 +212 @@
                 
         #ifdef _HAVE_PETSC_
-                _assert_(this->vector);
-                VecGetSize(this->vector,pM);
+                if(!this->vector)*pM=0;
+                else VecGetSize(this->vector,pM);
         #else
                 this->vector->GetSize(pM);
@@ -221 +224 @@
                 
         #ifdef _HAVE_PETSC_
-                _assert_(this->vector);
-                VecGetLocalSize(this->vector,pM);
+                if(!this->vector)*pM=0;
+                else VecGetLocalSize(this->vector,pM);
         #else
                 this->vector->GetLocalSize(pM);
@@ -236 +239 @@
         #ifdef _HAVE_PETSC_
                 Vec vec_output=NULL;
-                _assert_(this->vector);
-                VecDuplicate(this->vector,&vec_output);
-                output->vector=vec_output;
+                if(!this->vector){
+                        VecDuplicate(this->vector,&vec_output);
+                        output->vector=vec_output;
+                }
         #else
                 output->vector=this->vector->Duplicate();
@@ -300 +304 @@
 
         #ifdef _HAVE_PETSC_
-                VecCopy(this->vector,to->vector);
+                if(this->vector) VecCopy(this->vector,to->vector);
         #else
                 this->vector->Copy(to->vector);