dem_util.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/particles/particle_handler.h>
 
#include <meltpooldg/time_integration/solution_history.hpp>
 
#include <functional>
#include <vector>
 
 
namespace MeltPoolDG
{
  template <int dim>
  inline constexpr int axial_dim = [] {
    if constexpr (dim == 3)
      return 3;
    else
      return 1;
  }();
 
  template <int dim, typename number>
  dealii::Tensor<1, axial_dim<dim>, number>
  compute_torque(const dealii::Tensor<1, dim, number> &force,
                 const dealii::Tensor<1, dim, number> &lever_arm)
  {
    dealii::Tensor<1, axial_dim<dim>, number> torque;
    if constexpr (dim == 2)
      torque[0] = lever_arm[0] * force[1] - lever_arm[1] * force[0];
    if constexpr (dim == 3)
      torque = dealii::cross_product_3d(lever_arm, force);
    return torque;
  }
} // namespace MeltPoolDG