#include <floatmatrix.h>

Public Types
typedef double	Scalar

Public Member Functions
	FloatMatrix (std::size_t n, std::size_t m)
	FloatMatrix ()
	Creates zero sized matrix.
	FloatMatrix (const FloatMatrix &mat)
	Copy constructor.
	FloatMatrix (FloatMatrix &&mat) noexcept
	Copy constructor.
FloatMatrix &	operator= (const FloatMatrix &mat)
	Assignment operator, adjusts size of the receiver if necessary.
FloatMatrix &	operator= (FloatMatrix &&mat) noexcept
int	rows () const
int	cols () const
double &	operator() (std::size_t i, std::size_t j)
double	operator() (std::size_t i, std::size_t j) const
const double *	data () const
double *	data ()
template<std::size_t N, std::size_t M>
void	assignFloatMatrixF (const FloatMatrixF< N, M > &mat)
void	_resize_internal (int nr, int nc)
	FloatMatrix (std ::initializer_list< std ::initializer_list< double > > ini)
FloatMatrix &	operator= (std ::initializer_list< std ::initializer_list< double > >mat)
template<std::size_t M, std::size_t N>
	FloatMatrix (const FloatMatrixF< N, M > &src)
template<std::size_t N, std::size_t M>
FloatMatrix &	operator= (const FloatMatrixF< N, M > &mat)
void	checkBounds (Index i, Index j) const
int	giveNumberOfRows () const
	Returns number of rows of receiver.
int	giveNumberOfColumns () const
	Returns number of columns of receiver.
bool	isSquare () const
	Returns nonzero if receiver is square matrix.
bool	isNotEmpty () const
	Tests for empty matrix.
bool	isAllFinite () const
	Returns true if no element is NAN or infinite.
double	at (std::size_t i, std::size_t j) const
double &	at (std::size_t i, std::size_t j)
void	assemble (const FloatMatrix &src, const IntArray &loc)
void	assemble (const FloatMatrix &src, const IntArray &rowind, const IntArray &colind)
void	assembleT (const FloatMatrix &src, const IntArray &rowind, const IntArray &colind)
double	computeFrobeniusNorm () const
double	computeNorm (char p) const
double	computeReciprocalCondition (char p='1') const
void	beDiagonal (const FloatArray &diag)
double	giveTrace () const
double	giveDeterminant () const
void	zero ()
	Zeroes all coefficient of receiver.
void	beUnitMatrix ()
	Sets receiver to unity matrix.
void	bePinvID ()
	Sets receiver to the inverse of scaling matrix P multiplied by the deviatoric projector ID.
void	beTranspositionOf (const FloatMatrix &src)
void	beProductOf (const FloatMatrix &a, const FloatMatrix &b)
void	addProductOf (const FloatMatrix &a, const FloatMatrix &b)
void	addTProductOf (const FloatMatrix &a, const FloatMatrix &b)
void	beTProductOf (const FloatMatrix &a, const FloatMatrix &b)
void	beProductTOf (const FloatMatrix &a, const FloatMatrix &b)
void	beDyadicProductOf (const FloatArray &vec1, const FloatArray &vec2)
void	beNMatrixOf (const FloatArray &n, int nsd)
void	beLocalCoordSys (const FloatArray &normal)
void	setSubMatrix (const FloatMatrix &src, int sr, int sc)
void	beSubMatrixOf (const FloatMatrix &src, Index topRow, Index bottomRow, Index topCol, Index bottomCol)
void	beSubMatrixOf (const FloatMatrix &src, const IntArray &indxRow, const IntArray &indxCol)
void	addSubVectorRow (const FloatArray &src, int sr, int sc)
void	addSubVectorCol (const FloatArray &src, int sr, int sc)
void	copySubVectorRow (const FloatArray &src, int sr, int sc)
void	setColumn (const FloatArray &src, int c)
void	copyColumn (FloatArray &dest, int c) const
bool	beInverseOf (const FloatMatrix &src)
bool	solveForRhs (const FloatArray &b, FloatArray &answer, bool transpose=false)
bool	solveForRhs (const FloatMatrix &B, FloatMatrix &answer, bool transpose=false)
void	plusProductSymmUpper (const FloatMatrix &a, const FloatMatrix &b, double dV)
void	plusDyadSymmUpper (const FloatArray &a, double dV)
void	plusProductUnsym (const FloatMatrix &a, const FloatMatrix &b, double dV)
void	plusDyadUnsym (const FloatArray &a, const FloatArray &b, double dV)
void	plus_Nt_a_otimes_b_B (const FloatMatrix &N, const FloatArray &a, const FloatArray &b, const FloatMatrix &B, double dV=1)
void	add (const FloatMatrix &a)
void	add (double s, const FloatMatrix &a)
void	subtract (const FloatMatrix &a)
void	times (double f)
void	negated ()
void	symmetrized ()
void	rotatedWith (const FloatMatrix &r, char mode='n')
void	resize (Index rows, Index cols)
void	resizeWithData (Index, Index)
*Sets size of receiver to be an empty matrix It will have zero rows and zero columns size void	clear ()
bool	jaco_ (FloatArray &eval, FloatMatrix &v, int nf)
*Prints matrix to stdout Useful for debugging void	printYourself () const
void	printYourself (const std::string &name) const
void	printYourselfToFile (const std::string filename, const bool showDimensions=true) const
*Higher accuracy than printYourself void	pY () const
void	writeCSV (const std ::string &name) const
const double *	givePointer () const
double *	givePointer ()
void	beMatrixFormOfStress (const FloatArray &aArray)
void	beMatrixForm (const FloatArray &aArray)
void	changeComponentOrder ()
contextIOResultType	storeYourself (DataStream &stream) const
contextIOResultType	restoreYourself (DataStream &stream)
int	givePackSize (DataStream &buff) const

