itk::fem::LinearSystemWrapperVNL Class Reference

#include <itkFEMLinearSystemWrapperVNL.h>

Inheritance diagram for itk::fem::LinearSystemWrapperVNL:

Collaboration diagram for itk::fem::LinearSystemWrapperVNL:

Detailed Description

LinearSystemWrapper class that uses VNL numeric library functions to define a sparse linear system of equations.

See Also

LinearSystemWrapper

Definition at line 40 of file itkFEMLinearSystemWrapperVNL.h.

Public Types
typedef LinearSystemWrapper::Float	Float
typedef std::vector < MatrixRepresentation * >	MatrixHolder
typedef vnl_sparse_matrix< Float >	MatrixRepresentation
typedef LinearSystemWrapper	SuperClass
Public Types inherited from itk::fem::LinearSystemWrapper
typedef std::vector< unsigned int >	ColumnArray
typedef const Self *	ConstPointer
typedef Self *	Pointer
typedef LinearSystemWrapper	Self
typedef Solution	Superclass
Public Types inherited from itk::fem::Solution
typedef const Self *	ConstPointer
typedef double	Float
typedef Self *	Pointer
typedef Solution	Self
typedef Solution	Superclass

Public Member Functions
virtual void	AddMatrixValue (unsigned int i, unsigned int j, Float value, unsigned int matrixIndex) override
virtual void	AddSolutionValue (unsigned int i, Float value, unsigned int solutionIndex) override
virtual void	AddVectorValue (unsigned int i, Float value, unsigned int vectorIndex) override
virtual void	CopySolution2Vector (unsigned solutionIndex, unsigned int vectorIndex) override
virtual void	CopyVector2Solution (unsigned int vectorIndex, unsigned int solutionIndex) override
virtual void	DestroyMatrix (unsigned int matrixIndex) override
virtual void	DestroySolution (unsigned int solutionIndex) override
virtual void	DestroyVector (unsigned int vectorIndex) override
virtual Float	GetMatrixValue (unsigned int i, unsigned int j, unsigned int matrixIndex) const override
virtual Float	GetSolutionValue (unsigned int i, unsigned int solutionIndex) const override
virtual Float	GetVectorValue (unsigned int i, unsigned int vectorIndex) const override
virtual void	InitializeMatrix (unsigned int matrixIndex) override
virtual void	InitializeSolution (unsigned int solutionIndex) override
virtual void	InitializeVector (unsigned int vectorIndex) override
virtual bool	IsMatrixInitialized (unsigned int matrixIndex) override
virtual bool	IsSolutionInitialized (unsigned int solutionIndex) override
virtual bool	IsVectorInitialized (unsigned int vectorIndex) override
	LinearSystemWrapperVNL ()
virtual void	MultiplyMatrixMatrix (unsigned int resultMatrixIndex, unsigned int leftMatrixIndex, unsigned int rightMatrixIndex) override
virtual void	MultiplyMatrixVector (unsigned int resultVectorIndex, unsigned int matrixIndex, unsigned int vectorIndex) override
virtual void	ScaleMatrix (Float scale, unsigned int matrixIndex) override
virtual void	SetMatrixValue (unsigned int i, unsigned int j, Float value, unsigned int matrixIndex) override
virtual void	SetMaximumNonZeroValuesInMatrix (unsigned int, unsigned int)
virtual void	SetSolutionValue (unsigned int i, Float value, unsigned int solutionIndex) override
virtual void	SetVectorValue (unsigned int i, Float value, unsigned int vectorIndex) override
virtual void	Solve (void) override
virtual void	SwapMatrices (unsigned int matrixIndex1, unsigned int matrixIndex2) override
virtual void	SwapSolutions (unsigned int solutionIndex1, unsigned int solutionIndex2) override
virtual void	SwapVectors (unsigned int vectorIndex1, unsigned int vectorIndex2) override
virtual	~LinearSystemWrapperVNL ()
Public Member Functions inherited from itk::fem::LinearSystemWrapper
virtual void	AddMatrixMatrix (unsigned int matrixIndex1, unsigned int matrixIndex2)
virtual void	AddVectorVector (unsigned int vectorIndex1, unsigned int vectorIndex2)
virtual void	Clean ()
virtual void	CopyMatrix (unsigned int matrixIndex1, unsigned int matrixIndex2)
virtual void	CopySolution2Vector (unsigned int solutionIndex, unsigned int vectorIndex)=0
virtual void	CopyVector (unsigned int vectorSource, unsigned int vectorDestination)
virtual void	GetColumnsOfNonZeroMatrixElementsInRow (unsigned int row, ColumnArray &cols, unsigned int matrixIndex=0)
unsigned int	GetNumberOfMatrices () const
unsigned int	GetNumberOfSolutions () const
unsigned int	GetNumberOfVectors () const
unsigned int	GetSystemOrder () const
virtual void	MultiplyMatrixSolution (unsigned int resultVectorIndex, unsigned int matrixIndex, unsigned int solutionIndex)
virtual void	OptimizeMatrixStorage (unsigned int matrixIndex, unsigned int tempMatrixIndex)
virtual void	ReverseCuthillMckeeOrdering (ColumnArray &newNumbering, unsigned int matrixIndex=0)
void	ScaleSolution (Float scale, unsigned int solutionIndex=0)
void	ScaleVector (Float scale, unsigned int vectorIndex=0)
void	SetNumberOfMatrices (unsigned int nMatrices)
void	SetNumberOfSolutions (unsigned int nSolutions)
void	SetNumberOfVectors (unsigned int nVectors)
void	SetSystemOrder (unsigned int N)
virtual	~LinearSystemWrapper ()
	LinearSystemWrapper ()
