Nektar::NekMeshUtils::CADSystem Class Reference

Base class for CAD interface system. More...

#include <CADSystem.h>

Collaboration diagram for Nektar::NekMeshUtils::CADSystem:

Public Member Functions
	CADSystem (const std::string &name)
	Default constructor. More...
std::string	GetName ()
	Return the name of the CAD system. More...
bool	LoadCAD ()
	Initialises CAD and makes surface, curve and vertex maps. More...
void	Report ()
	Reports basic properties to screen. More...
Array< OneD, NekDouble >	GetBoundingBox ()
	Returns bounding box of the domain. More...
int	GetNumSurf ()
	Get the number of surfaces. More...
int	GetNumCurve ()
	Get the number of curves. More...
CADCurveSharedPtr	GetCurve (int i)
	Gets a curve from the map. More...
CADSurfSharedPtr	GetSurf (int i)
	Gets a surface from the map. More...
std::map< int, CADVertSharedPtr >	GetVerts ()
	Gets map of all vertices. More...
int	GetNumVerts ()
	Gets number of vertices. More...
bool	InsideShape (Array< OneD, NekDouble > loc)
	based on location in space, uses opencascade routines to determin if the point is within the domain. This routine is slow and should be used sparingly, it is smart enough to take and form of geometry More...

Private Member Functions
void	AddVert (int i, TopoDS_Shape in)
	Function to add curve to CADSystem::m_verts. More...
void	AddCurve (int i, TopoDS_Shape in, int fv, int lv)
	Function to add curve to CADSystem::m_curves. More...
void	AddSurf (int i, TopoDS_Shape in, std::vector< EdgeLoop > ein)
	Function to add surface to CADSystem::m_surfs. More...

Private Attributes
std::string	m_name
	Name of cad file to be opened, including file extension. More...
std::map< int, CADCurveSharedPtr >	m_curves
	Map of curves. More...
std::map< int, CADSurfSharedPtr >	m_surfs
	Map of surfaces. More...
std::map< int, CADVertSharedPtr >	m_verts
	Map of vertices. More...
TopoDS_Shape	shape
	OCC master object. More...

Friends
class	MemoryManager< CADSystem >

Detailed Description

Base class for CAD interface system.

A class which can load and interact with CAD for Nektar++ using OpenCascade. This class contains maps to subclasses surface and curves.

Definition at line 65 of file CADSystem.h.

Constructor & Destructor Documentation

Nektar::NekMeshUtils::CADSystem::CADSystem ( const std::string &  name )

inline

Default constructor.

Definition at line 73 of file CADSystem.h.

                                     : m_name(name)
    {
    }

Member Function Documentation

void Nektar::NekMeshUtils::CADSystem::AddCurve ( int  i, TopoDS_Shape  in, int  fv, int  lv  )

private

Function to add curve to CADSystem::m_curves.

Definition at line 358 of file CADSystem.cpp.

{
    CADCurveSharedPtr newCurve =
        MemoryManager<CADCurve>::AllocateSharedPtr(i, in);

    vector<CADVertSharedPtr> vs;
    vs.push_back(m_verts[fv]);
    vs.push_back(m_verts[lv]);
    m_curves[i] = newCurve;
    m_curves[i]->SetVert(vs);
}

void Nektar::NekMeshUtils::CADSystem::AddSurf ( int  i, TopoDS_Shape  in, std::vector< EdgeLoop >  ein  )

private

Function to add surface to CADSystem::m_surfs.

Definition at line 370 of file CADSystem.cpp.

References ASSERTL0.

{
    CADSurfSharedPtr newSurf =
        MemoryManager<CADSurf>::AllocateSharedPtr(i, in, ein);
    m_surfs[i] = newSurf;

    if (in.Orientation() == 0)
    {
        m_surfs[i]->SetReverseNomral();
    }

    int tote = 0;
    for (int i = 0; i < ein.size(); i++)
    {
        tote += ein[i].edges.size();
    }

    ASSERTL0(tote != 1, "cannot handle periodic curves");

    if (tote == 2)
    {
        m_surfs[i]->SetTwoC();
    }
}

