stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Dynare++: drop unused codepaths

time-shift
Sébastien Villemot 2019-03-29 15:26:04 +01:00
parent 44b40bbf09
commit 39896d8d0a
No known key found for this signature in database
GPG Key ID: 2CECE9350ECEBE4A
4 changed files with 0 additions and 64 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
dynare++/sylv/cc/BlockDiagonal.cc
View File

@ -30,12 +30,6 @@ BlockDiagonal::BlockDiagonal(const QuasiTriangular &t)
}
}
BlockDiagonal::BlockDiagonal(int p, const BlockDiagonal &b)
: QuasiTriangular(p, b),
row_len(b.row_len), col_len(b.col_len)
{
}
/* Put zeroes to right upper submatrix whose first column is defined
by ‘edge’ */
void

1
dynare++/sylv/cc/BlockDiagonal.hh
View File

@ -15,7 +15,6 @@ class BlockDiagonal : public QuasiTriangular
std::vector<int> row_len, col_len;
public:
BlockDiagonal(ConstVector d, int d_size);
BlockDiagonal(int p, const BlockDiagonal &b);
BlockDiagonal(const BlockDiagonal &b) = default;
BlockDiagonal(const QuasiTriangular &t);
BlockDiagonal &operator=(const QuasiTriangular &t)

54
dynare++/sylv/cc/QuasiTriangular.cc
View File

@ -462,66 +462,12 @@ QuasiTriangular::setMatrix(double r, const QuasiTriangular &t)
getData().add(r, t.getData());
}
void
QuasiTriangular::setMatrixViaIter(double r, const QuasiTriangular &t)
{
register double rr = r;
diag_iter dil = diag_begin();
const_diag_iter dir = t.diag_begin();
for (; dil != diag_end(); ++dil, ++dir)
{
dil->getAlpha() = rr*(*dir->getAlpha());
if (!dil->isReal())
{
dil->getBeta1() = rr*dir->getBeta1();
dil->getBeta2() = rr*dir->getBeta2();
}
col_iter cil = col_begin(*dil);
const_col_iter cir = t.col_begin(*dir);
for (; cil != col_end(*dil); ++cil, ++cir)
if (dil->isReal())
*cil = rr*(*cir);
else
{
cil.a() = rr*cir.a();
cil.b() = rr*cir.b();
}
}
}
void
QuasiTriangular::addMatrix(double r, const QuasiTriangular &t)
{
getData().add(r, t.getData());
}
void
QuasiTriangular::addMatrixViaIter(double r, const QuasiTriangular &t)
{
register double rr = r;
diag_iter dil = diag_begin();
const_diag_iter dir = t.diag_begin();
for (; dil != diag_end(); ++dil, ++dir)
{
dil->getAlpha() = (*dil->getAlpha()) + rr*(*dir->getAlpha());
if (!dil->isReal())
{
dil->getBeta1() += rr*dir->getBeta1();
dil->getBeta2() += rr*dir->getBeta2();
}
col_iter cil = col_begin(*dil);
const_col_iter cir = t.col_begin(*dir);
for (; cil != col_end(*dil); ++cil, ++cir)
if (dil->isReal())
*cil += rr*(*cir);
else
{
cil.a() += rr*cir.a();
cil.b() += rr*cir.b();
}
}
}
void
QuasiTriangular::addUnit()
{

3
dynare++/sylv/cc/QuasiTriangular.hh
View File

@ -441,9 +441,6 @@ private:
void multaKron(KronVector &x, const ConstKronVector &b) const;
/* x = x + (thisᵀ⊗I)·b */
void multaKronTrans(KronVector &x, const ConstKronVector &b) const;
/* implementation via iterators, useful for large matrices */
void setMatrixViaIter(double r, const QuasiTriangular &t);
void addMatrixViaIter(double r, const QuasiTriangular &t);
/* hide noneffective implementations of parents */
void multsVec(Vector &x, const ConstVector &d) const;
void multsVecTrans(Vector &x, const ConstVector &d) const;