Public Member Functions inherited from itk::fem::Solution
virtual	~Solution ()

Private Attributes
MatrixHolder *	m_Matrices
std::vector< vnl_vector< Float > * > *	m_Solutions
std::vector< vnl_vector< Float > * > *	m_Vectors

Additional Inherited Members
Protected Attributes inherited from itk::fem::LinearSystemWrapper
unsigned int	m_NumberOfMatrices
unsigned int	m_NumberOfSolutions
unsigned int	m_NumberOfVectors
unsigned int	m_Order

Member Typedef Documentation

typedef LinearSystemWrapper::Float itk::fem::LinearSystemWrapperVNL::Float

Definition at line 45 of file itkFEMLinearSystemWrapperVNL.h.

typedef std::vector<MatrixRepresentation *> itk::fem::LinearSystemWrapperVNL::MatrixHolder

Definition at line 54 of file itkFEMLinearSystemWrapperVNL.h.

typedef vnl_sparse_matrix<Float> itk::fem::LinearSystemWrapperVNL::MatrixRepresentation

Definition at line 51 of file itkFEMLinearSystemWrapperVNL.h.

typedef LinearSystemWrapper itk::fem::LinearSystemWrapperVNL::SuperClass

Definition at line 48 of file itkFEMLinearSystemWrapperVNL.h.

Constructor & Destructor Documentation

itk::fem::LinearSystemWrapperVNL::LinearSystemWrapperVNL ( )

inline

Definition at line 57 of file itkFEMLinearSystemWrapperVNL.h.

virtual itk::fem::LinearSystemWrapperVNL::~LinearSystemWrapperVNL ( )

virtual

Member Function Documentation

virtual void itk::fem::LinearSystemWrapperVNL::AddMatrixValue ( unsigned int  i, unsigned int  j, Float  value, unsigned int  matrixIndex  )

inlineoverridevirtual

Virtual function to add a value to a specific element of the A matrix.

Parameters

i	row of the element
j	column of the element
value	value to add to the existing element
matrixIndex	index of matrix to add value to

Implements itk::fem::LinearSystemWrapper.

Definition at line 96 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::AddSolutionValue ( unsigned int  i, Float  value, unsigned int  solutionIndex  )

inlineoverridevirtual

Virtual function to add a value of specific element of the solution vector.

Parameters

i	element Index in solution vector
value	new value of the element
solutionIndex	index of solution vector to add value to

Implements itk::fem::LinearSystemWrapper.

Definition at line 123 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::AddVectorValue ( unsigned int  i, Float  value, unsigned int  vectorIndex  )

inlineoverridevirtual

Virtual function to add a value to a specific element of the B vector.

Parameters

i	row of the element
value	value to add to the existing element
vectorIndex	index of vector to add value to