void Nektar::NekMeshUtils::CADSystem::AddVert ( int  i, TopoDS_Shape  in  )

private

Function to add curve to CADSystem::m_verts.

Definition at line 351 of file CADSystem.cpp.

{
    CADVertSharedPtr newVert = MemoryManager<CADVert>::AllocateSharedPtr(i, in);

    m_verts[i] = newVert;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::CADSystem::GetBoundingBox

(

)

Returns bounding box of the domain.

Gets the bounding box of the domain by considering the start and end points of each curve in the geometry.

Returns

Array with 6 entries: xmin, xmax, ymin, ymax, zmin and zmax.

Definition at line 64 of file CADSystem.cpp.

{
    Array<OneD, NekDouble> bound(6);
    bound[0] = numeric_limits<double>::max(); // xmin
    bound[1] = numeric_limits<double>::min(); // xmax
    bound[2] = numeric_limits<double>::max(); // ymin
    bound[3] = numeric_limits<double>::min(); // ymax
    bound[4] = numeric_limits<double>::max(); // zmin
    bound[5] = numeric_limits<double>::min(); // zmax

    for (int i = 1; i <= m_curves.size(); i++)
    {
        gp_Pnt start, end;
        CADCurveSharedPtr c         = GetCurve(i);
        Array<OneD, NekDouble> ends = c->GetMinMax();

        bound[0] = min(bound[0], min(ends[0], ends[3]));
        bound[1] = max(bound[1], max(ends[0], ends[3]));

        bound[2] = min(bound[2], min(ends[1], ends[4]));
        bound[3] = max(bound[3], max(ends[1], ends[4]));

        bound[4] = min(bound[4], min(ends[2], ends[5]));
        bound[5] = max(bound[5], max(ends[2], ends[5]));
    }

    return bound;
}

CADCurveSharedPtr Nektar::NekMeshUtils::CADSystem::GetCurve ( int  i )

inline

Gets a curve from the map.

Definition at line 123 of file CADSystem.h.

References ASSERTL0, Nektar::iterator, and m_curves.

    {
        std::map<int, CADCurveSharedPtr>::iterator search = m_curves.find(i);
        ASSERTL0(search != m_curves.end(), "curve does not exist");

        return search->second;
    }

string Nektar::NekMeshUtils::CADSystem::GetName

(

)

Return the name of the CAD system.

Definition at line 52 of file CADSystem.cpp.

{
    return m_name;
}

int Nektar::NekMeshUtils::CADSystem::GetNumCurve ( )

inline

Get the number of curves.

Definition at line 115 of file CADSystem.h.

References m_curves.

    {
        return m_curves.size();
    }

int Nektar::NekMeshUtils::CADSystem::GetNumSurf ( )

inline

Get the number of surfaces.

Definition at line 107 of file CADSystem.h.

References m_surfs.

    {
        return m_surfs.size();
    }

int Nektar::NekMeshUtils::CADSystem::GetNumVerts ( )

inline

Gets number of vertices.

Definition at line 153 of file CADSystem.h.

References m_verts.

    {
        return m_verts.size();
    }

CADSurfSharedPtr Nektar::NekMeshUtils::CADSystem::GetSurf ( int  i )

inline

Gets a surface from the map.

Definition at line 134 of file CADSystem.h.

References ASSERTL0, Nektar::iterator, and m_surfs.

    {
        std::map<int, CADSurfSharedPtr>::iterator search = m_surfs.find(i);
        ASSERTL0(search != m_surfs.end(), "surface does not exist");

        return search->second;
    }

std::map<int, CADVertSharedPtr> Nektar::NekMeshUtils::CADSystem::GetVerts ( )

inline

Gets map of all vertices.

Definition at line 145 of file CADSystem.h.

References m_verts.

    {
        return m_verts;
    }

bool Nektar::NekMeshUtils::CADSystem::InsideShape

(

Array< OneD, NekDouble >

loc

)

based on location in space, uses opencascade routines to determin if the point is within the domain. This routine is slow and should be used sparingly, it is smart enough to take and form of geometry

Definition at line 395 of file CADSystem.cpp.

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

    BRepClass3d_SolidClassifier test(shape, p, 1E-7);
    if (test.State() == TopAbs_IN || test.State() == TopAbs_ON)
    {
        return true;
    }
    else
    {
        return false;
    }
}

