Nektar::NekMeshUtils::CADCurveOCE Class Reference

#include <CADCurveOCE.h>

Inheritance diagram for Nektar::NekMeshUtils::CADCurveOCE:

Collaboration diagram for Nektar::NekMeshUtils::CADCurveOCE:

Public Member Functions
	CADCurveOCE ()
	~CADCurveOCE ()
virtual Array< OneD, NekDouble >	GetBounds ()
	Returns the minimum and maximum parametric coords t of the curve. More...
virtual NekDouble	Length (NekDouble ti, NekDouble tf)
	Calculates the arclength between the two paremetric points ti and tf. ti must be less than tf. More...
virtual Array< OneD, NekDouble >	P (NekDouble t)
	Gets the location (x,y,z) in an array out of the curve at point t. More...
virtual Array< OneD, NekDouble >	D2 (NekDouble t)
	Gets the second derivatives at t. More...
virtual NekDouble	tAtArcLength (NekDouble s)
	Calculates the parametric coordinate and arclength location defined by s. More...
virtual Array< OneD, NekDouble >	GetMinMax ()
	Gets the start and end of the curve. More...
virtual NekDouble	loct (Array< OneD, NekDouble > xyz)
virtual NekDouble	Curvature (NekDouble t)
virtual Array< OneD, NekDouble >	NormalWRT (NekDouble t, int surf)
virtual Array< OneD, NekDouble >	N (NekDouble t)
void	Initialise (int i, TopoDS_Shape in)
Public Member Functions inherited from Nektar::NekMeshUtils::CADCurve
	CADCurve ()
	Default constructor. More...
	~CADCurve ()
void	SetAdjSurf (std::pair< CADSurfSharedPtr, CADOrientation::Orientation > i)
	set the ids of the surfaces either side of the curve More...
std::vector< std::pair < CADSurfSharedPtr, CADOrientation::Orientation > >	GetAdjSurf ()
NekDouble	GetTotLength ()
void	SetVert (std::vector< CADVertSharedPtr > &falVert)
std::vector< CADVertSharedPtr >	GetVertex ()
CADOrientation::Orientation	GetOrienationWRT (int surf)
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 ()
virtual CADOrientation::Orientation	Orientation ()

Static Public Member Functions
static CADCurveSharedPtr	create ()

Static Public Attributes
static std::string	key

Private Member Functions
	Handle (Geom_Curve) m_c
	Alternate object used for reverse lookups. More...

Private Attributes
BRepAdaptor_Curve	m_occCurve
	OpenCascade object of the curve. More...
TopoDS_Edge	m_occEdge
	OpenCascade edge. More...

Additional Inherited Members
Protected Attributes inherited from Nektar::NekMeshUtils::CADCurve
NekDouble	m_length
	Length of edge. More...
std::vector< std::pair < CADSurfSharedPtr, CADOrientation::Orientation > >	m_adjSurfs
	List of surfaces which this curve belongs to. More...
std::vector< CADVertSharedPtr >	m_mainVerts
	list of end vertices 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 CADCurveOCE.h.

Constructor & Destructor Documentation

Nektar::NekMeshUtils::CADCurveOCE::CADCurveOCE ( )

inline

Definition at line 58 of file CADCurveOCE.h.

    {
    }

Nektar::NekMeshUtils::CADCurveOCE::~CADCurveOCE ( )

inline

Definition at line 62 of file CADCurveOCE.h.

    {
    }

Member Function Documentation

static CADCurveSharedPtr Nektar::NekMeshUtils::CADCurveOCE::create ( )

inlinestatic

Definition at line 51 of file CADCurveOCE.h.

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

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

NekDouble Nektar::NekMeshUtils::CADCurveOCE::Curvature ( NekDouble  t )

virtual

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 191 of file CADCurveOCE.cpp.

