LSProblem< dim > Class Template Reference

#include <ComputeCutCellMoments.H>

Public Member Functions
	~LSProblem ()
	LSProblem (const int &a_degreeP, const bool &a_useConstraints)
	LSProblem (const int &a_orderAccuracy, const int &a_degreeP, const bool &a_useConstraints, const IndexTM< Real, dim > &a_normal)
int	invertNormalEq (const Vector< Real > &a_rhs, Vector< Real > &a_residual)
int	getDegree ()
int	getOrderAccuracy ()
void	getMatrix (Real **a_matrix)
void	getRhs (Vector< Real > &a_rhs)
void	getUnknowns (Vector< Real > &a_unknown)
void	monoMaxMin (Real &a_maxVal, Real &a_minVal, const IndexTM< int, dim > &a_mono, const IFData< dim > &a_IFData)
void	computeBounds (const IndexTM< Real, dim > &a_dx, const CutCellMoments< dim > &a_ccm)
void	print (ostream &a_out) const
void	outputMatrix () const
void	outputRhs () const
void	outputUnknowns () const
void	outputBounds () const
const LocPthMoment &	getLocMonomialMapDegreeP () const
const PthMomentLoc &	getMonomialLocMapDegreeP () const
const PthMomentLoc &	getMonomialLocMapDegreePLess1 () const
Real	getUnknown (int loc)
int	getNumberDegP ()
int	getNumberDegPLess1 ()
int	numActiveBounds () const
void	setRhs (const Vector< Real > &a_rhs)
void	setMatrix ()
void	momentBounds (Real &a_lobnd, Real &a_hibnd, const IvDim &a_mono, const IFData< dim > &a_ifData)
void	differentiate (int &a_coeff, IvDim &a_Dmono, int &a_idir, const IvDim &a_mono)
int	nChooseR (int a_n, int a_r)
void	fillMap (PthMomentLoc &a_monoLoc, LocPthMoment &a_locMono, const int &a_degree)
int	numMonomials (const int &a_monoDegree)
int	factorial (const int &a_n, const int &a_m=0)
void	allocArray (const int &a_rows, const int &a_cols, Real **&a_A)
void	freeArray (const int &a_rows, const int &a_cols, Real **&a_A)

Public Attributes
int	m_order
int	m_degreeP
int	m_numActiveBounds
bool	m_useConstraints
IndexTM< Real, dim >	m_normal
PthMomentLoc	m_monoLocP
LocPthMoment	m_locMonoP
int	m_numP
PthMomentLoc	m_monoLocPLess1
LocPthMoment	m_locMonoPLess1
int	m_numPLess1
Real **	m_matrix
Vector< Real >	m_unknowns
Vector< Real >	m_rhs
Vector< Real >	m_lowerBound
Vector< Real >	m_upperBound

Private Types
typedef IndexTM< int, dim >	IvDim
typedef map< IvDim, int >	PthMomentLoc
typedef map< int, IvDim >	LocPthMoment

Member Typedef Documentation

IvDim

template<int dim>

typedef IndexTM<int,dim> LSProblem< dim >::IvDim

private

PthMomentLoc

template<int dim>

typedef map<IvDim,int > LSProblem< dim >::PthMomentLoc

private

LocPthMoment

template<int dim>

typedef map<int,IvDim> LSProblem< dim >::LocPthMoment

private

Constructor & Destructor Documentation

~LSProblem()

template<int dim>

LSProblem< dim >::~LSProblem

(

)

References dim.

LSProblem() [1/2]

template<int dim>

LSProblem< dim >::LSProblem	(	const int &	a_degreeP,
		const bool &	a_useConstraints
	)

References LSProblem< dim >::fillMap(), LSProblem< dim >::m_degreeP, LSProblem< dim >::m_locMonoP, LSProblem< dim >::m_locMonoPLess1, LSProblem< dim >::m_matrix, LSProblem< dim >::m_monoLocP, LSProblem< dim >::m_monoLocPLess1, LSProblem< dim >::m_numP, LSProblem< dim >::m_numPLess1, and LSProblem< dim >::numMonomials().

LSProblem() [2/2]

template<int dim>

