CompoundSymMatrix::CompoundSymMatrix(const CompoundSymMatrixSpace* owner_space)
:
SymMatrix(owner_space),
owner_space_(owner_space),
matrices_valid_(false)
{
for (Index irow=0; irow<NComps_Dim(); irow++) {
std::vector< SmartPtr<Matrix> > row(irow+1);
std::vector< SmartPtr<const Matrix> > const_row(irow+1);
comps_.push_back(row);
const_comps_.push_back(const_row);
}
}
CompoundMatrix::CompoundMatrix(const CompoundMatrixSpace* owner_space)
:
Matrix(owner_space),
owner_space_(owner_space),
matrices_valid_(false)
{
std::vector<SmartPtr<Matrix> > row(NComps_Cols());
std::vector<SmartPtr<const Matrix> > const_row(NComps_Cols());
for (Index irow=0; irow<NComps_Rows(); irow++) {
const_comps_.push_back(const_row);
comps_.push_back(row);
}
}