Nektar::LibUtilities::TimeIntegrationSolution Class Reference

#include <TimeIntegrationScheme.h>

Public Types
typedef Array< OneD, Array< OneD, Array< OneD, NekDouble > > >	TripleArray
typedef Array< OneD, Array< OneD, NekDouble > >	DoubleArray

Public Member Functions
	TimeIntegrationSolution (const TimeIntegrationSchemeKey &key, const DoubleArray &y, const NekDouble time, const NekDouble timestep)
	TimeIntegrationSolution (const TimeIntegrationSchemeKey &key, const TripleArray &y, const Array< OneD, NekDouble > &t)
	TimeIntegrationSolution (const TimeIntegrationSchemeKey &key, unsigned int nvar, unsigned int npoints)
	TimeIntegrationSolution (const TimeIntegrationSchemeKey &key)
const TimeIntegrationSchemeSharedPtr &	GetIntegrationScheme () const
const TimeIntegrationSchemeKey &	GetIntegrationSchemeKey () const
TimeIntegrationMethod	GetIntegrationMethod () const
const TripleArray &	GetSolutionVector () const
TripleArray &	UpdateSolutionVector ()
const DoubleArray &	GetSolution () const
DoubleArray &	UpdateSolution ()
const Array< OneD, const NekDouble > &	GetTimeVector () const
Array< OneD, NekDouble > &	UpdateTimeVector ()
NekDouble	GetTime () const
int	GetNsteps ()
int	GetFirstDim () const
int	GetSecondDim () const
unsigned int	GetNvalues (void) const
unsigned int	GetNderivs (void) const
const Array< OneD, const unsigned int > &	GetTimeLevelOffset ()
DoubleArray &	GetValue (const unsigned int timeLevelOffset)
DoubleArray &	GetDerivative (const unsigned int timeLevelOffset)
NekDouble	GetValueTime (const unsigned int timeLevelOffset)
void	SetValue (const unsigned int timeLevelOffset, const DoubleArray &y, const NekDouble t)
void	SetDerivative (const unsigned int timeLevelOffset, const DoubleArray &y, const NekDouble timestep)
void	SetSolVector (const int Offset, const DoubleArray &y)
void	RotateSolutionVector ()

Private Attributes
TimeIntegrationSchemeSharedPtr	m_scheme
TripleArray	m_solVector
Array< OneD, NekDouble >	m_t

Detailed Description

Definition at line 637 of file TimeIntegrationScheme.h.

Member Typedef Documentation

DoubleArray

typedef Array<OneD, Array<OneD, NekDouble> > Nektar::LibUtilities::TimeIntegrationSolution::DoubleArray

Definition at line 641 of file TimeIntegrationScheme.h.

TripleArray

typedef Array<OneD, Array<OneD, Array<OneD, NekDouble> > > Nektar::LibUtilities::TimeIntegrationSolution::TripleArray

Definition at line 640 of file TimeIntegrationScheme.h.

Constructor & Destructor Documentation

TimeIntegrationSolution() [1/4]

Nektar::LibUtilities::TimeIntegrationSolution::TimeIntegrationSolution	(	const TimeIntegrationSchemeKey &	key,
		const DoubleArray &	y,
		const NekDouble	time,
		const NekDouble	timestep
	)

Definition at line 74 of file TimeIntegrationScheme.cpp.

References GetTimeLevelOffset(), m_scheme, m_solVector, and m_t.

                                                                                  :
            m_scheme(TimeIntegrationSchemeManager()[key]),
            m_solVector(m_scheme->GetNsteps()),
            m_t(m_scheme->GetNsteps())
        {
            m_solVector[0] = y;
            m_t[0] = time;

            int nsteps         = m_scheme->GetNsteps();
            int nvar           = y.num_elements();
            int npoints        = y[0].num_elements();
            int nMultiStepVals = m_scheme->GetNmultiStepValues();
            const Array<OneD, const unsigned int>& timeLevels = GetTimeLevelOffset(); 
            for(int i = 1; i < nsteps; i++)
            {
                m_solVector[i] = Array<OneD, Array<OneD, NekDouble> >(nvar);
                for(int j = 0; j < nvar; j++)
                {
                    m_solVector[i][j] = Array<OneD,NekDouble>(npoints,0.0);
                }
                if(i < nMultiStepVals)
                {
                    m_t[i] = time - i*timestep*timeLevels[i];
                }
                else
                {
                    m_t[i] = timestep;
                }
            }
        }