bool Nektar::NekMeshUtils::CADSystem::LoadCAD

(

)

Initialises CAD and makes surface, curve and vertex maps.

Returns

true if completed successfully

TODO fix this test

Definition at line 93 of file CADSystem.cpp.

References ASSERTL0, Nektar::NekMeshUtils::EdgeLoop::center, Nektar::NekMeshUtils::curve, Nektar::NekMeshUtils::EdgeLoop::edgeo, Nektar::NekMeshUtils::EdgeLoop::edges, and Nektar::iterator.

{
    if (!boost::filesystem::exists(m_name.c_str()))
    {
        return false;
    }

    string ext;
    size_t pos = m_name.find(".");
    ext        = m_name.substr(pos);

    if (boost::iequals(ext, ".STEP") || boost::iequals(ext, ".STP"))
    {
        // Takes step file and makes OpenCascade shape
        STEPControl_Reader reader;
        reader = STEPControl_Reader();
        reader.ReadFile(m_name.c_str());
        reader.NbRootsForTransfer();
        reader.TransferRoots();
        shape = reader.OneShape();
        if (shape.IsNull())
        {
            return false;
        }
    }
    else if (boost::iequals(ext, ".IGES") || boost::iequals(ext, ".IGS"))
    {
        // Takes IGES file and makes OpenCascade shape
        IGESControl_Reader reader;
        reader = IGESControl_Reader();
        reader.ReadFile(m_name.c_str());
        reader.NbRootsForTransfer();
        reader.TransferRoots();
        shape = reader.OneShape();
        if (shape.IsNull())
        {
            return false;
        }
    }
    else
    {
        return false;
    }

    // faces and verts can be extracted straight from shape
    TopTools_IndexedMapOfShape mapOfVerts, mapOfFaces;
    TopExp::MapShapes(shape, TopAbs_VERTEX, mapOfVerts);
    TopExp::MapShapes(shape, TopAbs_FACE, mapOfFaces);

    // edges need to be built from loops around faces to elimiate degen and
    // hanging edges
    TopTools_IndexedMapOfShape mapOfEdges;

    // build map of verticies
    for (int i = 1; i <= mapOfVerts.Extent(); i++)
    {
        TopoDS_Shape v = mapOfVerts.FindKey(i);
        AddVert(i, v);
    }

    // For each face of the geometry, get the local edges which bound it. If
    // they are valid (their type != 7), then add them to an edge map. This
    // filters out the dummy edges which OCC uses.
    for (int i = 1; i <= mapOfFaces.Extent(); i++)
    {
        TopoDS_Shape face = mapOfFaces.FindKey(i);

        TopTools_IndexedMapOfShape localEdges;
        TopExp::MapShapes(face, TopAbs_EDGE, localEdges);

        for (int j = 1; j <= localEdges.Extent(); j++)
        {
            TopoDS_Shape edge       = localEdges.FindKey(j);
            BRepAdaptor_Curve curve = BRepAdaptor_Curve(TopoDS::Edge(edge));
            if (curve.GetType() != 7)
            {
                if (!(mapOfEdges.Contains(edge)))
                {
                    mapOfEdges.Add(edge);
                }
            }
        }
    }

    map<int, vector<int> > adjsurfmap; // from id of curve to list of ids of
                                       // surfs

    // Adds edges to our type and map
    for (int i = 1; i <= mapOfEdges.Extent(); i++)
    {
        TopoDS_Shape edge = mapOfEdges.FindKey(i);
        TopoDS_Vertex fv =
            TopExp::FirstVertex(TopoDS::Edge(edge), Standard_True);
        TopoDS_Vertex lv =
            TopExp::LastVertex(TopoDS::Edge(edge), Standard_True);

        if (edge.Orientation() == 0)
        {
            AddCurve(i, edge, mapOfVerts.FindIndex(fv),
                     mapOfVerts.FindIndex(lv));
        }
        else
        {
            AddCurve(i, edge, mapOfVerts.FindIndex(lv),
                     mapOfVerts.FindIndex(fv));
        }
    }

    // For each face, examine all the wires (i.e. bounding loops) and
    // investigates the loop. Using this information, connectivity is determined
    // and edges are associated with surfaces.
    for (int i = 1; i <= mapOfFaces.Extent(); i++)
    {
        TopoDS_Shape face = mapOfFaces.FindKey(i);

        TopTools_IndexedMapOfShape mapOfWires;
        TopExp::MapShapes(face, TopAbs_WIRE, mapOfWires);

        // this pice of code does an idiot check on the loops to make sure
        // they dont cross or touch
        if (mapOfWires.Extent() > 1)
        {
            TopoDS_Wire ow = BRepTools::OuterWire(TopoDS::Face(face));

            vector<TopoDS_Shape> wirefacecuts;
            vector<gp_Pnt> centersofcutfaces;

            for (int j = 1; j <= mapOfWires.Extent(); j++)
            {
                TopoDS_Shape wire = mapOfWires.FindKey(j);

                if (wire != ow)
                {
                    BRepBuilderAPI_MakeFace build(
                        BRep_Tool::Surface(TopoDS::Face(face)), 1e-7);
                    build.Add(TopoDS::Wire(wire));
                    TopoDS_Shape newface = build.Shape();
                    wirefacecuts.push_back(newface);
                    BRepAdaptor_Surface b =
                        BRepAdaptor_Surface(TopoDS::Face(newface));
                    NekDouble u, v;
                    u = (b.LastUParameter() - b.FirstUParameter()) / 2.0;
                    v = (b.LastVParameter() - b.FirstVParameter()) / 2.0;
                    centersofcutfaces.push_back(b.Value(u, v));
                }
            }
            for (int j = 0; j < wirefacecuts.size(); j++)
            {
                for (int k = 0; k < wirefacecuts.size(); k++)
                {
                    if (j == k)
                        continue;

                    /// TODO fix this test
                    BRepClass_FaceClassifier fc(TopoDS::Face(wirefacecuts[j]),
                                                centersofcutfaces[k], 1e-7);
                    // ASSERTL0(fc.State() == 1, "Internal face loops make this
                    // cad impossible to mesh");
                }
            }
        }

        vector<EdgeLoop> edgeloops;

        // now we acutally analyse the loops for cad building
        for (int j = 1; j <= mapOfWires.Extent(); j++)
        {
            EdgeLoop edgeloop;

            TopoDS_Shape wire = mapOfWires.FindKey(j);

            ShapeAnalysis_Wire wiretest(TopoDS::Wire(wire), TopoDS::Face(face),
                                        1E-6);

            // calculate the center of the wire
            GProp_GProps massProps;
            BRepGProp::SurfaceProperties(wire, massProps);
            gp_Pnt gPt = massProps.CentreOfMass();
            // alternative locuv methods
            ShapeAnalysis_Surface sas(BRep_Tool::Surface(TopoDS::Face(face)));
            sas.SetDomain(
                BRepAdaptor_Surface(TopoDS::Face(face)).FirstUParameter(),
                BRepAdaptor_Surface(TopoDS::Face(face)).LastUParameter(),
                BRepAdaptor_Surface(TopoDS::Face(face)).FirstVParameter(),
                BRepAdaptor_Surface(TopoDS::Face(face)).LastVParameter());
            gp_Pnt2d p2 = sas.ValueOfUV(gPt, 1e-7);
            Array<OneD, NekDouble> cen(2);
            cen[0]          = p2.X();
            cen[1]          = p2.Y();
            edgeloop.center = cen;

            BRepTools_WireExplorer exp;

            exp.Init(TopoDS::Wire(wire));

            while (exp.More())
            {
                TopoDS_Shape edge = exp.Current();

                if (mapOfEdges.Contains(edge))
                {
                    int e = mapOfEdges.FindIndex(edge);
                    edgeloop.edges.push_back(m_curves[e]);
                    edgeloop.edgeo.push_back(exp.Orientation());
                    adjsurfmap[e].push_back(i);
                }

                exp.Next();
            }

            edgeloops.push_back(edgeloop);
        }

        AddSurf(i, face, edgeloops);
    }

    // attempts to identify properties of the vertex on the degen edge
    for (int i = 1; i <= mapOfFaces.Extent(); i++)
    {
        TopoDS_Shape face = mapOfFaces.FindKey(i);

        TopTools_IndexedMapOfShape localEdges;
        TopExp::MapShapes(face, TopAbs_EDGE, localEdges);

        for (int j = 1; j <= localEdges.Extent(); j++)
        {
            TopoDS_Shape edge = localEdges.FindKey(j);
            if (BRep_Tool::Degenerated(TopoDS::Edge(edge)))
            {
                cout << "degen edge on face " << i << endl;
                gp_Pnt2d p1, p2;

                BRep_Tool::UVPoints(TopoDS::Edge(edge), TopoDS::Face(face), p1,
                                    p2);

                m_verts[mapOfVerts.FindIndex(TopExp::FirstVertex(
                            TopoDS::Edge(edge), Standard_True))]
                    ->SetDegen(i, m_surfs[i], (p1.X() + p2.X()) / 2.0,
                               (p1.Y() + p2.Y()) / 2.0);
            }
        }
    }

    // This checks that all edges are bound by two surfaces, sanity check.
    for (map<int, vector<int> >::iterator it = adjsurfmap.begin();
         it != adjsurfmap.end(); it++)
    {
        ASSERTL0(it->second.size() == 2, "no three curve surfaces");
        vector<CADSurfSharedPtr> sfs;
        for (int i = 0; i < it->second.size(); i++)
        {
            sfs.push_back(m_surfs[it->second[i]]);
        }
        m_curves[it->first]->SetAdjSurf(sfs);
    }
    return true;
}

