This branch includes all the requirements to add the option of Riemannien subspace in RASR:
nonlinear_eigh() returns the subspace decomposition using pymanopt trustregion.
unit test of nonlinear_eigh() : for now it fails
Please review the code, cost function, gradien and hessian mainly because pymanopt thrustregion fails to converge to correct solution. This issue has been mentioned in #14 .
Please don't merge this branch until we succeeded to make it work ^^(the output in the test_pymanopt.py should be a permutation matrix).
This branch includes all the requirements to add the option of Riemannien subspace in RASR:
Please review the code, cost function, gradien and hessian mainly because pymanopt thrustregion fails to converge to correct solution. This issue has been mentioned in #14 .
Please don't merge this branch until we succeeded to make it work ^^(the output in the test_pymanopt.py should be a permutation matrix).