Nektar::Collections::CoalescedGeomData Class Reference

#include <CoalescedGeomData.h>

Public Member Functions
	CoalescedGeomData (void)
virtual	~CoalescedGeomData (void)
const Array< OneD, const NekDouble > &	GetJac (std::vector< StdRegions::StdExpansionSharedPtr > &pColLExp)
const Array< OneD, const NekDouble > &	GetJacWithStdWeights (std::vector< StdRegions::StdExpansionSharedPtr > &pColLExp)
const Array< TwoD, const NekDouble > &	GetDerivFactors (std::vector< StdRegions::StdExpansionSharedPtr > &pColLExp)

Private Attributes
std::map< GeomData, Array< OneD, NekDouble > >	m_oneDGeomData
std::map< GeomData, Array< TwoD, NekDouble > >	m_twoDGeomData

Detailed Description

Definition at line 53 of file CoalescedGeomData.h.

Constructor & Destructor Documentation

CoalescedGeomData()

Nektar::Collections::CoalescedGeomData::CoalescedGeomData

(

void

)

Definition at line 45 of file CoalescedGeomData.cpp.

{
}

~CoalescedGeomData()

Nektar::Collections::CoalescedGeomData::~CoalescedGeomData ( void  )

virtual

Definition at line 49 of file CoalescedGeomData.cpp.

{
}

Member Function Documentation

GetDerivFactors()

const Array< TwoD, const NekDouble > & Nektar::Collections::CoalescedGeomData::GetDerivFactors

(

std::vector< StdRegions::StdExpansionSharedPtr > &

pColLExp

)

Definition at line 150 of file CoalescedGeomData.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::Collections::eDerivFactors, Vmath::Fill(), Nektar::LocalRegions::Expansion::GetMetricInfo(), and Vmath::Vcopy().

{
    if(m_twoDGeomData.count(eDerivFactors) == 0)
    {
        LibUtilities::PointsKeyVector ptsKeys = pCollExp[0]->GetPointsKeys();

        int nElmts = pCollExp.size();
        const int coordim = pCollExp[0]->GetCoordim();
        int dim = ptsKeys.size();

        // set up Cached Jacobians to be continuous
        int npts = 1;
        for (int i = 0; i < dim; ++i)
        {
            npts   *= ptsKeys[i].GetNumPoints();
        }


        Array<TwoD, NekDouble> newDFac(dim*coordim,npts*nElmts);

        //copy Jacobians into a continuous list and set new chatched value
        int cnt = 0;
        for(int i = 0; i < nElmts; ++i)
        {
            const StdRegions::StdExpansion * sep = &(*pCollExp[i]);
            const LocalRegions::Expansion  * lep = dynamic_cast<const LocalRegions::Expansion*>( sep );

            const Array<TwoD, const NekDouble> Dfac = lep->GetMetricInfo()->GetDerivFactors( ptsKeys );

            if( lep->GetMetricInfo()->GetGtype() == SpatialDomains::eDeformed)
            {
                for (int j = 0; j < dim*coordim; ++j)
                {
                    Vmath::Vcopy(npts, &Dfac[j][0], 1, &newDFac[j][cnt], 1);
                }
            }
            else
            {
                for (int j = 0; j < dim*coordim; ++j)
                {
                    Vmath::Fill(npts, Dfac[j][0], &newDFac[j][cnt], 1);
                }
            }
            cnt += npts;
        }

        m_twoDGeomData[eDerivFactors] = newDFac;
    }

    return m_twoDGeomData[eDerivFactors];
}

GetJac()

const Array< OneD, const NekDouble > & Nektar::Collections::CoalescedGeomData::GetJac

(

std::vector< StdRegions::StdExpansionSharedPtr > &

pColLExp

)

Definition at line 53 of file CoalescedGeomData.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::Collections::eJac, Vmath::Fill(), Nektar::LocalRegions::Expansion::GetMetricInfo(), and Vmath::Vcopy().

