Changeset 24630 for issm/trunk-jpl/src/c/modules/ConstraintsStatex/RiftConstraintsState.cpp

Timestamp:

03/03/20 12:19:15 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: speeding up ConstraintState if no penalties and no rift found

File:

• issm/trunk-jpl/src/c/modules/ConstraintsStatex/RiftConstraintsState.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/ConstraintsStatex/RiftConstraintsState.cpp

@@ -6 +6 @@
 
 /*current module: */
-/*RiftIsPresent(Loads* loads,int configuration_type){{{*/
-int RiftIsPresent(Loads* loads,int configuration_type){
-
-        int i;
-
-        int found=0;
-        int mpi_found=0;
-
-        /*go though loads, and figure out if one of the loads is a Riftfront: */
-        for (i=0;i<loads->Size();i++){
-                Load* load=(Load*)loads->GetObjectByOffset(i);
-                if(RiftfrontEnum==loads->GetEnum(i)){
-                        found=1;
-                        break;
-                }
-        }
-
-        ISSM_MPI_Reduce (&found,&mpi_found,1,ISSM_MPI_INT,ISSM_MPI_SUM,0,IssmComm::GetComm() );
-        ISSM_MPI_Bcast(&mpi_found,1,ISSM_MPI_INT,0,IssmComm::GetComm());                
-        found=mpi_found;
-
-        return found;
-}
-/*}}}*/
 /*RiftConstraintsState(int* pconverged, int* pnum_unstable_constraints,Loads* loads,int min_mechanical_constraints,int configuration_type){{{*/
 void RiftConstraintsState(int* pconverged, int* pnum_unstable_constraints,Loads* loads,int min_mechanical_constraints,int configuration_type){