Implementation of matrix containing floating point numbers. More...

#include <floatmatrix.h>

Public Member Functions
	FloatMatrix (int n, int m)
	Creates matrix of given size. More...

	FloatMatrix ()
	Creates zero sized matrix. More...

	FloatMatrix (const FloatArray &vector, bool transpose=false)
	Constructor. More...

	FloatMatrix (const FloatMatrix &mat)
	Copy constructor. More...

	FloatMatrix (FloatMatrix &&mat)
	Copy constructor. More...

	FloatMatrix (std::initializer_list< std::initializer_list< double > >mat)
	Initializer list constructor. More...

FloatMatrix &	operator= (std::initializer_list< std::initializer_list< double > >mat)
	Assignment operator. More...

FloatMatrix &	operator= (std::initializer_list< FloatArray >mat)
	Assignment operator. More...

FloatMatrix &	operator= (const FloatMatrix &mat)
	Assignment operator, adjusts size of the receiver if necessary. More...

FloatMatrix &	operator= (FloatMatrix &&mat)

	~FloatMatrix ()
	Destructor. More...

void	checkBounds (int i, int j) const
	Checks size of receiver towards requested bounds. More...

int	giveNumberOfRows () const
	Returns number of rows of receiver. More...

int	giveNumberOfColumns () const
	Returns number of columns of receiver. More...

bool	isSquare () const
	Returns nonzero if receiver is square matrix. More...

bool	isNotEmpty () const
	Tests for empty matrix. More...

bool	isFinite () const
	Returns true if no element is NAN or infinite. More...

double	at (int i, int j) const
	Coefficient access function. More...

double &	at (int i, int j)
	Coefficient access function. More...

double &	operator() (int i, int j)
	Coefficient access function. More...

double	operator() (int i, int j) const
	Coefficient access function. More...

void	assemble (const FloatMatrix &src, const IntArray &loc)
	Assembles the contribution using localization array into receiver. More...

void	assemble (const FloatMatrix &src, const IntArray &rowind, const IntArray &colind)
	Assembles the contribution using localization array into receiver. More...

void	assemble (const FloatMatrix &src, const int rowind, const int colind)
	Assembles the contribution using localization array into receiver. More...

double	computeFrobeniusNorm () const
	Computes the Frobenius norm of the receiver. More...

double	computeNorm (char p) const
	Computes the operator norm of the receiver. More...

double	computeReciprocalCondition (char p= '1') const
	Computes the conditioning of the receiver. More...

void	beDiagonal (const FloatArray &diag)
	Modifies receiver to be a diagonal matrix with the components specified in diag. More...

double	giveTrace () const
	Returns the trace of the receiver. More...

double	giveDeterminant () const
	Returns the trace (sum of diagonal components) of the receiver. More...

void	zero ()
	Zeroes all coefficient of receiver. More...

void	beUnitMatrix ()
	Sets receiver to unity matrix. More...

void	bePinvID ()
	Sets receiver to the inverse of scaling matrix P multiplied by the deviatoric projector ID. More...

void	beTranspositionOf (const FloatMatrix &src)
	Assigns to the receiver the transposition of parameter. More...

void	beProductOf (const FloatMatrix &a, const FloatMatrix &b)
	Assigns to the receiver product of $a \cdot b$ . More...

void	addProductOf (const FloatMatrix &a, const FloatMatrix &b)
	Adds to the receiver product of $a \cdot b$ . More...

void	addTProductOf (const FloatMatrix &a, const FloatMatrix &b)
	Adds to the receiver product of $a^{\mathrm{T}} \cdot b$ . More...

void	beTProductOf (const FloatMatrix &a, const FloatMatrix &b)
	Assigns to the receiver product of $a^{\mathrm{T}} \cdot b$ . More...

void	beProductTOf (const FloatMatrix &a, const FloatMatrix &b)
	Assigns to the receiver product of $a \cdot b^{\mathrm{T}}$ . More...

void	beDyadicProductOf (const FloatArray &vec1, const FloatArray &vec2)
	Assigns to the receiver the dyadic product $v_1 \cdot v_2^{\mathrm{T}}$ . More...

void	beNMatrixOf (const FloatArray &n, int nsd)
	Assigns the receiver to be a repeated diagonal matrix. More...

void	beLocalCoordSys (const FloatArray &normal)
	Makes receiver the local coordinate for the given normal. More...

void	setSubMatrix (const FloatMatrix &src, int sr, int sc)
	Adds the given matrix as sub-matrix to receiver. More...

void	setTSubMatrix (const FloatMatrix &src, int sr, int sc)
	Adds the given matrix as sub-matrix to receiver. More...

void	beSubMatrixOf (const FloatMatrix &src, int topRow, int bottomRow, int topCol, int bottomCol)
	Assigns to the receiver the sub-matrix of another matrix. More...

void	beSubMatrixOf (const FloatMatrix &src, const IntArray &indxRow, const IntArray &indxCol)
	Modifies receiver to be a sub-matrix of another matrix. More...

void	addSubVectorRow (const FloatArray &src, int sr, int sc)
	Adds given vector to receiver starting at given position. More...

void	addSubVectorCol (const FloatArray &src, int sr, int sc)
	Adds given vector to receiver starting at given position. More...

void	copySubVectorRow (const FloatArray &src, int sr, int sc)
	Copy (set) given vector to receiver row sr, starting at column sc. More...

void	setColumn (const FloatArray &src, int c)
	Sets the values of the matrix in specified column. More...

void	copyColumn (FloatArray &dest, int c) const
	Fetches the values from the specified column. More...

void	beInverseOf (const FloatMatrix &src)
	Modifies receiver to become inverse of given parameter. More...

bool	solveForRhs (const FloatArray &b, FloatArray &answer, bool transpose=false)
	Solves the system of linear equations $K\cdot a = b$ . More...

void	solveForRhs (const FloatMatrix &B, FloatMatrix &answer, bool transpose=false)
	Solves the system of linear equations $K\cdot A = B$ . More...

void	plusProductSymmUpper (const FloatMatrix &a, const FloatMatrix &b, double dV)
	Adds to the receiver the product $a^{\mathrm{T}}\cdot b \mathrm{d}V$ . More...

void	plusDyadSymmUpper (const FloatArray &a, double dV)
	Adds to the receiver the dyadic product $a \otimes a \mathrm{d}V$ . More...

void	plusProductUnsym (const FloatMatrix &a, const FloatMatrix &b, double dV)
	Adds to the receiver the product $a^{\mathrm{T}} \cdot b \mathrm{d}V$ . More...

void	plusDyadUnsym (const FloatArray &a, const FloatArray &b, double dV)
	Adds to the receiver the product $a \otimes b \mathrm{d}V$ . More...

void	add (const FloatMatrix &a)
	Adds matrix to the receiver. More...

void	add (double s, const FloatMatrix &a)
	Adds matrix to the receiver scaled by s. More...

void	subtract (const FloatMatrix &a)
	Subtracts matrix from the receiver. More...

void	times (double f)
	Multiplies receiver by factor f. More...

void	negated ()
	Changes sign of receiver values. More...

void	initFromVector (const FloatArray &vector, bool transposed)
	Assigns to receiver one column or one row matrix containing vector. More...

void	symmetrized ()
	Initializes the lower half of the receiver according to the upper half. More...

void	rotatedWith (const FloatMatrix &r, char mode= 'n')
	Returns the receiver 'a' transformed using give transformation matrix r. More...

void	resize (int rows, int cols)
	Checks size of receiver towards requested bounds. More...

void	resizeWithData (int, int)
	Checks size of receiver towards requested bounds. More...

void	hardResize (int r, int c)
	Resizing that enforces reallocation of memory. More...

void	clear ()
	Sets size of receiver to be an empty matrix. It will have zero rows and zero columns size. More...

bool	jaco_ (FloatArray &eval, FloatMatrix &v, int nf)
	Computes eigenvalues and eigenvectors of receiver (must be symmetric) The receiver is preserved. More...

void	printYourself () const
	Prints matrix to stdout. Useful for debugging. More...

void	printYourself (const std::string &name) const
	Print receiver on stdout with custom name. More...

void	printYourselfToFile (const std::string filename, const bool showDimensions=true) const
	Print matrix to file. More...

void	pY () const
	Higher accuracy than printYourself. More...

void	writeCSV (const std::string &name) const
	Writes receiver as CSV (comma seperated values) More...

const double *	givePointer () const
	Exposes the internal values of the matrix. More...

double *	givePointer ()

void	beMatrixFormOfStress (const FloatArray &aArray)
	Reciever will be a 3x3 matrix formed from a vector with either 9 or 6 components. More...

void	beMatrixForm (const FloatArray &aArray)

void	changeComponentOrder ()
	Swaps the indices in the 6x6 matrix such that it converts between OOFEM's and Abaqus' way of writing matrices. More...

contextIOResultType	storeYourself (DataStream &stream) const

contextIOResultType	restoreYourself (DataStream &stream)

int	givePackSize (DataStream &buff) const

Protected Attributes
int	nRows
	Number of rows. More...

