It looks like your code is wrong, and numpy 1.12 happened to let you get
away with it
This line:
evals = evals[evals > tolerance]
Reduces the eigenvalues to only those which are greater than the tolerance
When you do U[:, evals > tolerance], evals > tolerance is just going
to be [True,
True,
Dear ALL,
I am trying to port an eigenalysis function that runs smoothly on Numpy
1.12 but fail miserably on Numpy 1.13 or higher with the dreadful error
"boolean index did not match indexed array along dimension 1".
Here is a fragment of the code, where the error occurrs:
evals, evecs =
Draft agenda: https://hackmd.io/f_e6dnssTkuIC0Pa0TEV6w?view
There is a section for community suggested topics, feel free to join the
conversation and add in topics that need attention.
BIDS team
___
NumPy-Discussion mailing list