Nektar::Collections::CoalescedGeomData Class Reference

#include <CoalescedGeomData.h>

Collaboration diagram for Nektar::Collections::CoalescedGeomData:

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 54 of file CoalescedGeomData.h.

Constructor & Destructor Documentation

Nektar::Collections::CoalescedGeomData::CoalescedGeomData

(

void

)

Definition at line 44 of file CoalescedGeomData.cpp.

{
}

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

virtual

Definition at line 48 of file CoalescedGeomData.cpp.

{
}

Member Function Documentation

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

(

std::vector< StdRegions::StdExpansionSharedPtr > &

pColLExp

)

Definition at line 147 of file CoalescedGeomData.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::Collections::eDerivFactors, Vmath::Fill(), npts, 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 Array<TwoD, const NekDouble> Dfac =
                    pCollExp[i]->GetMetricInfo()->GetDerivFactors(ptsKeys);

            if (pCollExp[i]->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];
}

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

(

std::vector< StdRegions::StdExpansionSharedPtr > &

pColLExp

)

Definition at line 52 of file CoalescedGeomData.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::Collections::eJac, Vmath::Fill(), npts, 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 Array<OneD, const NekDouble> jac =
                                pCollExp[i]->GetMetricInfo()->GetJac(ptsKeys);

            if (pCollExp[i]->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];
}

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

(

std::vector< StdRegions::StdExpansionSharedPtr > &

pColLExp

)

Definition at line 100 of file CoalescedGeomData.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::Collections::eJacWithStdWeights, Vmath::Fill(), npts, 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 Array<OneD, const NekDouble> jac =
                            pCollExp[i]->GetMetricInfo()->GetJac(ptsKeys);

            if (pCollExp[i]->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

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

private

Definition at line 71 of file CoalescedGeomData.h.

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

private

Definition at line 72 of file CoalescedGeomData.h.