pwld_00_constrdestr.cc

Go to the documentation of this file.

#include "PieceWiseLinearDiscontinuous.h"
 
#include "math/UnknownManager/unknown_manager.h"
 
#include "chi_runtime.h"
#include "chi_log.h"
 
#include "utils/chi_timer.h"
 
namespace chi_math::spatial_discretization
{
 
// ###################################################################
/**Constructor.*/
PieceWiseLinearDiscontinuous::PieceWiseLinearDiscontinuous(
  const chi_mesh::MeshContinuum& grid,
  chi_math::QuadratureOrder q_order,
  chi_math::CoordinateSystemType cs_type)
  : PieceWiseLinearBase(
      grid, q_order, SDMType::PIECEWISE_LINEAR_DISCONTINUOUS, cs_type)
{
  CreateCellMappings();
 
  OrderNodes();
}
 
// ###################################################################
/**Construct a shared object using the protected constructor.*/
std::shared_ptr<PieceWiseLinearDiscontinuous> PieceWiseLinearDiscontinuous::New(
  const chi_mesh::MeshContinuum& grid,
  QuadratureOrder q_order /*=QuadratureOrder::SECOND*/,
  CoordinateSystemType cs_type /*=CoordinateSystemType::CARTESIAN*/)
 
{
  const auto PWLD = SpatialDiscretizationType::PIECEWISE_LINEAR_DISCONTINUOUS;
  // First try to find an existing spatial discretization that matches the
  // one requested.
  for (auto& sdm : Chi::sdm_stack)
    if (sdm->Type() == PWLD and
        std::addressof(sdm->Grid()) == std::addressof(grid) and
        sdm->GetCoordinateSystemType() == cs_type)
    {
      auto fe_ptr = std::dynamic_pointer_cast<FiniteElementBase>(sdm);
 
      ChiLogicalErrorIf(not fe_ptr, "Casting failure to FE");
 
      if (fe_ptr->GetQuadratureOrder() != q_order) break;
 
      auto sdm_ptr =
        std::dynamic_pointer_cast<PieceWiseLinearDiscontinuous>(fe_ptr);
 
      ChiLogicalErrorIf(not sdm_ptr, "Casting failure");
 
      return sdm_ptr;
    }
 
  auto new_sdm = std::shared_ptr<PieceWiseLinearDiscontinuous>(
    new PieceWiseLinearDiscontinuous(grid, q_order, cs_type));
 
  Chi::sdm_stack.push_back(new_sdm);
 
  return new_sdm;
}
 
} // namespace chi_math::spatial_discretization