Static Public Member Functions
static FloatMatrix	fromIniList (std ::initializer_list< std ::initializer_list< double > >)
static FloatMatrix	fromCols (std ::initializer_list< FloatArray >mat)
static FloatMatrix	fromArray (const FloatArray &vector, bool transpose=false)

Static Public Attributes
static constexpr int	Dim = 2

Protected Attributes
std::size_t	nRows
	Number of rows.
std::size_t	nColumns
	Number of columns.
std ::vector< double >	values
	Values of matrix stored column wise.

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 83 of file floatmatrix.h.

Member Typedef Documentation

◆ Scalar

typedef double oofem::FloatMatrix::Scalar

Definition at line 103 of file floatmatrix.h.

Constructor & Destructor Documentation

◆ FloatMatrix() [1/6]

oofem::FloatMatrix::FloatMatrix	(	std::size_t	n,
		std::size_t	m )

inline

Creates matrix of given size.

Parameters

n	Number of rows.
m	Requested number of columns.

Definition at line 109 of file floatmatrix.h.

References nColumns, nRows, and values.

Referenced by add(), add(), addProductOf(), addTProductOf(), assemble(), assemble(), assembleT(), beInverseOf(), bePinvID(), beProductOf(), beProductTOf(), beSubMatrixOf(), beSubMatrixOf(), beTProductOf(), beTranspositionOf(), changeComponentOrder(), computeReciprocalCondition(), FloatMatrix(), FloatMatrix(), fromArray(), fromCols(), fromIniList(), giveDeterminant(), jaco_(), operator<<, operator=(), operator=(), operator=(), operator=(), plus_Nt_a_otimes_b_B(), plusProductSymmUpper(), plusProductUnsym(), resizeWithData(), rotatedWith(), setSubMatrix(), solveForRhs(), solveForRhs(), and subtract().

◆ FloatMatrix() [2/6]

oofem::FloatMatrix::FloatMatrix ( )

inline

Creates zero sized matrix.

Definition at line 111 of file floatmatrix.h.

References nColumns, nRows, and values.

◆ FloatMatrix() [3/6]

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

inline

Copy constructor.

Definition at line 113 of file floatmatrix.h.

References FloatMatrix(), nColumns, nRows, and values.

◆ FloatMatrix() [4/6]

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

inlinenoexcept

Copy constructor.

Definition at line 115 of file floatmatrix.h.

References FloatMatrix(), nColumns, nRows, and values.

◆ FloatMatrix() [5/6]

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

inline

Definition at line 168 of file floatmatrix.h.

References fromIniList().

◆ FloatMatrix() [6/6]

template<std::size_t M, std::size_t N>

oofem::FloatMatrix::FloatMatrix ( const FloatMatrixF< N, M > & src )

inline

Definition at line 173 of file floatmatrix.h.

References assignFloatMatrixF().

Member Function Documentation

◆ _resize_internal()

void oofem::FloatMatrix::_resize_internal	(	int	nr,
		int	nc )

Definition at line 68 of file floatmatrix.C.

References nColumns, nRows, and values.

Referenced by beDyadicProductOf(), beInverseOf(), beProductOf(), beProductTOf(), beTProductOf(), beTranspositionOf(), and fromIniList().

◆ add() [1/2]

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 962 of file floatmatrix.C.

References _LOOP_MATRIX, cols(), FloatMatrix(), givePointer(), isNotEmpty(), OOFEM_ERROR, and rows().

◆ add() [2/2]

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 991 of file floatmatrix.C.

References _LOOP_MATRIX, cols(), FloatMatrix(), givePointer(), isNotEmpty(), OOFEM_ERROR, rows(), and times().

◆ addProductOf()

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

Adds to the receiver product of \( a \cdot b \) .

Definition at line 464 of file floatmatrix.C.

References at(), cols(), FloatMatrix(), givePointer(), OOFEM_ERROR, and rows().

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

◆ addSubVectorCol()

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 622 of file floatmatrix.C.

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

Referenced by oofem::Shell7Base::giveTractionBC().

◆ addSubVectorRow()

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 600 of file floatmatrix.C.

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

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

◆ addTProductOf()

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

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

Definition at line 499 of file floatmatrix.C.

References at(), cols(), FloatMatrix(), givePointer(), OOFEM_ERROR, and rows().

◆ assemble() [1/2]

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 181 of file floatmatrix.C.

References at(), oofem::IntArray::at(), FloatMatrix(), isSquare(), OOFEM_ERROR, rows(), and oofem::IntArray::size().

◆ assemble() [2/2]

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 222 of file floatmatrix.C.

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

◆ assembleT()

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

Assembles the transposed 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 255 of file floatmatrix.C.

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

Referenced by oofem::MPElement::assembleTermContributionT().

◆ assignFloatMatrixF()

template<std::size_t N, std::size_t M>

void oofem::FloatMatrix::assignFloatMatrixF ( const FloatMatrixF< N, M > & mat )

inline

Definition at line 162 of file floatmatrix.h.

References cols(), N, resize(), and rows().

Referenced by FloatMatrix(), and operator=().

◆ at() [1/2]

double & oofem::FloatMatrix::at	(	std::size_t	i,
		std::size_t	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 231 of file floatmatrix.h.

◆ at() [2/2]

double oofem::FloatMatrix::at	(	std::size_t	i,
		std::size_t	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 221 of file floatmatrix.h.

◆ beDiagonal()

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 1111 of file floatmatrix.C.

References oofem::diag(), and resize().

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

◆ beDyadicProductOf()

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 535 of file floatmatrix.C.

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

Referenced by oofem::IsotropicDamageMaterial1::computeEta(), oofem::PrescribedGradientBCWeak::computeField(), oofem::tet21ghostsolid::computeStiffnessMatrix(), oofem::IntElLine1PF::computeStiffnessMatrix_ud(), oofem::FEI2dLineHermite::evaldNdx(), oofem::BTamNTerm::evaluate_lin(), oofem::NTcN::evaluate_lin(), oofem::NTN::evaluate_lin(), oofem::IsotropicDamageMaterial::give1dStressStiffMtrx(), oofem::SADGBoundaryElement::giveCharacteristicMatrix(), oofem::RankineMatGrad::giveGradientDamageStiffnessMatrix_du_BB(), oofem::MisesMatGrad::givePlaneStrainStiffMtrx(), oofem::IsotropicDamageMaterial::givePlaneStressStiffMtrx(), oofem::MisesMat::givePlaneStressStiffMtrx(), oofem::AnisotropicDamageMaterial::obtainAlpha3(), and oofem::IntElLine1PF::solveForLocalDamage().

◆ beInverseOf()

bool 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.

Returns: False if K is singular, otherwise true.

Definition at line 768 of file floatmatrix.C.

References _resize_internal(), at(), cols(), oofem::det(), FloatMatrix(), oofem::FloatArray::givePointer(), givePointer(), oofem::IntArray::givePointer(), isSquare(), nRows, OOFEM_ERROR, OOFEM_WARNING, rows(), and oofem::zero().

◆ beLocalCoordSys()

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 1646 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::Shell7BaseXFEM::computeCohesiveForces(), oofem::IntElSurfQuad1::computeTransformationMatrixAt(), oofem::IntElSurfTr1::computeTransformationMatrixAt(), oofem::Shell7BaseXFEM::giveFailedInterfaceNumber(), oofem::Shell7BaseXFEM::giveRecoveredTransverseInterfaceStress(), and oofem::Shell7BaseXFEM::recoverShearStress().

◆ beMatrixForm()

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

Definition at line 1689 of file floatmatrix.C.

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

Referenced by oofem::SolidShell::computeBEmatrixAt(), oofem::tet21ghostsolid::computeBoundarySurfaceLoadVector(), oofem::NLStructuralElement::computeCurrentVolume(), oofem::WeakPeriodicBoundaryCondition::computeDeformationGradient(), oofem::SolidShell::computeEVector(), oofem::tet21ghostsolid::computeLoadVector(), oofem::MaterialForceEvaluator::computeMaterialForce(), oofem::tet21ghostsolid::computeStiffnessMatrix(), oofem::LinearElasticMaterial::giveEshelbyStressVector_PlaneStrain(), oofem::LargeStrainMasterMaterialGrad::giveFirstPKStressVectorGrad(), oofem::Element::giveGlobalIPValue(), oofem::tet21ghostsolid::giveInternalForcesVectorGivenSolution(), oofem::tet21ghostsolid::giveInternalForcesVectorGivenSolutionDebug(), and oofem::AnisotropicDamageMaterial::giveRealStressVector().

◆ beMatrixFormOfStress()

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 1806 of file floatmatrix.C.

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

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

◆ beNMatrixOf()

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 1633 of file floatmatrix.C.

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

◆ bePinvID()

void oofem::FloatMatrix::bePinvID ( )

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

Definition at line 1600 of file floatmatrix.C.

References FloatMatrix(), and resize().

◆ beProductOf()

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 365 of file floatmatrix.C.

References _resize_internal(), at(), cols(), FloatMatrix(), givePointer(), OOFEM_ERROR, and rows().

◆ beProductTOf()

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 431 of file floatmatrix.C.

References _resize_internal(), at(), cols(), FloatMatrix(), givePointer(), OOFEM_ERROR, and rows().

Referenced by oofem::StructuralInterfaceElementPhF::computeBd_matrixAt(), oofem::Shell7Base::computeCauchyStressVector(), oofem::XfemStructuralElementInterface::computeCohesiveTangent(), oofem::SolidShell::computeGeometricStiffness(), oofem::Truss3dnl2::computeInitialStressStiffness(), oofem::LIBeam3dNL2::computeStiffnessMatrix(), oofem::LIBeam3dNL::computeStiffnessMatrix(), oofem::GradientDamageElement::computeStiffnessMatrix_ud(), oofem::NURBSSpace3dElement::drawScalar(), oofem::FEI2dTrQuad::evaldNdx(), oofem::PoissonTerm::evaluate_lin(), oofem::LargeStrainMasterMaterialGrad::give3dKappaMatrix(), oofem::RefinedElement::giveCompatibleBcDofArray(), oofem::MPlasticMaterial::giveConsistentStiffnessMatrix(), oofem::MPlasticMaterial::giveElastoPlasticStiffnessMatrix(), oofem::Shell7Base::giveLayerContributionToSR(), oofem::Shell7Base::nodalLeastSquareFitFromIP(), and rotatedWith().

◆ beSubMatrixOf() [1/2]

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 937 of file floatmatrix.C.

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

◆ beSubMatrixOf() [2/2]

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

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 val`ue to be added.
sc	Determines the column position (receiver's 1-based index) of first src value to be added. / void setTSubMatrix(const FloatMatrix &src, int sr, int sc); /* Assigns to the receiver the sub-matrix of another matrix.
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 901 of file floatmatrix.C.

References at(), cols(), FloatMatrix(), OOFEM_ERROR, resize(), and rows().

◆ beTProductOf()

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 398 of file floatmatrix.C.

References _resize_internal(), at(), cols(), FloatMatrix(), givePointer(), OOFEM_ERROR, and rows().

◆ beTranspositionOf()

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

Assigns to the receiver the transposition of parameter. Grows or shrinks if necessary.

Definition at line 351 of file floatmatrix.C.

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

◆ beUnitMatrix()

void oofem::FloatMatrix::beUnitMatrix ( )

Sets receiver to unity matrix.

Definition at line 1070 of file floatmatrix.C.

References at(), cols(), isSquare(), OOFEM_ERROR, rows(), and oofem::zero().

◆ changeComponentOrder()

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 1723 of file floatmatrix.C.

References at(), cols(), FloatMatrix(), and rows().

◆ checkBounds()

void oofem::FloatMatrix::checkBounds	(	Index	i,
		Index	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 161 of file floatmatrix.C.

References cols(), OOFEM_ERROR, and rows().

Referenced by operator()(), and operator()().

◆ clear()

*Sets size of receiver to be an empty matrix It will have zero rows and zero columns size void oofem::FloatMatrix::clear ( )

inline

Definition at line 540 of file floatmatrix.h.

References resize().

◆ cols()

int oofem::FloatMatrix::cols ( ) const

inline

◆ computeFrobeniusNorm()

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 1170 of file floatmatrix.C.

References _LOOP_MATRIX.

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

◆ computeNorm()

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 1179 of file floatmatrix.C.

References at(), cols(), oofem::FloatArray::givePointer(), givePointer(), oofem::max(), oofem::norm(), OOFEM_ERROR, and rows().

Referenced by computeReciprocalCondition().

◆ 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 1769 of file floatmatrix.C.

References computeNorm(), FloatMatrix(), giveDeterminant(), oofem::FloatArray::givePointer(), givePointer(), oofem::IntArray::givePointer(), oofem::inv(), isSquare(), OOFEM_ERROR, and oofem::rcond().

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

◆ copyColumn()

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 654 of file floatmatrix.C.

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

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

◆ copySubVectorRow()

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 1222 of file floatmatrix.C.

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

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

◆ data() [1/2]

double * oofem::FloatMatrix::data ( )

inline

Definition at line 158 of file floatmatrix.h.

References values.

◆ data() [2/2]

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

inline

Definition at line 157 of file floatmatrix.h.

References values.

Referenced by givePointer(), and givePointer().

◆ fromArray()

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

static

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 1046 of file floatmatrix.C.

References cols(), FloatMatrix(), resize(), rows(), and oofem::FloatArray::size().

◆ fromCols()

FloatMatrix oofem::FloatMatrix::fromCols ( std ::initializer_list< FloatArray > mat )

static

Definition at line 144 of file floatmatrix.C.

References FloatMatrix(), OOFEM_ERROR, and rows().

Referenced by oofem::Shell7BaseXFEM::giveLocalCZNodeCoordsForExport(), and oofem::Shell7BaseXFEM::giveLocalNodeCoordsForExport().

◆ fromIniList()

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

static

Definition at line 126 of file floatmatrix.C.

References _resize_internal(), FloatMatrix(), OOFEM_ERROR, and rows().

Referenced by oofem::Structural3DElement::computeConstitutiveMatrix_dPdF_At(), oofem::Structural3DElement::computeConstitutiveMatrixAt(), FloatMatrix(), and operator=().

◆ giveDeterminant()

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 1085 of file floatmatrix.C.

References cols(), FloatMatrix(), isSquare(), OOFEM_ERROR, and rows().

Referenced by oofem::Quad1MindlinShell3D::computeBmatrixAt(), oofem::tet21ghostsolid::computeBoundarySurfaceLoadVector(), oofem::NLStructuralElement::computeCurrentVolume(), oofem::SolidShell::computeEASBmatrixAt(), oofem::PolygonLine::computeIntersectionPoints(), oofem::MeshQualityErrorEstimator::computeJacobianError(), oofem::tet21ghostsolid::computeLoadVector(), oofem::TR21_2D_SUPG::computeQuadraticFunct(), oofem::TR21_2D_SUPG::computeQuadraticFunct(), computeReciprocalCondition(), oofem::tet21ghostsolid::computeStiffnessMatrix(), oofem::BSplineInterpolation::evaldNdx(), oofem::FEI3dHexaLin::evaldNdx(), oofem::FEI3dHexaQuad::evaldNdx(), oofem::FEI3dTetQuad::evaldNdx(), oofem::FEI3dWedgeLin::evaldNdx(), oofem::FEI3dWedgeQuad::evaldNdx(), oofem::NURBSInterpolation::evaldNdx(), oofem::TSplineInterpolation::evaldNdx(), oofem::tet21ghostsolid::giveInternalForcesVectorGivenSolution(), oofem::tet21ghostsolid::giveInternalForcesVectorGivenSolutionDebug(), oofem::FEI3dWedgeLin::giveTransformationJacobian(), oofem::FEI3dWedgeQuad::giveTransformationJacobian(), oofem::FEInterpolation::giveTransformationJacobian(), and oofem::Triangle::isOrientedAnticlockwise().

◆ giveNumberOfColumns()

int oofem::FloatMatrix::giveNumberOfColumns ( ) const

inline

Returns number of columns of receiver.

Definition at line 203 of file floatmatrix.h.

References cols().

◆ giveNumberOfRows()

int oofem::FloatMatrix::giveNumberOfRows ( ) const

inline

Returns number of rows of receiver.

Definition at line 200 of file floatmatrix.h.

References rows().

◆ givePackSize()

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

Definition at line 1938 of file floatmatrix.C.

References cols(), oofem::DataStream::givePackSizeOfDouble(), oofem::DataStream::givePackSizeOfSizet(), and rows().

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

◆ givePointer() [1/2]

double * oofem::FloatMatrix::givePointer ( )

inline

Definition at line 580 of file floatmatrix.h.

References data().

◆ givePointer() [2/2]

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 579 of file floatmatrix.h.

References data().

Referenced by add(), add(), addProductOf(), addTProductOf(), oofem::PetscSparseMtrx::assemble(), oofem::PetscSparseMtrx::assemble(), beInverseOf(), oofem::FloatArray::beProductOf(), beProductOf(), beProductTOf(), oofem::FloatArray::beTProductOf(), beTProductOf(), computeNorm(), computeReciprocalCondition(), oofem::AbaqusUserElement::giveInternalForcesVector(), oofem::FloatArray::plusProduct(), plusProductSymmUpper(), plusProductUnsym(), solveForRhs(), solveForRhs(), subtract(), oofem::PetscSparseMtrx::times(), oofem::PetscSparseMtrx::timesT(), and oofem::PetscSparseMtrx::toFloatMatrix().

◆ giveTrace()

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 1121 of file floatmatrix.C.

References cols(), isSquare(), OOFEM_ERROR, and rows().

Referenced by oofem::AnisotropicDamageMaterial::computeKappa().

◆ isAllFinite()

bool oofem::FloatMatrix::isAllFinite ( ) const

Returns true if no element is NAN or infinite.

Definition at line 339 of file floatmatrix.C.

References cols(), and rows().

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

◆ isNotEmpty()

bool oofem::FloatMatrix::isNotEmpty ( ) const

inline

Tests for empty matrix.

Definition at line 210 of file floatmatrix.h.

References cols(), and rows().

Referenced by add(), add(), oofem::EngngModel::assemble(), oofem::EngngModel::assemble(), oofem::NonStationaryTransportProblem::assembleAlgorithmicPartOfRhs(), oofem::DGProblem::assembleDirichletBcRhsVector(), oofem::DIIDynamic::assembleDirichletBcRhsVector(), oofem::EngngModel::assembleExtrapolatedForces(), oofem::EngngModel::assemblePrescribedExtrapolatedForces(), oofem::TransportElement::computeBCMtrxAt(), oofem::LIBeam2dNL::computeInitialStressMatrix(), oofem::NlDEIDynamic::computeMassMtrx(), oofem::TransportGradientDirichlet::computeTangent(), oofem::TransportGradientPeriodic::computeTangent(), oofem::SUPGElement::giveCharacteristicMatrix(), oofem::TR_SHELL01::giveCharacteristicMatrix(), oofem::TR_SHELL02::giveCharacteristicMatrix(), oofem::Element::giveRotationMatrix(), oofem::ScalarAdvectionLhsAssembler::matrixFromElement(), oofem::SUPGTangentAssembler::matrixFromElement(), oofem::TMLhsAssembler::matrixFromElement(), plus_Nt_a_otimes_b_B(), plusDyadSymmUpper(), plusDyadUnsym(), plusProductSymmUpper(), plusProductUnsym(), oofem::DEIDynamic::solveYourselfAt(), subtract(), oofem::DIIDynamic::timesMtrx(), oofem::NonLinearDynamic::timesMtrx(), oofem::ScalarAdvectionRhsAssembler::vectorFromElement(), and oofem::SUPGInternalForceAssembler::vectorFromElement().

◆ isSquare()

bool oofem::FloatMatrix::isSquare ( ) const

inline

Returns nonzero if receiver is square matrix.

Definition at line 208 of file floatmatrix.h.

References cols(), and rows().

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

◆ jaco_()

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 1953 of file floatmatrix.C.

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

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().

◆ negated()

void oofem::FloatMatrix::negated ( )

Changes sign of receiver values.

Definition at line 1164 of file floatmatrix.C.

References _LOOP_MATRIX.

Referenced by oofem::MixedGradientPressureNeumann::assemble(), oofem::MixedGradientPressureWeakPeriodic::assemble(), oofem::PrescribedGradientBCNeumann::assemble(), oofem::TransverseReinfConstraint::assemble(), oofem::PrescribedDispSlipBCNeumannRC::assembleOnReinfStress(), oofem::PrescribedDispSlipBCNeumannRC::assembleOnStress(), oofem::PrescribedDispSlipBCNeumannRC::assembleOnTransferStress(), oofem::TR1_2D_SUPG::computeOutFlowBCTerm_MB(), oofem::TR1_2D_SUPG_AXI::computeOutFlowBCTerm_MB(), oofem::PrescribedGradientBCNeumann::computeTangent(), oofem::TransportGradientNeumann::computeTangent(), oofem::TransportGradientDirichlet::computeXi(), oofem::MPlasticMaterial2::giveConsistentStiffnessMatrix(), oofem::MPlasticMaterial2::giveElastoPlasticStiffnessMatrix(), oofem::MPlasticMaterial::giveElastoPlasticStiffnessMatrix(), oofem::LinearElasticMaterial::giveEshelbyStressVector_PlaneStrain(), oofem::AbaqusUserElement::giveInternalForcesVector(), oofem::TransverseReinfConstraint::integrateTangentOnConcrete(), and subtract().

◆ operator()() [1/2]

double & oofem::FloatMatrix::operator()	(	std::size_t	i,
		std::size_t	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 138 of file floatmatrix.h.

References checkBounds(), nRows, and values.

◆ operator()() [2/2]

double oofem::FloatMatrix::operator()	(	std::size_t	i,
		std::size_t	j ) const

inline

Coefficient access function. Implements 0-based indexing.

Parameters

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

Definition at line 150 of file floatmatrix.h.

References checkBounds(), nRows, and values.

◆ operator=() [1/4]

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

inline

Assignment operator, adjusts size of the receiver if necessary.

Definition at line 117 of file floatmatrix.h.

References FloatMatrix(), nColumns, nRows, and values.

◆ operator=() [2/4]

template<std::size_t N, std::size_t M>

FloatMatrix & oofem::FloatMatrix::operator= ( const FloatMatrixF< N, M > & mat )

inline

Definition at line 186 of file floatmatrix.h.

References assignFloatMatrixF(), and FloatMatrix().

◆ operator=() [3/4]

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

inlinenoexcept

Definition at line 123 of file floatmatrix.h.

References FloatMatrix(), nColumns, nRows, and values.

◆ operator=() [4/4]

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

inline

Definition at line 169 of file floatmatrix.h.

References FloatMatrix(), and fromIniList().

◆ plus_Nt_a_otimes_b_B()

void oofem::FloatMatrix::plus_Nt_a_otimes_b_B	(	const FloatMatrix &	N,
		const FloatArray &	a,
		const FloatArray &	b,
		const FloatMatrix &	B,
		double	dV = 1 )

Adds to the receiver the product \(Nt a \otimes b B \mathrm{d}V\). If the receiver has zero size, it is expanded.

Parameters

N	matrix N in equation.
a	Array a in equation.
b	Array b in equation.
B	matrix N in equation.
dV	Scaling factor.

Definition at line 745 of file floatmatrix.C.

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

Referenced by oofem::StructuralPenaltyContactBoundaryCondition::computeTangentFromContact().

◆ plusDyadSymmUpper()

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 667 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), cols(), oofem::FloatArray::givePointer(), oofem::FloatArray::giveSize(), isNotEmpty(), resize(), and rows().

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::MITC4Shell::computeStiffnessMatrix(), oofem::Quad1MindlinShell3D::computeStiffnessMatrix(), oofem::SurfaceTensionBoundaryCondition::computeTangentFromElement(), oofem::TransportGradientDirichlet::computeXi(), oofem::Shell7Base::NodalRecoveryMI_computeNNMatrix(), oofem::Tr_Warp::ZZNodalRecoveryMI_computeNNMatrix(), and oofem::ZZNodalRecoveryModelInterface::ZZNodalRecoveryMI_computeNNMatrix().

◆ plusDyadUnsym()

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 723 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), cols(), oofem::FloatArray::givePointer(), oofem::FloatArray::giveSize(), isNotEmpty(), resize(), and rows().

Referenced by oofem::tet21ghostsolid::computeLoadVector(), oofem::TR1_2D_SUPG_AXI::computeOutFlowBCTerm_MB(), oofem::tet21ghostsolid::computeStiffnessMatrix(), oofem::MixedGradientPressureWeakPeriodic::computeStress(), oofem::SurfaceTensionBoundaryCondition::computeTangentFromElement(), oofem::TransportGradientDirichlet::computeXi(), oofem::NTaTmTe::evaluate_lin(), oofem::TransportGradientNeumann::integrateTangent(), oofem::Shell7Base::NodalRecoveryMI_computeNValProduct(), oofem::IDNLMaterial::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::MisesMatNl::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::RankineMatNl::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::TrabBoneNL3D::NonlocalMaterialStiffnessInterface_addIPContribution(), oofem::Tr_Warp::ZZNodalRecoveryMI_computeNValProduct(), and oofem::ZZNodalRecoveryModelInterface::ZZNodalRecoveryMI_computeNValProduct().

◆ plusProductSymmUpper()

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 1243 of file floatmatrix.C.

References at(), cols(), FloatMatrix(), givePointer(), isNotEmpty(), resize(), and rows().

◆ plusProductUnsym()

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 691 of file floatmatrix.C.

References at(), cols(), FloatMatrix(), givePointer(), isNotEmpty(), resize(), and rows().

◆ printYourself() [1/2]

void oofem::FloatMatrix::printYourself ( ) const

Definition at line 1501 of file floatmatrix.C.

References at(), cols(), and rows().

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

◆ printYourself() [2/2]

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 1542 of file floatmatrix.C.

References at(), cols(), and rows().

◆ printYourselfToFile()

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

Print matrix to file.

Parameters

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

Definition at line 1520 of file floatmatrix.C.

References at(), cols(), OOFEM_ERROR, and rows().

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

◆ pY()

void oofem::FloatMatrix::pY ( ) const

Definition at line 1567 of file floatmatrix.C.

References at(), cols(), and rows().

◆ resize()

void oofem::FloatMatrix::resize	(	Index	rows,
		Index	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 79 of file floatmatrix.C.

References nColumns, nRows, rows(), and values.

◆ resizeWithData()

void oofem::FloatMatrix::resizeWithData	(	Index	rows,
		Index	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 91 of file floatmatrix.C.

References at(), cols(), FloatMatrix(), oofem::min(), nColumns, nRows, resize(), rows(), and values.

◆ restoreYourself()

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

Definition at line 1909 of file floatmatrix.C.

References oofem::CIO_IOERR, oofem::CIO_OK, oofem::DataStream::read(), and resize().

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

◆ rotatedWith()

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 1612 of file floatmatrix.C.

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

◆ rows()

int oofem::FloatMatrix::rows ( ) const

inline

◆ setColumn()

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 641 of file floatmatrix.C.

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

Referenced by oofem::XfemStructuralElementInterface::computeCohesiveTangent(), oofem::tet21ghostsolid::computeNumericStiffnessMatrix(), oofem::tet21ghostsolid::computeNumericStiffnessMatrixDebug(), oofem::MixedGradientPressureWeakPeriodic::computeTangents(), oofem::FE2FluidMaterial::computeTangents3D(), oofem::DKTPlate::computeVertexBendingMoments(), oofem::FEI3dHexaQuad::edgeEvaldNdx(), oofem::FEI3dHexaLin::evaldNdx(), oofem::FEI3dHexaQuad::evaldNdx(), oofem::FEI3dTetQuad::evaldNdx(), oofem::FEI3dWedgeLin::evaldNdx(), oofem::FEI3dWedgeQuad::evaldNdx(), oofem::AbaqusUserMaterial::give3dMaterialStiffnessMatrix(), oofem::MFrontUserMaterial::give3dMaterialStiffnessMatrix(), oofem::StructuralFE2Material::give3dMaterialStiffnessMatrix(), oofem::StructuralPythonMaterial::give3dMaterialStiffnessMatrix(), oofem::FEI3dHexaLin::giveJacobianMatrixAt(), oofem::FEI3dHexaQuad::giveJacobianMatrixAt(), oofem::FEI3dTetQuad::giveJacobianMatrixAt(), oofem::FEI3dWedgeLin::giveJacobianMatrixAt(), oofem::FEI3dWedgeQuad::giveJacobianMatrixAt(), oofem::StructuralSlipFE2Material::giveSensitivities(), oofem::MixedGradientPressureWeakPeriodic::integrateTractionXTangent(), oofem::PLHoopStressCirc::propagateInterface(), oofem::InverseIteration::solve(), oofem::SparseLinearSystemNM::solve(), and oofem::SubspaceIteration::solve().

◆ setSubMatrix()

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 549 of file floatmatrix.C.

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

Referenced by oofem::ContactPair::compute_dNdxi_matrices(), oofem::FEContactPoint::compute_dNdxi_matrix(), oofem::XfemElementInterface::ComputeBOrBHMatrix(), oofem::TransportElement::computeFlow(), oofem::MITC4Shell::computeGtoLRotationMatrix(), oofem::Quad1MindlinShell3D::computeGtoLRotationMatrix(), oofem::TR_SHELL11::computeGtoLRotationMatrix(), oofem::ContactPair::computeNmatrix(), oofem::Shell7Base::giveLayerContributionToSR(), and oofem::Homogenize::herveZaoui().

◆ solveForRhs() [1/2]

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 1303 of file floatmatrix.C.

References oofem::FloatArray::at(), at(), beTranspositionOf(), cols(), FloatMatrix(), oofem::FloatArray::givePointer(), givePointer(), oofem::IntArray::givePointer(), isSquare(), nRows, OOFEM_ERROR, rows(), oofem::FloatArray::size(), and oofem::transpose().

Referenced by oofem::MMALeastSquareProjection::__mapVariable(), oofem::MPlasticMaterial2::closestPointReturn(), oofem::MPlasticMaterial::closestPointReturn(), oofem::KelvinChainMaterial::computeCharCoefficients(), oofem::MaxwellChainMaterial::computeCharCoefficients(), oofem::TransportGradientNeumann::computeEta(), oofem::SPRNodalRecoveryModel::computePatch(), oofem::MPlasticMaterial2::cuttingPlaneReturn(), oofem::MPlasticMaterial::cuttingPlaneReturn(), oofem::LEPlic::doCellDLS(), oofem::StructuralMaterial::giveFirstPKStressVector_StressControl(), oofem::RCM2Material::giveRealPrincipalStressVector3d(), oofem::FCMMaterial::giveRealStressVector(), oofem::StructuralMaterial::giveRealStressVector_ShellStressControl(), oofem::FEI3dHexaConst::global2local(), oofem::FEI3dHexaLin::global2local(), oofem::FEI3dHexaQuad::global2local(), oofem::FEI3dTetQuad::global2local(), oofem::FEI3dWedgeLin::global2local(), oofem::FEI3dWedgeQuad::global2local(), oofem::FEInterpolation2d::global2local(), oofem::LSpace::NodalAveragingRecoveryMI_computeNodalValue(), oofem::FluidMaterialEvaluator::solveYourself(), and oofem::StructuralMaterialEvaluator::solveYourself().

◆ solveForRhs() [2/2]

bool 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 1397 of file floatmatrix.C.

References at(), beTranspositionOf(), cols(), FloatMatrix(), giveNumberOfColumns(), givePointer(), oofem::IntArray::givePointer(), isSquare(), nRows, OOFEM_ERROR, rows(), and oofem::transpose().

◆ storeYourself()

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

Definition at line 1884 of file floatmatrix.C.

References oofem::CIO_IOERR, oofem::CIO_OK, cols(), and rows().

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

◆ subtract()

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

Subtracts matrix from the receiver.

Parameters

a	Matrix to be subtracted.

Definition at line 1020 of file floatmatrix.C.

References _LOOP_MATRIX, cols(), FloatMatrix(), givePointer(), isNotEmpty(), negated(), OOFEM_ERROR, and rows().

Referenced by oofem::Shell7BaseXFEM::computeCohesiveNmatrixAt(), oofem::AnisotropicDamageMaterial::computeDamageTensor(), oofem::Truss3dnl2::computeInitialStressStiffness(), oofem::Shell7Base::computePressureTangentMatrix(), oofem::AnisotropicDamageMaterial::computeSecantOperator(), oofem::PrescribedDispSlipBCDirichletRC::computeTangent(), oofem::PrescribedGenStrainShell7::computeTangent(), oofem::PrescribedGradient::computeTangent(), oofem::PrescribedGradientBCPeriodic::computeTangent(), oofem::TransportGradientDirichlet::computeTangent(), oofem::TransportGradientPeriodic::computeTangent(), oofem::MixedGradientPressureDirichlet::computeTangents(), oofem::FE2FluidMaterial::computeTangents3D(), oofem::StructuralInterfaceMaterial::give2dStiffnessMatrix_Eng(), oofem::StructuralFE2Material::give3dMaterialStiffnessMatrix(), oofem::MPlasticMaterial::giveConsistentStiffnessMatrix(), oofem::FCMMaterial::giveMaterialStiffnessMatrix(), oofem::AnisotropicDamageMaterial::giveRealStressVector(), oofem::operator-(), and oofem::operator-=().

◆ symmetrized()

void oofem::FloatMatrix::symmetrized ( )

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

Definition at line 1137 of file floatmatrix.C.

References at(), cols(), OOFEM_ERROR, and rows().

◆ times()

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

Multiplies receiver by factor f.

Parameters

f	Factor to multiply by.

Definition at line 1155 of file floatmatrix.C.

References _LOOP_MATRIX.

◆ writeCSV()

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

Writes receiver as CSV (comma seperated values)

Parameters

name Filename

Definition at line 1586 of file floatmatrix.C.

References at(), cols(), and rows().

◆ zero()

void oofem::FloatMatrix::zero ( )

Zeroes all coefficient of receiver.

Definition at line 1064 of file floatmatrix.C.

References _LOOP_MATRIX.

◆ operator<<

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

friend

Definition at line 2126 of file floatmatrix.C.

References cols(), FloatMatrix(), and rows().

Member Data Documentation

◆ Dim

int oofem::FloatMatrix::Dim = 2

staticconstexpr

Definition at line 102 of file floatmatrix.h.

◆ nColumns

std::size_t oofem::FloatMatrix::nColumns

protected

Number of columns.

Definition at line 97 of file floatmatrix.h.

Referenced by _resize_internal(), cols(), FloatMatrix(), FloatMatrix(), FloatMatrix(), FloatMatrix(), operator=(), operator=(), plus_Nt_a_otimes_b_B(), resize(), and resizeWithData().

◆ nRows

std::size_t oofem::FloatMatrix::nRows

protected

Number of rows.

Definition at line 95 of file floatmatrix.h.

Referenced by _resize_internal(), beInverseOf(), FloatMatrix(), FloatMatrix(), FloatMatrix(), FloatMatrix(), operator()(), operator()(), operator=(), operator=(), plus_Nt_a_otimes_b_B(), resize(), resizeWithData(), rows(), solveForRhs(), and solveForRhs().

◆ values

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

protected

Values of matrix stored column wise.

Definition at line 99 of file floatmatrix.h.

Referenced by _resize_internal(), data(), data(), FloatMatrix(), FloatMatrix(), FloatMatrix(), FloatMatrix(), operator()(), operator()(), operator=(), operator=(), plus_Nt_a_otimes_b_B(), resize(), and resizeWithData().

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

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ Scalar

Constructor & Destructor Documentation

◆ FloatMatrix() [1/6]

◆ FloatMatrix() [2/6]

◆ FloatMatrix() [3/6]

◆ FloatMatrix() [4/6]

◆ FloatMatrix() [5/6]

◆ FloatMatrix() [6/6]

Member Function Documentation

◆ _resize_internal()

◆ add() [1/2]

◆ add() [2/2]

◆ addProductOf()

◆ addSubVectorCol()

◆ addSubVectorRow()

◆ addTProductOf()

◆ assemble() [1/2]

◆ assemble() [2/2]

◆ assembleT()

◆ assignFloatMatrixF()

◆ at() [1/2]

◆ at() [2/2]

◆ beDiagonal()

◆ beDyadicProductOf()

◆ beInverseOf()

◆ beLocalCoordSys()

◆ beMatrixForm()

◆ beMatrixFormOfStress()

◆ beNMatrixOf()

◆ bePinvID()

◆ beProductOf()

◆ beProductTOf()

◆ beSubMatrixOf() [1/2]

◆ beSubMatrixOf() [2/2]

◆ beTProductOf()

◆ beTranspositionOf()

◆ beUnitMatrix()

◆ changeComponentOrder()

◆ checkBounds()

◆ clear()

◆ cols()

◆ computeFrobeniusNorm()

◆ computeNorm()

◆ computeReciprocalCondition()

◆ copyColumn()

◆ copySubVectorRow()

◆ data() [1/2]

◆ data() [2/2]

◆ fromArray()

◆ fromCols()

◆ fromIniList()

◆ giveDeterminant()

◆ giveNumberOfColumns()

◆ giveNumberOfRows()

◆ givePackSize()

◆ givePointer() [1/2]

◆ givePointer() [2/2]

◆ giveTrace()

◆ isAllFinite()

◆ isNotEmpty()

◆ isSquare()

◆ jaco_()

◆ negated()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ operator=() [1/4]

◆ operator=() [2/4]

◆ operator=() [3/4]

◆ operator=() [4/4]

◆ plus_Nt_a_otimes_b_B()

◆ plusDyadSymmUpper()

◆ plusDyadUnsym()