Nektar::MultiRegions::Preconditioner Class Reference

#include <Preconditioner.h>

Inheritance diagram for Nektar::MultiRegions::Preconditioner:

Public Member Functions
	Preconditioner (const std::shared_ptr< GlobalLinSys > &plinsys, const AssemblyMapSharedPtr &pLocToGloMap)
virtual	~Preconditioner ()
void	DoPreconditioner (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput)
void	DoPreconditionerWithNonVertOutput (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const Array< OneD, NekDouble > &pNonVertOutput, Array< OneD, NekDouble > &pVertForce=NullNekDouble1DArray)
void	DoTransformBasisToLowEnergy (Array< OneD, NekDouble > &pInOut)
void	DoTransformCoeffsFromLowEnergy (Array< OneD, NekDouble > &pInOut)
void	DoTransformCoeffsToLowEnergy (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput)
void	DoTransformBasisFromLowEnergy (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput)
void	BuildPreconditioner ()
void	InitObject ()
Array< OneD, NekDouble >	AssembleStaticCondGlobalDiagonals ()
	Performs global assembly of diagonal entries to global Schur complement matrix. More...
const DNekScalBlkMatSharedPtr &	GetBlockTransformedSchurCompl () const
const DNekScalBlkMatSharedPtr &	GetBlockCMatrix () const
const DNekScalBlkMatSharedPtr &	GetBlockInvDMatrix () const
const DNekScalBlkMatSharedPtr &	GetBlockSchurCompl () const
const DNekScalBlkMatSharedPtr &	GetBlockTransformationMatrix () const
const DNekScalBlkMatSharedPtr &	GetBlockTransposedTransformationMatrix () const
DNekScalMatSharedPtr	TransformedSchurCompl (int offset, int bndoffset, const std::shared_ptr< DNekScalMat > &loc_mat)

Protected Member Functions
virtual DNekScalMatSharedPtr	v_TransformedSchurCompl (int offset, int bndoffset, const std::shared_ptr< DNekScalMat > &loc_mat)
	Get block elemental transposed transformation matrix \(\mathbf{R}^{T}\). More...
virtual void	v_InitObject ()
virtual void	v_DoPreconditioner (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput)
	Apply a preconditioner to the conjugate gradient method. More...
virtual void	v_DoPreconditionerWithNonVertOutput (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const Array< OneD, NekDouble > &pNonVertOutput, Array< OneD, NekDouble > &pVertForce)
	Apply a preconditioner to the conjugate gradient method with an output for non-vertex degrees of freedom. More...
virtual void	v_DoTransformBasisToLowEnergy (Array< OneD, NekDouble > &pInOut)
	Transform from original basis to low energy basis. More...
virtual void	v_DoTransformCoeffsFromLowEnergy (Array< OneD, NekDouble > &pInOut)
	Transform from low energy coeffs to orignal basis. More...
virtual void	v_DoTransformCoeffsToLowEnergy (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput)
	Multiply by the block inverse transformation matrix. More...
virtual void	v_DoTransformBasisFromLowEnergy (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput)
	Multiply by the block transposed inverse transformation matrix. More...
virtual void	v_BuildPreconditioner ()

Protected Attributes
const std::weak_ptr< GlobalLinSys >	m_linsys
PreconditionerType	m_preconType
DNekMatSharedPtr	m_preconditioner
std::weak_ptr< AssemblyMap >	m_locToGloMap
LibUtilities::CommSharedPtr	m_comm

Private Member Functions
void	NullPreconditioner (void)

Static Private Attributes
static std::string	lookupIds []
static std::string	def

Detailed Description

This class implements preconditioning for the conjugate gradient matrix solver.

This class implements diagonal preconditioning for the conjugate gradient matrix solver.

Definition at line 65 of file Preconditioner.h.

Constructor & Destructor Documentation

Preconditioner()

Nektar::MultiRegions::Preconditioner::Preconditioner	(	const std::shared_ptr< GlobalLinSys > &	plinsys,
		const AssemblyMapSharedPtr &	pLocToGloMap
	)

Definition at line 76 of file Preconditioner.cpp.

    : m_linsys(plinsys), m_preconType(pLocToGloMap->GetPreconType()),
      m_locToGloMap(pLocToGloMap)
{
}

~Preconditioner()

virtual Nektar::MultiRegions::Preconditioner::~Preconditioner ( )

inlinevirtual

Definition at line 73 of file Preconditioner.h.

    {
    }

Member Function Documentation

AssembleStaticCondGlobalDiagonals()

Array< OneD, NekDouble > Nektar::MultiRegions::Preconditioner::AssembleStaticCondGlobalDiagonals

