Masscon.h

Go to the documentation of this file.

 
#ifndef _MASSCON_H_
#define _MASSCON_H_
 
/*Headers:*/
/*{{{*/
#include "./Definition.h"
#include "../datastructures/datastructures.h"
#include "./Elements/Element.h"
#include "./Elements/Elements.h"
#include "./FemModel.h"
#include "../classes/Params/Parameters.h"
/*}}}*/
 
class Masscon: public Object, public Definition{
 
   public: 
 
      int         definitionenum;
      char*       name;
      IssmDouble* levelset;
      int         M;
      
      /*Masscon constructors, destructors :*/
      Masscon(){/*{{{*/
 
         this->definitionenum = -1;
         this->name = NULL;
         this->levelset=NULL;
         this->M=0;
 
      }
      /*}}}*/
      Masscon(char* in_name, int in_definitionenum, IssmDouble* levelsetin, int Min){ /*{{{*/
 
         this->definitionenum=in_definitionenum;
         this->name   = xNew<char>(strlen(in_name)+1);
         xMemCpy<char>(this->name,in_name,strlen(in_name)+1);
 
         this->levelset   = xNew<IssmDouble>(Min);
         xMemCpy<IssmDouble>(this->levelset, levelsetin, Min);
 
         this->M=Min;
 
      }
      /*}}}*/
      ~Masscon(){/*{{{*/
         if(this->name)xDelete(this->name); 
         if(this->levelset)xDelete(this->levelset);
      }
      /*}}}*/
      /*Object virtual function resolutoin: */
      Object* copy() {/*{{{*/
         Masscon* mf = new Masscon(this->name,this->definitionenum,this->levelset,this->M);
         return (Object*) mf;
      }
      /*}}}*/
      void DeepEcho(void){/*{{{*/
         this->Echo();
      }
      /*}}}*/
      void Echo(void){/*{{{*/
         _printf_(" Masscon: " << this->name << " " << this->definitionenum << "\n");
         _printf_("    levelset: " << this->levelset << "\n");
         _printf_("    M: " << this->M << "\n");
      }
      /*}}}*/
      int Id(void){/*{{{*/
         return -1;
      }
      /*}}}*/
      void Marshall(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction){/*{{{*/
         _error_("not implemented yet!"); 
      } 
      /*}}}*/
      int ObjectEnum(void){/*{{{*/
         return MassconEnum;
      }
      /*}}}*/
      /*Definition virtual function resolutoin: */
      int DefinitionEnum(){/*{{{*/
 
         return this->definitionenum;
      }
      /*}}}*/
      char* Name(){/*{{{*/
 
         char* name2=xNew<char>(strlen(this->name)+1);
         xMemCpy(name2,this->name,strlen(this->name)+1);
 
         return name2;
      }
      /*}}}*/
       IssmDouble Response(FemModel* femmodel){/*{{{*/
 
          int i;
          IssmDouble mass_t=0.;
          IssmDouble all_mass_t=0.;
 
 
          for(i=0;i<femmodel->elements->Size();i++){
             Element* element=(Element*)femmodel->elements->GetObjectByOffset(i);
             mass_t+=element->Masscon(this->levelset);
          }
 
          ISSM_MPI_Allreduce ( (void*)&mass_t,(void*)&all_mass_t,1,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,IssmComm::GetComm());
          mass_t=all_mass_t;
          
          return mass_t;
       }
         /*}}}*/
};
 
#endif  /* _MASSCON_H_ */