source: issm/trunk/src/c/DataSet/Loads.cpp@ 4219

/*
 * \file Loads.c
 * \brief: implementation of the Loads class, derived from DataSet class
 */

/*Headers: {{{1*/
#ifdef HAVE_CONFIG_H
        #include "config.h"
#else
#error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
#endif

#include <vector>
#include <functional>
#include <algorithm>
#include <iostream>

#include "./DataSet.h"
#include "../shared/shared.h"
#include "../include/include.h"
#include "../EnumDefinitions/EnumDefinitions.h"

using namespace std;
/*}}}*/

/*Object constructors and destructor*/
/*FUNCTION Loads::Loads(){{{1*/
Loads::Loads(){
        return;
}
/*}}}*/
/*FUNCTION Loads::Loads(int in_enum){{{1*/
Loads::Loads(int in_enum): DataSet(in_enum){
        //do nothing;
        return;
}
/*}}}*/
/*FUNCTION Loads::~Loads(){{{1*/
Loads::~Loads(){
        return;
}
/*}}}*/

/*Numerics:*/
/*FUNCTION Loads::MeltingIsPresent{{{1*/
int   Loads::MeltingIsPresent(){

        int i;
        int found=0;
        int mpi_found=0;

        for(i=0;i<this->Size();i++){
                Object* object=(Object*)this->GetObjectByOffset(i);
                if (object->Enum()==PengridEnum){
                        found=1;
                        break;
                }
        }
        
        #ifdef _PARALLEL_
        MPI_Reduce (&found,&mpi_found,1,MPI_INT,MPI_SUM,0,MPI_COMM_WORLD );
        MPI_Bcast(&mpi_found,1,MPI_INT,0,MPI_COMM_WORLD);                
        found=mpi_found;
        #endif

        return found;
}
/*}}}*/
/*FUNCTION Loads::MeltingConstraints{{{1*/
void  Loads::MeltingConstraints(int* pconverged, int* pnum_unstable_constraints){

        int i;

        int unstable=0;
        int num_unstable_constraints=0;
        int converged=0;
        int sum_num_unstable_constraints=0;

        num_unstable_constraints=0;     

        /*Enforce constraints: */
        for(i=0;i<this->Size();i++){
                Object* object=(Object*)this->GetObjectByOffset(i);
                if(object->Enum()==PengridEnum){

                        Pengrid* pengrid=(Pengrid*)object;
                        pengrid->PenaltyConstrain(&unstable);
                        num_unstable_constraints+=unstable;
                }
        }

        #ifdef _PARALLEL_
        MPI_Reduce (&num_unstable_constraints,&sum_num_unstable_constraints,1,MPI_INT,MPI_SUM,0,MPI_COMM_WORLD );
        MPI_Bcast(&sum_num_unstable_constraints,1,MPI_INT,0,MPI_COMM_WORLD);                
        num_unstable_constraints=sum_num_unstable_constraints;
        #endif

        /*Have we converged? : */
        if (num_unstable_constraints==0) converged=1;
        else converged=0;

        /*Assign output pointers: */
        *pconverged=converged;
        *pnum_unstable_constraints=num_unstable_constraints;
}
/*}}}*/