TimeIntegrationSolution() [2/4]

Nektar::LibUtilities::TimeIntegrationSolution::TimeIntegrationSolution	(	const TimeIntegrationSchemeKey &	key,
		const TripleArray &	y,
		const Array< OneD, NekDouble > &	t
	)

Definition at line 108 of file TimeIntegrationScheme.cpp.

References ASSERTL1, and m_scheme.

                                                                                         :
            m_scheme(TimeIntegrationSchemeManager()[key]),
            m_solVector(y),
            m_t(t)
        {
            ASSERTL1(y.num_elements()==m_scheme->GetNsteps(),"Amount of Entries does not match number of (multi-) steps");
        }

TimeIntegrationSolution() [3/4]

Nektar::LibUtilities::TimeIntegrationSolution::TimeIntegrationSolution	(	const TimeIntegrationSchemeKey &	key,
		unsigned int	nvar,
		unsigned int	npoints
	)

Definition at line 118 of file TimeIntegrationScheme.cpp.

References m_scheme, and m_solVector.

                                                                              :
            m_scheme(TimeIntegrationSchemeManager()[key]),
            m_solVector(m_scheme->GetNsteps()),
            m_t(m_scheme->GetNsteps())
        {
            for(int i = 0; i < m_scheme->GetNsteps(); i++)
            {
                m_solVector[i] = Array<OneD, Array<OneD, NekDouble> >(nvar);
                for(int j = 0; j < nvar; j++)
                {
                    m_solVector[i][j] = Array<OneD,NekDouble>(npoints);
                }
            }         
        }

TimeIntegrationSolution() [4/4]

Nektar::LibUtilities::TimeIntegrationSolution::TimeIntegrationSolution

(

const TimeIntegrationSchemeKey &

key

)

Definition at line 135 of file TimeIntegrationScheme.cpp.

                                                                                           :
            m_scheme(TimeIntegrationSchemeManager()[key]),
            m_solVector(m_scheme->GetNsteps()),
            m_t(m_scheme->GetNsteps())
        {      
        }

Member Function Documentation

GetDerivative()

DoubleArray& Nektar::LibUtilities::TimeIntegrationSolution::GetDerivative ( const unsigned int  timeLevelOffset )

inline

Definition at line 763 of file TimeIntegrationScheme.h.

References ASSERTL1.

            {
                int nMultiStepVals = m_scheme->GetNmultiStepValues();
                int size           = m_scheme->GetNsteps();
                const Array<OneD, const unsigned int>& offsetvec = GetTimeLevelOffset();

                for(int i = nMultiStepVals; i < size; i++)
                {
                    if( timeLevelOffset == offsetvec[i] )
                    {
                        return m_solVector[i];
                    }
                }
                ASSERTL1(false,"The solution vector of this scheme does not contain a derivative at the requested time-level");
                return m_solVector[0];
            }

GetFirstDim()

int Nektar::LibUtilities::TimeIntegrationSolution::GetFirstDim ( ) const

inline

Definition at line 712 of file TimeIntegrationScheme.h.

            {
                return m_solVector[0].num_elements();
            }

GetIntegrationMethod()

TimeIntegrationMethod Nektar::LibUtilities::TimeIntegrationSolution::GetIntegrationMethod ( ) const

inline

Definition at line 667 of file TimeIntegrationScheme.h.

            {
                return m_scheme->GetIntegrationMethod();
            }

GetIntegrationScheme()

const TimeIntegrationSchemeSharedPtr& Nektar::LibUtilities::TimeIntegrationSolution::GetIntegrationScheme ( ) const

