HeatSolver.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 "pde/solver/timed/HeatSolver.hpp"
 
// #include "ode/solver/RungeKutta.hpp"
 
#include "pde/problem/timed/HeatProblem.hpp"
#include "pde/assemblers/FE/FEHeatAssembler.hpp"
#include "pde/assemblers/FE/FEDivKGradBCAssembler.hpp"
 
HeatSolver::HeatSolver(HeatProblem *problem) :
  FETimedSolver(problem)
{
//  m_scaleSystem = 1.0 / problem->getTimeStep();
}
 
HeatSolver::~HeatSolver()
{
}
 
void
HeatSolver::initializeAssemblers()
{
  auto hpb = this->getHeatProblem();
 
 
  FEHeatAssembler* opAsb = ceps::getNew<FEHeatAssembler>(m_fe);
  for (Unknown* u: m_problem->getUnknowns())
    if (u->isSpatial())
      opAsb->setKForUnknown(u,hpb->getDiffusionTensor(u));
  m_opAsb = opAsb; 
 
  m_bcAsb  = ceps::getNew<FEDivKGradBCAssembler>(m_fe);  
  m_bcAsb->setMatrix(m_lhs.get());
}
 
HeatProblem*
HeatSolver::getHeatProblem() const
{
  return ceps::runtimeCast<HeatProblem*>(m_problem);
}