Implements itk::fem::LinearSystemWrapper.

Definition at line 111 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::CopySolution2Vector ( unsigned  solutionIndex, unsigned int  vectorIndex  )

overridevirtual

virtual void itk::fem::LinearSystemWrapperVNL::CopyVector2Solution ( unsigned int  vectorIndex, unsigned int  solutionIndex  )

overridevirtual

Copy a vector to a solution vector

Parameters

vectorIndex	index of a vector to copy
solutionIndex	index of a solution to copy the solution to

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::DestroyMatrix ( unsigned int  matrixIndex )

overridevirtual

Free the memory from a matrix

Parameters

matrixIndex

index of matrix to destroy

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::DestroySolution ( unsigned int  solutionIndex )

overridevirtual

Free the memory from a solution vector

Parameters

solutionIndex

index of solution vector to destroy

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::DestroyVector ( unsigned int  vectorIndex )

overridevirtual

Free the memory from a vector

Parameters

vectorIndex

index of vector to destroy

Implements itk::fem::LinearSystemWrapper.

virtual Float itk::fem::LinearSystemWrapperVNL::GetMatrixValue ( unsigned int  i, unsigned int  j, unsigned int  matrixIndex  ) const

inlineoverridevirtual

Virtual function to get a value of a specific element of a matrix.

Parameters

i	row of the element
j	column of the element
matrixIndex	index of matrix to get value from

Implements itk::fem::LinearSystemWrapper.

Definition at line 86 of file itkFEMLinearSystemWrapperVNL.h.

virtual Float itk::fem::LinearSystemWrapperVNL::GetSolutionValue ( unsigned int  i, unsigned int  solutionIndex  ) const

overridevirtual

Returns value of i-th element in a solution vector. This value is calculated generalized displacement of the i-th degree of freedom in a FEM problem. Note that in general there may be several solution vectors. You can select which one do you want by passing the second parameter.

Parameters

i	element index in solution vector
solutionIndex	index of solution vector to get value from

Note

If the solution vector doesn't exist (problem was not yet solved), or the index i is out of range, the function returns 0.0.

Implements itk::fem::Solution.

virtual Float itk::fem::LinearSystemWrapperVNL::GetVectorValue ( unsigned int  i, unsigned int  vectorIndex  ) const

inlineoverridevirtual

Virtual function to get a value of a specific element of the B vector.

Parameters

i	row of the element
vectorIndex	index of vector to get value from

Implements itk::fem::LinearSystemWrapper.

Definition at line 101 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::InitializeMatrix ( unsigned int  matrixIndex )

overridevirtual

Initialization of the A matrix. First any existing data for matrix A must be be destroyed, and then a new matrix is created in the memory. All elements in A must be set to zero.

Parameters

matrixIndex

index of matrix to initialize

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::InitializeSolution ( unsigned int  solutionIndex )

overridevirtual

Initialization of a solution vector. Existing memory must be destroyed and the new solution vector is created in memory. All values should be set to zero.

Parameters

solutionIndex

index of solution vector to initialize

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::InitializeVector ( unsigned int  vectorIndex )

overridevirtual

Initialization of the a vector. First any existing data for vector B must be destroyed, then new vector is created in the memory. All elements in B must be set to zero.

Implements itk::fem::LinearSystemWrapper.

virtual bool itk::fem::LinearSystemWrapperVNL::IsMatrixInitialized ( unsigned int  matrixIndex )

overridevirtual

Check to see if matrix is initialized

Parameters

matrixIndex

index of matrix to examine

Implements itk::fem::LinearSystemWrapper.

virtual bool itk::fem::LinearSystemWrapperVNL::IsSolutionInitialized ( unsigned int  solutionIndex )

overridevirtual

Check to see if solution vector is initialized

Parameters

solutionIndex

index of solution vector to examine

Implements itk::fem::LinearSystemWrapper.

virtual bool itk::fem::LinearSystemWrapperVNL::IsVectorInitialized ( unsigned int  vectorIndex )

overridevirtual

Check to see if vector is initialized

Parameters

vectorIndex

vector of index to examine

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::MultiplyMatrixMatrix ( unsigned int  resultMatrixIndex, unsigned int  leftMatrixIndex, unsigned int  rightMatrixIndex  )

overridevirtual

