Currently, DynamicSparseMatrix allows us to write to a coefficient outside the range, and leads to an unhelpful error later.

eg:

Code: Select all

DynamicSparseMatrix<float> A(5,5);
A.coeffRef(10,4) = 1;
cout << "A = "<<A<<"\n";


This leads to an error in the cout statement. Now, either A should reject the assignment, or change its dimension. The former would be better than the latter, I think.