source: issm/trunk/src/c/modules/InputConvergencex/InputConvergencex.cpp@ 13975

Last change on this file since 13975 was 13975, checked in by Mathieu Morlighem, 12 years ago

merged trunk-jpl and trunk for revision 13974

File size: 1.5 KB
Line 
1/*!\file InputConvergencex
2 */
3
4#include "./InputConvergencex.h"
5#include "../../shared/shared.h"
6#include "../../include/include.h"
7#include "../../io/io.h"
8#include "../../toolkits/toolkits.h"
9#include "../../EnumDefinitions/EnumDefinitions.h"
10
11bool InputConvergencex(Elements* elements,Nodes* nodes,Vertices* vertices,Loads* loads,Materials* materials,Parameters* parameters,int* enums, int num_enums, int* criterionenums, IssmDouble* criterionvalues,int num_criterionenums){
12
13 /*intermediary:*/
14 int i;
15 bool converged;
16 int num_notconverged = 0;
17 int total_notconverged;
18 IssmDouble *eps = NULL;
19 Element *element = NULL;
20
21 /*allocate dynamic memory: */
22 eps=xNew<IssmDouble>(num_criterionenums);
23
24 /*Go through elements, and ask them to do the job: */
25 for(i=0;i<elements->Size();i++){
26 element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
27 converged=element->InputConvergence(eps,enums,num_enums,criterionenums,criterionvalues,num_criterionenums);
28 if(!converged) num_notconverged++;
29 }
30
31 /*In parallel, we need to gather the converged status: */
32 #ifdef _HAVE_MPI_
33 MPI_Allreduce ( (void*)&num_notconverged,(void*)&total_notconverged,1,MPI_INT,MPI_SUM,IssmComm::GetComm());
34 num_notconverged=total_notconverged;
35 #endif
36 if(VerboseConvergence()) _pprintLine_(" #elements above convergence criterion = " << num_notconverged);
37
38 /*Free ressources:*/
39 xDelete<IssmDouble>(eps);
40
41 /*return: */
42 if (num_notconverged){
43 return false;}
44 else{
45 return true;}
46}
Note: See TracBrowser for help on using the repository browser.