Changeset 26324

Timestamp:

06/15/21 08:30:57 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: use Pid instead of Sid to improve scaling

File:

• issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp

@@ -553 +553 @@
 
                                         if((surface_crevasse-surface>0. || crevassedepth-thickness>0.) && bed<0.){
-                                                vec_constraint_nodes->SetValue(node->Sid(),1.0,INS_VAL);
+                                                vec_constraint_nodes->SetValue(node->Pid(),1.0,INS_VAL);
                                         }
                                 }
@@ -582 +582 @@
                                 for(int in=0;in<numnodes;in++){
                                         Node* node=element->GetNode(in);
-                                        if(constraint_nodes[node->Sid()]==1.){
+                                        if(constraint_nodes[node->Pid()]==1.){
                                                 isconnected = true;
                                                 break;
@@ -600 +600 @@
                                                 surface_input->GetInputValue(&surface,gauss);
 
-                                                if((surface_crevasse-surface>0. || crevassedepth-thickness>0.) && bed<0. && levelset>-300. && levelset<0. && constraint_nodes[node->Sid()]==0.){
+                                                if((surface_crevasse-surface>0. || crevassedepth-thickness>0.) && bed<0. && levelset>-300. && levelset<0. && constraint_nodes[node->Pid()]==0.){
                                                         local_nflipped++;
-                                                        vec_constraint_nodes->SetValue(node->Sid(),1.0,INS_VAL);
+                                                        vec_constraint_nodes->SetValue(node->Pid(),1.0,INS_VAL);
                                                 }
                                         }
@@ -631 +631 @@
                                 if(!node->IsActive()) continue;
 
-                                if(constraint_nodes[node->Sid()]>0.){
+                                if(constraint_nodes[node->Pid()]>0.){
                                         node->ApplyConstraint(0,+1.);
                                 }