leadfield.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 _LEADFIELD_H
#define _LEADFIELD_H
 
#include <string>
#include <vector>
 
#include "fem.h"
#include "stimulate.h"
 
namespace opencarp {
 
class Electrics;
 
class Leadfield
{
 public:
  Leadfield();
  ~Leadfield();
 
  int run(Electrics & elec);
 
 private:
  std::vector<std::string>        lead_names;
  std::vector<std::vector<float>> trace_data;
  std::vector<sf_vec*>            Z;         
  sf_vec*                         Im         = nullptr;
  std::string                     vm_file;
  float                           ecg_timedt = 0.0f;
  std::string                     lf_dir;
 
  void construct(Electrics & elec);
  void load(Electrics & elec);
  void apply(sf_mat & Ki, sf_vec & vm);
  void apply_unit_curr_lead(const std::string & lead, Electrics & elec);
  void write_leadfield_igb(const std::string & lead_name, sf_vec * phie_i);
  void collect_lead_names(const SF::vector<stimulus> & stimuli);
};
 
} // namespace opencarp
 
#endif