Nektar::NekMeshUtils::CADSurfOCE Class Reference

#include <CADSurfOCE.h>

Inheritance diagram for Nektar::NekMeshUtils::CADSurfOCE:

Collaboration diagram for Nektar::NekMeshUtils::CADSurfOCE:

Public Member Functions
	CADSurfOCE ()
	~CADSurfOCE ()
void	Initialise (int i, TopoDS_Shape in)
virtual Array< OneD, NekDouble >	GetBounds ()
	Get the limits of the parametric space for the surface. More...
virtual Array< OneD, NekDouble >	N (Array< OneD, NekDouble > uv)
	Get the normal vector at parametric point u,v. More...
virtual Array< OneD, NekDouble >	D1 (Array< OneD, NekDouble > uv)
	Get the set of first derivatives at parametric point u,v. More...
virtual Array< OneD, NekDouble >	D2 (Array< OneD, NekDouble > uv)
	Get the set of second derivatives at parametric point u,v. More...
virtual Array< OneD, NekDouble >	P (Array< OneD, NekDouble > uv)
	Get the x,y,z at parametric point u,v. More...
virtual Array< OneD, NekDouble >	locuv (Array< OneD, NekDouble > p)
	Performs a reverse look up to find u,v and x,y,z. More...
virtual NekDouble	DistanceTo (Array< OneD, NekDouble > p)
	does unconstrained locuv to project point from anywhere and calculate the distance between the orthonormal projection to the surface and the point More...
virtual void	ProjectTo (Array< OneD, NekDouble > &tp, Array< OneD, NekDouble > &uv)
	takes a point from anywhere find the nearest surface point and its uv More...
virtual NekDouble	Curvature (Array< OneD, NekDouble > uv)
	returns curvature at point uv More...
Public Member Functions inherited from Nektar::NekMeshUtils::CADSurf
	CADSurf ()
	Default constructor. More...
	~CADSurf ()
std::vector < CADSystem::EdgeLoopSharedPtr >	GetEdges ()
	Get the loop structures which bound the cad surface. More...
void	SetEdges (std::vector< CADSystem::EdgeLoopSharedPtr > ein)
CADOrientation::Orientation	Orientation ()
	query reversed normal More...
Public Member Functions inherited from Nektar::NekMeshUtils::CADObject
	CADObject ()
	Default constructor. More...
virtual	~CADObject ()
int	GetId ()
	Return ID of the vertex. More...
CADType::cadType	GetType ()

Static Public Member Functions
static CADSurfSharedPtr	create ()
Static Public Member Functions inherited from Nektar::NekMeshUtils::CADSurf
static void	OrientateEdges (CADSurfSharedPtr surf, std::vector< CADSystem::EdgeLoopSharedPtr > &ein)

Static Public Attributes
static std::string	key

Private Member Functions
void	Test (Array< OneD, NekDouble > uv)
	Function which tests the the value of uv used is within the surface. More...
	Handle (Geom_Surface) m_s
	Alternate OpenCascade object for surface. Used by reverse lookup. More...

Private Attributes
BRepAdaptor_Surface	m_occSurface
	OpenCascade object for surface. More...
Array< OneD, NekDouble >	m_bounds
	parametric bounds More...
ShapeAnalysis_Surface *	m_sas
	locuv object (stored because it gets faster with stored information) More...

Additional Inherited Members
Protected Attributes inherited from Nektar::NekMeshUtils::CADSurf
std::vector < CADSystem::EdgeLoopSharedPtr >	m_edges
	List of bounding edges in loops with orientation. More...
Protected Attributes inherited from Nektar::NekMeshUtils::CADObject
int	m_id
	ID of the vert. More...
CADType::cadType	m_type
	type of the cad object More...
CADOrientation::Orientation	m_orientation
	orientation of the CADObject More...

Detailed Description

Definition at line 47 of file CADSurfOCE.h.

Constructor & Destructor Documentation

Nektar::NekMeshUtils::CADSurfOCE::CADSurfOCE ( )

inline

Definition at line 58 of file CADSurfOCE.h.

    {
    }

Nektar::NekMeshUtils::CADSurfOCE::~CADSurfOCE ( )

inline

Definition at line 62 of file CADSurfOCE.h.

    {
    }

Member Function Documentation

static CADSurfSharedPtr Nektar::NekMeshUtils::CADSurfOCE::create ( )

inlinestatic

Definition at line 51 of file CADSurfOCE.h.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

    {
        return MemoryManager<CADSurfOCE>::AllocateSharedPtr();
    }