LSProblem< dim >::LSProblem	(	const int &	a_orderAccuracy,
		const int &	a_degreeP,
		const bool &	a_useConstraints,
		const IndexTM< Real, dim > &	a_normal
	)

References dim, LSProblem< dim >::fillMap(), LSProblem< dim >::m_degreeP, LSProblem< dim >::m_locMonoP, LSProblem< dim >::m_locMonoPLess1, LSProblem< dim >::m_lowerBound, LSProblem< dim >::m_monoLocP, LSProblem< dim >::m_monoLocPLess1, LSProblem< dim >::m_numP, LSProblem< dim >::m_numPLess1, LSProblem< dim >::m_rhs, LSProblem< dim >::m_unknowns, LSProblem< dim >::m_upperBound, LSProblem< dim >::numMonomials(), Vector< T >::resize(), and LSProblem< dim >::setMatrix().

Member Function Documentation

invertNormalEq()

template<int dim>

int LSProblem< dim >::invertNormalEq	(	const Vector< Real > &	a_rhs,
		Vector< Real > &	a_residual
	)

References Abs(), dim, MayDay::Error(), ConstrainedLS::INCONSISTENT_BOUNDS, LSquares::LeastSquares(), LSProblem< dim >::m_lowerBound, LSProblem< dim >::m_matrix, LSProblem< dim >::m_numActiveBounds, LSProblem< dim >::m_numP, LSProblem< dim >::m_numPLess1, LSProblem< dim >::m_rhs, LSProblem< dim >::m_unknowns, LSProblem< dim >::m_upperBound, LSProblem< dim >::m_useConstraints, ConstrainedLS::numberActiveConstraints(), pout(), Vector< T >::resize(), ConstrainedLS::SINGULAR, ConstrainedLS::solveBoundConstrained(), ConstrainedLS::SUCCESS, ConstrainedLS::UNCONVERGED, and ConstrainedLS::UNDERDETERMINED.

Referenced by ComputeCutCellMoments< dim >::computeMoments(), and MinimalCCCM< dim >::computeMomentsRecursively().

getDegree()

template<int dim>

int LSProblem< dim >::getDegree ( )

inline

References LSProblem< dim >::m_degreeP.

getOrderAccuracy()

template<int dim>

int LSProblem< dim >::getOrderAccuracy ( )

inline

References LSProblem< dim >::computeBounds(), LSProblem< dim >::getMatrix(), LSProblem< dim >::getRhs(), LSProblem< dim >::getUnknowns(), LSProblem< dim >::m_order, LSProblem< dim >::monoMaxMin(), LSProblem< dim >::outputBounds(), LSProblem< dim >::outputMatrix(), LSProblem< dim >::outputRhs(), LSProblem< dim >::outputUnknowns(), and LSProblem< dim >::print().

getMatrix()

template<int dim>

void LSProblem< dim >::getMatrix

(

Real **

a_matrix

)

Referenced by LSProblem< dim >::getOrderAccuracy().

getRhs()

template<int dim>

void LSProblem< dim >::getRhs

(

Vector< Real > &

a_rhs

)

Referenced by LSProblem< dim >::getOrderAccuracy().

getUnknowns()

template<int dim>

void LSProblem< dim >::getUnknowns

(

Vector< Real > &

a_unknown

)

Referenced by LSProblem< dim >::getOrderAccuracy().

monoMaxMin()

template<int dim>

void LSProblem< dim >::monoMaxMin	(	Real &	a_maxVal,
		Real &	a_minVal,
		const IndexTM< int, dim > &	a_mono,
		const IFData< dim > &	a_IFData
	)

References dim, IFData< dim >::m_cellCenterCoord, IFData< dim >::m_cornerSigns, IFData< dim >::m_localCoord, and POW().

Referenced by LSProblem< dim >::getOrderAccuracy().

computeBounds()

template<int dim>

void LSProblem< dim >::computeBounds	(	const IndexTM< Real, dim > &	a_dx,
		const CutCellMoments< dim > &	a_ccm
	)

References dim, CutCellMoments< dim >::m_EBmoments, CutCellMoments< dim >::m_IFData, CutCellMoments< dim >::m_moments, Max(), and IndexTM< T, N >::sum().

