qpOASES
3.2.1
An Implementation of the Online Active Set Strategy
|
Interfaces matrix-vector operations tailored to general dense matrices. More...
#include <Matrices.hpp>
Public Member Functions | |
DenseMatrix () | |
DenseMatrix (int_t m, int_t n, int_t lD, real_t *v) | |
virtual | ~DenseMatrix () |
virtual void | free () |
virtual Matrix * | duplicate () const |
virtual real_t | diag (int_t i) const |
virtual BooleanType | isDiag () const |
virtual real_t | getNorm (int_t type=2) const |
virtual real_t | getRowNorm (int_t rNum, int_t type=2) const |
virtual returnValue | getRowNorm (real_t *norm, int_t type=2) const |
virtual returnValue | getRow (int_t rNum, const Indexlist *const icols, real_t alpha, real_t *row) const |
virtual returnValue | getCol (int_t cNum, const Indexlist *const irows, real_t alpha, real_t *col) const |
virtual returnValue | getSparseSubmatrix (int_t irowsLength, const int_t *const irowsNumber, int_t icolsLength, const int_t *const icolsNumber, int_t rowoffset, int_t coloffset, int_t &numNonzeros, int_t *irn, int_t *jcn, real_t *avals, BooleanType only_lower_triangular=BT_FALSE) const |
virtual returnValue | times (int_t xN, real_t alpha, const real_t *x, int_t xLD, real_t beta, real_t *y, int_t yLD) const |
virtual returnValue | transTimes (int_t xN, real_t alpha, const real_t *x, int_t xLD, real_t beta, real_t *y, int_t yLD) const |
virtual returnValue | times (const Indexlist *const irows, const Indexlist *const icols, int_t xN, real_t alpha, const real_t *x, int_t xLD, real_t beta, real_t *y, int_t yLD, BooleanType yCompr=BT_TRUE) const |
virtual returnValue | transTimes (const Indexlist *const irows, const Indexlist *const icols, int_t xN, real_t alpha, const real_t *x, int_t xLD, real_t beta, real_t *y, int_t yLD) const |
virtual returnValue | addToDiag (real_t alpha) |
virtual real_t * | full () const |
virtual returnValue | print (const char *name=0) const |
virtual returnValue | writeToFile (FILE *output_file, const char *prefix) const |
virtual returnValue | getSparseSubmatrix (const Indexlist *const irows, const Indexlist *const icols, int_t rowoffset, int_t coloffset, int_t &numNonzeros, int_t *irn, int_t *jcn, real_t *avals, BooleanType only_lower_triangular=BT_FALSE) const |
virtual returnValue | getSparseSubmatrix (const Indexlist *const irows, int_t idx_icol, int_t rowoffset, int_t coloffset, int_t &numNonzeros, int_t *irn, int_t *jcn, real_t *avals, BooleanType only_lower_triangular=BT_FALSE) const |
virtual returnValue | getSparseSubmatrix (int_t idx_row, const Indexlist *const icols, int_t rowoffset, int_t coloffset, int_t &numNonzeros, int_t *irn, int_t *jcn, real_t *avals, BooleanType only_lower_triangular=BT_FALSE) const |
BooleanType | needToFreeMemory () const |
void | doFreeMemory () |
void | doNotFreeMemory () |
Protected Attributes | |
int_t | nRows |
int_t | nCols |
int_t | leaDim |
real_t * | val |
BooleanType | freeMemory |
Dense matrix class (row major format).
DenseMatrix::DenseMatrix | ( | ) | [inline] |
Default constructor.
Referenced by duplicate().
DenseMatrix::DenseMatrix | ( | int_t | m, |
int_t | n, | ||
int_t | lD, | ||
real_t * | v | ||
) | [inline] |
Constructor from vector of values. Caution: Data pointer must be valid throughout lifetime
m | Number of rows. |
n | Number of columns. |
lD | Leading dimension. |
v | Values. |
DenseMatrix::~DenseMatrix | ( | ) | [virtual] |
Destructor.
References BT_TRUE, free(), and Matrix::needToFreeMemory().
returnValue DenseMatrix::addToDiag | ( | real_t | alpha | ) | [virtual] |
real_t DenseMatrix::diag | ( | int_t | i | ) | const [virtual] |
void Matrix::doFreeMemory | ( | ) | [inline, inherited] |
Enables de-allocation of internal memory.
References BT_TRUE, and Matrix::freeMemory.
Referenced by QProblemB::createDiagSparseMat(), duplicate(), SparseMatrix::duplicate(), SparseMatrixRow::duplicate(), SymDenseMat::duplicateSym(), SymSparseMat::duplicateSym(), QProblemB::setupQPdataFromFile(), QProblem::setupQPdataFromFile(), solveOqpBenchmark(), SparseMatrix::SparseMatrix(), and SparseMatrixRow::SparseMatrixRow().
void Matrix::doNotFreeMemory | ( | ) | [inline, inherited] |
Disables de-allocation of internal memory.
Referenced by SparseMatrix::free(), SparseMatrixRow::free(), Matrix::Matrix(), SparseMatrix::SparseMatrix(), and SparseMatrixRow::SparseMatrixRow().
Matrix * DenseMatrix::duplicate | ( | ) | const [virtual] |
Returns a deep-copy of the Matrix object.
Implements Matrix.
Reimplemented in SymDenseMat.
References BT_TRUE, DenseMatrix(), Matrix::doFreeMemory(), nCols, Matrix::needToFreeMemory(), nRows, real_t, and val.
void DenseMatrix::free | ( | ) | [virtual] |
real_t * DenseMatrix::full | ( | ) | const [virtual] |
returnValue DenseMatrix::getCol | ( | int_t | cNum, |
const Indexlist *const | irows, | ||
real_t | alpha, | ||
real_t * | col | ||
) | const [virtual] |
Retrieve indexed entries of matrix column multiplied by alpha.
cNum | Column number. |
irows | Index list specifying rows. |
alpha | Scalar factor. |
col | Output column vector. |
Implements Matrix.
References BT_TRUE, isEqual(), leaDim, Indexlist::length, Indexlist::number, SUCCESSFUL_RETURN, and val.
real_t DenseMatrix::getNorm | ( | int_t | type = 2 | ) | const [virtual] |
Get the N-norm of the matrix
type | Norm type, 1: one-norm, 2: Euclidean norm. |
Implements Matrix.
References nCols, nRows, REFER_NAMESPACE_QPOASES, and val.
Referenced by getRowNorm().
returnValue DenseMatrix::getRow | ( | int_t | rNum, |
const Indexlist *const | icols, | ||
real_t | alpha, | ||
real_t * | row | ||
) | const [virtual] |
Retrieve indexed entries of matrix row multiplied by alpha.
rNum | Row number. |
icols | Index list specifying columns. |
alpha | Scalar factor. |
row | Output row vector. |
Implements Matrix.
References BT_TRUE, isEqual(), leaDim, Indexlist::length, nCols, Indexlist::number, SUCCESSFUL_RETURN, and val.
real_t DenseMatrix::getRowNorm | ( | int_t | rNum, |
int_t | type = 2 |
||
) | const [virtual] |
returnValue DenseMatrix::getRowNorm | ( | real_t * | norm, |
int_t | type = 2 |
||
) | const [virtual] |
Get the N-norm of all rows
norm | Norm of each row. |
type | Norm type, 1: one-norm, 2: Euclidean norm. |
Implements Matrix.
References getNorm(), leaDim, nCols, nRows, REFER_NAMESPACE_QPOASES, SUCCESSFUL_RETURN, and val.
BEGIN_NAMESPACE_QPOASES returnValue Matrix::getSparseSubmatrix | ( | const Indexlist *const | irows, |
const Indexlist *const | icols, | ||
int_t | rowoffset, | ||
int_t | coloffset, | ||
int_t & | numNonzeros, | ||
int_t * | irn, | ||
int_t * | jcn, | ||
real_t * | avals, | ||
BooleanType | only_lower_triangular = BT_FALSE |
||
) | const [virtual, inherited] |
Retrieve entries of submatrix in Harwell-Boeing sparse format. If irn, jcn, and avals are null, this only counts the number of nonzeros. Otherwise, numNonzeros containts the size of irn, jcn, and avals on entry, and the written number of entries on return.
irows | Index list specifying rows. |
icols | Index list specifying columns. |
rowoffset | Offset for row entries. |
coloffset | Offset for row entries. |
numNonzeros | Number of nonzeros in submatrix. |
irn | Row position of entries (as position in irows) plus rowoffset. |
jcn | Column position of entries (as position in irows) plus coloffset. |
avals | Numerical values of the entries. |
only_lower_triangular | if true, only the lower triangular portion is returned. This can only be true for symmetric matrices and if irows==jcols. |
References Indexlist::getLength(), and Indexlist::getNumberArray().
Referenced by SQProblemSchur::addConstraint(), Matrix::getSparseSubmatrix(), SQProblemSchur::removeBound(), and SQProblemSchur::resetSchurComplement().
returnValue Matrix::getSparseSubmatrix | ( | const Indexlist *const | irows, |
int_t | idx_icol, | ||
int_t | rowoffset, | ||
int_t | coloffset, | ||
int_t & | numNonzeros, | ||
int_t * | irn, | ||
int_t * | jcn, | ||
real_t * | avals, | ||
BooleanType | only_lower_triangular = BT_FALSE |
||
) | const [virtual, inherited] |
Retrieve entries of submatrix in Harwell-Boeing sparse format. If irn, jcn, and avals are null, this only counts the number of nonzeros. Otherwise, numNonzeros containts the size of irn, jcn, and avals on entry, and the written number of entries on return. This version retrieves one column.
irows | Index list specifying rows. |
idx_icol | Index list specifying columns. |
rowoffset | Offset for row entries. |
coloffset | Offset for row entries. |
numNonzeros | Number of nonzeros in submatrix. |
irn | Row position of entries (as position in irows) plus rowoffset. |
jcn | Column position of entries (as position in irows) plus coloffset. |
avals | Numerical values of the entries. |
only_lower_triangular | if true, only the lower triangular portion is returned. This can only be true for symmetric matrices and if irows==jcols. |
References Indexlist::getLength(), Indexlist::getNumberArray(), and Matrix::getSparseSubmatrix().
returnValue Matrix::getSparseSubmatrix | ( | int_t | idx_row, |
const Indexlist *const | icols, | ||
int_t | rowoffset, | ||
int_t | coloffset, | ||
int_t & | numNonzeros, | ||
int_t * | irn, | ||
int_t * | jcn, | ||
real_t * | avals, | ||
BooleanType | only_lower_triangular = BT_FALSE |
||
) | const [virtual, inherited] |
Retrieve entries of submatrix in Harwell-Boeing sparse format. If irn, jcn, and avals are null, this only counts the number of nonzeros. Otherwise, numNonzeros containts the size of irn, jcn, and avals on entry, and the written number of entries on return. This version retrieves one row.
idx_row | Row number. |
icols | Index list specifying columns. |
rowoffset | Offset for row entries. |
coloffset | Offset for row entries. |
numNonzeros | Number of nonzeros in submatrix. |
irn | Row position of entries (as position in irows) plus rowoffset. |
jcn | Column position of entries (as position in irows) plus coloffset. |
avals | Numerical values of the entries. |
only_lower_triangular | if true, only the lower triangular portion is returned. This can only be true for symmetric matrices and if irows==jcols. |
References Indexlist::getLength(), Indexlist::getNumberArray(), and Matrix::getSparseSubmatrix().
returnValue DenseMatrix::getSparseSubmatrix | ( | int_t | irowsLength, |
const int_t *const | irowsNumber, | ||
int_t | icolsLength, | ||
const int_t *const | icolsNumber, | ||
int_t | rowoffset, | ||
int_t | coloffset, | ||
int_t & | numNonzeros, | ||
int_t * | irn, | ||
int_t * | jcn, | ||
real_t * | avals, | ||
BooleanType | only_lower_triangular = BT_FALSE |
||
) | const [virtual] |
Retrieve entries of submatrix in Harwell-Boeing sparse format. If irn, jcn, and avals are null, this only counts the number of nonzeros. Otherwise, numNonzeros containts the size of irn, jcn, and avals on entry, and the written number of entries on return.
irowsLength | Number of rows. |
irowsNumber | Array with row numbers. |
icolsLength | Number of columns. |
icolsNumber | Array with column numbers. |
rowoffset | Offset for row entries. |
coloffset | Offset for row entries. |
numNonzeros | Number of nonzeros in submatrix. |
irn | Row position of entries (as position in irows) plus rowoffset. |
jcn | Column position of entries (as position in irows) plus coloffset. |
avals | Numerical values of the entries. |
only_lower_triangular | if true, only the lower triangular portion is returned. This can only be true for symmetric matrices and if irows==jcols. |
Implements Matrix.
References BT_FALSE, isZero(), leaDim, real_t, RET_INVALID_ARGUMENTS, SUCCESSFUL_RETURN, THROWERROR, and val.
BooleanType DenseMatrix::isDiag | ( | ) | const [virtual] |
BooleanType Matrix::needToFreeMemory | ( | ) | const [inline, inherited] |
Returns whether internal memory needs to be de-allocated.
References Matrix::freeMemory.
Referenced by duplicate(), SymDenseMat::duplicateSym(), ~DenseMatrix(), SparseMatrix::~SparseMatrix(), and SparseMatrixRow::~SparseMatrixRow().
returnValue DenseMatrix::print | ( | const char * | name = 0 | ) | const [virtual] |
Prints matrix to screen.
name | Name of matrix. |
Implements Matrix.
References nCols, nRows, REFER_NAMESPACE_QPOASES, and val.
returnValue DenseMatrix::times | ( | int_t | xN, |
real_t | alpha, | ||
const real_t * | x, | ||
int_t | xLD, | ||
real_t | beta, | ||
real_t * | y, | ||
int_t | yLD | ||
) | const [virtual] |
Evaluate Y=alpha*A*X + beta*Y.
xN | Number of vectors to multiply. |
alpha | Scalar factor for matrix vector product. |
x | Input vector to be multiplied. |
xLD | Leading dimension of input x. |
beta | Scalar factor for y. |
y | Output vector of results. |
yLD | Leading dimension of output y. |
Implements Matrix.
References GEMM, getMax(), nCols, nRows, SUCCESSFUL_RETURN, and val.
returnValue DenseMatrix::times | ( | const Indexlist *const | irows, |
const Indexlist *const | icols, | ||
int_t | xN, | ||
real_t | alpha, | ||
const real_t * | x, | ||
int_t | xLD, | ||
real_t | beta, | ||
real_t * | y, | ||
int_t | yLD, | ||
BooleanType | yCompr = BT_TRUE |
||
) | const [virtual] |
Evaluate matrix vector product with submatrix given by Indexlist.
irows | Index list specifying rows. |
icols | Index list specifying columns. |
xN | Number of vectors to multiply. |
alpha | Scalar factor for matrix vector product. |
x | Input vector to be multiplied. |
xLD | Leading dimension of input x. |
beta | Scalar factor for y. |
y | Output vector of results. |
yLD | Leading dimension of output y. |
yCompr | Compressed storage for y. |
Implements Matrix.
References BT_FALSE, BT_TRUE, isEqual(), Indexlist::iSort, isZero(), leaDim, Indexlist::length, nCols, Indexlist::number, SUCCESSFUL_RETURN, and val.
returnValue DenseMatrix::transTimes | ( | int_t | xN, |
real_t | alpha, | ||
const real_t * | x, | ||
int_t | xLD, | ||
real_t | beta, | ||
real_t * | y, | ||
int_t | yLD | ||
) | const [virtual] |
Evaluate Y=alpha*A'*X + beta*Y.
xN | Number of vectors to multiply. |
alpha | Scalar factor for matrix vector product. |
x | Input vector to be multiplied. |
xLD | Leading dimension of input x. |
beta | Scalar factor for y. |
y | Output vector of results. |
yLD | Leading dimension of output y. |
Implements Matrix.
References GEMM, getMax(), nCols, nRows, SUCCESSFUL_RETURN, and val.
returnValue DenseMatrix::transTimes | ( | const Indexlist *const | irows, |
const Indexlist *const | icols, | ||
int_t | xN, | ||
real_t | alpha, | ||
const real_t * | x, | ||
int_t | xLD, | ||
real_t | beta, | ||
real_t * | y, | ||
int_t | yLD | ||
) | const [virtual] |
Evaluate matrix transpose vector product.
irows | Index list specifying rows. |
icols | Index list specifying columns. |
xN | Number of vectors to multiply. |
alpha | Scalar factor for matrix vector product. |
x | Input vector to be multiplied. |
xLD | Leading dimension of input x. |
beta | Scalar factor for y. |
y | Output vector of results. |
yLD | Leading dimension of output y. |
Implements Matrix.
References BT_FALSE, BT_TRUE, isEqual(), Indexlist::iSort, isZero(), leaDim, Indexlist::length, Indexlist::number, SUCCESSFUL_RETURN, and val.
returnValue DenseMatrix::writeToFile | ( | FILE * | output_file, |
const char * | prefix | ||
) | const [virtual] |
Write matrix to file.
Implements Matrix.
References RET_NOT_YET_IMPLEMENTED, and THROWERROR.
BooleanType Matrix::freeMemory [protected, inherited] |
Indicating whether internal memory needs to be de-allocated.
Referenced by Matrix::doFreeMemory(), and Matrix::needToFreeMemory().
int_t DenseMatrix::leaDim [protected] |
Leading dimension.
Referenced by addToDiag(), SymDenseMat::bilinear(), diag(), getCol(), getRow(), getRowNorm(), getSparseSubmatrix(), isDiag(), times(), and transTimes().
int_t DenseMatrix::nCols [protected] |
Number of columns.
Referenced by addToDiag(), duplicate(), SymDenseMat::duplicateSym(), full(), getNorm(), getRow(), getRowNorm(), isDiag(), print(), times(), and transTimes().
int_t DenseMatrix::nRows [protected] |
Number of rows.
Referenced by addToDiag(), duplicate(), SymDenseMat::duplicateSym(), full(), getNorm(), getRowNorm(), isDiag(), print(), times(), and transTimes().
real_t* DenseMatrix::val [protected] |
Vector of entries.
Referenced by addToDiag(), SymDenseMat::bilinear(), diag(), duplicate(), SymDenseMat::duplicateSym(), free(), full(), getCol(), getNorm(), getRow(), getRowNorm(), getSparseSubmatrix(), isDiag(), print(), times(), and transTimes().