cristicbz / scid

Scientific library for the D programming language
Boost Software License 1.0
23 stars 8 forks source link

Transposing Triangular Matrices Modifies the Zero Element #86

Open dsimcha opened 12 years ago

dsimcha commented 12 years ago
void main() {
     import scid.matvec, std.stdio;
     auto mat = TriangularMatrix!double([[1.0, 2.0], [0.0, 3.0]]);
     eval(mat.t);
}

core.exception.AssertError@scid\storage\triangular.d(145): Modification of zero element in triangle matrix.

Looks like the fallback copying routine's copying the matrix vector-by-vector doesn't work in the triangular-to-triangular case.