Changeset 25267

Timestamp:

07/10/20 23:01:28 (5 years ago)

Author:

Eric.Larour

Message:

CHG: added assert and fixed pretty big bug.

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -6058 +6058 @@
         else if(type==TransientInput2Enum){
 
-
                 IssmDouble* steps=NULL;
                 int nsteps;
@@ -6065 +6064 @@
 
                 /*retrieve transient input:*/
-                transientinput= this->inputs2->GetTransientInput(name);
-                transientinput2= this->inputs2->GetTransientInput(DummyEnum);
+                transientinput= this->inputs2->GetTransientInput(name); _assert_(transientinput);
+                transientinput2= this->inputs2->GetTransientInput(DummyEnum); _assert_(transientinput2);
 
                 /*retrieve time steps: */
@@ -6101 +6100 @@
                                 for (int j=0;j<NUMVERTICES;j++){
                                         values[j]=triainput->element_values[j];
-                                        if(partition[this->vertices[i]->Sid()]!=-1){
+                                        if(partition[this->vertices[j]->Sid()]!=-1){
                                                 if(nt==1) values[j]*=distributed_values[(int)partition[this->vertices[j]->Sid()]];//we scale all the time steps  with the same distributed_value 
                                                 else values[j]*=distributed_values[(int)partition[this->vertices[j]->Sid()]*nsteps+i];//we scale all the time steps with distributed value for each step