inline

Definition at line 657 of file TimeIntegrationScheme.h.

            {
                return m_scheme;
            }

GetIntegrationSchemeKey()

const TimeIntegrationSchemeKey& Nektar::LibUtilities::TimeIntegrationSolution::GetIntegrationSchemeKey ( ) const

inline

Definition at line 662 of file TimeIntegrationScheme.h.

            {
                return m_scheme->GetIntegrationSchemeKey();
            }

GetNderivs()

unsigned int Nektar::LibUtilities::TimeIntegrationSolution::GetNderivs ( void  ) const

inline

Definition at line 731 of file TimeIntegrationScheme.h.

            {
                return m_scheme->GetNmultiStepDerivs();
            }

GetNsteps()

int Nektar::LibUtilities::TimeIntegrationSolution::GetNsteps ( void  )

inline

Definition at line 707 of file TimeIntegrationScheme.h.

            {
                return m_scheme->GetNsteps();
            }

GetNvalues()

unsigned int Nektar::LibUtilities::TimeIntegrationSolution::GetNvalues ( void  ) const

inline

Definition at line 724 of file TimeIntegrationScheme.h.

            {
                return m_scheme->GetNmultiStepValues();
            }

GetSecondDim()

int Nektar::LibUtilities::TimeIntegrationSolution::GetSecondDim ( ) const

inline

Definition at line 717 of file TimeIntegrationScheme.h.

            {
                return m_solVector[0][0].num_elements();
            }

GetSolution()

const DoubleArray& Nektar::LibUtilities::TimeIntegrationSolution::GetSolution ( ) const

inline

Definition at line 682 of file TimeIntegrationScheme.h.

            {
                return m_solVector[0];
            }

GetSolutionVector()

const TripleArray& Nektar::LibUtilities::TimeIntegrationSolution::GetSolutionVector ( ) const

inline

Definition at line 672 of file TimeIntegrationScheme.h.

            {
                return m_solVector;
            }

GetTime()

NekDouble Nektar::LibUtilities::TimeIntegrationSolution::GetTime ( ) const

inline

Definition at line 702 of file TimeIntegrationScheme.h.

            {
                return m_t[0];
            }

GetTimeLevelOffset()

const Array<OneD, const unsigned int>& Nektar::LibUtilities::TimeIntegrationSolution::GetTimeLevelOffset ( )

inline

Definition at line 738 of file TimeIntegrationScheme.h.

Referenced by TimeIntegrationSolution().

            {
                return m_scheme->GetTimeLevelOffset();
            }

GetTimeVector()

const Array<OneD, const NekDouble>& Nektar::LibUtilities::TimeIntegrationSolution::GetTimeVector ( ) const

inline

Definition at line 692 of file TimeIntegrationScheme.h.

            {
                return m_t;
            }

GetValue()

DoubleArray& Nektar::LibUtilities::TimeIntegrationSolution::GetValue ( const unsigned int  timeLevelOffset )

inline

Definition at line 745 of file TimeIntegrationScheme.h.

References ASSERTL1.

            {
                int nMultiStepVals = m_scheme->GetNmultiStepValues();
                const Array<OneD, const unsigned int>& offsetvec = GetTimeLevelOffset();

                for(int i = 0; i < nMultiStepVals; i++)
                {
                    if( timeLevelOffset == offsetvec[i] )
                    {
                        return m_solVector[i];
                    }
                }
                ASSERTL1(false,"The solution vector of this scheme does not contain a value at the requested time-level");
                return m_solVector[0];
            }

GetValueTime()

NekDouble Nektar::LibUtilities::TimeIntegrationSolution::GetValueTime ( const unsigned int  timeLevelOffset )

inline

Definition at line 782 of file TimeIntegrationScheme.h.

