electric_integrators.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: LicenseRef-APL-1.1
 
#ifndef _ELEC_INTEGRATORS_H
#define _ELEC_INTEGRATORS_H
 
#include "SF_base.h"
#include "basics.h"
#include "fem_types.h"
#include "sf_interface.h"
 
namespace opencarp {
 
void read_el_scale_vec(const char* file, mesh_t mt, SF::vector<double>& el_scale, int& el_scale_dpn);
 
#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