CardiacSolver.hpp

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
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#pragma once
 
#include "pde/solver/timed/HeatSolver.hpp"
#include "ode/solver/IonicSolver.hpp"
#include "cardiac/problem/CardiacProblem.hpp"
#include "cardiac/common/FECardiacParameters.hpp"
#include "cardiac/common/ActivationTracker.hpp"
 
class CardiacSolver : public HeatSolver
{
 
  public:
  
    virtual ~CardiacSolver();
 
    CardiacSolver(CardiacProblem* problem);
 
    void 
    solve() override;
 
    CepsString 
    getResultsLocationMessage() const;
 
    CardiacProblem* 
    getCardiacProblem() const;
 
    ActivationTracker*
    getActivationTracker() const;
 
  protected:
 
    void
    assembleAndSolve() override;
 
    CepsBool 
    finished() const override;
 
    // --------------------====
    // LHS computation
 
    // /// @brief Computes the stiffness matrix including conductivity
    // void computeCardiacStiffnessMatrix();
 
    // /// @brief Computes the stiffness matrix including conductivity
    // void computeCardiacStiffnessMatrixForSources();
 
    // --------------------====
    // RHS compute and solve
 
    void fillSourceTermsVector(CepsReal tn) override;
 
    // /// @brief Updates matrix and RHS for RK numerical schemes (actually performs all RK substeps)
    // void updateAndSolveRK();
 
    // -------------------------==
    // I/O
 
    void output(DHVecPtr dummy, CepsBool immediatWriting = true) override;
 
 
  protected:
 
    // Shorcuts
    ActivationTracker* m_activationTracker; 
 
    // Ionic solvers
    CepsVector<IonicSolver*> m_ionSolvers;  
    CepsString               m_ionSolvOpts; 
 
    // Output stuff
    CepsBool               m_writeIIon;                
    CepsBool               m_writeIapp;                
    CepsReal               m_postProcessPeriod;        
    std::set<CepsUInt>     m_saveDepolarizationTime;   
 
    // Run control
    CepsBool m_stopAtCompleteActivation;     
    CepsBool m_stopAtCompleteRepolarization; 
 
};