FEDivKGradAssembler.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/assemblers/FE/FEAssembler.hpp"
#include "pde/common/functions/SAFunc.hpp"
 
class FEDivKGradAssembler : public FEAssembler
{
 
  public:
 
    FEDivKGradAssembler();
 
    explicit FEDivKGradAssembler(
      FiniteElements*               fe, 
      CepsBool                      isBoundary = false, 
      const CepsSet<CepsAttribute>& attributes = {CepsUniversal},
      const CepsVector<Unknown*>&   unknowns   = {}
    );
 
    virtual
    ~FEDivKGradAssembler();
 
    void
    setKForUnknown(Unknown* u, CepsMathScalar k);
  
    void
    setKForUnknown(Unknown* u, SAFunc<CepsMathScalar>* k);
 
    void
    setKForUnknown(Unknown* u, CepsMathTensor k);
  
    void
    setKForUnknown(Unknown* u, SAFunc<CepsMathTensor>* k);
 
    CepsBool
    isChangingBetweenTimes(CepsReal t1, CepsReal t2) const override;
 
  protected:
 
    CepsVector<CepsMathTensor>
    getDiffusionTensorsForUnknowns(
      const CepsVector<Unknown*>& us, 
      FEBase*                     elem,
      CepsReal3D                  x, 
      CepsReal                    t
    ) const;
  
    void
    computeBlocksOnElementAtQuadPoint(
      FEBase*                  element,
      CepsReal3D               xQ,
      CepsReal                 t,
      const CepsMathDynamic1D& phi,
      const CepsMathDynamic2D& gradPhi
    ) override;
 
 
 
  protected:
    CepsMap<CepsUnknownIndex, CepsMathScalar>          m_cstScalars;  
    CepsMap<CepsUnknownIndex, SAFunc<CepsMathScalar>*> m_fctScalars;  
    CepsMap<CepsUnknownIndex, CepsMathTensor>          m_cstTensors;  
    CepsMap<CepsUnknownIndex, SAFunc<CepsMathTensor>*> m_fctTensors;  
    CepsBool m_isTimeDependant;                                       
};