FileInterpolatorSAFunc.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/common/functions/SAFunc.hpp"
 
#include <vtkSmartPointer.h>
#include <vtkUnstructuredGridReader.h>
 
class AbstractPdeProblem;
 
class FileInterpolatorSAFunc: public SAFunc<CepsReal>
{
 
  public:
 
    explicit FileInterpolatorSAFunc(
      AbstractPdeProblem* pb,
      CepsReal            snapTimeRef,
      const CepsString&   filesBase,
      CepsUInt            tOrder = 1
    );
 
    FileInterpolatorSAFunc(const FileInterpolatorSAFunc&) = default;
 
    FileInterpolatorSAFunc&
    operator=(const FileInterpolatorSAFunc&) = default;
 
    CepsReal
    virtual eval(CepsStandardArgs args) final;
 
    CepsEnum
    getFlags() const final;
 
 
  protected:
 
    void
    computeInterpolationWeights();
 
    CepsString
    getVtkFileName(CepsUInt i) const;
 
    void
    loadFilesForTime(CepsReal t);
 
    AbstractPdeProblem* m_problem;     
    CepsReal            m_snapTimeRef; 
    CepsUInt            m_tOrder;      
    CepsString          m_filesBase;   
    CepsUInt            m_nDigitsFile; 
 
    CepsReal m_loadedTMin;      
    CepsReal m_loadedTMax;      
    CepsUInt m_loadedSnapIdMin; 
    CepsUInt m_loadedSnapIdMax; 
    CepsVector<CepsMap<CepsUnknownIndex,CepsVector<CepsReal>>> m_loadedData;
    CepsVector<CepsMap<CepsUnknownIndex,CepsReal>> m_loadedData0D;
 
    CepsMap<CepsDofGlobalIndex,CepsMap<vtkIdType,CepsReal>> m_itpWeights;
 
    CepsReal m_ptLocatorTolerance;
 
};
 
namespace ceps {
  CepsEnum
  getFlagsOf(FileInterpolatorSAFunc* f);
}