int	nColumns
	Number of columns. More...

std::vector< double >	values
	Values of matrix stored column wise. More...

Friends
std::ostream &	operator<< (std::ostream &out, const FloatMatrix &r)

Detailed Description

Implementation of matrix containing floating point numbers.

FloatMatrix can grow and shrink to requested size if required. Implementation includes many computing and manipulation methods. Rows and Columns are indexed starting from 1.

The matrix stores its nRows*nColumns coefficients in 'values'. The coefficients are stored column by column, like in Fortran.

Tasks:

Storing and retrieving a coefficient (method 'at') ;
Performing standard operations : addition, multiplication, transposition, inversion, lumping, rotation, etc ;
Introduced allocatedSize variable to allow dynamic rescaling of matrix size possibly without memory reallocation. At startup matrix occupies space given by allocatedSpace = nRows*nColums. Then there can be further request for resizing matrix to smaller dimension then we only change nRows and nColumns variables , but allocatedSize variable remain untouched - expecting possible matrix grow and then re-using previously allocated space. If further request for growing then is necessary memory reallocation. This process is controlled in resize member function.

Author: Mikael Öhman; Jim Brouzoulis; many others (please add yourselves)

Definition at line 94 of file floatmatrix.h.

Constructor & Destructor Documentation

oofem::FloatMatrix::FloatMatrix	(	int	n,
		int	m
	)

inline

Creates matrix of given size.

Parameters

n	Number of rows.
m	Requested number of columns.

Definition at line 110 of file floatmatrix.h.

oofem::FloatMatrix::FloatMatrix ( )

inline

Creates zero sized matrix.

Definition at line 112 of file floatmatrix.h.

oofem::FloatMatrix::FloatMatrix	(	const FloatArray &	vector,
		bool	transpose = `false`
	)

Constructor.

Creates float matrix from float vector. Vector may be stored row wise or column wise, depending on second parameter.

Parameters

vector	Float vector from which matrix is constructed
transpose	If false (default) then a matrix of size (vector->giveSize(),1) will be created and initialized, if true then a matrix of size (1,vector->giveSize()) will be created.

Definition at line 113 of file floatmatrix.C.

References oofem::FloatArray::giveSize(), nColumns, nRows, oofem::FloatArray::values, and values.

oofem::FloatMatrix::FloatMatrix ( const FloatMatrix & mat )

inline

Copy constructor.

Definition at line 123 of file floatmatrix.h.

oofem::FloatMatrix::FloatMatrix ( FloatMatrix && mat )

inline

Copy constructor.

Definition at line 125 of file floatmatrix.h.

oofem::FloatMatrix::FloatMatrix ( std::initializer_list< std::initializer_list< double > > mat )

Initializer list constructor.

Definition at line 131 of file floatmatrix.C.

References nRows, OOFEM_ERROR, RESIZE, and values.

oofem::FloatMatrix::~FloatMatrix ( )

inline

Destructor.

Definition at line 146 of file floatmatrix.h.

Member Function Documentation

void oofem::FloatMatrix::add ( const FloatMatrix & a )

Adds matrix to the receiver.

If receiver has zero size, size is accordingly adjusted.

Parameters

a	Matrix to be added.

Definition at line 1023 of file floatmatrix.C.

References givePointer(), isNotEmpty(), nColumns, nRows, OOFEM_ERROR, operator=(), and values.

void oofem::FloatMatrix::add	(	double	s,
		const FloatMatrix &	a
	)

Adds matrix to the receiver scaled by s.

If receiver has zero size, size is accordingly adjusted.

Parameters

s	Scaling factor.
a	Matrix to be added.

Definition at line 1054 of file floatmatrix.C.

References givePointer(), isNotEmpty(), nColumns, nRows, OOFEM_ERROR, operator=(), times(), and values.

void oofem::FloatMatrix::addProductOf	(	const FloatMatrix &	a,
		const FloatMatrix &	b
	)

Adds to the receiver product of $a \cdot b$ .

Definition at line 427 of file floatmatrix.C.

References at(), givePointer(), nColumns, nRows, and OOFEM_ERROR.

Referenced by oofem::Shell7BaseXFEM::discComputeStiffness().

void oofem::FloatMatrix::addSubVectorCol	(	const FloatArray &	src,
		int	sr,
		int	sc
	)

Adds given vector to receiver starting at given position.

Parameters

src	Source matrix.
sr	Starting row position.
sc	Starting column position.

Definition at line 627 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), giveNumberOfColumns(), giveNumberOfRows(), oofem::FloatArray::giveSize(), oofem::max(), and resizeWithData().

Referenced by oofem::Node2NodeContactL::computeContactTangent(), oofem::ContactElement::ContactElement(), and oofem::Shell7Base::giveTractionBC().

void oofem::FloatMatrix::addSubVectorRow	(	const FloatArray &	src,
		int	sr,
		int	sc
	)

Adds given vector to receiver starting at given position.

Parameters

src	Source matrix.
sr	Starting row position.
sc	Starting column position.

Definition at line 607 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), giveNumberOfColumns(), giveNumberOfRows(), oofem::FloatArray::giveSize(), oofem::max(), and resizeWithData().

Referenced by oofem::Node2NodeContactL::computeContactTangent(), oofem::ContactElement::ContactElement(), oofem::Shell7Base::giveL2contribution(), oofem::Shell7Base::giveLayerContributionToSR(), oofem::Shell7Base::giveSPRcontribution(), and oofem::Shell7Base::nodalLeastSquareFitFromIP().

void oofem::FloatMatrix::addTProductOf	(	const FloatMatrix &	a,
		const FloatMatrix &	b
	)

Adds to the receiver product of $a^{\mathrm{T}} \cdot b$ .

Definition at line 459 of file floatmatrix.C.

References at(), givePointer(), nColumns, nRows, and OOFEM_ERROR.

void oofem::FloatMatrix::assemble	(	const FloatMatrix &	src,
		const IntArray &	loc
	)

Assembles the contribution using localization array into receiver.

The receiver must have dimensions large enough to localize contribution.

Parameters

src	Source to be assembled.
loc	Localization indices.

Definition at line 251 of file floatmatrix.C.

References oofem::IntArray::at(), at(), giveNumberOfRows(), oofem::IntArray::giveSize(), isSquare(), and OOFEM_ERROR.

void oofem::FloatMatrix::assemble	(	const FloatMatrix &	src,
		const IntArray &	rowind,
		const IntArray &	colind
	)

Assembles the contribution using localization array into receiver.

The receiver must have dimensions large enough to localize contribution.

Parameters

src	Source to be assembled.
rowind	Row localization indices.
colind	Column localization indices.

Definition at line 277 of file floatmatrix.C.

References oofem::IntArray::at(), at(), giveNumberOfColumns(), giveNumberOfRows(), oofem::IntArray::giveSize(), and OOFEM_ERROR.

void oofem::FloatMatrix::assemble	(	const FloatMatrix &	src,
		const int *	rowind,
		const int *	colind
	)

Assembles the contribution using localization array into receiver.

The receiver must have dimensions large enough to localize contribution.

Parameters

src	Source to be assembled.
rowind	Row localization indices.
colind	Column localization indices.

Definition at line 305 of file floatmatrix.C.

References at(), giveNumberOfColumns(), and giveNumberOfRows().

double oofem::FloatMatrix::at	(	int	i,
		int	j
	)		const

inline

Coefficient access function.

Returns value of coefficient at given position of the receiver. Implements 1-based indexing.

Parameters

i	Row position of coefficient.
j	Column position of coefficient.

Definition at line 176 of file floatmatrix.h.

double& oofem::FloatMatrix::at	(	int	i,
		int	j
	)

inline

Coefficient access function.

Returns value of coefficient at given position of the receiver. Implements 1-based indexing.

Parameters

i	Row position of coefficient.
j	Column position of coefficient.

Definition at line 187 of file floatmatrix.h.

void oofem::FloatMatrix::beDiagonal ( const FloatArray & diag )

Modifies receiver to be a diagonal matrix with the components specified in diag.

Returns: Determinant of receiver.

Definition at line 1433 of file floatmatrix.C.

References oofem::FloatArray::giveSize(), and resize().

Referenced by oofem::LineDistributedSpring::computeConstitutiveMatrixAt(), oofem::LumpedMassElement::computeLumpedMassMatrix(), oofem::NodalSpringElement::computeLumpedMassMatrix(), oofem::NodalSpringElement::computeStiffnessMatrix(), and oofem::WinklerMaterial::give3dBeamSubSoilStiffMtrx().

void oofem::FloatMatrix::beDyadicProductOf	(	const FloatArray &	vec1,
		const FloatArray &	vec2
	)

Assigns to the receiver the dyadic product $v_1 \cdot v_2^{\mathrm{T}}$ .

Grows or shrinks if necessary.

Definition at line 492 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), oofem::FloatArray::giveSize(), and RESIZE.

