FECardiacBCAssembler.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 "cardiac/assembler/FE/FECardiacBCAssembler.hpp"
#include "cardiac/problem/CardiacProblem.hpp"
 
FECardiacBCAssembler::FECardiacBCAssembler(
  CardiacProblem*               pb,
  FiniteElements*               fe, 
  const CepsVector<Unknown*>&   unknowns
) : FECardiacBCAssembler(pb, pb->getTissueAttributes(), fe, unknowns)
{
}
 
FECardiacBCAssembler::FECardiacBCAssembler(
  CardiacProblem*               pb,
  const CepsSet<CepsAttribute>  &attributes,
  FiniteElements*               fe, 
  const CepsVector<Unknown*>&   unknowns
) :
    FEDivKGradBCAssembler(fe,attributes,unknowns),
    FECardiacParameters(pb)
{
}
 
void
FECardiacBCAssembler::computeBlocksOnElementAtQuadPoint(
  FEBase*                  elem,
  CepsReal3D               xQ,
  CepsReal                 t,
  const CepsMathDynamic1D& phi,
  const CepsMathDynamic2D& gradPhi
)
{
  // a * sigma * grad_n u  + b * u = g
  // => sigma * grad_n u = g / a - (b / a) * u
 
  // Fetch the unknown for which we assemble
  CepsUnknownIndex uid = m_bc->getUnknownId();
  Unknown*         u   = m_problem->getUnknown(uid);
 
  // Check that we can assemble on that element for this unknown.
  // If not, the vectors following structures are empty
  // (we still need to call setupBlocks though !)
  auto us = extractUnknownsFrom({u},elem);
  setupBlocksOnElementForUnknown(elem,u);
  if (us.empty())
    return;
 
  // Evaluate the BC coefficients at the quadrature point
  CepsStandardArgs args;
  args.t         = t;
  args.x         = xQ;
  args.cellId    = elem->getGeomObject()->getGlobalIndex();
  args.attr      = elem->getGeomObject()->getAttributes();
  args.unknownId = uid;
  CepsReal alpha = m_bc->getAlpha(args);
  CepsReal beta  = m_bc->getBeta(args);
  CepsReal gamma = m_bc->eval(args);
 
  // Get Am and Cm at position xQ
  CepsReal Am = m_Am->eval(args);
  CepsReal Cm = m_Cm->eval(args);
 
  // Now computing the submatrix/subvector 
  m_bMat[0][0] += 1. / (Am * Cm) * beta  / alpha * phi * phi.transpose();
  m_bVec[0]    += 1. / (Am * Cm) * gamma / alpha * phi;
 
  return;
}