LegendrePolynomial.cpp

Go to the documentation of this file.

/** \file LegendrePolynomial.cpp
 * \brief Implementation of base for Legendre polynomials
 */
 
namespace MoFEM {
 
MoFEMErrorCode LegendrePolynomialCtx::query_interface(
    boost::typeindex::type_index type_index,
    UnknownInterface **iface) const {
  *iface = const_cast<LegendrePolynomialCtx *>(this);
  return 0;
}
 
MoFEMErrorCode LegendrePolynomial::query_interface(
    boost::typeindex::type_index type_index,
    UnknownInterface **iface) const {
  *iface = const_cast<LegendrePolynomial *>(this);
  return 9;
}
 
MoFEMErrorCode
LegendrePolynomial::getValue(MatrixDouble &pts,
                             boost::shared_ptr<BaseFunctionCtx> ctx_ptr) {
 
  MoFEMFunctionBeginHot;
  auto ctx = ctx_ptr->getInterface<LegendrePolynomialCtx>();
  ctx->baseFunPtr->resize(pts.size2(), ctx->P + 1, false);
  ctx->baseDiffFunPtr->resize(pts.size2(), ctx->dIm * (ctx->P + 1), false);
  double *l = NULL;
  double *diff_l = NULL;
  for (unsigned int gg = 0; gg < pts.size2(); gg++) {
    if (ctx->baseFunPtr)
      l = &((*ctx->baseFunPtr)(gg, 0));
    if (ctx->baseDiffFunPtr)
      diff_l = &((*ctx->baseDiffFunPtr)(gg, 0));
    ierr = (ctx->basePolynomialsType0)(ctx->P, pts(0, gg), ctx->diffS, l,
                                       diff_l, ctx->dIm);
    CHKERRG(ierr);
  }
  MoFEMFunctionReturnHot(0);
}
 
} // namespace MoFEM