Referenced by oofem::Node2NodeContact::computeContactTangent(), oofem::IsotropicDamageMaterial1::computeEta(), oofem::PrescribedGradientBCWeak::computeField(), oofem::tet21ghostsolid::computeStiffnessMatrix(), oofem::IntElLine1PF::computeStiffnessMatrix_ud(), oofem::TrabBone3D::constructDerivativeOfPlasFlowDirec(), oofem::ContactElement::ContactElement(), oofem::FEI2dLineHermite::evaldNdx(), oofem::MisesMat::give3dLSMaterialStiffnessMatrix(), oofem::MisesMatGrad::give3dMaterialStiffnessMatrix(), oofem::TrabBoneGrad3D::give3dMaterialStiffnessMatrix(), oofem::MisesMat::give3dMaterialStiffnessMatrix(), oofem::TrabBoneNL3D::give3dMaterialStiffnessMatrix(), oofem::TrabBone3D::give3dMaterialStiffnessMatrix(), oofem::IntMatBilinearCZFagerstrom::give3dStiffnessMatrix_dTdj(), oofem::IntMatBilinearCZJansson::give3dStiffnessMatrix_dTdj(), oofem::NonlinearMassTransferMaterial::giveCharacteristicMatrix(), oofem::NonlinearFluidMaterial::giveDeviatoricStiffnessMatrix(), oofem::MisesMatGrad::givePlaneStrainStiffMtrx(), oofem::IsotropicDamageMaterial::givePlaneStressStiffMtrx(), oofem::AnisotropicDamageMaterial::obtainAlpha3(), and oofem::IntElLine1PF::solveForLocalDamage().

void oofem::FloatMatrix::beInverseOf ( const FloatMatrix & src )

Modifies receiver to become inverse of given parameter.

Size of receiver will be adjusted.

Parameters

src	Matrix to be inverted.

Definition at line 835 of file floatmatrix.C.

References at(), oofem::IntArray::givePointer(), oofem::FloatArray::givePointer(), givePointer(), isSquare(), nColumns, nRows, OOFEM_ERROR, RESIZE, and zero().

void oofem::FloatMatrix::beLocalCoordSys ( const FloatArray & normal )

Makes receiver the local coordinate for the given normal.

Implemented for 2D and 3D.

Parameters

normal Normal (normalized).

Definition at line 515 of file floatmatrix.C.

References oofem::FloatArray::add(), oofem::FloatArray::at(), at(), oofem::FloatArray::beVectorProductOf(), oofem::FloatArray::dotProduct(), oofem::FloatArray::giveSize(), oofem::FloatArray::normalize(), OOFEM_ERROR, and resize().

Referenced by oofem::ParticleTopologyDescription::calculateShortestDistance(), oofem::ParticleTopologyDescription::collectNeighbors(), oofem::Shell7BaseXFEM::computeCohesiveTangentAt(), oofem::IntElSurfTr1::computeTransformationMatrixAt(), oofem::ParticleTopologyDescription::generatePSLG(), oofem::Shell7BaseXFEM::giveFailedInterfaceNumber(), oofem::Shell7BaseXFEM::giveRecoveredTransverseInterfaceStress(), and oofem::Shell7BaseXFEM::recoverShearStress().

void oofem::FloatMatrix::beMatrixForm ( const FloatArray & aArray )

Definition at line 1657 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), oofem::FloatArray::giveSize(), OOFEM_ERROR, and resize().

void oofem::FloatMatrix::beMatrixFormOfStress ( const FloatArray & aArray )

Reciever will be a 3x3 matrix formed from a vector with either 9 or 6 components.

Order of matrix components in vector: 11, 22, 33, 23, 13, 12, 32, 31, 21 If size(aArray) = 6, a symmetric matrix will be created.

Parameters

aArray Array to transform.

Definition at line 1774 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), oofem::IntArray::givePointer(), oofem::FloatArray::givePointer(), givePointer(), oofem::FloatArray::giveSize(), nRows, OOFEM_ERROR, oofem::FloatArray::resize(), and resize().

Referenced by oofem::Shell7Base::computeCauchyStressVector(), oofem::SolidShell::computeGeometricStiffness(), oofem::Shell7Base::computeInterLaminarStressesAt(), oofem::tet21ghostsolid::computeStiffnessMatrix(), oofem::LinearElasticMaterial::giveEshelbyStressVector_PlaneStrain(), oofem::tet21ghostsolid::giveInternalForcesVectorGivenSolution(), and oofem::tet21ghostsolid::giveInternalForcesVectorGivenSolutionDebug().

void oofem::FloatMatrix::beNMatrixOf	(	const FloatArray &	n,
		int	nsd
	)

Assigns the receiver to be a repeated diagonal matrix.

Parameters

n	Vector with components which will appear in respective diagonal.
nsd	Number of spatial dimensions

Definition at line 505 of file floatmatrix.C.

References oofem::FloatArray::giveSize(), and resize().

Referenced by oofem::PrescribedGradientBCWeak::assembleTangentGPContributionNew(), oofem::PrescribedGradientBCWeak::compute_x_times_N_1(), oofem::PrescribedGradientBCWeak::compute_x_times_N_2(), oofem::TransportElement::computeBoundaryEdgeLoadVector(), oofem::TransportElement::computeBoundarySurfaceLoadVector(), oofem::tet21ghostsolid::computeBoundarySurfaceLoadVector(), oofem::MITC4Shell::computeEdgeNMatrix(), oofem::StructuralElement::computeEdgeNMatrix(), oofem::PFEMElement2d::computeEdgeNMatrixAt(), oofem::PrescribedGradientBCWeak::computeExtForceElContrib(), oofem::PrescribedGradientBCWeak::computeField(), oofem::TransportElement::computeLoadVector(), oofem::XfemElementInterface::computeNCohesive(), oofem::PhaseFieldElement::computeNd_matrixAt(), oofem::StructuralInterfaceElementPhF::computeNd_matrixAt(), oofem::PlaneStressStructuralElementEvaluator::computeNMatrixAt(), oofem::Space3dStructuralElementEvaluator::computeNMatrixAt(), oofem::Tr_Warp::computeNmatrixAt(), oofem::QDKTPlate::computeNmatrixAt(), oofem::DKTPlate::computeNmatrixAt(), oofem::Truss3d::computeNmatrixAt(), oofem::TransportElement::computeNmatrixAt(), oofem::StructuralElement::computeNmatrixAt(), oofem::Tet1_3D_SUPG::computeNuMatrix(), oofem::TR21_2D_SUPG::computeNuMatrix(), oofem::Quad10_2D_SUPG::computeNuMatrix(), oofem::QDKTPlate::computeSurfaceNMatrix(), oofem::MITC4Shell::computeSurfaceNMatrix(), oofem::StructuralElement::computeSurfaceNMatrix(), oofem::Structural3DElement::computeSurfaceNMatrixAt(), oofem::TransportElement::computeTangentFromEdgeLoad(), oofem::TransportElement::computeTangentFromSurfaceLoad(), oofem::MixedGradientPressureNeumann::integrateDevTangent(), oofem::StokesFlowVelocityHomogenization::integrateNMatrix(), oofem::PrescribedGradientBCNeumann::integrateTangent(), oofem::MixedGradientPressureWeakPeriodic::integrateTractionVelocityTangent(), oofem::MixedGradientPressureNeumann::integrateVolTangent(), and oofem::XfemElementInterface::XfemElementInterface_createEnrNmatrixAt().

void oofem::FloatMatrix::bePinvID ( )

Sets receiver to the inverse of scaling matrix P multiplied by the deviatoric projector ID.

Definition at line 1347 of file floatmatrix.C.

References resize(), and values.

Referenced by oofem::MisesMatGrad::give3dMaterialStiffnessMatrix(), and oofem::MisesMat::give3dMaterialStiffnessMatrix().

void oofem::FloatMatrix::beProductOf	(	const FloatMatrix &	a,
		const FloatMatrix &	b
	)

Assigns to the receiver product of $a \cdot b$ .

Grows or shrinks if necessary.

Definition at line 337 of file floatmatrix.C.

References at(), givePointer(), nColumns, nRows, OOFEM_ERROR, and RESIZE.

void oofem::FloatMatrix::beProductTOf	(	const FloatMatrix &	a,
		const FloatMatrix &	b
	)

Assigns to the receiver product of $a \cdot b^{\mathrm{T}}$ .

Grows or shrinks if necessary.

Definition at line 397 of file floatmatrix.C.

References at(), givePointer(), nColumns, nRows, OOFEM_ERROR, and RESIZE.