Referenced by ComputeCutCellMoments< dim >::computeMoments(), and LSProblem< dim >::getOrderAccuracy().

print()

template<int dim>

void LSProblem< dim >::print

(

ostream &

a_out

)

const

References dim, LSProblem< dim >::m_degreeP, LSProblem< dim >::m_locMonoP, LSProblem< dim >::m_locMonoPLess1, LSProblem< dim >::m_monoLocP, LSProblem< dim >::m_monoLocPLess1, LSProblem< dim >::outputBounds(), LSProblem< dim >::outputMatrix(), LSProblem< dim >::outputRhs(), and LSProblem< dim >::outputUnknowns().

Referenced by ComputeCutCellMoments< dim >::computeMoments(), MinimalCCCM< dim >::computeMomentsRecursively(), LSProblem< dim >::getOrderAccuracy(), and operator<<().

outputMatrix()

template<int dim>

void LSProblem< dim >::outputMatrix

(

)

const

References dim, LSProblem< dim >::m_locMonoP, LSProblem< dim >::m_locMonoPLess1, LSProblem< dim >::m_matrix, and pout().

Referenced by LSProblem< dim >::getOrderAccuracy(), and LSProblem< dim >::print().

outputRhs()

template<int dim>

void LSProblem< dim >::outputRhs

(

)

const

References LSProblem< dim >::m_rhs, pout(), and Vector< T >::size().

Referenced by LSProblem< dim >::getOrderAccuracy(), and LSProblem< dim >::print().

outputUnknowns()

template<int dim>

void LSProblem< dim >::outputUnknowns

(

)

const

References LSProblem< dim >::m_unknowns, pout(), and Vector< T >::size().

Referenced by LSProblem< dim >::getOrderAccuracy(), and LSProblem< dim >::print().

outputBounds()

template<int dim>

void LSProblem< dim >::outputBounds

(

)

const

References LSProblem< dim >::m_lowerBound, LSProblem< dim >::m_unknowns, LSProblem< dim >::m_upperBound, LSProblem< dim >::m_useConstraints, pout(), and Vector< T >::size().

Referenced by LSProblem< dim >::getOrderAccuracy(), and LSProblem< dim >::print().

getLocMonomialMapDegreeP()

template<int dim>

const LocPthMoment& LSProblem< dim >::getLocMonomialMapDegreeP ( ) const

inline

References LSProblem< dim >::m_locMonoP.

Referenced by MinimalCCCM< dim >::computeRhs(), and ComputeCutCellMoments< dim >::computeRhs().

getMonomialLocMapDegreeP()

template<int dim>

const PthMomentLoc& LSProblem< dim >::getMonomialLocMapDegreeP ( ) const

inline

References LSProblem< dim >::m_monoLocP.

Referenced by CutCellMoments< dim >::addBdMoments(), ComputeCutCellMoments< dim >::computeMoments(), and MinimalCCCM< dim >::computeMomentsRecursively().

getMonomialLocMapDegreePLess1()

template<int dim>

const PthMomentLoc& LSProblem< dim >::getMonomialLocMapDegreePLess1 ( ) const

inline

References LSProblem< dim >::m_monoLocPLess1.

Referenced by CutCellMoments< dim >::addBdMoments(), ComputeCutCellMoments< dim >::computeMoments(), and MinimalCCCM< dim >::computeMomentsRecursively().

getUnknown()

template<int dim>

Real LSProblem< dim >::getUnknown ( int  loc )

inline

References LSProblem< dim >::m_unknowns.

Referenced by ComputeCutCellMoments< dim >::computeMoments(), and MinimalCCCM< dim >::computeMomentsRecursively().

getNumberDegP()

template<int dim>

int LSProblem< dim >::getNumberDegP ( )

inline

References LSProblem< dim >::m_numP.

Referenced by ComputeCutCellMoments< dim >::computeMoments(), MinimalCCCM< dim >::computeMomentsRecursively(), MinimalCCCM< dim >::computeRhs(), and ComputeCutCellMoments< dim >::computeRhs().

getNumberDegPLess1()

template<int dim>

int LSProblem< dim >::getNumberDegPLess1 ( )

inline

