Fix bug 736: LDLT isPositive returns false for a positive semidefinite matrix
Add unit test covering this case.