Nektar::LibUtilities::CADSurf Class Reference

#include <CADSurf.h>

Collaboration diagram for Nektar::LibUtilities::CADSurf:

Public Member Functions
	CADSurf (int i, TopoDS_Shape in, std::vector< std::vector< std::pair< int, int > > > ein)
	Default constructor. More...
std::vector< std::vector< std::pair< int, int > > >	GetEdges ()
	Get the IDs of the edges which bound the surface. More...
Array< OneD, NekDouble >	GetBounds ()
	Get the limits of the parametric space for the surface. More...
Array< OneD, NekDouble >	N (Array< OneD, NekDouble > uv)
	Get the normal vector at parametric point u,v. More...
Array< OneD, NekDouble >	D1 (Array< OneD, NekDouble > uv)
	Get the set of first derivatives at parametric point u,v. More...
Array< OneD, NekDouble >	D2 (Array< OneD, NekDouble > uv)
	Get the set of second derivatives at parametric point u,v. More...
Array< OneD, NekDouble >	P (Array< OneD, NekDouble > uv)
	Get the x,y,z at parametric point u,v. More...
Array< OneD, NekDouble >	locuv (Array< OneD, NekDouble > p)
	Performs a reverse look up to find u,v and x,y,z. More...

Private Member Functions
	Handle (Geom_Surface) m_s
	Alternate OpenCascade object for surface. Used by reverse lookup. More...

Private Attributes
int	m_ID
	ID of surface. More...
BRepAdaptor_Surface	m_occSurface
	OpenCascade object for surface. More...
std::vector< std::vector< std::pair< int, int > > >	m_edges
	List of bounding edges in loops with orientation. More...

Friends
class	MemoryManager< CADSurf >

Detailed Description

Definition at line 50 of file CADSurf.h.

Constructor & Destructor Documentation

Nektar::LibUtilities::CADSurf::CADSurf	(	int	i,
		TopoDS_Shape	in,
		std::vector< std::vector< std::pair< int, int > > >	ein
	)

Default constructor.

Definition at line 43 of file CADSurf.cpp.

References m_occSurface.

                                                      : m_ID(i), m_edges(ein)
{
    // this bit of code changes the units of the cad from mm opencascade
    // defualt to m
    gp_Trsf transform;
    gp_Pnt ori(0.0, 0.0, 0.0);
    transform.SetScale(ori, 1.0 / 1000.0);
    TopLoc_Location mv(transform);
    m_s = BRep_Tool::Surface(TopoDS::Face(in));
    in.Move(mv);
    m_occSurface = BRepAdaptor_Surface(TopoDS::Face(in));
}

Member Function Documentation

Array< OneD, NekDouble > Nektar::LibUtilities::CADSurf::D1

(

Array< OneD, NekDouble >

uv

)

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.

Definition at line 151 of file CADSurf.cpp.

References m_occSurface.

{
    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::LibUtilities::CADSurf::D2

(

Array< OneD, NekDouble >

uv

)

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

Definition at line 178 of file CADSurf.cpp.

References m_occSurface.

{
    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;
}

Array<OneD, NekDouble> Nektar::LibUtilities::CADSurf::GetBounds ( )

inline

Get the limits of the parametric space for the surface.

Returns

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

Definition at line 80 of file CADSurf.h.

References m_occSurface.

        {
            Array<OneD,NekDouble> b(4);

            b[0] = m_occSurface.FirstUParameter();
            b[1] = m_occSurface.LastUParameter();
            b[2] = m_occSurface.FirstVParameter();
            b[3] = m_occSurface.LastVParameter();

            return b;
        }

std::vector<std::vector<std::pair<int,int> > > Nektar::LibUtilities::CADSurf::GetEdges ( )

inline

Get the IDs of the edges which bound the surface.

The edges are organsised into two vectors, which are grouped into the continuous loops of the bounding edges, then the edges, which are a pair of integers. The first item is the edge ID and the second is an integer that indicates whether this edge is orientated forwards or backwards on this surface to form the loop.

Definition at line 70 of file CADSurf.h.

References m_edges.

        {
            return m_edges;
        }

Nektar::LibUtilities::CADSurf::Handle ( Geom_Surface  )

private

Alternate OpenCascade object for surface. Used by reverse lookup.

Array< OneD, NekDouble > Nektar::LibUtilities::CADSurf::locuv

(

Array< OneD, NekDouble >

p

)

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

Todo:

create a check so that if the calculated uv is out of bounds

Definition at line 63 of file CADSurf.cpp.

References ASSERTL0, ASSERTL1, Nektar::NekConstants::GeomTol, and m_occSurface.

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

    GeomAPI_ProjectPointOnSurf projection(loc, m_s,
                        m_occSurface.FirstUParameter(),
                        m_occSurface.LastUParameter(),
                        m_occSurface.FirstVParameter(),
                        m_occSurface.LastVParameter(),
                        Extrema_ExtAlgo_Tree);

    ASSERTL0(projection.NbPoints() > 0, "locuv failed");

    Quantity_Parameter ui;
    Quantity_Parameter vi;

    projection.Parameters(1,ui,vi);

    /// @todo create a check so that if the calculated uv is out of bounds

    Array<OneD, NekDouble> uvr(2);
    uvr[0] = ui;
    uvr[1] = vi;

    ASSERTL1(projection.Distance(1) < NekConstants::GeomTol,
                "large locuv distance");

    return uvr;
}

Array< OneD, NekDouble > Nektar::LibUtilities::CADSurf::N

(

Array< OneD, NekDouble >

uv

)

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.

Definition at line 118 of file CADSurf.cpp.

References m_occSurface.

{
    Array<OneD, NekDouble> normal(3);
    gp_Pnt Loc;
    gp_Vec D1U, D1V;
    m_occSurface.D1(uv[0], uv[1], Loc, D1U, D1V);
    gp_Vec n = D1U.Crossed(D1V);

    if (n.X() == 0 && n.Y() == 0 && n.Z() == 0)
    {
        // Return bad normal
        normal[0] = 0.0;
        normal[1] = 0.0;
        normal[2] = 0.0;
    }
    else
    {
        n.Normalize();
        normal[0] = n.X();
        normal[1] = n.Y();
        normal[2] = n.Z();
    }

    return normal;
}

Array< OneD, NekDouble > Nektar::LibUtilities::CADSurf::P

(

Array< OneD, NekDouble >

uv

)

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

Parameters

uv	Array of u and v parametric coords.

Returns

Array of xyz location.

Todo:

create bound checking

Definition at line 100 of file CADSurf.cpp.

References m_occSurface.

{
    /// @todo create bound checking
    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;
}

Friends And Related Function Documentation

friend class MemoryManager< CADSurf >

friend

Definition at line 53 of file CADSurf.h.

Member Data Documentation

std::vector<std::vector<std::pair<int,int> > > Nektar::LibUtilities::CADSurf::m_edges

private

List of bounding edges in loops with orientation.

Definition at line 106 of file CADSurf.h.

Referenced by GetEdges().

int Nektar::LibUtilities::CADSurf::m_ID

private

ID of surface.

Definition at line 100 of file CADSurf.h.

BRepAdaptor_Surface Nektar::LibUtilities::CADSurf::m_occSurface

private

OpenCascade object for surface.

Definition at line 102 of file CADSurf.h.

Referenced by CADSurf(), D1(), D2(), GetBounds(), locuv(), N(), and P().