Changeset 17229

Timestamp:

02/07/14 09:39:17 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed spawn input in case interpolation is P0

Location:

issm/trunk-jpl/src/c/classes/Inputs

Files:

• PentaInput.cpp (modified) (1 diff)
• TriaInput.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp

@@ -94 +94 @@
                 IssmDouble newvalues[3]; 
 
-                /*Create arrow of indices depending on location (0=base 1=surface)*/
+                /*Create array of indices depending on location (0=base 1=surface)*/
                 int indices[3];
                 switch(location){

issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp

@@ -100 +100 @@
         /*output*/
         SegInput* outinput=NULL;
-        IssmDouble newvalues[2]; //Assume P1 interpolation only for now
-
-        /*Create arrow of indices depending on location (0=base 1=surface)*/
-
-        newvalues[0]=this->values[index1];
-        newvalues[1]=this->values[index2];
-
-        /*Create new Seg input*/
-        outinput=new SegInput(this->enum_type,&newvalues[0],P1Enum);
+
+        if(this->element_type==P0Enum){
+                outinput=new SegInput(this->enum_type,&this->values[0],P0Enum);
+        }
+        else{
+                /*Assume P1 interpolation only for now*/
+                IssmDouble newvalues[2];
+
+                /*Create array of indices depending on location (0=base 1=surface)*/
+                newvalues[0]=this->values[index1];
+                newvalues[1]=this->values[index2];
+
+                /*Create new Seg input*/
+                outinput=new SegInput(this->enum_type,&newvalues[0],P1Enum);
+        }
 
         /*Assign output*/