Referenced by oofem::StructuralInterfaceElementPhF::computeBd_matrixAt(), oofem::Shell7Base::computeCauchyStressVector(), oofem::XfemStructuralElementInterface::computeCohesiveTangent(), oofem::Shell7Base::computeFAt(), oofem::SolidShell::computeGeometricStiffness(), oofem::LIBeam3dNL2::computeStiffnessMatrix(), oofem::LIBeam3dNL::computeStiffnessMatrix(), oofem::NURBSSpace3dElement::drawScalar(), oofem::FEI2dTrQuad::evaldNdx(), oofem::FEI2dQuadLin::evaldNdx(), oofem::FEI2dQuadQuad::evaldNdx(), oofem::LargeStrainMasterMaterialGrad::give3dKappaMatrix(), oofem::MisesMat::give3dLSMaterialStiffnessMatrix(), oofem::MooneyRivlinMaterial::give3dMaterialStiffnessMatrix_dPdF(), oofem::RefinedElement::giveCompatibleBcDofArray(), oofem::MPlasticMaterial::giveConsistentStiffnessMatrix(), oofem::MITC4Shell::givedNdx(), oofem::Shell7Base::giveDualBase(), oofem::MPlasticMaterial::giveElastoPlasticStiffnessMatrix(), oofem::MisesMat::giveFirstPKStressVector_3d(), oofem::AbaqusUserMaterial::giveFirstPKStressVector_3d(), oofem::Shell7Base::giveLayerContributionToSR(), oofem::MooneyRivlinMaterial::givePlaneStrainStiffMtrx_dPdF(), oofem::Shell7Base::nodalLeastSquareFitFromIP(), rotatedWith(), oofem::MDM::rotateTensor4(), and oofem::PrescribedGenStrainShell7::setDeformationGradient().

void oofem::FloatMatrix::beSubMatrixOf	(	const FloatMatrix &	src,
		int	topRow,
		int	bottomRow,
		int	topCol,
		int	bottomCol
	)

Assigns to the receiver the sub-matrix of another matrix.

Parameters

src	Matrix from which sub-matrix is taken
topRow	Index of top row, where sub matrix row index starts
bottomRow	Index of bottom row, where sub matrix ends (including this row).
topCol	Index of top column of sub-matrix.
bottomCol	index of bottom column of sub-matrix.

Definition at line 962 of file floatmatrix.C.

References at(), nColumns, nRows, OOFEM_ERROR, and resize().

Referenced by oofem::Shell7BaseXFEM::computeCohesiveTangent(), oofem::NLStructuralElement::computeInitialStressMatrix(), oofem::MaterialForceEvaluator::computeMaterialForce(), oofem::Shell7BaseXFEM::computeStiffnessMatrix(), oofem::NLStructuralElement::computeStiffnessMatrix(), oofem::StructuralElement::computeStiffnessMatrix(), oofem::IntElLine1PF::computeStiffnessMatrix_dd(), oofem::TransportGradientPeriodic::computeTangent(), oofem::StructuralInterfaceMaterial::give2dStiffnessMatrix_Eng(), oofem::PlasticMaterial::giveConsistentStiffnessMatrix(), oofem::MPlasticMaterial::giveConsistentStiffnessMatrix(), oofem::MITC4Shell::givedNdx(), oofem::StructuralMaterial::giveFirstPKStressVector_1d(), oofem::StructuralMaterial::giveFirstPKStressVector_PlaneStress(), oofem::IntMatBilinearCZFagerstromRate::giveFirstPKTraction_3d(), oofem::IntMatBilinearCZJansson::giveFirstPKTraction_3d(), oofem::IntMatBilinearCZFagerstrom::giveFirstPKTraction_3d(), oofem::StructuralMaterial::giveRealStressVector_ShellStressControl(), oofem::StructuralMaterial::giveRealStressVector_StressControl(), oofem::StructuralMaterial::giveReducedMatrixForm(), oofem::StructuralMaterial::giveReducedSymMatrixForm(), oofem::Shell7BaseXFEM::OLDcomputeStiffnessMatrix(), oofem::Shell7Base::recoverShearStress(), oofem::StructuralMaterialEvaluator::solveYourself(), and oofem::FluidMaterialEvaluator::solveYourself().

void oofem::FloatMatrix::beSubMatrixOf	(	const FloatMatrix &	src,
		const IntArray &	indxRow,
		const IntArray &	indxCol
	)

Modifies receiver to be a sub-matrix of another matrix.

Parameters

src	Matrix from which sub-matrix is taken
indxRow	Array of the row indices to be extracted
indxCol	Array of the column indices to be extracted

Definition at line 998 of file floatmatrix.C.

References oofem::IntArray::at(), at(), giveNumberOfColumns(), giveNumberOfRows(), oofem::IntArray::giveSize(), oofem::IntArray::maximum(), oofem::IntArray::minimum(), OOFEM_ERROR, and resize().

void oofem::FloatMatrix::beTProductOf	(	const FloatMatrix &	a,
		const FloatMatrix &	b
	)

Assigns to the receiver product of $a^{\mathrm{T}} \cdot b$ .

Grows or shrinks if necessary.

Definition at line 367 of file floatmatrix.C.

References at(), givePointer(), nColumns, nRows, OOFEM_ERROR, and RESIZE.

void oofem::FloatMatrix::beTranspositionOf ( const FloatMatrix & src )

Assigns to the receiver the transposition of parameter.

Grows or shrinks if necessary.

Definition at line 323 of file floatmatrix.C.

References at(), giveNumberOfColumns(), giveNumberOfRows(), and RESIZE.

void oofem::FloatMatrix::beUnitMatrix ( )

Sets receiver to unity matrix.

Definition at line 1332 of file floatmatrix.C.

References at(), isSquare(), nColumns, nRows, OOFEM_ERROR, and zero().

Referenced by oofem::PrescribedGradientBCWeak::assemble(), oofem::XfemStructuralElementInterface::computeCohesiveTangent(), oofem::Element::computeDofTransformationMatrix(), oofem::XfemStructuralElementInterface::computeGlobalCohesiveTractionVector(), oofem::MisesMat::computeGLPlasticStrain(), oofem::RigidArmNode::computeMasterContribution(), oofem::PrescribedGradientBCNeumann::computeTangent(), oofem::StructuralInterfaceElement::computeTraction(), oofem::ElementSide::computeTransformation(), oofem::StructuralMaterial::convert_dSdE_2_dPdF(), oofem::VTKXMLExportModule::getCellVariableFromIS(), oofem::IntMatElastic::give3dStiffnessMatrix_dTdj(), oofem::NonlinearMassTransferMaterial::giveCharacteristicMatrix(), oofem::StructuralInterfaceElement::giveEngTraction(), oofem::StructuralInterfaceMaterial::giveEngTraction_3d(), oofem::AbaqusUserMaterial::giveFirstPKStressVector_3d(), oofem::IntMatCoulombContact::giveGeneralStiffnessMatrix(), oofem::CompoDamageMat::giveRealStressVector(), oofem::AbaqusUserMaterial::giveRealStressVector_3d(), oofem::tet21ghostsolid::giveRowTransformationMatrix(), oofem::StructuralMaterial::giveStiffnessMatrix(), oofem::OrthotropicLinearElasticMaterial::giveTensorRotationMatrix(), oofem::OrthotropicLinearElasticMaterial::initializeFrom(), oofem::IntMatBilinearCZFagerstromStatus::initTempStatus(), oofem::IntMatBilinearCZFagerstromStatus::IntMatBilinearCZFagerstromStatus(), oofem::IntMatBilinearCZJanssonStatus::IntMatBilinearCZJanssonStatus(), oofem::LargeStrainMasterMaterialStatus::LargeStrainMasterMaterialStatus(), oofem::LSPrimaryVariableMapper::mapPrimaryVariables(), oofem::GJacobi::solve(), and oofem::StructuralInterfaceMaterialStatus::StructuralInterfaceMaterialStatus().

void oofem::FloatMatrix::changeComponentOrder ( )

Swaps the indices in the 6x6 matrix such that it converts between OOFEM's and Abaqus' way of writing matrices.

Currently used to convert the 6x6 Jacobian from Abaqus UMAT to OOFEM.

Definition at line 1691 of file floatmatrix.C.

References at(), nColumns, and nRows.

Referenced by oofem::AbaqusUserMaterial::giveFirstPKStressVector_3d(), and oofem::AbaqusUserMaterial::giveRealStressVector_3d().

void oofem::FloatMatrix::checkBounds	(	int	i,
		int	j
	)		const

Checks size of receiver towards requested bounds.

Current implementation will call exit(1), if positions are outside bounds.

Parameters

i	Required number of rows.
j	Required number of columns.

Definition at line 189 of file floatmatrix.C.

References nColumns, nRows, and OOFEM_ERROR.

Referenced by isFinite().

void oofem::FloatMatrix::clear ( )

inline

Sets size of receiver to be an empty matrix. It will have zero rows and zero columns size.

Definition at line 516 of file floatmatrix.h.

double oofem::FloatMatrix::computeFrobeniusNorm ( ) const

Computes the Frobenius norm of the receiver.

The Frobenius norm is defined as the square root of the sum of the absolute squares of its elements.

Returns: Frobenius norm.

Definition at line 1613 of file floatmatrix.C.

References values.

Referenced by oofem::TransportGradientNeumann::computeTangent(), oofem::StructuralFE2Material::give3dMaterialStiffnessMatrix(), oofem::RVEStokesFlow::giveIPValue(), oofem::FE2FluidMaterial::giveStiffnessMatrices(), oofem::FluidMaterialEvaluator::solveYourself(), oofem::TR1_2D_SUPG2_AXI::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG2::updateStabilizationCoeffs(), oofem::TR1_2D_SUPG::updateStabilizationCoeffs(), oofem::TR21_2D_SUPG::updateStabilizationCoeffs(), and oofem::Quad10_2D_SUPG::updateStabilizationCoeffs().