{

    if(m_oneDGeomData.count(eJac) == 0)
    {

        LibUtilities::PointsKeyVector ptsKeys = pCollExp[0]->GetPointsKeys();
        int nElmts = pCollExp.size();

        // set up Cached Jacobians to be continuous
        int npts = 1;
        for (int i = 0; i < ptsKeys.size(); ++i)
        {
            npts   *= ptsKeys[i].GetNumPoints();
        }


        Array<OneD, NekDouble> newjac(npts*nElmts);

        //copy Jacobians into a continuous list and set new chatched value
        int cnt = 0;
        for(int i = 0; i < nElmts; ++i)
        {
            const StdRegions::StdExpansion * sep = &(*pCollExp[i]);
            const LocalRegions::Expansion  * lep = dynamic_cast<const LocalRegions::Expansion*>( sep );
          
            const Array<OneD, const NekDouble> jac = lep->GetMetricInfo()->GetJac( ptsKeys );

            if( lep->GetMetricInfo()->GetGtype() == SpatialDomains::eDeformed )
            {
                Vmath::Vcopy(npts, &jac[0], 1, &newjac[cnt], 1);
            }
            else
            {
                Vmath::Fill(npts, jac[0], &newjac[cnt], 1);
            }

            cnt += npts;
        }

        m_oneDGeomData[eJac] = newjac;
    }

    return m_oneDGeomData[eJac];
}

GetJacWithStdWeights()

const Array< OneD, const NekDouble > & Nektar::Collections::CoalescedGeomData::GetJacWithStdWeights

(

std::vector< StdRegions::StdExpansionSharedPtr > &

pColLExp

)

Definition at line 102 of file CoalescedGeomData.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::Collections::eJacWithStdWeights, Vmath::Fill(), Nektar::LocalRegions::Expansion::GetMetricInfo(), and Vmath::Vcopy().

{
    if(m_oneDGeomData.count(eJacWithStdWeights) == 0)
    {
        LibUtilities::PointsKeyVector ptsKeys = pCollExp[0]->GetPointsKeys();
        int nElmts = pCollExp.size();

        // set up Cached Jacobians to be continuous
        int npts = 1;
        for (int i = 0; i < ptsKeys.size(); ++i)
        {
            npts   *= ptsKeys[i].GetNumPoints();
        }


        Array<OneD, NekDouble> newjac(npts*nElmts), tmp;

        //copy Jacobians into a continuous list and set new chatched value
        int cnt = 0;
        for(int i = 0; i < nElmts; ++i)
        {
            const StdRegions::StdExpansion * sep = &(*pCollExp[i]);
            const LocalRegions::Expansion  * lep = dynamic_cast<const LocalRegions::Expansion*>( sep );

            const Array<OneD, const NekDouble> jac = lep->GetMetricInfo()->GetJac(ptsKeys);

            if( lep->GetMetricInfo()->GetGtype() == SpatialDomains::eDeformed )
            {
                Vmath::Vcopy(npts, &jac[0], 1, &newjac[cnt], 1);
            }
            else
            {
                Vmath::Fill(npts, jac[0], &newjac[cnt], 1);
            }

            pCollExp[0]->MultiplyByStdQuadratureMetric(newjac + cnt,
                                                       tmp = newjac + cnt);
            cnt += npts;
        }

        m_oneDGeomData[eJacWithStdWeights] = newjac;
    }

    return m_oneDGeomData[eJacWithStdWeights];
}

Member Data Documentation

m_oneDGeomData

std::map<GeomData,Array<OneD, NekDouble> > Nektar::Collections::CoalescedGeomData::m_oneDGeomData

private

Definition at line 70 of file CoalescedGeomData.h.

m_twoDGeomData

std::map<GeomData,Array<TwoD, NekDouble> > Nektar::Collections::CoalescedGeomData::m_twoDGeomData

private

Definition at line 71 of file CoalescedGeomData.h.