node.cc

Go to the documentation of this file.

#include "mesh/node.h"
 
#include <vector>
 
#include "banks/source_bank.h"
 
namespace hummingbird {
double DistanceBetweenNodes(const Node& node_1, const Node& node_2) {
  double x_part = node_1.x - node_2.x;
  double y_part = node_1.y - node_2.y;
  double z_part = node_1.z - node_2.z;
  return std::sqrt(x_part * x_part + y_part * y_part + z_part * z_part);
}
 
void UpdateScalarFlux(Node& node,
                      const QuadratureBase<Ordinate> angular_quadrature) {
  std::vector<QuadraturePair> quad_pairs;
  quad_pairs.reserve(angular_quadrature.n_points());
  for (auto i = 0; i < angular_quadrature.n_points(); i++)
    quad_pairs.push_back(QuadraturePair(i, node.angular_fluxes.at(i)));
  node.scalar_flux = angular_quadrature.Integrate(quad_pairs);
}
 
std::vector<double> ComputeNodeSourceFluxes(
    const Node& node, const int material_id, const int source_id,
    const MaterialBank& material_bank, const SourceBank& source_bank,
    const QuadratureBase<Ordinate>& angular_quadrature) {
  std::vector<double> source_fluxes;
  source_fluxes.reserve(angular_quadrature.n_points());
  for (auto n = 0; n < angular_quadrature.n_points(); n++) {
    double ind_source = source_bank.GetByID(source_id)->EvaluateAtNode(
        node, angular_quadrature.GetAbscissa(n));
    double inscattering = material_bank.GetByID(material_id).scattering_xs /
                          (4.0 * M_PI) * node.scalar_flux;
    source_fluxes.push_back(inscattering + ind_source);
  }
  return source_fluxes;
}
}  // namespace hummingbird