double oofem::FloatMatrix::computeNorm ( char p ) const

Computes the operator norm of the receiver.

Parameters

p	Norm type, '1' for 1 norm, '2' for 2 norm.

Returns: Norm of receiver.

Todo:: Use this when obtaining eigen values is implemented.

Definition at line 1618 of file floatmatrix.C.

References at(), giveNumberOfRows(), oofem::FloatArray::givePointer(), givePointer(), oofem::max(), nColumns, oofem::norm(), nRows, and OOFEM_ERROR.

Referenced by computeReciprocalCondition().

double oofem::FloatMatrix::computeReciprocalCondition ( char p = '1' ) const

Computes the conditioning of the receiver.

From 0 to 1, where 0 is singular and 1 best. The receiver must be square. Works identically as MATLAB/Octaves rcond().

Parameters

p	Norm type, '1' for 1 norm, '2' for 2 norm.

Returns: Conditioning of receiver.

Definition at line 1737 of file floatmatrix.C.

References beInverseOf(), computeNorm(), giveDeterminant(), oofem::IntArray::givePointer(), oofem::FloatArray::givePointer(), givePointer(), isSquare(), nColumns, nRows, and OOFEM_ERROR.

Referenced by oofem::MeshQualityErrorEstimator::computeJacobianError().

void oofem::FloatMatrix::copyColumn	(	FloatArray &	dest,
		int	c
	)		const

Fetches the values from the specified column.

Output array is resized to fit column.

Parameters

dest	Array to copy values to.
c	Column to copy.

Definition at line 662 of file floatmatrix.C.

References oofem::FloatArray::begin(), giveNumberOfColumns(), giveNumberOfRows(), OOFEM_ERROR, oofem::FloatArray::resize(), and values.

Referenced by oofem::FloatArray::beColumnOf(), oofem::AbaqusUserElement::giveInternalForcesVector(), oofem::IntMatBilinearCZJanssonStatus::IntMatBilinearCZJanssonStatus(), oofem::SparseLinearSystemNM::solve(), and oofem::PetscSparseMtrx::times().

void oofem::FloatMatrix::copySubVectorRow	(	const FloatArray &	src,
		int	sr,
		int	sc
	)

Copy (set) given vector to receiver row sr, starting at column sc.

Parameters

src	Source matrix.
sr	Starting row position.
sc	Starting column position.

Definition at line 677 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), giveNumberOfColumns(), giveNumberOfRows(), oofem::FloatArray::giveSize(), oofem::max(), and resizeWithData().

Referenced by oofem::DofManager::computeM2LTransformation().

double oofem::FloatMatrix::giveDeterminant ( ) const

Returns the trace (sum of diagonal components) of the receiver.

The receiver should be a square matrix.

Returns: Trace of receiver.

Definition at line 1408 of file floatmatrix.C.

References isSquare(), nColumns, nRows, OOFEM_ERROR, and values.

int oofem::FloatMatrix::giveNumberOfColumns ( ) const

inline

Returns number of columns of receiver.

Definition at line 158 of file floatmatrix.h.

int oofem::FloatMatrix::giveNumberOfRows ( ) const

inline

Returns number of rows of receiver.

Definition at line 156 of file floatmatrix.h.

int oofem::FloatMatrix::givePackSize ( DataStream & buff ) const

Definition at line 1905 of file floatmatrix.C.

References oofem::DataStream::givePackSizeOfDouble(), oofem::DataStream::givePackSizeOfInt(), nColumns, and nRows.

Referenced by oofem::MDM::estimatePackSize().

const double* oofem::FloatMatrix::givePointer ( ) const

inline

Exposes the internal values of the matrix.

Should typically not be used outside of matrix classes.

Returns: Pointer to the values of the matrix.

Definition at line 558 of file floatmatrix.h.

Referenced by add(), addProductOf(), addTProductOf(), oofem::PetscSparseMtrx::assemble(), beInverseOf(), beMatrixFormOfStress(), oofem::FloatArray::beProductOf(), beProductOf(), beProductTOf(), oofem::FloatArray::beTProductOf(), beTProductOf(), computeNorm(), computeReciprocalCondition(), oofem::AbaqusUserMaterial::giveFirstPKStressVector_3d(), oofem::AbaqusUserElement::giveInternalForcesVector(), oofem::AbaqusUserMaterial::giveRealStressVector_3d(), oofem::PetscSolver::petsc_solve(), plusDyadSymmUpper(), plusDyadUnsym(), oofem::FloatArray::plusProduct(), plusProductSymmUpper(), plusProductUnsym(), solveForRhs(), subtract(), oofem::PetscSparseMtrx::times(), oofem::PetscSparseMtrx::timesT(), and oofem::PetscSparseMtrx::toFloatMatrix().

double* oofem::FloatMatrix::givePointer ( )

inline

Definition at line 559 of file floatmatrix.h.

References oofem::operator<<().

double oofem::FloatMatrix::giveTrace ( ) const

Returns the trace of the receiver.

The receiver should be a square matrix.

Returns: Trace of receiver.

Definition at line 1443 of file floatmatrix.C.

References isSquare(), nColumns, nRows, OOFEM_ERROR, and values.

Referenced by oofem::AnisotropicDamageMaterial::computeKappa(), oofem::AnisotropicDamageMaterial::computePlaneStressSigmaZ(), and oofem::AnisotropicDamageMaterialStatus::printOutputAt().

void oofem::FloatMatrix::hardResize	(	int	r,
		int	c
	)

Resizing that enforces reallocation of memory.

Data is zeroed.

Parameters

r	Number of rows.
c	Number of columns.

Definition at line 1396 of file floatmatrix.C.

References nColumns, nRows, and values.

void oofem::FloatMatrix::initFromVector	(	const FloatArray &	vector,
		bool	transposed
	)

Assigns to receiver one column or one row matrix containing vector.

Parameters

vector	Source vector.
transposed	If false then (vector->giveSize(),1) FloatMatrix is assigned else (1,vector->giveSize()) FloatMatrix is assigned.

Definition at line 1308 of file floatmatrix.C.

References oofem::FloatArray::giveSize(), nColumns, nRows, oofem::FloatArray::values, and values.

Referenced by oofem::PlasticMaterial::giveRealStressVector().

bool oofem::FloatMatrix::isFinite ( ) const

Returns true if no element is NAN or infinite.

Definition at line 209 of file floatmatrix.C.

References at(), checkBounds(), nRows, operator()(), and values.

Referenced by oofem::LSPrimaryVariableMapper::mapPrimaryVariables().

bool oofem::FloatMatrix::isNotEmpty ( ) const

inline

Tests for empty matrix.

Definition at line 162 of file floatmatrix.h.

Referenced by add(), oofem::EngngModel::assemble(), oofem::NonStationaryTransportProblem::assembleAlgorithmicPartOfRhs(), oofem::DIIDynamic::assembleDirichletBcRhsVector(), oofem::EngngModel::assembleExtrapolatedForces(), oofem::EngngModel::assemblePrescribedExtrapolatedForces(), oofem::TransportElement::computeBCMtrxAt(), oofem::TrPlaneStrRot3d::computeGtoLRotationMatrix(), oofem::CCTPlate3d::computeGtoLRotationMatrix(), oofem::DKTPlate3d::computeGtoLRotationMatrix(), oofem::MITC4Shell::computeGtoLRotationMatrix(), oofem::LIBeam2dNL::computeInitialStressMatrix(), oofem::TrPlaneStrRot3d::computeLoadGToLRotationMtrx(), oofem::NlDEIDynamic::computeMassMtrx(), oofem::TransportGradientPeriodic::computeTangent(), oofem::TransportGradientDirichlet::computeTangent(), oofem::SUPGElement::giveCharacteristicMatrix(), oofem::TR_SHELL02::giveCharacteristicMatrix(), oofem::TR_SHELL01::giveCharacteristicMatrix(), oofem::Element::giveGlobalIPValue(), oofem::TrPlaneStrRot3d::giveLocalCoordinates(), oofem::Element::giveLocalCoordinateSystem(), oofem::Structural3DElement::giveMaterialOrientationAt(), oofem::Structural2DElement::giveMaterialOrientationAt(), oofem::Element::giveRotationMatrix(), oofem::SUPGTangentAssembler::matrixFromElement(), plusDyadSymmUpper(), plusDyadUnsym(), plusProductSymmUpper(), plusProductUnsym(), oofem::AbaqusUserElement::postInitialize(), oofem::DEIDynamic::solveYourselfAt(), subtract(), oofem::DIIDynamic::timesMtrx(), and oofem::SUPGInternalForceAssembler::vectorFromElement().

bool oofem::FloatMatrix::isSquare ( ) const

inline

Returns nonzero if receiver is square matrix.

Definition at line 160 of file floatmatrix.h.

