ActivationTracker.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 "common/CepsCommon.hpp"
 
#include "pde/common/writers/VtkWriter.hpp"
//  #include "pde/common/writers/MeditSolutionWriter.hpp"
#include "pde/common/writers/TimeWriter.hpp"
 
#include "linearAlgebra/DistributedHaloVector.hpp"
#include "geometry/Geometry.hpp"
 
class CardiacSolver;
class CardiacProblem;
 
class ActivationTracker : public CepsObject
{
  public:
 
    ActivationTracker(CardiacSolver* solver);
 
    ~ActivationTracker();
 
    void update(
      CepsInt iter, 
      CepsReal time, 
      DHVecPtr solution, 
      DHVecPtr prevSolution
    );
 
    void writeActivationMap();
 
    CepsBool allPointsHaveBeenActivated() const;
 
    CepsBool allPointsHaveBeenRepolarized() const;
 
    CepsString getByeByeMessage() const; 
 
    CepsVector<CepsReal>&
    getActivatedVolume();
 
 
  protected:
 
    void
    setupWithParameters(InputParameters* params);
 
    void
    initializeActivationMap();
 
    CepsBool
    allSeen(const CepsMap<CepsUnknownIndex,CepsSet<CepsDofGlobalIndex>>&) const;
 
 
  protected:
 
    CardiacProblem* m_problem;    
    Geometry*       m_geom;       
    TimeWriter*     m_timeWriter; 
 
    CepsUInt m_nbIterPostProcess; 
 
    CepsMap<CepsUnknownIndex,CepsReal3D> m_activationTimeData; 
    static constexpr CepsUInt _threshold = 0; 
    static constexpr CepsUInt _apdpctg   = 1; 
    static constexpr CepsUInt _minV4Peak = 2; 
 
    DHVecPtr m_activationTimes; 
    DHVecPtr m_activated;       
    DHVecPtr m_peakValues;      
    DHVecPtr m_peakTimes;       
    DHVecPtr m_APD50;           
    DHVecPtr m_APD;             
 
    CepsMap<CepsUnknownIndex,CepsSet<CepsDofGlobalIndex>> m_activationNotYetSeen; 
    CepsMap<CepsUnknownIndex,CepsSet<CepsDofGlobalIndex>> m_activationSeen;       
    CepsMap<CepsUnknownIndex,CepsSet<CepsDofGlobalIndex>> m_foundPeak;            
    CepsMap<CepsUnknownIndex,CepsSet<CepsDofGlobalIndex>> m_foundAPD50;           
    CepsMap<CepsUnknownIndex,CepsSet<CepsDofGlobalIndex>> m_foundAPD;             
 
    CepsVector<CepsReal> m_activatedVolume; 
    CepsReal             m_tissueVolume;    
 
    CepsUInt m_nToBeSeen; 
 
};