(

)

Performs global assembly of diagonal entries to global Schur complement matrix.

Definition at line 176 of file Preconditioner.cpp.

{
    auto asmMap = m_locToGloMap.lock();
 
    int nGlobalBnd = asmMap->GetNumGlobalBndCoeffs();
    int nDirBnd    = asmMap->GetNumGlobalDirBndCoeffs();
    int rows       = nGlobalBnd - nDirBnd;
 
    DNekScalBlkMatSharedPtr loc_mat;
    DNekScalMatSharedPtr bnd_mat;
    int sign1, sign2, gid1, gid2, i, j, n, cnt;
    Array<OneD, NekDouble> diagonals(rows, 0.0);
 
    // Extract diagonal contributions of globally assembled
    // schur complement matrix
    for (cnt = n = 0; n < m_linsys.lock()->GetNumBlocks(); ++n)
    {
        // Get statically condensed local matrix.
        loc_mat = (m_linsys.lock())->GetStaticCondBlock(n);
 
        // Extract boundary block.
        bnd_mat = loc_mat->GetBlock(0, 0);
 
        // Offset by number of rows.
        int bnd_row = bnd_mat->GetRows();
 
        for (i = 0; i < bnd_row; ++i)
        {
            gid1  = asmMap->GetLocalToGlobalBndMap(cnt + i) - nDirBnd;
            sign1 = asmMap->GetLocalToGlobalBndSign(cnt + i);
 
            if (gid1 < 0)
            {
                continue;
            }
 
            for (j = 0; j < bnd_row; ++j)
            {
                gid2  = asmMap->GetLocalToGlobalBndMap(cnt + j) - nDirBnd;
                sign2 = asmMap->GetLocalToGlobalBndSign(cnt + j);
 
                if (gid2 == gid1)
                {
                    diagonals[gid1] += sign1 * sign2 * (*bnd_mat)(i, j);
                }
            }
        }
        cnt += bnd_row;
    }
 
    return diagonals;
}

References m_linsys, and m_locToGloMap.

Referenced by Nektar::MultiRegions::PreconditionerDiagonal::StaticCondDiagonalPreconditionerSum().

BuildPreconditioner()

void Nektar::MultiRegions::Preconditioner::BuildPreconditioner ( )

inline

Definition at line 235 of file Preconditioner.h.

{
    v_BuildPreconditioner();
}

References v_BuildPreconditioner().

DoPreconditioner()

void Nektar::MultiRegions::Preconditioner::DoPreconditioner ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput  )

inline

Definition at line 178 of file Preconditioner.h.

{
    v_DoPreconditioner(pInput, pOutput);
}

References v_DoPreconditioner().

DoPreconditionerWithNonVertOutput()

void Nektar::MultiRegions::Preconditioner::DoPreconditionerWithNonVertOutput ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput, const Array< OneD, NekDouble > &  pNonVertOutput, Array< OneD, NekDouble > &  pVertForce = NullNekDouble1DArray  )

inline

Definition at line 187 of file Preconditioner.h.

{
    v_DoPreconditionerWithNonVertOutput(pInput, pOutput, pNonVertOutput,
                                        pVertForce);
}

References v_DoPreconditionerWithNonVertOutput().

DoTransformBasisFromLowEnergy()

void Nektar::MultiRegions::Preconditioner::DoTransformBasisFromLowEnergy ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput  )

inline

Definition at line 226 of file Preconditioner.h.

{
    v_DoTransformBasisFromLowEnergy(pInput, pOutput);
}

References v_DoTransformBasisFromLowEnergy().

DoTransformBasisToLowEnergy()

void Nektar::MultiRegions::Preconditioner::DoTransformBasisToLowEnergy ( Array< OneD, NekDouble > &  pInOut )

inline

Definition at line 199 of file Preconditioner.h.

{
    v_DoTransformBasisToLowEnergy(pInOut);
}

References v_DoTransformBasisToLowEnergy().

DoTransformCoeffsFromLowEnergy()

void Nektar::MultiRegions::Preconditioner::DoTransformCoeffsFromLowEnergy ( Array< OneD, NekDouble > &  pInOut )

inline

Definition at line 208 of file Preconditioner.h.

{
    v_DoTransformCoeffsFromLowEnergy(pInput);
}

References v_DoTransformCoeffsFromLowEnergy().

DoTransformCoeffsToLowEnergy()

void Nektar::MultiRegions::Preconditioner::DoTransformCoeffsToLowEnergy ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput  )

inline

Definition at line 217 of file Preconditioner.h.

{
    v_DoTransformCoeffsToLowEnergy(pInput, pOutput);
}

