Ignore:
Timestamp:
06/18/10 20:48:42 (15 years ago)
Author:
Eric.Larour
Message:

Diagnostic solution now compiles in parallel! yeehah!

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Elements/Penta.cpp

    r4055 r4057  
    53165316/*}}}*/
    53175317/*FUNCTION Penta::InputConvergence(int* pconverged, double* eps, int* enums,int num_enums,int* criterionenums,double* criterionvalues,int num_criterionenums){{{1*/
    5318 void  Penta::InputConvergence(int* pconverged,double* eps, int* enums,int num_enums,int* criterionenums,double* criterionvalues,int num_criterionenums){a
    5319 
     5318void  Penta::InputConvergence(int* pconverged,double* eps, int* enums,int num_enums,int* criterionenums,double* criterionvalues,int num_criterionenums){
     5319
     5320        int     i;
    53205321        Input** new_inputs=NULL;
    53215322        Input** old_inputs=NULL;
     
    53265327       
    53275328        for(i=0;i<num_enums/2;i++){
    5328                 new_inputs[i]=(Input*)this->inputs->GetInput(new_enums[2*i+0]);
    5329                 old_inputs[i]=(Input*)this->inputs->GetInput(old_enums[2*i+1]);
    5330                 if(!new_inputs[i])ISSMERROR("%s%s"," could not find input with enum ",EnumAsString(new_enums[2*i+0]));
    5331                 if(!old_inputs[i])ISSMERROR("%s%s"," could not find input with enum ",EnumAsString(old_enums[2*i+0]));
     5329                new_inputs[i]=(Input*)this->inputs->GetInput(enums[2*i+0]);
     5330                old_inputs[i]=(Input*)this->inputs->GetInput(enums[2*i+1]);
     5331                if(!new_inputs[i])ISSMERROR("%s%s"," could not find input with enum ",EnumAsString(enums[2*i+0]));
     5332                if(!old_inputs[i])ISSMERROR("%s%s"," could not find input with enum ",EnumAsString(enums[2*i+0]));
    53325333        }
    53335334
    53345335        /*ok, we've got the inputs (new and old), now loop throught the number of criterions and fill the eps array:*/
    53355336        for(i=0;i<num_criterionenums;i++){
    5336                 IsInputConverged(eps+i,new_inputs,old_inputs,num_enums/2,,criterionenums[i]);
     5337                IsInputConverged(eps+i,new_inputs,old_inputs,num_enums/2,criterionenums[i]);
    53375338                if(eps[i]>criterionvalues[i]) converged=0;
    53385339        }
Note: See TracChangeset for help on using the changeset viewer.