void GetSubMatrix::operator<<(const BaseMatrix& bmx)
{
REPORT
Tracer tr("SubMatrix(<<)"); GeneralMatrix* gmx = 0;
Try
{
SetUpLHS(); gmx = ((BaseMatrix&)bmx).Evaluate();
if (row_number != gmx->Nrows() || col_number != gmx->Ncols())
Throw(IncompatibleDimensionsException());
MatrixRow mrx(gmx, LoadOnEntry);
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--)
{
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Copy(mrx); mr.Next(); mrx.Next();
}
gmx->tDelete();
}
CatchAll
{
if (gmx) gmx->tDelete();
ReThrow;
}
}
void GetSubMatrix::operator+=(const BaseMatrix& bmx)
{
REPORT
Tracer tr("SubMatrix(+=)"); GeneralMatrix* gmx = 0;
// MatrixConversionCheck mcc; // Check for loss of info
Try
{
SetUpLHS(); gmx = ((BaseMatrix&)bmx).Evaluate();
if (row_number != gmx->Nrows() || col_number != gmx->Ncols())
Throw(IncompatibleDimensionsException());
if (gm->type().is_symmetric() &&
( ! gmx->type().is_symmetric() || row_skip != col_skip) )
Throw(ProgramException("Illegal operation on symmetric"));
MatrixRow mrx(gmx, LoadOnEntry);
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--)
{
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Check(mrx); // check for loss of info
sub.Add(mrx); mr.Next(); mrx.Next();
}
gmx->tDelete();
}
CatchAll
{
if (gmx) gmx->tDelete();
ReThrow;
}
}
void GetSubMatrix::operator*=(Real r) {
REPORT
Tracer tr("SubMatrix(*= or /= Real)");
MatrixConversionCheck mcc; // Check for loss of info
Try {
SetUpLHS();
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--) {
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Multiply(r); mr.Next();
}
#ifdef TEMPS_DESTROYED_QUICKLY
delete this;
#endif
}
CatchAll {
#ifdef TEMPS_DESTROYED_QUICKLY
delete this;
#endif
ReThrow;
}
}
void GetSubMatrix::operator-=(const BaseMatrix& bmx) {
REPORT
Tracer tr("SubMatrix(-=)"); GeneralMatrix* gmx = 0;
MatrixConversionCheck mcc; // Check for loss of info
Try {
SetUpLHS(); gmx = ((BaseMatrix&)bmx).Evaluate();
if (row_number != gmx->Nrows() || col_number != gmx->Ncols())
Throw(IncompatibleDimensionsException());
MatrixRow mrx(gmx, LoadOnEntry);
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--) {
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Check(mrx); // check for loss of info
sub.Sub(mrx); mr.Next(); mrx.Next();
}
gmx->tDelete();
#ifdef TEMPS_DESTROYED_QUICKLY
delete this;
#endif
}
CatchAll {
if (gmx) gmx->tDelete();
#ifdef TEMPS_DESTROYED_QUICKLY
delete this;
#endif
ReThrow;
}
}
void GetSubMatrix::operator=(Real r)
{
REPORT
Tracer tr("SubMatrix(=Real)");
SetUpLHS();
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--)
{
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Copy(r); mr.Next();
}
}
void GetSubMatrix::operator<<(const int* r)
{
REPORT
Tracer tr("SubMatrix(<<int*)");
SetUpLHS();
if (row_skip+row_number > gm->Nrows() || col_skip+col_number > gm->Ncols())
Throw(SubMatrixDimensionException());
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--)
{
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Copy(r); mr.Next();
}
}
void GetSubMatrix::inject(const GeneralMatrix& gmx)
{
REPORT
Tracer tr("SubMatrix(inject)");
SetUpLHS();
if (row_number != gmx.Nrows() || col_number != gmx.Ncols())
Throw(IncompatibleDimensionsException());
MatrixRow mrx((GeneralMatrix*)(&gmx), LoadOnEntry);
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--)
{
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Inject(mrx); mr.Next(); mrx.Next();
}
}
void GetSubMatrix::operator+=(Real r)
{
REPORT
Tracer tr("SubMatrix(+= or -= Real)");
// MatrixConversionCheck mcc; // Check for loss of info
Try
{
SetUpLHS();
MatrixRow mr(gm, LoadOnEntry+StoreOnExit+DirectPart, row_skip);
// do need LoadOnEntry
MatrixRowCol sub; int i = row_number;
while (i--)
{
mr.SubRowCol(sub, col_skip, col_number); // put values in sub
sub.Check(); // check for loss of info
sub.Add(r); mr.Next();
}
}
CatchAll
{
ReThrow;
}
}