References v_DoTransformCoeffsToLowEnergy().

GetBlockCMatrix()

const DNekScalBlkMatSharedPtr& Nektar::MultiRegions::Preconditioner::GetBlockCMatrix ( ) const

inline

GetBlockInvDMatrix()

const DNekScalBlkMatSharedPtr& Nektar::MultiRegions::Preconditioner::GetBlockInvDMatrix ( ) const

inline

GetBlockSchurCompl()

const DNekScalBlkMatSharedPtr& Nektar::MultiRegions::Preconditioner::GetBlockSchurCompl ( ) const

inline

GetBlockTransformationMatrix()

const DNekScalBlkMatSharedPtr& Nektar::MultiRegions::Preconditioner::GetBlockTransformationMatrix ( ) const

inline

GetBlockTransformedSchurCompl()

const DNekScalBlkMatSharedPtr& Nektar::MultiRegions::Preconditioner::GetBlockTransformedSchurCompl ( ) const

inline

GetBlockTransposedTransformationMatrix()

const DNekScalBlkMatSharedPtr& Nektar::MultiRegions::Preconditioner::GetBlockTransposedTransformationMatrix ( ) const

inline

InitObject()

void Nektar::MultiRegions::Preconditioner::InitObject ( )

inline

Definition at line 161 of file Preconditioner.h.

{
    v_InitObject();
}

References v_InitObject().

NullPreconditioner()

void Nektar::MultiRegions::Preconditioner::NullPreconditioner ( void  )

private

TransformedSchurCompl()

DNekScalMatSharedPtr Nektar::MultiRegions::Preconditioner::TransformedSchurCompl ( int  offset, int  bndoffset, const std::shared_ptr< DNekScalMat > &  loc_mat  )

inline

Definition at line 169 of file Preconditioner.h.

{
    return v_TransformedSchurCompl(offset, bndoffset, loc_mat);
}

References v_TransformedSchurCompl().

v_BuildPreconditioner()

void Nektar::MultiRegions::Preconditioner::v_BuildPreconditioner ( )

protectedvirtual

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, Nektar::MultiRegions::PreconditionerLinearWithLowEnergy, Nektar::MultiRegions::PreconditionerLinearWithDiag, Nektar::MultiRegions::PreconditionerLinearWithBlock, Nektar::MultiRegions::PreconditionerLinear, Nektar::MultiRegions::PreconditionerNull, Nektar::MultiRegions::PreconditionerDiagonal, and Nektar::MultiRegions::PreconditionerBlock.

Definition at line 157 of file Preconditioner.cpp.

{
}

Referenced by BuildPreconditioner().

v_DoPreconditioner()

void Nektar::MultiRegions::Preconditioner::v_DoPreconditioner ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput  )

protectedvirtual

Apply a preconditioner to the conjugate gradient method.

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, Nektar::MultiRegions::PreconditionerLinearWithLowEnergy, Nektar::MultiRegions::PreconditionerLinearWithDiag, Nektar::MultiRegions::PreconditionerLinearWithBlock, Nektar::MultiRegions::PreconditionerLinear, Nektar::MultiRegions::PreconditionerNull, Nektar::MultiRegions::PreconditionerDiagonal, and Nektar::MultiRegions::PreconditionerBlock.

Definition at line 100 of file Preconditioner.cpp.

{
    boost::ignore_unused(pInput, pOutput);
    NEKERROR(ErrorUtil::efatal, "Method does not exist");
}

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by DoPreconditioner().

v_DoPreconditionerWithNonVertOutput()

void Nektar::MultiRegions::Preconditioner::v_DoPreconditionerWithNonVertOutput ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput, const Array< OneD, NekDouble > &  pNonVertOutput, Array< OneD, NekDouble > &  pVertForce  )

protectedvirtual

Apply a preconditioner to the conjugate gradient method with an output for non-vertex degrees of freedom.

Reimplemented in Nektar::MultiRegions::PreconditionerLinear.

Definition at line 111 of file Preconditioner.cpp.

{
    boost::ignore_unused(pInput, pOutput, pNonVertOutput, pVertForce);
    NEKERROR(ErrorUtil::efatal, "Method does not exist");
}

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by DoPreconditionerWithNonVertOutput().

v_DoTransformBasisFromLowEnergy()

void Nektar::MultiRegions::Preconditioner::v_DoTransformBasisFromLowEnergy ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput  )

protectedvirtual

Multiply by the block transposed inverse transformation matrix.

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, and Nektar::MultiRegions::PreconditionerLinearWithLowEnergy.

Definition at line 150 of file Preconditioner.cpp.

