Changeset 4385

Timestamp:

07/01/10 10:10:41 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added check in GetSolutionFromInputsx that gssize>0

Location:

issm/trunk/src/c

Files:

• modules/GetSolutionFromInputsx/GetSolutionFromInputsx.cpp (modified) (1 diff)
• objects/Elements/Tria.cpp (modified) (1 diff)
• shared/Elements/elements.h (modified) (2 diffs)

issm/trunk/src/c/modules/GetSolutionFromInputsx/GetSolutionFromInputsx.cpp

@@ -32 +32 @@
         /*Get size of vector: */
         gsize=nodes->NumberOfDofs(analysis_type);
+        if (gsize==0) ISSMERROR("Allocating a Vec of size 0 as gsize=0 for analysis: %s",EnumAsString(analysis_type));
         
         /*Initialize solution: */

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -5212 +5212 @@
 
         /*Add value to global vector*/
+        printarray(&doflist[0],numdof);
         VecSetValues(solution,numdof,doflist,(const double*)values,INSERT_VALUES);
 

issm/trunk/src/c/shared/Elements/elements.h

@@ -12 +12 @@
 double Paterson(double temperature);
 int GetVerticesCoordinates(double* xyz,  Node** nodes, int numgrids);
-inline void printarray(double* array,int lines,int cols){
+inline void printarray(double* array,int lines,int cols=1){
         int i,j;
         printf("\n");
@@ -24 +24 @@
         printf("\n");
 }
-inline void printarray(int* array,int lines,int cols){
+inline void printarray(int* array,int lines,int cols=1){
         int i,j;
         printf("\n");