GetSubMatrix& BaseMatrix::Column(int first_col) const
#else
GetSubMatrix BaseMatrix::Column(int first_col) const
#endif
{
REPORT
Tracer tr("SubMatrix(column)");
int c = first_col - 1;
if (c<0) Throw(SubMatrixDimensionException());
#ifdef TEMPS_DESTROYED_QUICKLY
GetSubMatrix* x = new GetSubMatrix(this, 0, -1, c, 1, false);
MatrixErrorNoSpace(x);
return *x;
#else
return GetSubMatrix(this, 0, -1, c, 1, false);
#endif
}
GetSubMatrix& BaseMatrix::Row(int first_row) const
#else
GetSubMatrix BaseMatrix::Row(int first_row) const
#endif
{
REPORT
Tracer tr("SubMatrix(row)");
int a = first_row - 1;
if (a<0) Throw(SubMatrixDimensionException());
#ifdef TEMPS_DESTROYED_QUICKLY
GetSubMatrix* x = new GetSubMatrix(this, a, 1, 0, -1, false);
MatrixErrorNoSpace(x);
return *x;
#else
return GetSubMatrix(this, a, 1, 0, -1, false);
#endif
}
GeneralMatrix* MatrixType::New(int nr, int nc, BaseMatrix* bm) const
{
// make a new matrix with the given attributes
Tracer tr("New"); GeneralMatrix* gm;
switch (attribute)
{
case Valid:
if (nc==1) { gm = new ColumnVector(nr); break; }
if (nr==1) { gm = new RowVector(nc); break; }
gm = new Matrix(nr, nc); break;
case Valid+Symmetric:
gm = new SymmetricMatrix(nr); break;
case Valid+Band:
{
MatrixBandWidth bw = bm->BandWidth();
gm = new BandMatrix(nr,bw.lower,bw.upper); break;
}
case Valid+Symmetric+Band:
gm = new SymmetricBandMatrix(nr,bm->BandWidth().lower); break;
case Valid+Upper:
gm = new UpperTriangularMatrix(nr); break;
case Valid+Diagonal+Symmetric+Band+Upper+Lower:
gm = new DiagonalMatrix(nr); break;
case Valid+Band+Upper:
gm = new UpperBandMatrix(nr,bm->BandWidth().upper); break;
case Valid+Lower:
gm = new LowerTriangularMatrix(nr); break;
case Valid+Band+Lower:
gm = new LowerBandMatrix(nr,bm->BandWidth().lower); break;
default:
Throw(ProgramException("Invalid matrix type"));
}
MatrixErrorNoSpace(gm); gm->Protect(); return gm;
}
GetSubMatrix& BaseMatrix::Rows(int first_row, int last_row) const
#else
GetSubMatrix BaseMatrix::Rows(int first_row, int last_row) const
#endif
{
REPORT
Tracer tr("SubMatrix(rows)");
int a = first_row - 1; int b = last_row - first_row + 1;
if (a<0 || b<0) Throw(SubMatrixDimensionException());
// allow zero rows or columns
#ifdef TEMPS_DESTROYED_QUICKLY
GetSubMatrix* x = new GetSubMatrix(this, a, b, 0, -1, false);
MatrixErrorNoSpace(x);
return *x;
#else
return GetSubMatrix(this, a, b, 0, -1, false);
#endif
}
GetSubMatrix& BaseMatrix::Columns(int first_col, int last_col) const
#else
GetSubMatrix BaseMatrix::Columns(int first_col, int last_col) const
#endif
{
REPORT
Tracer tr("SubMatrix(columns)");
int c = first_col - 1; int d = last_col - first_col + 1;
if (c<0 || d<0) Throw(SubMatrixDimensionException());
// allow zero rows or columns
#ifdef TEMPS_DESTROYED_QUICKLY
GetSubMatrix* x = new GetSubMatrix(this, 0, -1, c, d, false);
MatrixErrorNoSpace(x);
return *x;
#else
return GetSubMatrix(this, 0, -1, c, d, false);
#endif
}
GeneralMatrix* Matrix::MakeSolver()
{
REPORT
GeneralMatrix* gm = new CroutMatrix(*this);
MatrixErrorNoSpace(gm); gm->ReleaseAndDelete(); return gm;
}