{
    GeomLProp_CLProps d(m_c,2,1e-8);
    d.SetParameter(t);

    return d.Curvature() * 1000.0;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADCurveOCE::D2 ( NekDouble  t )

virtual

Gets the second derivatives at t.

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 111 of file CADCurveOCE.cpp.

{
    Array<OneD, NekDouble> out(9);
    gp_Pnt loc;
    gp_Vec d1, d2;
    m_occCurve.D2(t, loc, d1, d2);

    out[0] = loc.X();
    out[1] = loc.Y();
    out[2] = loc.Z();
    out[3] = d1.X();
    out[4] = d1.Y();
    out[5] = d1.Z();
    out[6] = d2.X();
    out[7] = d2.Y();
    out[8] = d2.Z();

    return out;
}

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

virtual

Returns the minimum and maximum parametric coords t of the curve.

Returns

Array of two entries, min and max parametric coordinate.

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 199 of file CADCurveOCE.cpp.

Referenced by Initialise().

{
    Array<OneD, NekDouble> t(2);
    t[0] = m_occCurve.FirstParameter();
    t[1] = m_occCurve.LastParameter();

    return t;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADCurveOCE::GetMinMax ( )

virtual

Gets the start and end of the curve.

Returns

Array with 6 entries of endpoints x1,y1,z1,x2,y2,z2.

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 208 of file CADCurveOCE.cpp.

{
    Array<OneD, NekDouble> locs(6);

    gp_Pnt start =
        BRep_Tool::Pnt(TopExp::FirstVertex(m_occEdge, Standard_True));
    gp_Pnt end = BRep_Tool::Pnt(TopExp::LastVertex(m_occEdge, Standard_True));

    locs[0] = start.X();
    locs[1] = start.Y();
    locs[2] = start.Z();
    locs[3] = end.X();
    locs[4] = end.Y();
    locs[5] = end.Z();

    return locs;
}

Nektar::NekMeshUtils::CADCurveOCE::Handle ( Geom_Curve  )

private

Alternate object used for reverse lookups.

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

inline

Definition at line 77 of file CADCurveOCE.h.

References GetBounds(), Nektar::NekMeshUtils::CADObject::m_id, Nektar::NekMeshUtils::CADCurve::m_length, m_occCurve, and m_occEdge.

    {
        gp_Trsf transform;
        gp_Pnt ori(0.0, 0.0, 0.0);
        transform.SetScale(ori, 1.0 / 1000.0);
        TopLoc_Location mv(transform);
        TopoDS_Shape cp = in;
        in.Move(mv);

        m_occEdge  = TopoDS::Edge(in);
        m_occCurve = BRepAdaptor_Curve(m_occEdge);

        GProp_GProps System;
        BRepGProp::LinearProperties(m_occEdge, System);
        m_length = System.Mass();

        Array<OneD, NekDouble> b = GetBounds();
        m_c = BRep_Tool::Curve(TopoDS::Edge(cp), b[0], b[1]);

        m_id   = i;
    }

NekDouble Nektar::NekMeshUtils::CADCurveOCE::Length ( NekDouble  ti, NekDouble  tf  )

virtual

Calculates the arclength between the two paremetric points ti and tf. ti must be less than tf.

Parameters

ti	First parametric coordinate.
tf	Second parametric coordinate.

Returns

Arc length between ti and tf.

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 71 of file CADCurveOCE.cpp.

{
    Array<OneD, NekDouble> b = GetBounds();
    Handle(Geom_Curve) NewCurve = new Geom_TrimmedCurve(m_c, ti, tf);
    TopoDS_Edge NewEdge = BRepBuilderAPI_MakeEdge(NewCurve);
    GProp_GProps System;
    BRepGProp::LinearProperties(NewEdge, System);
    return System.Mass() / 1000.0;
}

NekDouble Nektar::NekMeshUtils::CADCurveOCE::loct ( Array< OneD, NekDouble >  xyz )

virtual

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 81 of file CADCurveOCE.cpp.

References CellMLToNektar.cellml_metadata::p.

{
    NekDouble t = 0.0;
    Array<OneD, NekDouble> b = GetBounds();

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

    ShapeAnalysis_Curve sac;
    gp_Pnt p;
    sac.Project(m_c,loc,1e-8 ,p,t);

    if(p.Distance(loc) > 1e-5)
    {
        cerr << "large loct distance" << endl;
    }
    return t;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADCurveOCE::N ( NekDouble  t )

virtual

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 168 of file CADCurveOCE.cpp.

{
    GeomLProp_CLProps d(m_c,2,1e-8);
    d.SetParameter(t+1e-8);

    gp_Vec d2 = d.D2();
    if(d2.Magnitude() < 1e-8)
    {
        //no normal, stright line
        return Array<OneD, NekDouble>(3,0.0);
    }

    gp_Dir n;
    d.Normal(n);

    Array<OneD, NekDouble> N(3);
    N[0] = n.X();
    N[1] = n.Y();
    N[2] = n.Z();

    return N;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADCurveOCE::NormalWRT ( NekDouble  t, int  surf  )

virtual

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 131 of file CADCurveOCE.cpp.

References ASSERTL0, Nektar::StdRegions::eBackwards, class_topology::P, and CellMLToNektar.cellml_metadata::p.

{
    Array<OneD, NekDouble> p = P(t);
    pair<CADSurfSharedPtr, CADOrientation::Orientation> surface;
    ASSERTL0(m_adjSurfs.size() == 1, "This will only work in 2D for one surface at the moment");
    surface = m_adjSurfs[0];

    Array<OneD, NekDouble> uv = surface.first->locuv(p);
    Array<OneD, NekDouble> d1 = surface.first->D1(uv);

    NekDouble t1 = t - 1e-8;
    NekDouble t2 = t + 1e-8;

    if(surface.second == CADOrientation::eBackwards)
    {
        swap(t1, t2);
    }

    Array<OneD, NekDouble> uv1 = surface.first->locuv(P(t1));
    Array<OneD, NekDouble> uv2 = surface.first->locuv(P(t2));

    NekDouble du = uv2[1] - uv1[1];
    NekDouble dv = -1.0*(uv2[0] - uv1[0]);

    Array<OneD, NekDouble> N(3,0.0);
    N[0] = (d1[3] * du + d1[6] * dv) / 2.0;
    N[1] = (d1[4] * du + d1[7] * dv) / 2.0;
    N[2] = (d1[5] * du + d1[8] * dv) / 2.0;

    NekDouble mag = sqrt(N[0]*N[0] + N[1]*N[1] + N[2]*N[2]);
    N[0] /= mag;
    N[1] /= mag;
    N[2] /= mag;

    return N;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADCurveOCE::P ( NekDouble  t )

virtual

Gets the location (x,y,z) in an array out of the curve at point t.

Parameters

t	Parametric coordinate

Returns

Array of x,y,z

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 99 of file CADCurveOCE.cpp.

{
    Array<OneD, NekDouble> location(3);
    gp_Pnt loc = m_occCurve.Value(t);

    location[0] = loc.X();
    location[1] = loc.Y();
    location[2] = loc.Z();

    return location;
}

NekDouble Nektar::NekMeshUtils::CADCurveOCE::tAtArcLength ( NekDouble  s )

virtual

Calculates the parametric coordinate and arclength location defined by s.

Parameters

s	Arclength location.

Returns

Calculated parametric coordinate.

Todo:

This really needs improving for accuracy.

Implements Nektar::NekMeshUtils::CADCurve.

Definition at line 48 of file CADCurveOCE.cpp.

{
    NekDouble dt =
        (m_occCurve.LastParameter() - m_occCurve.FirstParameter()) / (1000);
    NekDouble t = m_occCurve.FirstParameter();

    NekDouble len = 0.0;

    while (len <= s)
    {
        gp_Pnt P1, P2;
        gp_Vec drdt1, drdt2;

        m_occCurve.D1(t, P1, drdt1);
        t += dt;
        m_occCurve.D1(t, P2, drdt2);

        len += (drdt1.Magnitude() + drdt2.Magnitude()) / 2.0 * dt;
    }

    return t - dt;
}

Member Data Documentation

std::string Nektar::NekMeshUtils::CADCurveOCE::key

static

Initial value:

= GetCADCurveFactory().RegisterCreatorFunction(
    "oce", CADCurveOCE::create, "CADCurveOCE")

Definition at line 56 of file CADCurveOCE.h.

BRepAdaptor_Curve Nektar::NekMeshUtils::CADCurveOCE::m_occCurve

private

OpenCascade object of the curve.

Definition at line 101 of file CADCurveOCE.h.

Referenced by Initialise().

TopoDS_Edge Nektar::NekMeshUtils::CADCurveOCE::m_occEdge

private

OpenCascade edge.

Definition at line 103 of file CADCurveOCE.h.

Referenced by Initialise().