Perform a matrix*matrix operation and store the result in the linear system

Parameters

leftMatrixIndex	index of left matrix
rightMatrixIndex	index of right matrix
resultMatrixIndex	index of matrix where solution is stored

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::MultiplyMatrixVector ( unsigned int  resultVectorIndex, unsigned int  matrixIndex, unsigned int  vectorIndex  )

overridevirtual

Perform a matrix*vector operation and store the result in the linear system

Parameters

matrixIndex	index of matrix to multiply
vectorIndex	index of vector to multiply
resultVectorIndex	index of vector where result is store

Reimplemented from itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::ScaleMatrix ( Float  scale, unsigned int  matrixIndex  )

overridevirtual

Multiplies all elements of a matrix by a scalar

Parameters

scale	scalar to multiply all matrix values by
matrixIndex	index of matrix to modify

Reimplemented from itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::SetMatrixValue ( unsigned int  i, unsigned int  j, Float  value, unsigned int  matrixIndex  )

inlineoverridevirtual

Virtual function to set a value of a specific element of the A matrix.

Parameters

i	row of the element
j	column of the element
value	new value of the element
matrixIndex	index of matrix to set value in

Implements itk::fem::LinearSystemWrapper.

Definition at line 91 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::SetMaximumNonZeroValuesInMatrix ( unsigned  int, unsigned  int  )

inlinevirtual

Definition at line 81 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::SetSolutionValue ( unsigned int  i, Float  value, unsigned int  solutionIndex  )

inlineoverridevirtual

Virtual function to set a value of specific element of the solution vector.

Parameters

i	element Index in solution vector
value	new value of the element
solutionIndex	index of solution vector to set value in

Implements itk::fem::LinearSystemWrapper.

Definition at line 118 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::SetVectorValue ( unsigned int  i, Float  value, unsigned int  vectorIndex  )

inlineoverridevirtual

Virtual function to set a value of a specific element of the B vector.

Parameters

i	row of the element
value	new value of the element
vectorIndex	index of vector to set value in

Implements itk::fem::LinearSystemWrapper.

Definition at line 106 of file itkFEMLinearSystemWrapperVNL.h.

virtual void itk::fem::LinearSystemWrapperVNL::Solve ( void  )

overridevirtual

Solves the linear system and creates the solution vector, which can later be accessed via GetSolutionValue(i,SolutionIndex) member function. Here all the major processing is done with calls to external numeric library.

Note

This function can only be called after the linear system was properly assembled.

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::SwapMatrices ( unsigned int  matrixIndex1, unsigned int  matrixIndex2  )

overridevirtual

Swaps access indices of any 2 matrices in the linear system

Parameters

matrixIndex1	index of a matrix to swap
matrixIndex2	index of matrix to swap with

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::SwapSolutions ( unsigned int  solutionIndex1, unsigned int  solutionIndex2  )

overridevirtual

Swaps access indices of any 2 solution vectors in the linear system

Parameters

solutionIndex1	index of a solution vector to swap
solutionIndex2	index of solution vector to swap with

Implements itk::fem::LinearSystemWrapper.

virtual void itk::fem::LinearSystemWrapperVNL::SwapVectors ( unsigned int  vectorIndex1, unsigned int  vectorIndex2  )

overridevirtual

Swaps access indices of any 2 vectors in the linear system

Parameters

vectorIndex1	index of a vector to swap
vectorIndex2	index of vector to swap with

Implements itk::fem::LinearSystemWrapper.

Member Data Documentation

MatrixHolder* itk::fem::LinearSystemWrapperVNL::m_Matrices

private

vector of pointers to VNL sparse matrices

Definition at line 152 of file itkFEMLinearSystemWrapperVNL.h.

std::vector<vnl_vector<Float> *>* itk::fem::LinearSystemWrapperVNL::m_Solutions

private

vector of pointers to VNL vectors

Definition at line 158 of file itkFEMLinearSystemWrapperVNL.h.

std::vector<vnl_vector<Float> *>* itk::fem::LinearSystemWrapperVNL::m_Vectors

private

vector of pointers to VNL vectors

Definition at line 155 of file itkFEMLinearSystemWrapperVNL.h.

---

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

• itkFEMLinearSystemWrapperVNL.h