NekDouble Nektar::NekMeshUtils::CADSurfOCE::Curvature ( Array< OneD, NekDouble >  uv )

virtual

returns curvature at point uv

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 129 of file CADSurfOCE.cpp.

{
#if defined(NEKTAR_DEBUG)
    Test(uv);
#endif

    Array<OneD, NekDouble> n = N(uv);

    // a zero normal occurs at a signularity, CurvaturePoint
    // cannot be sampled here
    if (n[0] == 0 && n[1] == 0 && n[2] == 0)
    {
        return 0.0;
    }

    Array<OneD, NekDouble> r = D2(uv);

    // metric and curvature tensors
    NekDouble E = r[3] * r[3] + r[4] * r[4] + r[5] * r[5];
    NekDouble F = r[3] * r[6] + r[4] * r[7] + r[5] * r[8];
    NekDouble G = r[6] * r[6] + r[7] * r[7] + r[8] * r[8];
    NekDouble e = n[0] * r[9] + n[1] * r[10] + n[2] * r[11];
    NekDouble f = n[0] * r[15] + n[1] * r[16] + n[2] * r[17];
    NekDouble g = n[0] * r[12] + n[1] * r[13] + n[2] * r[14];

    // if det is zero cannot invert matrix, R=0 so must skip
    if (E * G - F * F < 1E-30)
    {
        return 0.0;
    }

    NekDouble K, H;

    K = (e * g - f * f) / (E * G - F * F);
    H = 0.5 * (e * G - 2 * f * F + g * E) / (E * G - F * F);

    NekDouble kv[2];
    kv[0] = abs(H + sqrt(H * H - K));
    kv[1] = abs(H - sqrt(H * H - K));

    return kv[0] > kv[1] ? kv[0] : kv[1];
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADSurfOCE::D1 ( Array< OneD, NekDouble >  uv )

virtual

Get the set of first derivatives at parametric point u,v.

Parameters

uv	Array of u and v parametric coords.

Returns

Array of xyz copmonents of first derivatives.

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 253 of file CADSurfOCE.cpp.

{
#if defined(NEKTAR_DEBUG)
    Test(uv);
#endif

    Array<OneD, NekDouble> r(9);
    gp_Pnt Loc;
    gp_Vec D1U, D1V;
    m_occSurface.D1(uv[0], uv[1], Loc, D1U, D1V);

    r[0] = Loc.X(); // x
    r[1] = Loc.Y(); // y
    r[2] = Loc.Z(); // z
    r[3] = D1U.X(); // dx/du
    r[4] = D1U.Y(); // dy/du
    r[5] = D1U.Z(); // dz/du
    r[6] = D1V.X(); // dx/dv
    r[7] = D1V.Y(); // dy/dv
    r[8] = D1V.Z(); // dz/dv

    return r;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADSurfOCE::D2 ( Array< OneD, NekDouble >  uv )

virtual

Get the set of second derivatives at parametric point u,v.

Parameters

uv	array of u and v parametric coords

Returns

array of xyz copmonents of second derivatives

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 277 of file CADSurfOCE.cpp.

{
#if defined(NEKTAR_DEBUG)
    Test(uv);
#endif

    Array<OneD, NekDouble> r(18);
    gp_Pnt Loc;
    gp_Vec D1U, D1V, D2U, D2V, D2UV;
    m_occSurface.D2(uv[0], uv[1], Loc, D1U, D1V, D2U, D2V, D2UV);

    r[0]  = Loc.X();  // x
    r[1]  = Loc.Y();  // y
    r[2]  = Loc.Z();  // z
    r[3]  = D1U.X();  // dx/dx
    r[4]  = D1U.Y();  // dy/dy
    r[5]  = D1U.Z();  // dz/dz
    r[6]  = D1V.X();  // dx/dx
    r[7]  = D1V.Y();  // dy/dy
    r[8]  = D1V.Z();  // dz/dz
    r[9]  = D2U.X();  // d2x/du2
    r[10] = D2U.Y();  // d2y/du2
    r[11] = D2U.Z();  // d2z/du2
    r[12] = D2V.X();  // d2x/dv2
    r[13] = D2V.Y();  // d2y/dv2
    r[14] = D2V.Z();  // d2z/dv2
    r[15] = D2UV.X(); // d2x/dudv
    r[16] = D2UV.Y(); // d2y/dudv
    r[17] = D2UV.Z(); // d2z/dudv

    return r;
}

NekDouble Nektar::NekMeshUtils::CADSurfOCE::DistanceTo ( Array< OneD, NekDouble >  p )

virtual

does unconstrained locuv to project point from anywhere and calculate the distance between the orthonormal projection to the surface and the point

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 172 of file CADSurfOCE.cpp.

{
    gp_Pnt loc(p[0] * 1000.0, p[1] * 1000.0, p[2] * 1000.0);

    // alternative locuv methods
    ShapeAnalysis_Surface sas(m_s);
    sas.SetDomain(m_bounds[0], m_bounds[1], m_bounds[2], m_bounds[3]);

    gp_Pnt2d p2 = sas.ValueOfUV(loc, 1e-7);

    gp_Pnt p3 = sas.Value(p2);

    return p3.Distance(loc);
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADSurfOCE::GetBounds ( )

virtual

Get the limits of the parametric space for the surface.

Returns

Array of 4 entries with parametric umin,umax,vmin,vmax.

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 76 of file CADSurfOCE.cpp.

{
    return m_bounds;
}

Nektar::NekMeshUtils::CADSurfOCE::Handle ( Geom_Surface  )

private

Alternate OpenCascade object for surface. Used by reverse lookup.

void Nektar::NekMeshUtils::CADSurfOCE::Initialise	(	int	i,
		TopoDS_Shape	in
	)

Definition at line 48 of file CADSurfOCE.cpp.

References Nektar::StdRegions::eBackwards.

{
    // this bit of code changes the units of the cad from mm opencascade
    // defualt to m

    m_s = BRep_Tool::Surface(TopoDS::Face(in));

    if (in.Orientation() == 1)
    {
        m_orientation = CADOrientation::eBackwards;
    }

    gp_Trsf transform;
    gp_Pnt ori(0.0, 0.0, 0.0);
    transform.SetScale(ori, 1.0 / 1000.0);
    TopLoc_Location mv(transform);

    in.Move(mv);
    m_occSurface = BRepAdaptor_Surface(TopoDS::Face(in));
    m_id         = i;

    m_bounds = Array<OneD, NekDouble>(4);
    BRepTools::UVBounds(TopoDS::Face(in), m_bounds[0], m_bounds[1], m_bounds[2],
                        m_bounds[3]);
    m_sas = new ShapeAnalysis_Surface(m_s);
    m_sas->SetDomain(m_bounds[0], m_bounds[1], m_bounds[2], m_bounds[3]);
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADSurfOCE::locuv ( Array< OneD, NekDouble >  p )

virtual

Performs a reverse look up to find u,v and x,y,z.

Parameters

p	Array of xyz location

Returns

The parametric location of xyz on this surface

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 81 of file CADSurfOCE.cpp.

References ASSERTL0.

{
    // has to transfer back to mm
    gp_Pnt loc(p[0] * 1000.0, p[1] * 1000.0, p[2] * 1000.0);

    Array<OneD, NekDouble> uvr(2);

    gp_Pnt2d p2 = m_sas->ValueOfUV(loc, 1e-3);
    uvr[0]      = p2.X();
    uvr[1]      = p2.Y();

    gp_Pnt p3 = m_sas->Value(p2);
    if (p3.Distance(loc) > 1.0)
    {
        cout << "large locuv distance " << p3.Distance(loc) << " " << m_id
             << endl;
    }

    // if the uv returned is slightly off the surface
    //(which ShapeAnalysis_Surface can do sometimes)
    if (uvr[0] < m_bounds[0] || uvr[0] > m_bounds[1] || uvr[1] < m_bounds[2] ||
        uvr[1] > m_bounds[3])
    {
        if (uvr[0] < m_bounds[0])
        {
            uvr[0] = m_bounds[0];
        }
        else if (uvr[0] > m_bounds[1])
        {
            uvr[0] = m_bounds[1];
        }
        else if (uvr[1] < m_bounds[2])
        {
            uvr[1] = m_bounds[2];
        }
        else if (uvr[1] > m_bounds[3])
        {
            uvr[1] = m_bounds[3];
        }
        else
        {
            ASSERTL0(false, "Cannot correct locuv");
        }
    }

    return uvr;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADSurfOCE::N ( Array< OneD, NekDouble >  uv )

virtual

Get the normal vector at parametric point u,v.

Parameters

uv	Array of u and v parametric coords.

Returns

Array of xyz components of normal vector.

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 223 of file CADSurfOCE.cpp.

References Nektar::StdRegions::eBackwards.

{
#if defined(NEKTAR_DEBUG)
    Test(uv);
#endif

    BRepLProp_SLProps slp(m_occSurface, 2, 1e-6);
    slp.SetParameters(uv[0], uv[1]);

    if (!slp.IsNormalDefined())
    {
        return Array<OneD, NekDouble>(3, 0.0);
    }

    gp_Dir d = slp.Normal();

    Array<OneD, NekDouble> normal(3);

    if (m_orientation == CADOrientation::eBackwards)
    {
        d.Reverse();
    }

    normal[0] = d.X();
    normal[1] = d.Y();
    normal[2] = d.Z();

    return normal;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADSurfOCE::P ( Array< OneD, NekDouble >  uv )

virtual

Get the x,y,z at parametric point u,v.

Parameters

uv	Array of u and v parametric coords.

Returns

Array of xyz location.

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 208 of file CADSurfOCE.cpp.

{
#if defined(NEKTAR_DEBUG)
    Test(uv);
#endif

    Array<OneD, NekDouble> location(3);
    gp_Pnt loc;
    loc         = m_occSurface.Value(uv[0], uv[1]);
    location[0] = loc.X();
    location[1] = loc.Y();
    location[2] = loc.Z();
    return location;
}

void Nektar::NekMeshUtils::CADSurfOCE::ProjectTo ( Array< OneD, NekDouble > &  tp, Array< OneD, NekDouble > &  uv  )

virtual

takes a point from anywhere find the nearest surface point and its uv

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 187 of file CADSurfOCE.cpp.

{
    gp_Pnt loc(tp[0] * 1000.0, tp[1] * 1000.0, tp[2] * 1000.0);

    // alternative locuv methods
    ShapeAnalysis_Surface sas(m_s);
    sas.SetDomain(m_bounds[0], m_bounds[1], m_bounds[2], m_bounds[3]);

    gp_Pnt2d p2 = sas.ValueOfUV(loc, 1e-7);

    gp_Pnt p3 = sas.Value(p2);

    tp[0] = p3.X() / 1000.0;
    tp[1] = p3.Y() / 1000.0;
    tp[2] = p3.Z() / 1000.0;

    uv[0] = p2.X();
    uv[1] = p2.Y();
}

void Nektar::NekMeshUtils::CADSurfOCE::Test ( Array< OneD, NekDouble >  uv )

privatevirtual

Function which tests the the value of uv used is within the surface.

Implements Nektar::NekMeshUtils::CADSurf.

Definition at line 310 of file CADSurfOCE.cpp.

References ASSERTL1.

{
    stringstream error;

    error << "Point not within parameter plane: ";

    bool passed = true;

    if (uv[0] < m_bounds[0])
    {
        if (fabs(uv[0] - m_bounds[0]) > 1E-6)
        {
            error << "U(" << uv[0] << ") is less than Umin(" << m_bounds[0]
                  << ")";
            passed = false;
        }
    }
    else if (uv[0] > m_bounds[1])
    {
        if (fabs(uv[0] - m_bounds[1]) > 1E-6)
        {
            error << "U(" << uv[0] << ") is greater than Umax(" << m_bounds[1]
                  << ")";
            passed = false;
        }
    }
    else if (uv[1] < m_bounds[2])
    {
        if (fabs(uv[1] - m_bounds[2]) > 1E-6)
        {
            error << "V(" << uv[1] << ") is less than Vmin(" << m_bounds[2]
                  << ")";
            passed = false;
        }
    }
    else if (uv[1] > m_bounds[3])
    {
        if (fabs(uv[1] - m_bounds[3]) > 1E-6)
        {
            error << "V(" << uv[1] << ") is greater than Vmax(" << m_bounds[3]
                  << ")";
            passed = false;
        }
    }

    error << " On Surface: " << GetId();
    ASSERTL1(passed, "Warning: " + error.str());
}

Member Data Documentation

std::string Nektar::NekMeshUtils::CADSurfOCE::key

static

Initial value:

= GetCADSurfFactory().RegisterCreatorFunction(
    "oce", CADSurfOCE::create, "CADSurfOCE")

Definition at line 56 of file CADSurfOCE.h.

Array<OneD, NekDouble> Nektar::NekMeshUtils::CADSurfOCE::m_bounds

private

parametric bounds

Definition at line 87 of file CADSurfOCE.h.

BRepAdaptor_Surface Nektar::NekMeshUtils::CADSurfOCE::m_occSurface

private

OpenCascade object for surface.

Definition at line 83 of file CADSurfOCE.h.

ShapeAnalysis_Surface* Nektar::NekMeshUtils::CADSurfOCE::m_sas

private

locuv object (stored because it gets faster with stored information)

Definition at line 89 of file CADSurfOCE.h.