fem_types.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
// openCARP is an open cardiac electrophysiology simulator.
//
// Copyright (C) 2020 openCARP project
//
// This program is licensed under the openCARP Academic Public License (APL)
// v1.0: You can use and redistribute it and/or modify it in non-commercial
// academic environments under the terms of APL as published by the openCARP
// project v1.0, or (at your option) any later version. Commercial use requires
// a commercial license (info@opencarp.org).
//
// This program is distributed without any warranty; see the openCARP APL for
// more details.
//
// You should have received a copy of the openCARP APL along with this program
// and can find it online: http://www.opencarp.org/license
// ----------------------------------------------------------------------------
 
#ifndef _FEM_TYPES_H
#define _FEM_TYPES_H
 
#include <map>
#include "basics.h"
#include "SF_vector.h"
//#include "constitutive_model.h"
//#include "constitutive_model_library.h"
 
namespace opencarp {
 
enum physMat_t {PhysMat=0, ElecMat=1, DiffuseMat=2, MechMat=3, OpticsMat=4};
 
enum cond_t {
  intra_cond, extra_cond, sum_cond, para_cond
};
 
enum mechMat_t {
  HOOKE, NEOHOOKE, MOONEYRIVLIN, GUCCIONE, USYK
};
 
 
// typedef CondTensModifier bidmT;
 
struct physMaterial {
  physMat_t     material_type; 
  double        fluct;         
};
 
 
struct elecMaterial : physMaterial
{
  double  InVal[3];    
  double  ExVal[3];    
  double  BathVal[3];  
  cond_t  g;           
  int    *augment;     
  int     nAug;        
};
 
 
struct diffMaterial : physMaterial
{
  double InVal[3];     
  double ExVal[3];     
  double BathVal[3];   
  double clearance;    
  double basal_conc;   
};
 
struct mechMaterial : physMaterial
{
  mechMat_t mat_t;              
  std::map<std::string, mechMat_t> mMat_map = {
    {"Hooke", mechMat_t::HOOKE},
    {"NeoHooke", mechMat_t::NEOHOOKE},
    {"MooneyRivlin", mechMat_t::MOONEYRIVLIN},
    {"Guccione", mechMat_t::GUCCIONE},
    {"Usyk", mechMat_t::USYK}
  }; //necessary to map command line arguments to material type
  inline void set_constitutive_model(std::string model_name)
  {
    mat_t = mMat_map[model_name];
    switch (mat_t) {
      case NEOHOOKE:
        log_msg(0, 0, 0, "Neo-Hooke chosen.");
        break;
      case HOOKE:
        log_msg(0, 0, 0, "Hooke not yet implemented.");
        break;
      case MOONEYRIVLIN:
        log_msg(0, 0, 0, "Mooney-Rivlin not yet implemented.");
        break;
      case GUCCIONE:
        log_msg(0, 0, 0, "Guccione chosen.");
        break;
      case USYK:
        log_msg(0, 0, 0, "Usyk not yet implemented");
      default:
        log_msg(0,0,0, "Choose an existing material model pls.");
        break;
    }
  }
};
 
struct opticsMaterial : physMaterial
{
  int    mode;         
  double Dil[3];       
  double Dem[3];       
  double muAil;        
  double muAem;        
};
 
// TODO: Write destructor for Regions
 
struct RegionSpecs {
  char         *regname;      
  int           regID;        
  int           nsubregs;     
  int          *subregtags;   
  physMaterial *material;     
  void         *customize;    
};
 
 
// struct rndFluct {
//   long          rseed;        //!< random seed
//   double        var;          //!<
// };
 
struct MaterialType {
  // rndFluct      rfluct;            //!< random material variations
  // double       *fluct     = NULL;  //!< random fluctuation vector
 
  SF::vector<int>         regionIDs;  
  SF::vector<RegionSpecs> regions;    
  SF::vector<double>      el_scale;      
};
 
}  // namespace opencarp
 
#endif