PacemakerPoissonSolver.cpp

Go to the documentation of this file.

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   This file is part of CEPS.
 
   CEPS is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.
 
   CEPS is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with CEPS (see file LICENSE at root of project).
   If not, see <https://www.gnu.org/licenses/>.
 
 
   Copyright 2019-2024 Inria, Universite de Bordeaux
 
   Authors, in alphabetical order:
 
     Pierre-Elliott BECUE, Florian CARO, Yves COUDIERE(*), Andjela DAVIDOVIC, 
     Charlie DOUANLA-LONTSI, Marc FUENTES, Mehdi JUHOOR, Michael LEGUEBE(*), 
     Pauline MIGERDITICHAN, Valentin PANNETIER(*), Nejib ZEMZEMI.
     * : currently active authors
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#include "cardiac/pacemaker/solver/PacemakerPoissonSolver.hpp"
 
#include "cardiac/assembler/FE/FECardiacBCAssembler.hpp"
#include "common/toolbox/VariableTimeStepper.hpp"
#include "cardiac/pacemaker/assembler/FE/FEPacemakerPoissonAssembler.hpp"
#include "cardiac/pacemaker/problem/PacemakerPoissonProblem.hpp"
#include "pde/assemblers/FE/FEIntegrator.hpp"
 
PacemakerPoissonSolver::PacemakerPoissonSolver(PacemakerPoissonProblem* problem) :
    HeatSolver             (problem),
    AbstractPacemakerSolver(problem,problem->getParameters())
{
  // renew the time stepper with a variable one
  ceps::destroyObject(m_timeStepper);
  m_timeStepper = ceps::getNew<VariableTimeStepper>(
    problem->getStartTime(), problem->getEndTime(), problem->getTimeStep()
  );
 
  if (ceps::isValidPtr(m_timeWriter))
  {
    m_timeWriter->addCustomData("Current I (uA)", &m_current);
    m_timeWriter->addCustomData("Measured voltage (mV)", &m_Vcroco);
    m_timeWriter->addCustomData("Mean of u on anode (mV)", &m_meanOfPotentialOnAnode);
    m_timeWriter->addCustomData("Mean of u on cathode (mV)", &m_meanOfPotentialOnCathode);
    m_timeWriter->addCustomData("Phase", &m_phase);
  }
}
 
// void
// PacemakerPoissonSolver::setupWithParameters(InputParameters* params)
// {
// }
 
void
PacemakerPoissonSolver::initializeAssemblers()
{
  auto pb = getPacemakerPoissonProblem();
 
  auto opAsb = ceps::getNew<FEPacemakerPoissonAssembler>(pb, m_fe);
  opAsb->setKForUnknown(pb->getUFieldUnknown(), pb->getDiffusionTensor());
  m_opAsb = opAsb;
 
  m_bcAsb = ceps::getNew<FEDivKGradBCAssembler>(m_fe);
}
 
PacemakerPoissonProblem*
PacemakerPoissonSolver::getPacemakerPoissonProblem() const
{
  return ceps::runtimeCast<PacemakerPoissonProblem*>(getTimedProblem());
}
 
void
PacemakerPoissonSolver::assembleAndSolve()
{
  CepsReal time = m_timeStepper->getTime();
  // CepsReal previous = m_timeStepper->getTimeStep();
  CepsReal dt       = getRequiredTimeStep(time,getPacemakerPoissonProblem()->getTimeStep());
 
  if (not ceps::equals(m_timeStepper->getTimeStep(), dt))
    m_timeStepper->setTimeStep(dt);
 
  HeatSolver::assembleAndSolve();
 
  // Compute additionnal values
  auto pb                    = getPacemakerPoissonProblem();
  auto anode                 = pb->getAnode();
  auto cathode               = pb->getCathode();
  auto u                     = pb->getUFieldUnknown();
  m_current                  = getPacemakerCurrent(m_uNp1, m_fe, time + dt);
  m_Vcroco                   = getAlligatorVoltage(m_uNp1, m_fe, time + dt);
  m_meanOfPotentialOnAnode   = getMeanOnElectrodeOf(m_uNp1, u, m_fe, anode);
  m_meanOfPotentialOnCathode = getMeanOnElectrodeOf(m_uNp1, u, m_fe, cathode);
  m_phase                    = static_cast<CepsReal>(pb->getPhase(time + dt));
}