Referenced by assemble(), beInverseOf(), beUnitMatrix(), computeReciprocalCondition(), oofem::StructuralElement::condense(), giveDeterminant(), giveTrace(), jaco_(), oofem::GJacobi::solve(), solveForRhs(), and oofem::StructuralMaterial::sortPrincDirAndValCloseTo().

bool oofem::FloatMatrix::jaco_	(	FloatArray &	eval,
		FloatMatrix &	v,
		int	nf
	)

Computes eigenvalues and eigenvectors of receiver (must be symmetric) The receiver is preserved.

Parameters

eval	Requested eigenvalues.
v	Requested eigenvectors (stored colum wise).
nf	Number of significant figures.

Returns: True if ok,otherwise false.

Definition at line 1912 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), giveNumberOfRows(), isSquare(), OOFEM_ERROR, oofem::FloatArray::resize(), and resize().

Referenced by oofem::AnisotropicDamageMaterial::computeDamageTensor(), oofem::MDM::computePDC(), oofem::AnisotropicDamageMaterial::computePlaneStressSigmaZ(), oofem::AnisotropicDamageMaterial::computePlaneStressStrain(), oofem::StrainVector::computePrincipalValDir(), oofem::StressVector::computePrincipalValDir(), oofem::StructuralMaterial::computePrincipalValDir(), oofem::AnisotropicDamageMaterial::computeSecantOperator(), oofem::LargeStrainMasterMaterial::giveFirstPKStressVector_3d(), oofem::LargeStrainMasterMaterialGrad::giveFirstPKStressVectorGrad(), oofem::AnisotropicDamageMaterial::giveIPValue(), oofem::AnisotropicDamageMaterial::giveRealStressVector(), oofem::AnisotropicDamageMaterial::obtainAlpha1(), oofem::AnisotropicDamageMaterial::obtainAlpha2(), and oofem::AnisotropicDamageMaterial::obtainAlpha3().

void oofem::FloatMatrix::negated ( )

Changes sign of receiver values.

Definition at line 1605 of file floatmatrix.C.

References values.

Referenced by oofem::PrescribedGradientBCNeumann::assemble(), oofem::MixedGradientPressureWeakPeriodic::assemble(), oofem::MixedGradientPressureNeumann::assemble(), oofem::Node2NodeContact::computeContactTangent(), oofem::ContactDefinition::computeContactTangent(), oofem::TR1_2D_SUPG_AXI::computeOutFlowBCTerm_MB(), oofem::TR1_2D_SUPG::computeOutFlowBCTerm_MB(), oofem::PrescribedGradientBCNeumann::computeTangent(), oofem::TransportGradientNeumann::computeTangent(), oofem::TransportGradientDirichlet::computeXi(), oofem::ContactElement::ContactElement(), oofem::MPlasticMaterial2::giveConsistentStiffnessMatrix(), oofem::MPlasticMaterial::giveElastoPlasticStiffnessMatrix(), oofem::MPlasticMaterial2::giveElastoPlasticStiffnessMatrix(), oofem::LinearElasticMaterial::giveEshelbyStressVector_PlaneStrain(), oofem::AbaqusUserElement::giveInternalForcesVector(), oofem::FETISolver::solve(), and subtract().

double& oofem::FloatMatrix::operator()	(	int	i,
		int	j
	)

inline

Coefficient access function.

Returns l-value of coefficient at given position of the receiver. Implements 0-based indexing.

Parameters

i	Row position of coefficient.
j	Column position of coefficient.

Definition at line 199 of file floatmatrix.h.

Referenced by isFinite().

double oofem::FloatMatrix::operator()	(	int	i,
		int	j
	)		const

inline

Coefficient access function.

Implements 0-based indexing.

Parameters

i	Row position of coefficient.
j	Column position of coefficient.

Definition at line 209 of file floatmatrix.h.

FloatMatrix & oofem::FloatMatrix::operator= ( std::initializer_list< std::initializer_list< double > > mat )

Assignment operator.

Definition at line 149 of file floatmatrix.C.

References nRows, OOFEM_ERROR, RESIZE, and values.

Referenced by add(), and subtract().

FloatMatrix & oofem::FloatMatrix::operator= ( std::initializer_list< FloatArray > mat )

Assignment operator.

Definition at line 169 of file floatmatrix.C.

References nRows, OOFEM_ERROR, RESIZE, and values.

FloatMatrix& oofem::FloatMatrix::operator= ( const FloatMatrix & mat )

inline

Assignment operator, adjusts size of the receiver if necessary.

Definition at line 133 of file floatmatrix.h.

References nColumns, nRows, and values.

FloatMatrix& oofem::FloatMatrix::operator= ( FloatMatrix && mat )

inline

Definition at line 139 of file floatmatrix.h.

void oofem::FloatMatrix::plusDyadSymmUpper	(	const FloatArray &	a,
		double	dV
	)

Adds to the receiver the dyadic product $a \otimes a \mathrm{d}V$ .

If the receiver has zero size, it is expanded. Computes only the upper half of receiver.

Parameters

a	Array a in equation.
dV	Scaling factor.

Definition at line 756 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), oofem::FloatArray::givePointer(), givePointer(), oofem::FloatArray::giveSize(), isNotEmpty(), nColumns, nRows, and values.

Referenced by oofem::TransportElement::computeBCSubMtrxAt(), oofem::TransportElement::computeCapacitySubMatrix(), oofem::IsotropicDamageMaterial1::computeEta(), oofem::TransportElement::computeIntSourceLHSSubMatrix(), oofem::TR1_2D_SUPG_AXI::computePenetrationWithResistanceBCTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeSlipWithFrictionBCTerm_MB(), oofem::Hexa21Stokes::computeStiffnessMatrix(), oofem::Tet21Stokes::computeStiffnessMatrix(), oofem::Tr21Stokes::computeStiffnessMatrix(), oofem::Tr1BubbleStokes::computeStiffnessMatrix(), oofem::Tet1BubbleStokes::computeStiffnessMatrix(), oofem::Quad1MindlinShell3D::computeStiffnessMatrix(), oofem::SurfaceTensionBoundaryCondition::computeTangentFromElement(), oofem::TransportGradientDirichlet::computeXi(), oofem::Shell7Base::NodalRecoveryMI_computeNNMatrix(), oofem::Tr_Warp::ZZNodalRecoveryMI_computeNNMatrix(), and oofem::ZZNodalRecoveryModelInterface::ZZNodalRecoveryMI_computeNNMatrix().

void oofem::FloatMatrix::plusDyadUnsym	(	const FloatArray &	a,
		const FloatArray &	b,
		double	dV
	)

Adds to the receiver the product $a \otimes b \mathrm{d}V$ .

If the receiver has zero size, it is expanded.

Parameters

a	Array a in equation.
b	Array b in equation.
dV	Scaling factor.

Definition at line 811 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), oofem::FloatArray::givePointer(), givePointer(), oofem::FloatArray::giveSize(), isNotEmpty(), nColumns, nRows, and values.

Referenced by oofem::ConcreteDPM2::computeDDRhoDDStress(), oofem::tet21ghostsolid::computeLoadVector(), oofem::TR1_2D_SUPG_AXI::computeOutFlowBCTerm_MB(), oofem::tet21ghostsolid::computeStiffnessMatrix(), oofem::Hexa21Stokes::computeStiffnessMatrix(), oofem::Tet21Stokes::computeStiffnessMatrix(), oofem::Tr21Stokes::computeStiffnessMatrix(), oofem::Tr1BubbleStokes::computeStiffnessMatrix(), oofem::Tet1BubbleStokes::computeStiffnessMatrix(), oofem::MixedGradientPressureWeakPeriodic::computeStress(), oofem::SurfaceTensionBoundaryCondition::computeTangentFromElement(), oofem::TransportGradientDirichlet::computeXi(), oofem::TutorialMaterial::give3dMaterialStiffnessMatrix(), oofem::GradDpElement::giveInternalForcesVector(), oofem::TransportGradientNeumann::integrateTangent(), oofem::Shell7Base::NodalRecoveryMI_computeNValProduct(), oofem::TrabBoneNL3D::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::MisesMatNl::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::RankineMatNl::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::IDNLMaterial::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::Tr_Warp::ZZNodalRecoveryMI_computeNValProduct(), and oofem::ZZNodalRecoveryModelInterface::ZZNodalRecoveryMI_computeNValProduct().

void oofem::FloatMatrix::plusProductSymmUpper	(	const FloatMatrix &	a,
		const FloatMatrix &	b,
		double	dV
	)

Adds to the receiver the product $a^{\mathrm{T}}\cdot b \mathrm{d}V$ .

If the receiver has zero size, it is expanded. Assumes that receiver and product $a^{\mathrm{T}}\cdot b \mathrm{d}V$ are symmetric matrices. Computes only the upper half of receiver.

Parameters

a	Matrix a in equation.
b	Matrix b in equation.
dV	Scaling factor.

Definition at line 698 of file floatmatrix.C.

References at(), givePointer(), isNotEmpty(), nColumns, nRows, and values.