References ASSERTL1.

            {
                int nMultiStepVals = m_scheme->GetNmultiStepValues();
                const Array<OneD, const unsigned int>& offsetvec = GetTimeLevelOffset();

                for(int i = 0; i < nMultiStepVals; i++)
                {
                    if( timeLevelOffset == offsetvec[i] )
                    {
                        return m_t[i];
                    }
                }
                ASSERTL1(false,"The solution vector of this scheme does not contain a value at the requested time-level");
                return m_t[0];
            }

RotateSolutionVector()

void Nektar::LibUtilities::TimeIntegrationSolution::RotateSolutionVector ( )

inline

Definition at line 845 of file TimeIntegrationScheme.h.

            {
                int i;
                int nMultiStepVals = m_scheme->GetNmultiStepValues();
                int size           = m_scheme->GetNsteps();
                for(i = (nMultiStepVals-1); i > 0; i--)
                {
                    m_solVector[i] = m_solVector[i-1];
                }

                for(i = (size-1); i > nMultiStepVals; i--)
                {
                    m_solVector[i] = m_solVector[i-1];
                }
            }

SetDerivative()

void Nektar::LibUtilities::TimeIntegrationSolution::SetDerivative ( const unsigned int  timeLevelOffset, const DoubleArray &  y, const NekDouble  timestep  )

inline

Definition at line 820 of file TimeIntegrationScheme.h.

            {
                int nMultiStepVals = m_scheme->GetNmultiStepValues();
                int size           = m_scheme->GetNsteps();
                const Array<OneD, const unsigned int>& offsetvec = GetTimeLevelOffset();

                for(int i = nMultiStepVals; i < size; i++)
                {
                    if( timeLevelOffset == offsetvec[i] )
                    {
                        m_solVector[i] = y;
                        m_t[i] = timestep;
                        return;
                    }
                }
            }

SetSolVector()

void Nektar::LibUtilities::TimeIntegrationSolution::SetSolVector ( const int  Offset, const DoubleArray &  y  )

inline

Definition at line 838 of file TimeIntegrationScheme.h.

            {
                m_solVector[Offset] = y;
            }

SetValue()

void Nektar::LibUtilities::TimeIntegrationSolution::SetValue ( const unsigned int  timeLevelOffset, const DoubleArray &  y, const NekDouble  t  )

inline

Definition at line 801 of file TimeIntegrationScheme.h.

            {
                int nMultiStepVals = m_scheme->GetNmultiStepValues();
                const Array<OneD, const unsigned int>& offsetvec = GetTimeLevelOffset();

                for(int i = 0; i < nMultiStepVals; i++)
                {
                    if( timeLevelOffset == offsetvec[i] )
                    {
                        m_solVector[i] = y;
                        m_t[i] = t;
                        return;
                    }
                }
            }

UpdateSolution()

DoubleArray& Nektar::LibUtilities::TimeIntegrationSolution::UpdateSolution ( )

inline

Definition at line 687 of file TimeIntegrationScheme.h.

            {
                return m_solVector[0];
            }

UpdateSolutionVector()

TripleArray& Nektar::LibUtilities::TimeIntegrationSolution::UpdateSolutionVector ( )

inline

Definition at line 677 of file TimeIntegrationScheme.h.

            {
                return m_solVector;
            }

UpdateTimeVector()

Array<OneD, NekDouble>& Nektar::LibUtilities::TimeIntegrationSolution::UpdateTimeVector ( )

inline

Definition at line 697 of file TimeIntegrationScheme.h.

            {
                return m_t;
            }

Member Data Documentation

m_scheme

TimeIntegrationSchemeSharedPtr Nektar::LibUtilities::TimeIntegrationSolution::m_scheme

private

Definition at line 863 of file TimeIntegrationScheme.h.

Referenced by TimeIntegrationSolution().

m_solVector

TripleArray Nektar::LibUtilities::TimeIntegrationSolution::m_solVector

private

Definition at line 864 of file TimeIntegrationScheme.h.

Referenced by TimeIntegrationSolution().

m_t

Array<OneD,NekDouble> Nektar::LibUtilities::TimeIntegrationSolution::m_t

private

Definition at line 865 of file TimeIntegrationScheme.h.

Referenced by TimeIntegrationSolution().