Nektar::LibUtilities::CADSystem Class Reference

Base class for CAD interface system. More...

#include <CADSystem.h>

Collaboration diagram for Nektar::LibUtilities::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 and curve maps. More...
void	Report ()
	Reports basic properties to screen. More...
Array< OneD, NekDouble >	GetBoundingBox ()
	Returns bounding box of the domain. More...
int	GetNumSurf ()
	Return number of surfaces. More...
int	GetNumCurve ()
	Return number of curves. More...
const CADCurveSharedPtr	GetCurve (int i)
	Gets curve type from map. More...
CADSurfSharedPtr	GetSurf (int i)
	Gets suface from map. More...
int	GetEPC ()
	Return Euler-Poincare number. More...

Private Member Functions
void	AddCurve (int i, TopoDS_Shape in)
	Private function to add curve to CADSystem::m_curves. More...
void	AddSurf (int i, TopoDS_Shape in, std::vector< std::vector< std::pair< int, int > > > ein)
	Private 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...
int	m_epc
	Euler-Poincare number of the CAD. More...
std::map< int, CADCurveSharedPtr >	m_curves
	map of curves More...
std::map< int, CADSurfSharedPtr >	m_surfs
	map of surfaces 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 59 of file CADSystem.h.

Constructor & Destructor Documentation

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

inline

Default constructor.

Definition at line 67 of file CADSystem.h.

                                                              : m_name(name)
        {
        }

Member Function Documentation

void Nektar::LibUtilities::CADSystem::AddCurve ( int  i, TopoDS_Shape  in  )

private

Private function to add curve to CADSystem::m_curves.

Definition at line 261 of file CADSystem.cpp.

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

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

private

Private function to add surface to CADSystem::m_surfs.

Definition at line 268 of file CADSystem.cpp.

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

Array< OneD, NekDouble > Nektar::LibUtilities::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 77 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;
}

const CADCurveSharedPtr Nektar::LibUtilities::CADSystem::GetCurve ( int  i )

inline

Gets curve type from map.

Definition at line 95 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;
        }

int Nektar::LibUtilities::CADSystem::GetEPC ( )

inline

Return Euler-Poincare number.

Definition at line 119 of file CADSystem.h.

References m_epc.

        {
            return m_epc;
        }

string Nektar::LibUtilities::CADSystem::GetName

(

)

Return the name of the CAD system.

Definition at line 53 of file CADSystem.cpp.

{
    return m_name;
}

int Nektar::LibUtilities::CADSystem::GetNumCurve ( )

inline

Return number of curves.

Definition at line 87 of file CADSystem.h.

References m_curves.

        {
            return m_curves.size();
        }

int Nektar::LibUtilities::CADSystem::GetNumSurf ( )

inline

Return number of surfaces.

Definition at line 79 of file CADSystem.h.

References m_surfs.

        {
            return m_surfs.size();
        }

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

inline

Gets suface from map.

Definition at line 107 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;
        }

bool Nektar::LibUtilities::CADSystem::LoadCAD

(

)

Initialises CAD and makes surface and curve maps.

Returns

true if completed successfully

Definition at line 111 of file CADSystem.cpp.

References ASSERTL0, 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);

    TopoDS_Shape shape;

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

    // From OpenCascade maps calculate Euler-Poincare number.
    TopTools_IndexedMapOfShape fc, vc, ec;
    TopExp::MapShapes(shape, TopAbs_FACE, fc);
    TopExp::MapShapes(shape, TopAbs_EDGE, ec);
    TopExp::MapShapes(shape, TopAbs_VERTEX, vc);

    m_epc = vc.Extent() - ec.Extent() + fc.Extent();

    TopTools_IndexedMapOfShape mapOfFaces;
    TopTools_IndexedMapOfShape mapOfEdges;
    TopExp::MapShapes(shape, TopAbs_FACE, mapOfFaces);

    // 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;

    // Adds edges to our type and map
    for(int i = 1; i <= mapOfEdges.Extent(); i++)
    {
        TopoDS_Shape edge = mapOfEdges.FindKey(i);
        AddCurve(i, edge);
    }

    // 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++)
    {
        vector<vector<pair<int,int> > > edges;

        TopoDS_Shape face = mapOfFaces.FindKey(i);

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

        for(int j = 1; j <= mapOfWires.Extent(); j++)
        {
            vector<pair<int,int> > edgeloop;

            TopoDS_Shape wire = mapOfWires.FindKey(j);

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

            BRepTools_WireExplorer exp;

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

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

                if(mapOfEdges.Contains(edge))
                {
                    pair<int,int> e;
                    e.first = mapOfEdges.FindIndex(edge);
                    adjsurfmap[e.first].push_back(i);
                    e.second = exp.Orientation();
                    edgeloop.push_back(e);
                }

                exp.Next();
            }

            edges.push_back(edgeloop);
        }

        AddSurf(i, face, edges);
    }

    // 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");
        m_curves[it->first]->SetAdjSurf(it->second);
    }

    return true;
}

void Nektar::LibUtilities::CADSystem::Report

(

)

Reports basic properties to screen.

Definition at line 61 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;
    cout << "\tCAD Euler-Poincaré characteristic: " << m_epc << endl;
}

Friends And Related Function Documentation

friend class MemoryManager< CADSystem >

friend

Definition at line 62 of file CADSystem.h.

Member Data Documentation

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

private

map of curves

Definition at line 135 of file CADSystem.h.

Referenced by GetCurve(), and GetNumCurve().

int Nektar::LibUtilities::CADSystem::m_epc

private

Euler-Poincare number of the CAD.

Definition at line 133 of file CADSystem.h.

Referenced by GetEPC().

std::string Nektar::LibUtilities::CADSystem::m_name

private

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

Definition at line 131 of file CADSystem.h.

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

private

map of surfaces

Definition at line 137 of file CADSystem.h.

Referenced by GetNumSurf(), and GetSurf().