{
    boost::ignore_unused(pInput, pOutput);
    NEKERROR(ErrorUtil::efatal, "Method does not exist");
}

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by DoTransformBasisFromLowEnergy().

v_DoTransformBasisToLowEnergy()

void Nektar::MultiRegions::Preconditioner::v_DoTransformBasisToLowEnergy ( Array< OneD, NekDouble > &  pInOut )

protectedvirtual

Transform from original basis to low energy basis.

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, and Nektar::MultiRegions::PreconditionerLinearWithLowEnergy.

Definition at line 123 of file Preconditioner.cpp.

{
    boost::ignore_unused(pInOut);
}

Referenced by DoTransformBasisToLowEnergy().

v_DoTransformCoeffsFromLowEnergy()

void Nektar::MultiRegions::Preconditioner::v_DoTransformCoeffsFromLowEnergy ( Array< OneD, NekDouble > &  pInOut )

protectedvirtual

Transform from low energy coeffs to orignal basis.

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, and Nektar::MultiRegions::PreconditionerLinearWithLowEnergy.

Definition at line 132 of file Preconditioner.cpp.

{
    boost::ignore_unused(pInOut);
}

Referenced by DoTransformCoeffsFromLowEnergy().

v_DoTransformCoeffsToLowEnergy()

void Nektar::MultiRegions::Preconditioner::v_DoTransformCoeffsToLowEnergy ( const Array< OneD, NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput  )

protectedvirtual

Multiply by the block inverse transformation matrix.

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, and Nektar::MultiRegions::PreconditionerLinearWithLowEnergy.

Definition at line 141 of file Preconditioner.cpp.

{
    boost::ignore_unused(pInput, pOutput);
}

Referenced by DoTransformCoeffsToLowEnergy().

v_InitObject()

void Nektar::MultiRegions::Preconditioner::v_InitObject ( )

protectedvirtual

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, Nektar::MultiRegions::PreconditionerLinearWithLowEnergy, Nektar::MultiRegions::PreconditionerLinearWithDiag, Nektar::MultiRegions::PreconditionerLinearWithBlock, Nektar::MultiRegions::PreconditionerLinear, Nektar::MultiRegions::PreconditionerNull, Nektar::MultiRegions::PreconditionerDiagonal, and Nektar::MultiRegions::PreconditionerBlock.

Definition at line 92 of file Preconditioner.cpp.

{
    NEKERROR(ErrorUtil::efatal, "Method does not exist");
}

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by InitObject().

v_TransformedSchurCompl()

DNekScalMatSharedPtr Nektar::MultiRegions::Preconditioner::v_TransformedSchurCompl ( int  offset, int  bndoffset, const std::shared_ptr< DNekScalMat > &  loc_mat  )

protectedvirtual

Get block elemental transposed transformation matrix \(\mathbf{R}^{T}\).

Reimplemented in Nektar::MultiRegions::PreconditionerLowEnergy, and Nektar::MultiRegions::PreconditionerLinearWithLowEnergy.

Definition at line 165 of file Preconditioner.cpp.

{
    boost::ignore_unused(offset, bnd_offset);
    return loc_mat;
}

Referenced by TransformedSchurCompl().

Member Data Documentation

def

std::string Nektar::MultiRegions::Preconditioner::def

staticprivate

Initial value:

=
    LibUtilities::SessionReader::RegisterDefaultSolverInfo("Preconditioner",
                                                           "Diagonal")

Definition at line 154 of file Preconditioner.h.

lookupIds

std::string Nektar::MultiRegions::Preconditioner::lookupIds

staticprivate

Initial value:

= {
    LibUtilities::SessionReader::RegisterEnumValue("Preconditioner", "Null",
                                                   eNull),
    LibUtilities::SessionReader::RegisterEnumValue("Preconditioner", "Diagonal",
                                                   eDiagonal),
    LibUtilities::SessionReader::RegisterEnumValue(
        "Preconditioner", "FullLinearSpaceWithDiagonal", eLinearWithDiagonal),
    LibUtilities::SessionReader::RegisterEnumValue("Preconditioner",
                                                   "FullLinearSpace", eLinear),
    LibUtilities::SessionReader::RegisterEnumValue(
        "Preconditioner", "LowEnergyBlock", eLowEnergy),
    LibUtilities::SessionReader::RegisterEnumValue(
        "Preconditioner", "FullLinearSpaceWithLowEnergyBlock",
        eLinearWithLowEnergy),
    LibUtilities::SessionReader::RegisterEnumValue("Preconditioner", "Block",
                                                   eBlock),
    LibUtilities::SessionReader::RegisterEnumValue(
        "Preconditioner", "FullLinearSpaceWithBlock", eLinearWithBlock),
}

