ionics.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: LicenseRef-APL-1.1
 
#ifndef _IONICS_H
#define _IONICS_H
 
#include "physics_types.h"
#include "sim_utils.h"
#include "stimulate.h"
#include "IGBheader.h"
#include "timers.h"
 
#include "sf_interface.h"
#include "MULTI_ION_IF.h"
#include "sv_init.h"
 
namespace opencarp {
 
 
void region_mask(mesh_t meshspec, SF::vector<RegionSpecs> & regspec,
                 SF::vector<int> & regionIDs, bool mask_elem, const char* reglist,
                 bool warn_on_default_tags = true);
 
class Ionics : public Basic_physic
{
  public:
  limpet::MULTI_IF* miif = NULL;
  mesh_t    ion_domain;
 
  generic_timing_stats comp_stats;
 
  Ionics(mesh_t gid) : ion_domain(gid)
  {
    switch(ion_domain) {
      case intra_elec_msh:
        name = "Myocard Ionics"; break;
      default:
        name = "Ionics"; break;
    }
  }
 
  void initialize();
  void destroy();
  void compute_step();
  void output_step();
  ~Ionics() = default;
 
  double timer_val(const int timer_id);
  std::string timer_unit(const int timer_id);
 
  private:
  double setup_MIIF(limpet::node_count_t nnodes, int nreg, IMPregion* impreg, int* mask,
             const char *start_fn, int numadjust, IMPVariableAdjustment *adjust,
             double time_step, bool close);
 
};
 
void compute_IIF(limpet::IonIfBase& pIF, limpet::GlobalData_t** impdata, limpet::node_index_t n);
 
struct sv_data {
  char*        name;        
  char**       imps;        
  char**       svNames;     
  int*         svInds;      
  int*         svSizes;     
  int*         svOff;       
  FILE_SPEC    outf;        
  void**       getsv;       
  sf_vec* ordered;     
  char*        units;       
  float        bogus;       
};
 
struct gvec_data {
  unsigned int           nRegs;        
  bool                   inclPS;       
  SF::vector<sv_data>    vecs;         
  SF::vector<int>        plugin_idx;   
  bool                   rm_dumps;     
};
 
void init_sv_gvec(gvec_data& GVs, limpet::MULTI_IF* miif, sf_vec & tmpl,
                  igb_output_manager & output_manager);
void assemble_sv_gvec(gvec_data & gvecs, limpet::MULTI_IF *miif);
}
#endif