References LSProblem< dim >::m_numPLess1.

numActiveBounds()

template<int dim>

int LSProblem< dim >::numActiveBounds ( ) const

inline

References LSProblem< dim >::allocArray(), LSProblem< dim >::differentiate(), LSProblem< dim >::factorial(), LSProblem< dim >::fillMap(), LSProblem< dim >::freeArray(), LSProblem< dim >::m_numActiveBounds, LSProblem< dim >::momentBounds(), LSProblem< dim >::nChooseR(), LSProblem< dim >::numMonomials(), LSProblem< dim >::setMatrix(), and LSProblem< dim >::setRhs().

Referenced by ComputeCutCellMoments< dim >::computeMoments().

setRhs()

template<int dim>

void LSProblem< dim >::setRhs

(

const Vector< Real > &

a_rhs

)

Referenced by LSProblem< dim >::numActiveBounds().

setMatrix()

template<int dim>

void LSProblem< dim >::setMatrix

(

)

References MayDay::Abort(), LSProblem< dim >::allocArray(), LSProblem< dim >::differentiate(), dim, LARGEINTVAL, LSProblem< dim >::m_matrix, LSProblem< dim >::m_monoLocP, LSProblem< dim >::m_monoLocPLess1, LSProblem< dim >::m_normal, LSProblem< dim >::m_numP, and LSProblem< dim >::m_numPLess1.

Referenced by LSProblem< dim >::LSProblem(), and LSProblem< dim >::numActiveBounds().

momentBounds()

template<int dim>

void LSProblem< dim >::momentBounds	(	Real &	a_lobnd,
		Real &	a_hibnd,
		const IvDim &	a_mono,
		const IFData< dim > &	a_ifData
	)

References dim, IFData< dim >::m_cellCenterCoord, IFData< dim >::m_cornerSigns, IFData< dim >::m_localCoord, and POW().

Referenced by LSProblem< dim >::numActiveBounds().

differentiate()

template<int dim>

void LSProblem< dim >::differentiate	(	int &	a_coeff,
		IvDim &	a_Dmono,
		int &	a_idir,
		const IvDim &	a_mono
	)

References MayDay::Abort(), dim, and LARGEINTVAL.

Referenced by LSProblem< dim >::numActiveBounds(), and LSProblem< dim >::setMatrix().

nChooseR()

template<int dim>

int LSProblem< dim >::nChooseR	(	int	a_n,
		int	a_r
	)

Referenced by LSProblem< dim >::numActiveBounds().

fillMap()

template<int dim>

void LSProblem< dim >::fillMap	(	PthMomentLoc &	a_monoLoc,
		LocPthMoment &	a_locMono,
		const int &	a_degree
	)

References generateMultiIndices(), and Vector< T >::size().

Referenced by LSProblem< dim >::LSProblem(), and LSProblem< dim >::numActiveBounds().

numMonomials()

template<int dim>

int LSProblem< dim >::numMonomials

(

const int &

a_monoDegree

)

References dim, LSProblem< dim >::factorial(), and LARGEINTVAL.

Referenced by LSProblem< dim >::LSProblem(), and LSProblem< dim >::numActiveBounds().

factorial()

template<int dim>

int LSProblem< dim >::factorial	(	const int &	a_n,
		const int &	a_m = 0
	)

References MayDay::Abort().

Referenced by LSProblem< dim >::numActiveBounds(), and LSProblem< dim >::numMonomials().

allocArray()

template<int dim>

void LSProblem< dim >::allocArray	(	const int &	a_rows,
		const int &	a_cols,
		Real **&	a_A
	)

Referenced by LSProblem< dim >::numActiveBounds(), and LSProblem< dim >::setMatrix().

freeArray()

template<int dim>

void LSProblem< dim >::freeArray	(	const int &	a_rows,
		const int &	a_cols,
		Real **&	a_A
	)

Referenced by LSProblem< dim >::numActiveBounds().

Member Data Documentation

m_order

template<int dim>

int LSProblem< dim >::m_order

Referenced by LSProblem< dim >::getOrderAccuracy().

m_degreeP

template<int dim>

int LSProblem< dim >::m_degreeP

