geometry_interface.cc

/* @file geometry_interface.cc
 *  @brief implementation of functions operating on generic geometry objects
 */
 
#include "geometry_interface.h"
 
#include "print_info.h"
 
namespace lf::geometry {
 
double Volume(const Geometry& geo) {
  const lf::base::dim_t refdim = geo.DimLocal();
 
  Eigen::Matrix<double, Eigen::Dynamic, 1> refc(refdim, 1);
  double refvol = 1.0;
  switch (geo.RefEl()) {
    case lf::base::RefEl::kPoint(): {
      return 1.0;
    }
    case lf::base::RefEl::kSegment(): {
      refc << 0.5;
      break;
    }
    case lf::base::RefEl::kTria(): {
      refvol = 0.5;
      refc << 1.0 / 3, 1.0 / 3;
      break;
    }
    case lf::base::RefEl::kQuad(): {
      refc << 0.5, 0.5;
      break;
    }
    default: {
      LF_VERIFY_MSG(false,
                    "Volume not available for " << geo.RefEl().ToString());
    }
  }  // end switch RefEl
  return (refvol * ((geo.IntegrationElement(refc))[0]));
}
 
}  // namespace lf::geometry