Definition at line 153 of file Preconditioner.h.

m_comm

LibUtilities::CommSharedPtr Nektar::MultiRegions::Preconditioner::m_comm

protected

Definition at line 122 of file Preconditioner.h.

Referenced by Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerCG(), Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerHDG(), Nektar::MultiRegions::PreconditionerLowEnergy::SetUpReferenceElements(), Nektar::MultiRegions::PreconditionerLowEnergy::v_BuildPreconditioner(), and Nektar::MultiRegions::PreconditionerLowEnergy::v_InitObject().

m_linsys

const std::weak_ptr<GlobalLinSys> Nektar::MultiRegions::Preconditioner::m_linsys

protected

Definition at line 118 of file Preconditioner.h.

Referenced by AssembleStaticCondGlobalDiagonals(), Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerCG(), Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerHDG(), Nektar::MultiRegions::PreconditionerDiagonal::DiagonalPreconditionerSum(), Nektar::MultiRegions::PreconditionerLowEnergy::SetupBlockTransformationMatrix(), Nektar::MultiRegions::PreconditionerLowEnergy::SetUpReferenceElements(), Nektar::MultiRegions::PreconditionerBlock::v_BuildPreconditioner(), Nektar::MultiRegions::PreconditionerLinear::v_BuildPreconditioner(), Nektar::MultiRegions::PreconditionerLowEnergy::v_BuildPreconditioner(), Nektar::MultiRegions::PreconditionerLinearWithBlock::v_InitObject(), Nektar::MultiRegions::PreconditionerLinearWithDiag::v_InitObject(), Nektar::MultiRegions::PreconditionerLinearWithLowEnergy::v_InitObject(), Nektar::MultiRegions::PreconditionerLowEnergy::v_InitObject(), and Nektar::MultiRegions::PreconditionerLowEnergy::v_TransformedSchurCompl().

m_locToGloMap

std::weak_ptr<AssemblyMap> Nektar::MultiRegions::Preconditioner::m_locToGloMap

protected

Definition at line 121 of file Preconditioner.h.

Referenced by AssembleStaticCondGlobalDiagonals(), Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerCG(), Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerHDG(), Nektar::MultiRegions::PreconditionerLowEnergy::CreateVariablePMask(), Nektar::MultiRegions::PreconditionerDiagonal::DiagonalPreconditionerSum(), Nektar::MultiRegions::PreconditionerLowEnergy::SetupBlockTransformationMatrix(), Nektar::MultiRegions::PreconditionerDiagonal::StaticCondDiagonalPreconditionerSum(), Nektar::MultiRegions::PreconditionerDiagonal::v_BuildPreconditioner(), Nektar::MultiRegions::PreconditionerLinear::v_BuildPreconditioner(), Nektar::MultiRegions::PreconditionerLowEnergy::v_BuildPreconditioner(), Nektar::MultiRegions::PreconditionerBlock::v_DoPreconditioner(), Nektar::MultiRegions::PreconditionerDiagonal::v_DoPreconditioner(), Nektar::MultiRegions::PreconditionerLinearWithLowEnergy::v_DoPreconditioner(), Nektar::MultiRegions::PreconditionerLowEnergy::v_DoPreconditioner(), Nektar::MultiRegions::PreconditionerLinear::v_DoPreconditionerWithNonVertOutput(), Nektar::MultiRegions::PreconditionerLowEnergy::v_DoTransformBasisFromLowEnergy(), Nektar::MultiRegions::PreconditionerLowEnergy::v_DoTransformBasisToLowEnergy(), Nektar::MultiRegions::PreconditionerLowEnergy::v_DoTransformCoeffsFromLowEnergy(), Nektar::MultiRegions::PreconditionerLowEnergy::v_DoTransformCoeffsToLowEnergy(), Nektar::MultiRegions::PreconditionerBlock::v_InitObject(), Nektar::MultiRegions::PreconditionerLinearWithBlock::v_InitObject(), Nektar::MultiRegions::PreconditionerLinearWithDiag::v_InitObject(), Nektar::MultiRegions::PreconditionerLinearWithLowEnergy::v_InitObject(), and Nektar::MultiRegions::PreconditionerLowEnergy::v_InitObject().

m_preconditioner

DNekMatSharedPtr Nektar::MultiRegions::Preconditioner::m_preconditioner

protected

Definition at line 120 of file Preconditioner.h.

m_preconType

PreconditionerType Nektar::MultiRegions::Preconditioner::m_preconType

protected

Definition at line 119 of file Preconditioner.h.