leadfield.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: Apache-2.0
 
#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