Implement OnTheLeft/OnTheRight for SparseTriangularView::solve
Submitted by Henrik Friberg
Assigned to Nobody
Link to original bugzilla bug (#835)
Version: 3.2
Description
The following does not compile:
Matrix3d A;
RowVector3d bT, xT;
A << 1, 2, 3,
4, 5, 6,
7, 8, 9;
bT << 1, 2, 3;
xT = A.triangularView<Upper>().solve<OnTheRight>(bT); // WORKS FINE!
printf("(%g, %g, %g)", xT(0), xT(1), xT(2));
spA.triangularView<Upper>().solve<OnTheRight>(bT); // COMPILE ERROR!
SparseMatrix<double> spA = A.sparseView();SparseTriangularView::solve accepts neither <OnTheLeft> nor <OnTheRight> as template argument, but neglecting the argument compiles fine..