fem_types.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: Apache-2.0
 
#ifndef _FEM_TYPES_H
#define _FEM_TYPES_H
 
#include "basics.h"
#include "SF_vector.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
};
 
// 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 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;    
};
 
#if WITH_EMI_MODEL
struct RegionSpecs_EMI {
  char         *regname;      
  int           regID;        
  int           nsubregs;     
  std::string   *subregtags;   
  physMaterial *material;     
  void         *customize;    
};
#endif
 
// 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;      
  int                     el_scale_dpn = 0; 
};
 
#if WITH_EMI_MODEL
struct MaterialType_EMI {
  // rndFluct      rfluct;            //!< random material variations
  // double       *fluct     = NULL;  //!< random fluctuation vector
 
  SF::vector<int>         regionIDs;  
  SF::vector<RegionSpecs_EMI> regions;    
  SF::vector<double>      el_scale;      
  int                     el_scale_dpn = 0; 
};
#endif
 
}  // namespace opencarp
 
#endif