Commit bd464931 authored by Gilles MARCKMANN's avatar Gilles MARCKMANN

Remove last reference to Trellis_Util::LagrangeMapping in xFEM and change it...

Remove last reference to Trellis_Util::LagrangeMapping in xFEM and change it in xmapping::xLagrangeMapping.
parent d88736e8
......@@ -9,9 +9,9 @@
#include <cassert>
#include <iostream>
// xmapping
#include "xLagrangeMapping.h"
// trellis
#include "LagrangeMapping.h"
#include "mEdge.h"
#include "mVertex.h"
......@@ -194,7 +194,8 @@ void xScalarFunctionDerivDiscXFEMOctree<DATAMANAGER>::computeEnr(AOMD::mEntity*
if (pMesh_refined)
{
const xMesh& Mesh_refined = refined_elements.at(*e);
Trellis_Util::LagrangeMapping mapEdgeTop(qnc);
// xmapping::xLagrangeMapping mapEdgeTop(qnc);
xmapping::xLagrangeMapping mapEdgeTop(xinterface::aomd::getRefElementType(*qnc), xinterface::aomd::getPoints(*qnc));
double x = 0., y = 0., z = 0.;
mapEdgeTop.eval(chi, 0., 0., x, y, z);
// find the subedge on boundary of e
......@@ -205,7 +206,9 @@ void xScalarFunctionDerivDiscXFEMOctree<DATAMANAGER>::computeEnr(AOMD::mEntity*
double u_qnc_sub_target = 0.;
for (AOMD::mEdge* qnc_sub : qnc_subedges)
{
Trellis_Util::LagrangeMapping mapEdgeCurrent(qnc_sub);
// xmapping::xLagrangeMapping mapEdgeCurrent(qnc_sub);
xmapping::xLagrangeMapping mapEdgeCurrent(xinterface::aomd::getRefElementType(*qnc_sub),
xinterface::aomd::getPoints(*qnc_sub));
double u, v, w;
mapEdgeCurrent.invert(x, y, z, u, v, w);
if (mapEdgeCurrent.inReferenceElement(u, v, w))
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment