File kernel/src/utils/SiconosAlgebra/BlockMatrix.hpp#

Go to the source code of this file

Object to handle block-matrices.

class BlockMatrix : public SiconosMatrix
#include <BlockMatrix.hpp>

“Block” matrix, ie container of matrices

A BlockMatrix is a boost::ublas::compressed_matrix of SP::SiconosMatrix.

The blocks positions are given by two Index objects, tabRow and tabCol.

If block 1 is n1xm1, block2 n2xm2, block3 n3xm3 …, then:

tabRow = [ n1 n1+n2 n1+n2+n3 …]

tabCol = [ m1 m1+m2 m1+m2+m3 …]

Public Functions

BlockMatrix(const SiconosMatrix &m)

copy constructor

Parameters:

m – a SiconosMatrix

BlockMatrix(const BlockMatrix &m)

copy constructor

Parameters:

m – a BlockMatrix

BlockMatrix(const std::vector<SP::SiconosMatrix> &m, unsigned int row, unsigned int col)

constructor with a list of pointer to SiconosMatrix (!links with pointer, no copy!)

Parameters:
  • m – a vector of SiconosMatrix

  • row – number of blocks in a row

  • col – number of col in a row

BlockMatrix(SP::SiconosMatrix A, SP::SiconosMatrix B, SP::SiconosMatrix C, SP::SiconosMatrix D)

contructor with a list of 4 pointer to SiconosMatrix (!links with pointer, no copy!)

Parameters:
  • A – block (0,0)

  • B – block (0,1)

  • C – block (1,0)

  • D – block (1,1)

~BlockMatrix(void) noexcept

destructor

inline virtual bool checkSymmetry(double tol) const override

determines if the matrix is symmetric up to a given tolerance

Returns:

true if the matrix is inversed

virtual unsigned int numberOfBlocks(unsigned int i) const override

get the number of block (i=0, row, i=1 col)

Parameters:

i – unsigned int(i=0, row, i=1 col)

Returns:

an unsigned int

virtual const DenseMat getDense(unsigned int row = 0, unsigned int col = 0) const override

get DenseMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a DenseMat

virtual const TriangMat getTriang(unsigned int row = 0, unsigned int col = 0) const override

get TriangMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a TriangMat

virtual const SymMat getSym(unsigned int row = 0, unsigned int col = 0) const override

get SymMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a SymMat

virtual const BandedMat getBanded(unsigned int row = 0, unsigned int col = 0) const override

get BandedMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a BandedMat

virtual const SparseMat getSparse(unsigned int row = 0, unsigned int col = 0) const override

get SparseMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a SparseMat

virtual const SparseCoordinateMat getSparseCoordinate(unsigned int row = 0, unsigned int col = 0) const override

get SparseCoordinateMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a SparseCoordinateMat

virtual const ZeroMat getZero(unsigned int row = 0, unsigned int col = 0) const override

get ZeroMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a ZeroMat

virtual const IdentityMat getIdentity(unsigned int row = 0, unsigned int col = 0) const override

get getIdentity matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

an IdentityMat

virtual DenseMat *dense(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on DenseMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a DenseMat*

virtual TriangMat *triang(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on TriangMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a TriangMat*

virtual SymMat *sym(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on SymMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – `an unsigned int, position of the block (column)

Returns:

a SymMat*

virtual BandedMat *banded(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on BandedMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a BandedMat*

virtual SparseMat *sparse(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on SparseMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a SparseMat*

virtual SparseCoordinateMat *sparseCoordinate(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on SparseCoordinateMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a SparseCoordinateMat*

virtual ZeroMat *zero_mat(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on ZeroMat matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

a ZeroMat*

virtual IdentityMat *identity(unsigned int row = 0, unsigned int col = 0) const override

get a pointer on Identity matrix

Parameters:
  • row – an unsigned int, position of the block (row)

  • col – an unsigned int, position of the block (column)

Returns:

an IdentityMat*

virtual double *getArray(unsigned int row = 0, unsigned int col = 0) const override

return the address of the array of double values of the matrix

Parameters:
  • row – position for the required block ->useless for SimpleMatrix

  • col – position for the required block ->useless for SimpleMatrix

Returns:

double* : the pointer on the double array

virtual void zero() override

sets all the values of the matrix to 0.0

virtual void randomize() override

Initialize the matrix with random values.

virtual void randomize_sym() override

Initialize a symmetric matrix with random values.

virtual void eye() override

set an identity matrix

virtual unsigned int size(unsigned int index) const override

get the number of rows or columns of the matrix

Parameters:

index – 0 for rows, 1 for columns

Returns:

an int

virtual void resize(unsigned int nbrow, unsigned int nbcol, unsigned int lower = 0, unsigned int upper = 0, bool b = true) override

resize the matrix with nbrow rows and nbcol columns, lower and upper are useful only for SparseMat.The existing elements of the Block matrix are preseved when specified.

Parameters:
  • nbrow

  • nbcol

  • lower

  • upper

  • b

virtual double normInf() const override

compute the infinite norm of the Block matrix

Returns:

a double

virtual void display() const override

display data on standard output

virtual void displayExpert(bool brief = true) const override

display data on standard output

virtual std::string toString() const override

put data of the matrix into a std::string

Returns:

std::string

virtual double &operator()(unsigned int i, unsigned int j) override

get or set the element matrix[i,j]

Parameters:
  • i – an unsigned int

  • j – an unsigned int

Returns:

the element matrix[i,j]

virtual double operator()(unsigned int i, unsigned int j) const override

get or set the element matrix[i,j]

Parameters:
  • i – an unsigned int

  • j – an unsigned int

Returns:

the element matrix[i,j]

virtual double getValue(unsigned int i, unsigned int j) const override

return the element matrix[i,j]

Parameters:
  • i – an unsigned int

  • j – an unsigned int

Returns:

a double

virtual void setValue(unsigned int i, unsigned int j, double value) override

set the element matrix[i,j]

Parameters:
  • i – an unsigned int i

  • j – an unsigned int j

  • value

virtual void trans() override

transpose in place: x->trans() is x = transpose of x.

virtual void trans(const SiconosMatrix &m) override

transpose a matrix: x->trans(m) is x = transpose of m.

Parameters:

m – the matrix to be transposed.

inline Index getTabRow() const

get the vector tabRow

Returns:

a vector of int

inline Index getTabCol() const

get the vector tabCol

Returns:

a vector of int

inline virtual const SP::Index tabRow() const override

get the vector tabRow

Returns:

a pointer to vector of int

inline virtual const SP::Index tabCol() const override

get the vector tabCol

Returns:

a pointer to vector of int

virtual SP::SiconosMatrix block(unsigned int row = 0, unsigned int col = 0) override

get block at position row-col

Parameters:
  • row – unsigned int

  • col – unsigned int

Returns:

SP::SiconosMatrix the requested block

virtual SPC::SiconosMatrix block(unsigned int row = 0, unsigned int col = 0) const override

get block at position row-col

Parameters:
  • row – unsigned int

  • col – unsigned int

Returns:

SP::SiconosMatrix the requested block

virtual void getRow(unsigned int r, SiconosVector &v) const override

get row index of current matrix and save it in v

Parameters:
  • r – index of required line

  • v[out] a vector

virtual void setRow(unsigned int r, const SiconosVector &v) override

set line row of the current matrix with vector v

Parameters:
  • r – index of required line

  • v – a vector

virtual void getCol(unsigned int c, SiconosVector &v) const override

get column index of current matrix and save it into vOut

Parameters:
  • c – index of required column

  • v[out] a vector

virtual void setCol(unsigned int c, const SiconosVector &v) override

set column col of the current matrix with vector

Parameters:
  • c – index of required column

  • v – a vector

void addSimple(unsigned int &i, unsigned int &j, const SiconosMatrix &m)

add a part of the input matrix (starting from (i,j) pos) to the current matrix

Parameters:
  • i – an unsigned int i (in-out)

  • j – an unsigned int j (in-out)

  • m – a SiconosMatrix (in-out)

void subSimple(unsigned int &i, unsigned int &j, const SiconosMatrix &m)

subtract a part of the input matrix (starting from (i,j) pos) to the current matrix

Parameters:
  • i – an unsigned int i (in-out)

  • j – an unsigned int j (in-out)

  • m – a SiconosMatrix (in-out)

virtual BlockMatrix &operator=(const SiconosMatrix &m) override

assignment

Parameters:

m – the matrix to be copied

Returns:

BlockMatrix&

BlockMatrix &operator=(const BlockMatrix &m)

assignment

Parameters:

m – the matrix to be copied

Returns:

BlockMatrix&

virtual BlockMatrix &operator=(const DenseMat &m) override

assignment

Parameters:

m – the matrix to be copied

Returns:

BlockMatrix&

virtual BlockMatrix &operator+=(const SiconosMatrix &m) override

operator +=

Parameters:

m – the matrix to add

Returns:

BlockMatrix&

virtual BlockMatrix &operator-=(const SiconosMatrix &m) override

operator -=

Parameters:

m – the matrix to subtract

Returns:

BlockMatrix&

inline virtual void updateNumericsMatrix() override#
virtual void PLUFactorizationInPlace() override

computes an LU factorization of a general M-by-N matrix using partial pivoting with row interchanges.

The result is returned in this (InPlace). Based on Blas dgetrf function.

virtual void Factorize() override

computes a factorization of a general M-by-N matrix The implementation is based on an internal NumericsMatrix

virtual void PLUInverseInPlace() override

compute inverse of this thanks to LU factorization with Partial pivoting.

This method inverts U and then computes inv(A) by solving the system inv(A)*L = inv(U) for inv(A). The result is returned in this (InPlace). Based on Blas dgetri function.

virtual void PLUForwardBackwardInPlace(SiconosMatrix &B) override

solves a system of linear equations A * X = B (A=this) with a general N-by-N matrix A using the LU factorization computed by PLUFactorizationInPlace.

Based on Blas dgetrs function.

Parameters:

B[inout] on input the RHS matrix b; on output: the result x

virtual void Solve(SiconosMatrix &B) override

solves a system of linear equations A * X = B (A=this) for a general N-by-N matrix A using the LU factorization computed by PLUFactorize.

Parameters:

B[inout] on input the RHS matrix b; on output the result x

virtual void PLUForwardBackwardInPlace(SiconosVector &B) override

solves a system of linear equations A * X = B (A=this) with a general N-by-N matrix A using the LU factorization computed by PLUFactorizationInPlace.

Based on Blas dgetrs function.

Parameters:

B[inout] on input the RHS matrix b; on output: the result x

virtual void Solve(SiconosVector &B) override

solves a system of linear equations A * X = B (A=this) for a general N-by-N matrix A using the LU factorization computed by PLUFactorize.

Parameters:

B[inout] on input the RHS matrix b; on output the result x

ACCEPT_STD_VISITORS()#
virtual size_t nnz(double tol = 1.e-14) override

number of non-zero in the matrix

Parameters:

tol – the tolerance under which a number is considered zero

Private Functions

ACCEPT_SERIALIZATION(BlockMatrix)#
inline BlockMatrix()#

default constructor

Private Members

SP::BlocksMat _mat#

A container of pointers to SiconosMatrix.

SP::Index _tabRow#

list of blocks dimension - tabRow[i] = tabRow[i-1] + ni, ni being the number of rows of block i.

SP::Index _tabCol#

list of blocks dimension - tabCol[i] = tabCol[i-1] + ni, ni being the number of columns of block i.

unsigned int _dimRow#

Number of rows (Warning: total number of scalar elements, not number of blocks)

unsigned int _dimCol#

Number of columns (Warning: total number of scalar elements, not number of blocks)

Friends

friend class SimpleMatrix
friend std::ostream &operator<<(std::ostream &os, const BlockMatrix &bm)

send data of the matrix to an ostream

Parameters:
Returns:

The same output stream

friend void scal(double, const SiconosMatrix&, SiconosMatrix&, bool)

multiplication of a matrix by a scalar, B = a*A (init = true) or B += a*A (init = false)

Parameters:
friend SiconosMatrix &operator*=(SiconosMatrix &m, const double &s)#
friend SiconosMatrix &operator/=(SiconosMatrix &m, const double &s)#