Fix nested_eval<Product<> > which wrongly returned a Product<> expression
Tridiagonalization::diagonal() and ::subDiagonal() did not work. Added unit-test