void oofem::FloatMatrix::plusProductUnsym	(	const FloatMatrix &	a,
		const FloatMatrix &	b,
		double	dV
	)

Adds to the receiver the product $a^{\mathrm{T}} \cdot b \mathrm{d}V$ .

If the receiver has zero size, it is expanded.

Parameters

a	Matrix a in equation.
b	Matrix b in equation.
dV	Scaling factor.

Definition at line 780 of file floatmatrix.C.

References at(), givePointer(), isNotEmpty(), nColumns, nRows, and values.

void oofem::FloatMatrix::printYourself ( ) const

Prints matrix to stdout. Useful for debugging.

Definition at line 1458 of file floatmatrix.C.

References at(), nColumns, and nRows.

Referenced by oofem::WeakPeriodicBoundaryCondition::computeElementTangent(), oofem::StructuralInterfaceMaterial::give2dStiffnessMatrix_Eng(), oofem::AbaqusUserMaterial::give3dMaterialStiffnessMatrix(), oofem::Shell7Base::giveLayerContributionToSR(), oofem::FE2FluidMaterial::giveStiffnessMatrices(), oofem::SkylineUnsym::printYourself(), oofem::Skyline::printYourself(), oofem::Shell7BaseXFEM::recoverShearStress(), oofem::GJacobi::solve(), oofem::SubspaceIteration::solve(), and oofem::StructuralMaterialEvaluator::solveYourself().

void oofem::FloatMatrix::printYourself ( const std::string & name ) const

Print receiver on stdout with custom name.

Parameters

name	Display name of reciever.

Definition at line 1499 of file floatmatrix.C.

References at(), nColumns, and nRows.

void oofem::FloatMatrix::printYourselfToFile	(	const std::string	filename,
		const bool	showDimensions = `true`
	)		const

Print matrix to file.

Parameters

filename Output filename showDimensions Determines if dimensions should be included in output. Default is true.

Definition at line 1477 of file floatmatrix.C.

References at(), nColumns, nRows, and OOFEM_ERROR.

Referenced by oofem::tet21ghostsolid::computeStiffnessMatrix().

void oofem::FloatMatrix::pY ( ) const

Higher accuracy than printYourself.

Definition at line 1524 of file floatmatrix.C.

References at(), nColumns, and nRows.

void oofem::FloatMatrix::resize	(	int	rows,
		int	cols
	)

Checks size of receiver towards requested bounds.

If dimension mismatch, size is adjusted accordingly. Content is always zeroed.

Parameters

rows	New number of rows.
cols	New number of columns.

Definition at line 1358 of file floatmatrix.C.

References nColumns, nRows, and values.

void oofem::FloatMatrix::resizeWithData	(	int	rows,
		int	columns
	)

Checks size of receiver towards requested bounds.

If dimension mismatch, size is adjusted accordingly. Note: New coefficients are initialized to zero, old are kept.

Definition at line 1369 of file floatmatrix.C.

References at(), giveNumberOfColumns(), giveNumberOfRows(), oofem::min(), nColumns, nRows, and values.

Referenced by addSubVectorCol(), addSubVectorRow(), oofem::Element::computeDofTransformationMatrix(), oofem::PrescribedGradientBCWeak::computeExtForceElContrib(), copySubVectorRow(), oofem::PrescribedGradient::give(), oofem::FCMMaterial::giveMaterialStiffnessMatrix(), oofem::SimpleCrossSection::giveMembraneRotStiffMtrx(), and oofem::Tr_Warp::ZZNodalRecoveryMI_computeNValProduct().

contextIOResultType oofem::FloatMatrix::restoreYourself ( DataStream & stream )

Definition at line 1877 of file floatmatrix.C.

References oofem::CIO_IOERR, oofem::CIO_OK, nColumns, nRows, oofem::DataStream::read(), and values.

Referenced by oofem::EigenValueDynamic::restoreContext(), oofem::RCSDEMaterialStatus::restoreContext(), oofem::RCSDMaterialStatus::restoreContext(), oofem::LinearStability::restoreContext(), oofem::LIBeam3d2::restoreContext(), oofem::MDMStatus::restoreContext(), oofem::RCM2MaterialStatus::restoreContext(), oofem::AnisotropicDamageMaterialStatus::restoreContext(), oofem::Node::restoreContext(), oofem::FCMMaterialStatus::restoreContext(), and oofem::MDM::unpackAndUpdateUnknowns().

void oofem::FloatMatrix::rotatedWith	(	const FloatMatrix &	r,
		char	mode = `'n'`
	)

Returns the receiver 'a' transformed using give transformation matrix r.

The method performs the operation $a = r^{\mathrm{T}} \cdot a \cdot r$ .

Parameters

r	Transformation matrix.
mode	If set to 't' then the transpose of the rotation matrix is used instead.

Definition at line 1557 of file floatmatrix.C.

References beProductOf(), beProductTOf(), beTProductOf(), and OOFEM_ERROR.

Referenced by oofem::EngngModel::assemble(), oofem::Shell7BaseXFEM::computeCohesiveForces(), oofem::Shell7BaseXFEM::computeCohesiveTangentAt(), oofem::Structural3DElement::computeConstitutiveMatrixAt(), oofem::PlaneStressElement::computeConstitutiveMatrixAt(), oofem::PlaneStrainElement::computeConstitutiveMatrixAt(), oofem::AxisymElement::computeConstitutiveMatrixAt(), oofem::NlDEIDynamic::computeMassMtrx(), oofem::IntElLine1IntPen::computeStiffnessMatrix(), oofem::IntElLine2IntPen::computeStiffnessMatrix(), oofem::StructuralInterfaceElement::computeStiffnessMatrix(), oofem::StructuralInterfaceElementPhF::computeStiffnessMatrix_uu(), oofem::IntElLine1PF::computeStiffnessMatrix_uu(), oofem::TransportGradientPeriodic::computeTangent(), oofem::TransportGradientDirichlet::computeTangent(), oofem::StructuralInterfaceElement::computeTraction(), oofem::RankineMat::evaluatePlaneStressStiffMtrx(), oofem::LayeredCrossSection::give2dPlateStiffMtrx(), oofem::WinklerMaterial::give3dBeamSubSoilStiffMtrx(), oofem::OrthotropicLinearElasticMaterial::give3dMaterialStiffnessMatrix(), oofem::CompoDamageMat::give3dMaterialStiffnessMatrix(), oofem::LayeredCrossSection::give3dShellStiffMtrx(), oofem::TrPlaneStrRot3d::giveCharacteristicTensor(), oofem::RerShell::giveCharacteristicTensor(), oofem::TrPlanestressRotAllman3d::giveCharacteristicTensor(), oofem::DKTPlate3d::giveCharacteristicTensor(), oofem::CCTPlate3d::giveCharacteristicTensor(), oofem::LinQuad3DPlaneStress::giveCharacteristicTensor(), oofem::RCM2Material::giveEffectiveMaterialStiffnessMatrix(), oofem::Element::giveGlobalIPValue(), oofem::FCMMaterial::giveMaterialStiffnessMatrix(), oofem::Shell7BaseXFEM::giveRecoveredTransverseInterfaceStress(), oofem::LayeredCrossSection::giveStiffnessMatrix_3d(), oofem::Shell7BaseXFEM::recoverShearStress(), oofem::DEIDynamic::solveYourselfAt(), oofem::DIIDynamic::timesMtrx(), oofem::TR_SHELL01::ZZErrorEstimatorI_computeLocalStress(), and oofem::TR_SHELL02::ZZErrorEstimatorI_computeLocalStress().

void oofem::FloatMatrix::setColumn	(	const FloatArray &	src,
		int	c
	)

Sets the values of the matrix in specified column.

If matrix size is zero, the size is adjusted.

Parameters

src	Array to set at column c.
c	Column to set.

Definition at line 648 of file floatmatrix.C.

References oofem::FloatArray::begin(), oofem::FloatArray::end(), giveNumberOfColumns(), giveNumberOfRows(), oofem::FloatArray::giveSize(), OOFEM_ERROR, and values.

void oofem::FloatMatrix::setSubMatrix	(	const FloatMatrix &	src,
		int	sr,
		int	sc
	)

Adds the given matrix as sub-matrix to receiver.

The sub-matrix values will be added to receivers corresponding receiver's values at positions (ri...ri+src.nrows, ci....ci+src.ncolumns). The size of receiver will be adjusted, if necessary.

Parameters

