Commit 6415d988 authored by Chris Cantwell's avatar Chris Cantwell
Browse files

Switched ExpList to hold Expansions rather than StdExpansions.

parent 2184c1e8
......@@ -134,7 +134,7 @@ namespace Nektar
Array<OneD, unsigned int> interiorMap;
PeriodicMap::const_iterator pIt;
const StdRegions::StdExpansionVector &locExpVector = *(locExp.GetExp());
const LocalRegions::ExpansionVector &locExpVector = *(locExp.GetExp());
LibUtilities::CommSharedPtr vCommRow = m_comm->GetRowComm();
m_globalToUniversalMap = Nektar::Array<OneD, int>(m_numGlobalCoeffs, -1);
......@@ -145,7 +145,7 @@ namespace Nektar
// Loop over all the elements in the domain to gather mesh data
for(i = 0; i < locExpVector.size(); ++i)
{
locExpansion = LocalRegions::Expansion::FromStdExp(locExpVector[i]);
locExpansion = locExpVector[i];
nVert += locExpansion->GetNverts();
nEdge += locExpansion->GetNedges();
nFace += locExpansion->GetNfaces();
......@@ -176,7 +176,7 @@ namespace Nektar
// Assemble global to universal mapping for this process
for(i = 0; i < locExpVector.size(); ++i)
{
locExpansion = LocalRegions::Expansion::FromStdExp(locExpVector[i]);
locExpansion = locExpVector[i];
nDim = locExpansion->GetShapeDimension();
cnt = locExp.GetCoeff_Offset(i);
......@@ -314,7 +314,7 @@ namespace Nektar
int i, j;
int nverts = 0;
const boost::shared_ptr<StdRegions::StdExpansionVector> exp
const boost::shared_ptr<LocalRegions::ExpansionVector> exp
= locexp.GetExp();
int nelmts = exp->size();
......
......@@ -152,7 +152,7 @@ namespace Nektar
int graphVertId = 0;
int globalId;
const StdRegions::StdExpansionVector &locExpVector = *(locExp.GetExp());
const LocalRegions::ExpansionVector &locExpVector = *(locExp.GetExp());
LocalRegions::SegExpSharedPtr locSegExp;
......
......@@ -167,7 +167,7 @@ namespace Nektar
StdRegions::Orientation edgeOrient;
Array<OneD, unsigned int> edgeInteriorMap;
Array<OneD, int> edgeInteriorSign;
const StdRegions::StdExpansionVector &locExpVector = *(locExp.GetExp());
const LocalRegions::ExpansionVector &locExpVector = *(locExp.GetExp());
m_signChange = false;
m_systemSingular = false;
Array<OneD, map<int,int> > ReorderedGraphVertId(2);
......@@ -567,7 +567,7 @@ namespace Nektar
StdRegions::StdExpansion2DSharedPtr locExpansion;
LocalRegions::SegExpSharedPtr bndSegExp;
MultiRegions::ExpList0DSharedPtr bndVertExp;
const StdRegions::StdExpansionVector &locExpVector = *(locExp.GetExp());
const LocalRegions::ExpansionVector &locExpVector = *(locExp.GetExp());
map<int,int>::iterator mapIt;
map<int,int>::const_iterator mapConstIt;
bool systemSingular = true;
......
......@@ -203,7 +203,7 @@ namespace Nektar
Array<OneD, int> faceInteriorSign;
PeriodicMap::const_iterator pIt;
const StdRegions::StdExpansionVector &locExpVector = *(locExp.GetExp());
const LocalRegions::ExpansionVector &locExpVector = *(locExp.GetExp());
m_signChange = false;
//m_systemSingular = false;
......
......@@ -89,7 +89,7 @@ namespace Nektar
Array<OneD,unsigned int> vmap;
LocalRegions::SegExpSharedPtr locSegExp;
const boost::shared_ptr<StdRegions::StdExpansionVector> exp1D = locExp.GetExp();
const boost::shared_ptr<LocalRegions::ExpansionVector> exp1D = locExp.GetExp();
m_numGlobalBndCoeffs = exp1D->size()+1;
m_numGlobalCoeffs = m_numGlobalBndCoeffs;
......@@ -206,7 +206,7 @@ namespace Nektar
LocalRegions::TriExpSharedPtr locTriExp;
SpatialDomains::Geometry1DSharedPtr SegGeom;
const boost::shared_ptr<StdRegions::StdExpansionVector> exp2D = locExp.GetExp();
const boost::shared_ptr<LocalRegions::ExpansionVector> exp2D = locExp.GetExp();
int nel = exp2D->size();
map<int, int> MeshEdgeId;
......@@ -649,7 +649,7 @@ namespace Nektar
SpatialDomains::Geometry2DSharedPtr FaceGeom;
StdRegions::StdExpansionSharedPtr locBndExp;
const boost::shared_ptr<StdRegions::StdExpansionVector> exp3D =
const boost::shared_ptr<LocalRegions::ExpansionVector> exp3D =
locExp.GetExp();
int nel = exp3D->size();
......@@ -1028,7 +1028,7 @@ namespace Nektar
int nDim = 0;
int i,j,k;
const StdRegions::StdExpansionVector &locExpVector = *(locExp.GetExp());
const LocalRegions::ExpansionVector &locExpVector = *(locExp.GetExp());
// Initialise the global to universal maps.
m_globalToUniversalBndMap = Nektar::Array<OneD, int>(m_numGlobalBndCoeffs, -1);
......@@ -1208,7 +1208,7 @@ namespace Nektar
int i;
int maxQuad = 0, quad = 0, nDim = 0, eid = 0, offset = 0;
const StdRegions::StdExpansionVector &locExpVector = *(locExp.GetExp());
const LocalRegions::ExpansionVector &locExpVector = *(locExp.GetExp());
int nTracePhys = trace->GetTotPoints();
......@@ -1243,7 +1243,7 @@ namespace Nektar
{
for (int i = 0; i < trace->GetExpSize(); ++i)
{
eid = LocalRegions::Expansion::FromStdExp(trace->GetExp(i))->GetGeom()->GetGlobalID();
eid = trace->GetExp(i)->GetGeom()->GetGlobalID();
offset = trace->GetPhys_Offset(i);
m_traceToUniversalMap[offset] = eid*maxQuad+1;
}
......@@ -1252,7 +1252,7 @@ namespace Nektar
{
for (int i = 0; i < trace->GetExpSize(); ++i)
{
eid = LocalRegions::Expansion::FromStdExp(trace->GetExp(i))->GetGeom()->GetGlobalID();
eid = trace->GetExp(i)->GetGeom()->GetGlobalID();
offset = trace->GetPhys_Offset(i);
quad = trace->GetExp(i)->GetTotPoints();
......
......@@ -94,7 +94,7 @@ namespace Nektar
pSession, graph2D, variable, false,
false);
m_exp = MemoryManager<StdRegions::StdExpansionVector>
m_exp = MemoryManager<LocalRegions::ExpansionVector>
::AllocateSharedPtr();
for(n = 0; n < m_planes.num_elements(); ++n)
......
......@@ -82,7 +82,7 @@ namespace Nektar
m_lines[0] = line_zero = MemoryManager<ContField1D>::AllocateSharedPtr(pSession,graph1D,variable);
m_exp = MemoryManager<StdRegions::StdExpansionVector>::AllocateSharedPtr();
m_exp = MemoryManager<LocalRegions::ExpansionVector>::AllocateSharedPtr();
nel = m_lines[0]->GetExpSize();
for(i = 0; i < nel; ++i)
......
......@@ -1350,8 +1350,8 @@ namespace Nektar
LocalRegions::ExpansionSharedPtr exp;
for(cnt = n = 0; n < GetExpSize(); ++n)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[n]);
for(i = 0; i < (*m_exp)[n]->GetNverts(); ++i)
exp = (*m_exp)[n];
for(i = 0; i < exp->GetNverts(); ++i)
{
id = exp->GetGeom()->GetVid(i);
......
......@@ -1532,11 +1532,9 @@ namespace Nektar
// Populate global ID map (takes global geometry ID to local
// expansion list ID).
LocalRegions::ExpansionSharedPtr exp;
for (i = 0; i < GetExpSize(); ++i)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
globalIdMap[exp->GetGeom()->GetGlobalID()] = i;
globalIdMap[(*m_exp)[i]->GetGeom()->GetGlobalID()] = i;
}
// Determine number of boundary condition expansions.
......
......@@ -101,7 +101,7 @@ namespace Nektar
// note that nzplanes can be larger than nzmodes
m_planes[0] = plane_zero = MemoryManager<DisContField2D>::AllocateSharedPtr(pSession,graph2D,variable,True,False);
m_exp = MemoryManager<StdRegions::StdExpansionVector>::AllocateSharedPtr();
m_exp = MemoryManager<LocalRegions::ExpansionVector>::AllocateSharedPtr();
nel = m_planes[0]->GetExpSize();
for(i = 0; i < nel; ++i)
......@@ -177,7 +177,7 @@ namespace Nektar
if(boundaryCondition->GetBoundaryConditionType() != SpatialDomains::ePeriodic)
{
boost::shared_ptr<StdRegions::StdExpansionVector> exp = MemoryManager<StdRegions::StdExpansionVector>::AllocateSharedPtr();
boost::shared_ptr<LocalRegions::ExpansionVector> exp = MemoryManager<LocalRegions::ExpansionVector>::AllocateSharedPtr();
for(n = 0; n < nplanes; ++n)
{
......
......@@ -101,7 +101,7 @@ namespace Nektar
//
m_lines[0] = line_zero = MemoryManager<DisContField1D>::AllocateSharedPtr(pSession,graph1D,variable);
m_exp = MemoryManager<StdRegions::StdExpansionVector>::AllocateSharedPtr();
m_exp = MemoryManager<LocalRegions::ExpansionVector>::AllocateSharedPtr();
nel = m_lines[0]->GetExpSize();
for(i = 0; i < nel; ++i)
......
......@@ -87,7 +87,7 @@ namespace Nektar
m_coeffs(),
m_phys(),
m_physState(false),
m_exp(MemoryManager<StdRegions::StdExpansionVector>
m_exp(MemoryManager<LocalRegions::ExpansionVector>
::AllocateSharedPtr()),
m_coeff_offset(),
m_phys_offset(),
......@@ -113,7 +113,7 @@ namespace Nektar
m_coeffs(),
m_phys(),
m_physState(false),
m_exp(MemoryManager<StdRegions::StdExpansionVector>
m_exp(MemoryManager<LocalRegions::ExpansionVector>
::AllocateSharedPtr()),
m_coeff_offset(),
m_phys_offset(),
......@@ -140,7 +140,7 @@ namespace Nektar
m_coeffs(),
m_phys(),
m_physState(false),
m_exp(MemoryManager<StdRegions::StdExpansionVector>
m_exp(MemoryManager<LocalRegions::ExpansionVector>
::AllocateSharedPtr()),
m_coeff_offset(),
m_phys_offset(),
......@@ -1303,15 +1303,13 @@ namespace Nektar
}
}
StdRegions::StdExpansionSharedPtr& ExpList::GetExp(
LocalRegions::ExpansionSharedPtr& ExpList::GetExp(
const Array<OneD, const NekDouble> &gloCoord)
{
Array<OneD, NekDouble> stdCoord(GetCoordim(0),0.0);
for (int i = 0; i < (*m_exp).size(); ++i)
{
LocalRegions::ExpansionSharedPtr exp =
LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
if (exp->GetGeom()->ContainsPoint(gloCoord))
if ((*m_exp)[i]->GetGeom()->ContainsPoint(gloCoord))
{
return (*m_exp)[i];
}
......@@ -1341,13 +1339,10 @@ namespace Nektar
NekDouble tol)
{
static int start = 0;
LocalRegions::ExpansionSharedPtr exp;
// start search at previous element or 0
for (int i = start; i < (*m_exp).size(); ++i)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
if (exp->GetGeom()->ContainsPoint(gloCoords, locCoords, tol))
if ((*m_exp)[i]->GetGeom()->ContainsPoint(gloCoords, locCoords, tol))
{
start = i;
return i;
......@@ -1356,9 +1351,7 @@ namespace Nektar
for (int i = 0; i < start; ++i)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
if (exp->GetGeom()->ContainsPoint(gloCoords, locCoords, tol))
if ((*m_exp)[i]->GetGeom()->ContainsPoint(gloCoords, locCoords, tol))
{
start = i;
return i;
......@@ -1393,9 +1386,7 @@ namespace Nektar
Array<OneD,Array<OneD, NekDouble> > loctangent;
// Assume whole array is of same coordinate dimension
LocalRegions::ExpansionSharedPtr exp =
LocalRegions::Expansion::FromStdExp((*m_exp)[0]);
int coordim = exp->GetGeom()->GetCoordim();
int coordim = (*m_exp)[0]->GetGeom()->GetCoordim();
ASSERTL0(tangents.num_elements() > 0,
"Must have storage for at least one tangent");
......@@ -1407,13 +1398,12 @@ namespace Nektar
for(i = 0; i < m_exp->size(); ++i)
{
// Get the number of points and normals for this expansion.
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
e_npoints = exp->GetTotPoints();
e_npoints = (*m_exp)[i]->GetTotPoints();
offset = m_phys_offset[i];
for (j = 0; j < tangents.num_elements(); ++j)
{
loctangent = exp->GetMetricInfo()->GetTangent(j);
loctangent = (*m_exp)[i]->GetMetricInfo()->GetTangent(j);
// Get the physical data offset for this expansion.
for (k = 0; k < coordim; ++k)
......@@ -1441,13 +1431,11 @@ namespace Nektar
m_session->LoadGeometricInfo("TangentCentreX",coords[0],0.0);
m_session->LoadGeometricInfo("TangentCentreY",coords[1],0.0);
LocalRegions::ExpansionSharedPtr exp;
// Apply geometric info to each expansion.
for (int i = 0; i < m_exp->size(); ++i)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
exp->GetMetricInfo()->SetTangentOrientation(dir);
exp->GetMetricInfo()->SetTangentCircularCentre(coords);
(*m_exp)[i]->GetMetricInfo()->SetTangentOrientation(dir);
(*m_exp)[i]->GetMetricInfo()->SetTangentCircularCentre(coords);
}
}
......@@ -2191,7 +2179,6 @@ namespace Nektar
for(s = startenum; s <= endenum; ++s)
{
LocalRegions::ExpansionSharedPtr exp;
LibUtilities::ShapeType shape;
std::vector<unsigned int> elementIDs;
std::vector<LibUtilities::BasisType> basis;
......@@ -2204,10 +2191,9 @@ namespace Nektar
for(int i = 0; i < (*m_exp).size(); ++i)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
if(exp->GetGeom()->GetShapeType() == (LibUtilities::ShapeType) s)
if((*m_exp)[i]->GetGeom()->GetShapeType() == (LibUtilities::ShapeType) s)
{
elementIDs.push_back(exp->GetGeom()->GetGlobalID());
elementIDs.push_back((*m_exp)[i]->GetGeom()->GetGlobalID());
if(first)
{
shape = (LibUtilities::ShapeType) s;
......@@ -2280,7 +2266,6 @@ namespace Nektar
void ExpList::v_AppendFieldData(LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector<NekDouble> &fielddata, Array<OneD, NekDouble> &coeffs)
{
LocalRegions::ExpansionSharedPtr exp;
int i;
// Determine mapping from element ids to location in
// expansion list
......@@ -2289,8 +2274,7 @@ namespace Nektar
map<int, int> ElmtID_to_ExpID;
for(i = 0; i < (*m_exp).size(); ++i)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
ElmtID_to_ExpID[exp->GetGeom()->GetGlobalID()] = i;
ElmtID_to_ExpID[(*m_exp)[i]->GetGeom()->GetGlobalID()] = i;
}
for(i = 0; i < fielddef->m_elementIDs.size(); ++i)
......@@ -2355,11 +2339,9 @@ namespace Nektar
// Loop in reverse order so that in case where using a Homogeneous
// expansion it sets geometry ids to first part of m_exp
// list. Otherwise will set to second (complex) expansion
LocalRegions::ExpansionSharedPtr exp;
for(i = (*m_exp).size()-1; i >= 0; --i)
{
exp = LocalRegions::Expansion::FromStdExp((*m_exp)[i]);
elmtToExpId[exp->GetGeom()->GetGlobalID()] = i;
elmtToExpId[(*m_exp)[i]->GetGeom()->GetGlobalID()] = i;
}
for (i = 0; i < fielddef->m_elementIDs.size(); ++i)
......
......@@ -40,7 +40,7 @@
#include <LibUtilities/Communication/Comm.h>
#include <LibUtilities/BasicUtils/SessionReader.h>
#include <MultiRegions/MultiRegions.hpp>
#include <StdRegions/StdExpansion.h>
#include <LocalRegions/Expansion.h>
#include <MultiRegions/GlobalMatrix.h>
#include <MultiRegions/GlobalMatrixKey.h>
#include <SpatialDomains/MeshGraph.h>
......@@ -555,15 +555,15 @@ namespace Nektar
}
/// This function returns the vector of elements in the expansion.
inline const boost::shared_ptr<StdRegions::StdExpansionVector> GetExp() const;
inline const boost::shared_ptr<LocalRegions::ExpansionVector> GetExp() const;
/// This function returns (a shared pointer to) the local elemental
/// expansion of the \f$n^{\mathrm{th}}\f$ element.
inline StdRegions::StdExpansionSharedPtr& GetExp(int n) const;
inline LocalRegions::ExpansionSharedPtr& GetExp(int n) const;
/// This function returns (a shared pointer to) the local elemental
/// expansion containing the arbitrary point given by \a gloCoord.
MULTI_REGIONS_EXPORT StdRegions::StdExpansionSharedPtr& GetExp(
MULTI_REGIONS_EXPORT LocalRegions::ExpansionSharedPtr& GetExp(
const Array<OneD, const NekDouble> &gloCoord);
/// This function returns the index of the local elemental
......@@ -908,7 +908,7 @@ namespace Nektar
* where most of the routines for the derived classes are defined
* in the #ExpList base class.
*/
boost::shared_ptr<StdRegions::StdExpansionVector> m_exp;
boost::shared_ptr<LocalRegions::ExpansionVector> m_exp;
/// Offset of elemental data into the array #m_coeffs
Array<OneD, int> m_coeff_offset;
......@@ -1788,7 +1788,7 @@ namespace Nektar
* @return (A shared pointer to) the local expansion of the
* \f$n^{\mathrm{th}}\f$ element.
*/
inline StdRegions::StdExpansionSharedPtr& ExpList::GetExp(int n) const
inline LocalRegions::ExpansionSharedPtr& ExpList::GetExp(int n) const
{
return (*m_exp)[n];
}
......@@ -1796,7 +1796,7 @@ namespace Nektar
/**
* @return (A const shared pointer to) the local expansion vector #m_exp
*/
inline const boost::shared_ptr<StdRegions::StdExpansionVector> ExpList::GetExp(void) const
inline const boost::shared_ptr<LocalRegions::ExpansionVector> ExpList::GetExp(void) const
{
return m_exp;
}
......
......@@ -89,7 +89,7 @@ namespace Nektar
const Array<OneD, const ExpListSharedPtr> &bndConstraint,
const Array<OneD, const SpatialDomains::BoundaryConditionShPtr>
&bndCond,
const StdRegions::StdExpansionVector &locexp,
const LocalRegions::ExpansionVector &locexp,
const SpatialDomains::MeshGraphSharedPtr &graph1D,
const map<int,int> &periodicVertices,
const bool DeclareCoeffPhysArrays)
......
......@@ -80,7 +80,7 @@ namespace Nektar
const Array<OneD,const ExpListSharedPtr> &bndConstraint,
const Array<OneD,const SpatialDomains
::BoundaryConditionShPtr> &bndCond,
const StdRegions::StdExpansionVector &locexp,
const LocalRegions::ExpansionVector &locexp,
const SpatialDomains::MeshGraphSharedPtr &graph1D,
const map<int,int> &periodicVertices,
const bool DeclareCoeffPhysArrays = true);
......
......@@ -408,7 +408,7 @@ namespace Nektar
const Array<OneD,const ExpListSharedPtr> &bndConstraint,
const Array<OneD, const SpatialDomains
::BoundaryConditionShPtr> &bndCond,
const StdRegions::StdExpansionVector &locexp,
const LocalRegions::ExpansionVector &locexp,
const SpatialDomains::MeshGraphSharedPtr &graph2D,
const PeriodicMap &periodicEdges,
const bool DeclareCoeffPhysArrays,
......@@ -789,7 +789,7 @@ namespace Nektar
void ExpList1D::SetUpPhysTangents(
const StdRegions::StdExpansionVector &locexp)
const LocalRegions::ExpansionVector &locexp)
{
map<int, int> EdgeGID;
int i,cnt,n,id;
......@@ -1032,7 +1032,7 @@ namespace Nektar
// }
void ExpList1D::v_SetUpPhysTangents(
const StdRegions::StdExpansionVector &locexp)
const LocalRegions::ExpansionVector &locexp)
{
SetUpPhysTangents(locexp);
}
......
......@@ -99,7 +99,7 @@ namespace Nektar
const Array<OneD,const ExpListSharedPtr> &bndConstraint,
const Array<OneD,const SpatialDomains
::BoundaryConditionShPtr> &bndCond,
const StdRegions::StdExpansionVector &locexp,
const LocalRegions::ExpansionVector &locexp,
const SpatialDomains::MeshGraphSharedPtr &graph2D,
const PeriodicMap &periodicEdges,
const bool DeclareCoeffPhysArrays = true,
......@@ -127,7 +127,7 @@ namespace Nektar
//MULTI_REGIONS_EXPORT void SetUpPhysNormals();//const StdRegions::StdExpansionVector &locexp);
/// Set up the tangents on each expansion
MULTI_REGIONS_EXPORT void SetUpPhysTangents(const StdRegions::StdExpansionVector &locexp);
MULTI_REGIONS_EXPORT void SetUpPhysTangents(const LocalRegions::ExpansionVector &locexp);
// direction = 1: Upwind
// direction = -1: Downwind
......@@ -169,7 +169,7 @@ namespace Nektar
/// Set up the tangents on each expansion.
virtual void v_SetUpPhysTangents(
const StdRegions::StdExpansionVector &locexp);
const LocalRegions::ExpansionVector &locexp);
virtual void v_WriteVtkPieceHeader(std::ofstream &outfile, int expansion);
......
......@@ -433,7 +433,7 @@ namespace Nektar
ExpList2D::ExpList2D(
const Array<OneD,const ExpListSharedPtr> &bndConstraint,
const Array<OneD,const SpatialDomains::BoundaryConditionShPtr> &bndCond,
const StdRegions::StdExpansionVector &locexp,
const LocalRegions::ExpansionVector &locexp,
const SpatialDomains::MeshGraphSharedPtr &graph3D,
const PeriodicMap &periodicFaces,
const bool DeclareCoeffPhysArrays,
......
......@@ -96,7 +96,7 @@ namespace Nektar
MULTI_REGIONS_EXPORT ExpList2D(
const Array<OneD,const ExpListSharedPtr> &bndConstraint,
const Array<OneD,const SpatialDomains::BoundaryConditionShPtr> &bndCond,
const StdRegions::StdExpansionVector &locexp,
const LocalRegions::ExpansionVector &locexp,
const SpatialDomains::MeshGraphSharedPtr &graph3D,
const PeriodicMap &periodicFaces,
const bool DeclareCoeffPhysArrays = true,
......
......@@ -53,7 +53,7 @@ namespace Nektar
const NekDouble lhom,
const bool useFFT,
const bool dealiasing,
const boost::shared_ptr<StdRegions::StdExpansionVector> &exp,
const boost::shared_ptr<LocalRegions::ExpansionVector> &exp,
const Array<OneD, ExpListSharedPtr> &planes):
ExpListHomogeneous1D(pSession,HomoBasis,lhom,useFFT,dealiasing)
{
......@@ -92,7 +92,7 @@ namespace Nektar
// note that nzplanes can be larger than nzmodes
m_planes[0] = plane_zero = MemoryManager<ExpList1D>::AllocateSharedPtr(m_session,graph1D, False);
m_exp = MemoryManager<StdRegions::StdExpansionVector>::AllocateSharedPtr();
m_exp = MemoryManager<LocalRegions::ExpansionVector>::AllocateSharedPtr();
nel = m_planes[0]->GetExpSize();
for(j = 0; j < nel; ++j)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment