source: issm/trunk-jpl/src/c/modules/ModelProcessorx/Damage/UpdateElementsDamage.cpp@ 16188

/*
 * UpdateElementsDamage:
 */

#include "../../../toolkits/toolkits.h"
#include "../../../classes/classes.h"
#include "../../../shared/shared.h"
#include "../../MeshPartitionx/MeshPartitionx.h"
#include "../ModelProcessorx.h"

void    UpdateElementsDamage(Elements* elements, IoModel* iomodel,int analysis_counter,int analysis_type){

        /*Update elements: */
        int counter=0;
        for(int i=0;i<iomodel->numberofelements;i++){
                if(iomodel->my_elements[i]){
                        Element* element=(Element*)elements->GetObjectByOffset(counter);
                        element->Update(i,iomodel,analysis_counter,analysis_type,P1Enum);
                        counter++;
                }
        }

        /*What input do I need to run my damage evolution model?*/
        iomodel->FetchDataToInput(elements,VxEnum);
        iomodel->FetchDataToInput(elements,VyEnum);
        iomodel->FetchDataToInput(elements,VzEnum);
        iomodel->FetchDataToInput(elements,DamageDEnum);
        iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum);
        iomodel->FetchDataToInput(elements,PressureEnum);

        bool dakota_analysis;
        iomodel->Constant(&dakota_analysis,QmuIsdakotaEnum);
        if(dakota_analysis){
                elements->InputDuplicate(DamageDEnum, QmuDamageDEnum);
        }
}