Matice.h

Go to the documentation of this file.

 
#ifndef MATICE_H_
#define MATICE_H_
 
/*Headers:*/
/*{{{*/
#include "./Material.h"
#include "../Hook.h"
class IoModel;
class Elements;
class Element;
class Loads;
class Nodes;
class Vertices;
class Materials;
class Parameters;
class Gauss;
class Input;
/*}}}*/
 
class Matice: public Material{
 
   private: 
      int      mid;
      bool     isdamaged;
      bool     isenhanced;
      int      rheology_law;
      Hook    *helement;
      Element *element;
 
   public:
      /*Matice constructors, destructors: {{{*/
      Matice();
      Matice(int mid,int i, IoModel* iomodel);
      Matice(int mid,int i, int materialtype);
      ~Matice();
      void Init(int mid,int i, int materialtype);
      /*}}}*/
      /*Object virtual functions definitions:{{{ */
      Object* copy();
      void  DeepEcho();
      void  Echo();
      int   Id(); 
      void Marshall(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction);
      int   ObjectEnum();
      /*}}}*/
      /*Material virtual functions resolution: {{{*/
      void   Configure(Elements* elements);
      Material*  copy2(Element* element);
      void       GetViscosity(IssmDouble* pviscosity, IssmDouble eps_eff,Gauss* gauss);
      void       GetViscosityBar(IssmDouble* pviscosity, IssmDouble eps_eff,Gauss* gauss);
      void       GetViscosityComplement(IssmDouble* pviscosity_complement, IssmDouble* pepsilon,Gauss* gauss);
      void       GetViscosityDComplement(IssmDouble*, IssmDouble*,Gauss* gauss);
      void       GetViscosityDerivativeEpsSquare(IssmDouble* pmu_prime, IssmDouble* pepsilon,Gauss* gauss);
      void       GetViscosity_B(IssmDouble* pviscosity, IssmDouble eps_eff,Gauss* gauss);
      void       GetViscosity_D(IssmDouble* pviscosity, IssmDouble eps_eff,Gauss* gauss);
      void       GetViscosity2dDerivativeEpsSquare(IssmDouble* pmu_prime, IssmDouble* pepsilon,Gauss* gauss);
      IssmDouble GetA(Gauss* gauss);
      IssmDouble GetAbar(Gauss* gauss);
      IssmDouble GetB(Gauss* gauss);
      IssmDouble GetBbar(Gauss* gauss);
      IssmDouble GetD(Gauss* gauss);
      IssmDouble GetDbar(Gauss* gauss);
      IssmDouble GetE(Gauss* gauss);
      IssmDouble GetEbar(Gauss* gauss);
      IssmDouble GetN();
      bool       IsDamage();
      bool       IsEnhanced();
      void       ResetHooks();
      void       SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters);
 
      void       ViscosityFSDerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon,Gauss* gauss);
      void       ViscosityHODerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon,Gauss* gauss);
      void       ViscositySSADerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon,Gauss* gauss);
 
      void       ViscosityFS(IssmDouble* pviscosity,int dim,IssmDouble* xyz_list,Gauss* gauss,Input2* vx_input,Input2* vy_input,Input2* vz_input);
      void       ViscosityHO(IssmDouble* pviscosity,int dim,IssmDouble* xyz_list,Gauss* gauss,Input2* vx_input,Input2* vy_input);
      void       ViscosityL1L2(IssmDouble* pviscosity,IssmDouble* xyz_list,Gauss* gauss,Input2* vx_input,Input2* vy_input,Input2* surf);
      void       ViscositySSA(IssmDouble* pviscosity,int dim,IssmDouble* xyz_list,Gauss* gauss,Input2* vx_input,Input2* vy_input);
      void       ViscosityBFS(IssmDouble* pmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input2* vx_input,Input2* vy_input,Input2* vz_input,IssmDouble eps_eff){_error_("not supported");};
      void       ViscosityBHO(IssmDouble* pmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input2* vx_input,Input2* vy_input,IssmDouble eps_eff){_error_("not supported");};
      void       ViscosityBSSA(IssmDouble* pmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input2* vx_input,Input2* vy_input,IssmDouble eps_eff){_error_("not supported");};
      /*}}}*/
};
 
#endif  /* _MATICE_H_ */