This is a small PR which changes the starting point used for D in the cSFS2D algorithm to a starting point similar to that used for the other methods. As the cholesky decomposition is then performed on D, care had to be taken to ensure that D was positive definite (rather than non-negative definite) for this method, so that the np.linalg.cholesky function would not error.
This is a small PR which changes the starting point used for D in the
cSFS2Dalgorithm to a starting point similar to that used for the other methods. As the cholesky decomposition is then performed on D, care had to be taken to ensure that D was positive definite (rather than non-negative definite) for this method, so that thenp.linalg.choleskyfunction would not error.