Go to the source code of this file.
◆ KillIcebergsx()
void KillIcebergsx |
( |
FemModel * |
femmodel | ) |
|
Definition at line 11 of file KillIcebergsx.cpp.
18 IssmDouble *local_mask = xNewZeroInit<IssmDouble>(nbv_local);
26 element_flag[i] =
true;
31 for(
int v=0;v<numvertices;v++) local_mask[element->
vertices[v]->
Lid()] = 1.;
40 bool keepsyncing =
true;
47 bool keepgoing =
true;
61 for(
int j=0;j<numvertices;j++){
63 if(local_mask[lid]>0.){
69 element_flag[i] =
true;
70 for(
int j=0;j<numvertices;j++){
72 if(local_mask[lid]==0.){
95 xDelete<bool>(element_flag);
103 bool deactivate =
false;
104 for(
int j=0;j<numvertices;j++){
106 if(local_mask[lid]==0.){
114 IssmDouble* values = xNew<IssmDouble>(numvertices);
115 for(
int j=0;j<numvertices;j++) values[j] = 1.;
117 xDelete<IssmDouble>(values);
123 xDelete<IssmDouble>(local_mask);
int ISSM_MPI_Reduce(void *sendbuf, void *recvbuf, int count, ISSM_MPI_Datatype datatype, ISSM_MPI_Op op, int root, ISSM_MPI_Comm comm)