void Nektar::NekMeshUtils::CADSystem::Report

(

)

Reports basic properties to screen.

Definition at line 57 of file CADSystem.cpp.

{
    cout << endl << "CAD report:" << endl;
    cout << "\tCAD has: " << m_curves.size() << " curves." << endl;
    cout << "\tCAD has: " << m_surfs.size() << " surfaces." << endl;
}

Friends And Related Function Documentation

friend class MemoryManager< CADSystem >

friend

Definition at line 68 of file CADSystem.h.

Member Data Documentation

std::map<int, CADCurveSharedPtr> Nektar::NekMeshUtils::CADSystem::m_curves

private

Map of curves.

Definition at line 176 of file CADSystem.h.

Referenced by GetCurve(), and GetNumCurve().

std::string Nektar::NekMeshUtils::CADSystem::m_name

private

Name of cad file to be opened, including file extension.

Definition at line 174 of file CADSystem.h.

std::map<int, CADSurfSharedPtr> Nektar::NekMeshUtils::CADSystem::m_surfs

private

Map of surfaces.

Definition at line 178 of file CADSystem.h.

Referenced by GetNumSurf(), and GetSurf().

std::map<int, CADVertSharedPtr> Nektar::NekMeshUtils::CADSystem::m_verts

private

Map of vertices.

Definition at line 180 of file CADSystem.h.

Referenced by GetNumVerts(), and GetVerts().

TopoDS_Shape Nektar::NekMeshUtils::CADSystem::shape

private

OCC master object.

Definition at line 182 of file CADSystem.h.