Referenced by LSProblem< dim >::getDegree(), LSProblem< dim >::LSProblem(), and LSProblem< dim >::print().

m_numActiveBounds

template<int dim>

int LSProblem< dim >::m_numActiveBounds

Referenced by LSProblem< dim >::invertNormalEq(), and LSProblem< dim >::numActiveBounds().

m_useConstraints

template<int dim>

bool LSProblem< dim >::m_useConstraints

Referenced by LSProblem< dim >::invertNormalEq(), and LSProblem< dim >::outputBounds().

m_normal

template<int dim>

IndexTM<Real,dim> LSProblem< dim >::m_normal

Referenced by LSProblem< dim >::setMatrix().

m_monoLocP

template<int dim>

PthMomentLoc LSProblem< dim >::m_monoLocP

Referenced by ComputeCutCellMoments< dim >::computeResiduals(), LSProblem< dim >::getMonomialLocMapDegreeP(), LSProblem< dim >::LSProblem(), LSProblem< dim >::print(), and LSProblem< dim >::setMatrix().

m_locMonoP

template<int dim>

LocPthMoment LSProblem< dim >::m_locMonoP

Referenced by LSProblem< dim >::getLocMonomialMapDegreeP(), LSProblem< dim >::LSProblem(), LSProblem< dim >::outputMatrix(), and LSProblem< dim >::print().

m_numP

template<int dim>

int LSProblem< dim >::m_numP

Referenced by ComputeCutCellMoments< dim >::computeResiduals(), LSProblem< dim >::getNumberDegP(), LSProblem< dim >::invertNormalEq(), LSProblem< dim >::LSProblem(), and LSProblem< dim >::setMatrix().

m_monoLocPLess1

template<int dim>

PthMomentLoc LSProblem< dim >::m_monoLocPLess1

Referenced by ComputeCutCellMoments< dim >::computeResiduals(), LSProblem< dim >::getMonomialLocMapDegreePLess1(), LSProblem< dim >::LSProblem(), LSProblem< dim >::print(), and LSProblem< dim >::setMatrix().

m_locMonoPLess1

template<int dim>

LocPthMoment LSProblem< dim >::m_locMonoPLess1

Referenced by LSProblem< dim >::LSProblem(), LSProblem< dim >::outputMatrix(), and LSProblem< dim >::print().

m_numPLess1

template<int dim>

int LSProblem< dim >::m_numPLess1

Referenced by ComputeCutCellMoments< dim >::computeResiduals(), LSProblem< dim >::getNumberDegPLess1(), LSProblem< dim >::invertNormalEq(), LSProblem< dim >::LSProblem(), and LSProblem< dim >::setMatrix().

m_matrix

template<int dim>

Real** LSProblem< dim >::m_matrix

Referenced by ComputeCutCellMoments< dim >::computeResiduals(), LSProblem< dim >::invertNormalEq(), LSProblem< dim >::LSProblem(), LSProblem< dim >::outputMatrix(), and LSProblem< dim >::setMatrix().

m_unknowns

template<int dim>

Vector<Real> LSProblem< dim >::m_unknowns

Referenced by LSProblem< dim >::getUnknown(), LSProblem< dim >::invertNormalEq(), LSProblem< dim >::LSProblem(), LSProblem< dim >::outputBounds(), and LSProblem< dim >::outputUnknowns().

m_rhs

template<int dim>

Vector<Real> LSProblem< dim >::m_rhs

Referenced by LSProblem< dim >::invertNormalEq(), LSProblem< dim >::LSProblem(), and LSProblem< dim >::outputRhs().

m_lowerBound

template<int dim>

Vector<Real> LSProblem< dim >::m_lowerBound

Referenced by LSProblem< dim >::invertNormalEq(), LSProblem< dim >::LSProblem(), and LSProblem< dim >::outputBounds().

m_upperBound

template<int dim>

Vector<Real> LSProblem< dim >::m_upperBound

Referenced by LSProblem< dim >::invertNormalEq(), LSProblem< dim >::LSProblem(), and LSProblem< dim >::outputBounds().

---

The documentation for this class was generated from the following files:

• ComputeCutCellMoments.H
• LSProblem.H
• LSProblemImplem.H