src	Sub-matrix to be inserted.
sr	Determines the row position (receiver's 1-based index) of first src value to be added.
sc	Determines the column position (receiver's 1-based index) of first src value to be added.

Definition at line 557 of file floatmatrix.C.

References giveNumberOfColumns(), giveNumberOfRows(), and OOFEM_ERROR.

Referenced by oofem::XfemElementInterface::ComputeBOrBHMatrix(), oofem::TransportElement::computeFlow(), oofem::MITC4Shell::computeGtoLRotationMatrix(), oofem::Quad1MindlinShell3D::computeGtoLRotationMatrix(), oofem::Shell7Base::computeLambdaGMatrices(), oofem::Shell7BaseXFEM::computeLambdaGMatricesDis(), and oofem::Shell7Base::giveLayerContributionToSR().

void oofem::FloatMatrix::setTSubMatrix	(	const FloatMatrix &	src,
		int	sr,
		int	sc
	)

Adds the given matrix as sub-matrix to receiver.

The same as setSubMatrix but transposes the source matrix.

Parameters

src	Sub-matrix to be inserted transposed.
sr	Determines the row position (receiver's 1-based index) of first src value to be added.
sc	Determines the column position (receiver's 1-based index) of first src value to be added.

Definition at line 582 of file floatmatrix.C.

References giveNumberOfColumns(), giveNumberOfRows(), and OOFEM_ERROR.

bool oofem::FloatMatrix::solveForRhs	(	const FloatArray &	b,
		FloatArray &	answer,
		bool	transpose = `false`
	)

Solves the system of linear equations $K\cdot a = b$ .

Uses Gaussian elimination with pivoting directly on receiver.

Parameters

b	RHS of linear system.
answer	Solution of linear equations.
transpose	Solves for the transpose of K.

Returns: False if K is singular, otherwise true.

Definition at line 1112 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), beTranspositionOf(), oofem::IntArray::givePointer(), oofem::FloatArray::givePointer(), givePointer(), oofem::FloatArray::giveSize(), isSquare(), nColumns, nRows, and OOFEM_ERROR.

Referenced by oofem::MMALeastSquareProjection::__mapVariable(), oofem::MPlasticMaterial::closestPointReturn(), oofem::MPlasticMaterial2::closestPointReturn(), oofem::KelvinChainMaterial::computeCharCoefficients(), oofem::MaxwellChainMaterial::computeCharCoefficients(), oofem::TransportGradientNeumann::computeEta(), oofem::SPRNodalRecoveryModel::computePatch(), oofem::MPlasticMaterial::cuttingPlaneReturn(), oofem::MPlasticMaterial2::cuttingPlaneReturn(), oofem::LEPlic::doCellDLS(), oofem::StructuralMaterial::giveFirstPKStressVector_1d(), oofem::StructuralMaterial::giveFirstPKStressVector_PlaneStress(), oofem::RCM2Material::giveRealPrincipalStressVector3d(), oofem::FCMMaterial::giveRealStressVector(), oofem::StructuralMaterial::giveRealStressVector_ShellStressControl(), oofem::StructuralMaterial::giveRealStressVector_StressControl(), oofem::DruckerPragerPlasticitySM::giveRegAlgorithmicStiffMatrix(), oofem::FEI3dHexaLin::global2local(), oofem::FEI3dTetQuad::global2local(), oofem::FEI3dWedgeLin::global2local(), oofem::FEInterpolation2d::global2local(), oofem::FEI3dHexaQuad::global2local(), oofem::LSpace::NodalAveragingRecoveryMI_computeNodalValue(), oofem::MITC4Shell::NodalAveragingRecoveryMI_computeNodalValue(), oofem::ConcreteDPM2::performRegularReturn(), oofem::ConcreteDPMStatus::restoreConsistency(), oofem::StructuralMaterialEvaluator::solveYourself(), and oofem::FluidMaterialEvaluator::solveYourself().

void oofem::FloatMatrix::solveForRhs	(	const FloatMatrix &	B,
		FloatMatrix &	answer,
		bool	transpose = `false`
	)

Solves the system of linear equations $K\cdot A = B$ .

Uses Gaussian elimination with pivoting directly on receiver.

Parameters

B	RHS of linear system.
answer	Solution of linear equations, each column corresponding to columns in B.
transpose	Solves for the transpose of K.

Definition at line 1204 of file floatmatrix.C.

References at(), beTranspositionOf(), giveNumberOfColumns(), giveNumberOfRows(), oofem::IntArray::givePointer(), givePointer(), isSquare(), nColumns, nRows, and OOFEM_ERROR.

contextIOResultType oofem::FloatMatrix::storeYourself ( DataStream & stream ) const

Definition at line 1852 of file floatmatrix.C.

References oofem::CIO_IOERR, oofem::CIO_OK, nColumns, and nRows.

Referenced by oofem::MDM::packUnknowns(), oofem::EigenValueDynamic::saveContext(), oofem::RCSDEMaterialStatus::saveContext(), oofem::RCSDMaterialStatus::saveContext(), oofem::LinearStability::saveContext(), oofem::LIBeam3d2::saveContext(), oofem::MDMStatus::saveContext(), oofem::RCM2MaterialStatus::saveContext(), oofem::AnisotropicDamageMaterialStatus::saveContext(), oofem::FCMMaterialStatus::saveContext(), and oofem::Node::saveContext().

void oofem::FloatMatrix::subtract ( const FloatMatrix & a )

Subtracts matrix from the receiver.

Parameters

a	Matrix to be subtracted.

Definition at line 1084 of file floatmatrix.C.

References givePointer(), isNotEmpty(), nColumns, negated(), nRows, OOFEM_ERROR, operator=(), and values.

Referenced by oofem::Shell7BaseXFEM::computeCohesiveNmatrixAt(), oofem::AnisotropicDamageMaterial::computeDamageTensor(), oofem::MisesMat::computeGLPlasticStrain(), oofem::Shell7Base::computePressureTangentMatrix(), oofem::Shell7BaseXFEM::computePressureTangentMatrixDis(), oofem::AnisotropicDamageMaterial::computeSecantOperator(), oofem::PrescribedGradientBCPeriodic::computeTangent(), oofem::TransportGradientPeriodic::computeTangent(), oofem::PrescribedGradient::computeTangent(), oofem::PrescribedGenStrainShell7::computeTangent(), oofem::MixedGradientPressureDirichlet::computeTangents(), oofem::RankineMat::evaluatePlaneStressStiffMtrx(), oofem::StructuralInterfaceMaterial::give2dStiffnessMatrix_Eng(), oofem::TutorialMaterial::give3dMaterialStiffnessMatrix(), oofem::StructuralFE2Material::give3dMaterialStiffnessMatrix(), oofem::IntMatBilinearCZJansson::give3dStiffnessMatrix_dTdj(), oofem::IntMatBilinearCZFagerstrom::give3dStiffnessMatrix_dTdj(), oofem::MPlasticMaterial::giveConsistentStiffnessMatrix(), oofem::FCMMaterial::giveMaterialStiffnessMatrix(), oofem::AnisotropicDamageMaterial::giveRealStressVector(), and oofem::FE2FluidMaterial::giveStiffnessMatrices().

void oofem::FloatMatrix::symmetrized ( )

Initializes the lower half of the receiver according to the upper half.

Definition at line 1576 of file floatmatrix.C.

References at(), nColumns, nRows, and OOFEM_ERROR.

void oofem::FloatMatrix::times ( double f )

Multiplies receiver by factor f.

Parameters

f	Factor to multiply by.

Definition at line 1594 of file floatmatrix.C.

References values.

void oofem::FloatMatrix::writeCSV ( const std::string & name ) const

Writes receiver as CSV (comma seperated values)

Parameters

name Filename

Definition at line 1543 of file floatmatrix.C.

References at(), nColumns, and nRows.

void oofem::FloatMatrix::zero ( )

Zeroes all coefficient of receiver.

Definition at line 1326 of file floatmatrix.C.

References values.

Friends And Related Function Documentation

std :: ostream& operator<<	(	std::ostream &	out,
		const FloatMatrix &	r
	)

friend

Definition at line 2098 of file floatmatrix.C.

Member Data Documentation

int oofem::FloatMatrix::nColumns

protected

Number of columns.

Definition at line 100 of file floatmatrix.h.

Referenced by add(), addProductOf(), addTProductOf(), beInverseOf(), beProductOf(), beProductTOf(), beSubMatrixOf(), beTProductOf(), beUnitMatrix(), changeComponentOrder(), checkBounds(), computeNorm(), computeReciprocalCondition(), FloatMatrix(), giveDeterminant(), givePackSize(), giveTrace(), hardResize(), initFromVector(), oofem::operator<<(), operator=(), plusDyadSymmUpper(), plusDyadUnsym(), plusProductSymmUpper(), plusProductUnsym(), printYourself(), printYourselfToFile(), pY(), resize(), resizeWithData(), restoreYourself(), solveForRhs(), storeYourself(), subtract(), symmetrized(), and writeCSV().

int oofem::FloatMatrix::nRows

protected

std :: vector< double > oofem::FloatMatrix::values

protected

Values of matrix stored column wise.

Definition at line 102 of file floatmatrix.h.

Referenced by add(), bePinvID(), computeFrobeniusNorm(), copyColumn(), FloatMatrix(), giveDeterminant(), giveTrace(), hardResize(), initFromVector(), isFinite(), negated(), operator=(), plusDyadSymmUpper(), plusDyadUnsym(), plusProductSymmUpper(), plusProductUnsym(), resize(), resizeWithData(), restoreYourself(), setColumn(), subtract(), times(), and zero().

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

Public Member Functions

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation