electric_integrators.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 _ELEC_INTEGRATORS_H
#define _ELEC_INTEGRATORS_H

#include "SF_base.h"
#include "basics.h"
#include "fem_types.h"

namespace opencarp {

#define EP_MAX_IPOINTS 128
#define EP_MAX_LPOINTS 8

class elec_stiffness_integrator : public SF::matrix_integrator<mesh_int_t,mesh_real_t>
{
  MaterialType & material;   

  SF::dmat<double> cond;          
  SF::dmat<double> shape;         
  SF::dmat<double> gshape;        
  double J[9];                    
  SF::Point lpts[EP_MAX_LPOINTS]; 

  int p_order = 1;  

  SF::Point           ipts[EP_MAX_IPOINTS];    
  double              w   [EP_MAX_IPOINTS];    

  public:
  elec_stiffness_integrator(MaterialType & inp_mat) : material(inp_mat)
  {}

  void operator() (const SF::element_view<mesh_int_t,mesh_real_t> & elem, SF::dmat<double> & buff);
  void dpn(mesh_int_t & row_dpn, mesh_int_t & col_dpn);
};

class mass_integrator : public SF::matrix_integrator<mesh_int_t,mesh_real_t>
{
  SF::dmat<double> shape;         
  double J[9];                    
  SF::Point lpts[EP_MAX_LPOINTS]; 

  int p_order = 1;  

  SF::Point           ipts[EP_MAX_IPOINTS];    
  double              w   [EP_MAX_IPOINTS];    

  public:
  void operator() (const SF::element_view<mesh_int_t,mesh_real_t> & elem, SF::dmat<double> & buff);
  void dpn(mesh_int_t & row_dpn, mesh_int_t & col_dpn);
};

}  // namespace opencarp

#endif