lbs_01a_inputchecks.cc

Go to the documentation of this file.

#include "lbs_solver.h"
 
#include "mesh/MeshHandler/chi_meshhandler.h"
#include "mesh/MeshContinuum/chi_meshcontinuum.h"
 
#include "chi_runtime.h"
#include "chi_log.h"
 
//###################################################################
/**Performs general input checks before initialization continues.*/
void lbs::LBSSolver::PerformInputChecks()
{
  if (groups_.empty())
  {
    Chi::log.LogAllError()
      << "LinearBoltzmann::SteadyStateSolver: No groups added to solver.";
    Chi::Exit(EXIT_FAILURE);
  }
 
  num_groups_ = groups_.size();
 
  if (groupsets_.empty())
  {
    Chi::log.LogAllError()
      << "LinearBoltzmann::SteadyStateSolver: No group-sets added to solver.";
    Chi::Exit(EXIT_FAILURE);
  }
  int grpset_counter=0;
  for (auto& group_set : groupsets_)
  {
    if (group_set.groups_.empty())
    {
      Chi::log.LogAllError()
        << "LinearBoltzmann::SteadyStateSolver: No groups added to groupset "
        << grpset_counter << ".";
      Chi::Exit(EXIT_FAILURE);
    }
    ++grpset_counter;
  }
  if (options_.sd_type == chi_math::SpatialDiscretizationType::UNDEFINED)
  {
    Chi::log.LogAllError()
      << "LinearBoltzmann::SteadyStateSolver: No discretization_ method set.";
    Chi::Exit(EXIT_FAILURE);
  }
 
  grid_ptr_ = chi_mesh::GetCurrentHandler().GetGrid();
 
  if (grid_ptr_ == nullptr)
  {
    Chi::log.LogAllError()
      << "LinearBoltzmann::SteadyStateSolver: No grid_ptr_ available from region.";
    Chi::Exit(EXIT_FAILURE);
  }
 
  //======================================== Determine geometry type
  using namespace chi_mesh;
  const auto grid_attribs = grid_ptr_->Attributes();
  if (grid_attribs & DIMENSION_1)
    options_.geometry_type = GeometryType::ONED_SLAB;
  else if (grid_attribs & DIMENSION_2)
    options_.geometry_type = GeometryType::TWOD_CARTESIAN;
  else if (grid_attribs & DIMENSION_3)
    options_.geometry_type = GeometryType::THREED_CARTESIAN;
  else
    ChiLogicalError("Cannot deduce geometry type from mesh.");
 
 
}