Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

error in update_spatial_parallel #32

Open
fxpena opened this issue Oct 23, 2018 · 1 comment
Open

error in update_spatial_parallel #32

fxpena opened this issue Oct 23, 2018 · 1 comment

Comments

@fxpena
Copy link

fxpena commented Oct 23, 2018

Hi I'm using the demo for large_data_1p. After initializing the residual, the script does update_spatial_parallel and I got a warning:

Warning: Matrix is singular to working precision.

In com (line 24)
In determine_search_location (line 58)
In Sources2D/update_spatial_parallel (line 66)

Then an error:

Error using determine_search_location (line 66)
Input matrix contains NaN or Inf.

Error in Sources2D/update_spatial_parallel (line 66)
IND = sparse(logical(determine_search_location(obj.A, search_method, options)));

I'm having trouble understanding why this error occurs and what it means. Does it have to do with the step where neurons are picked from the residual? The seeds that were picked in that step hardly coincide with the seeds picked during initialization.

@edonzis
Copy link

edonzis commented Jan 17, 2019

I'm having the same issue, but I'm running msRunCNMFE_large. When I track down the error it seems that the coordinates (??) for some of the neurons end up being NaN. The error comes in the parafor loop starting at line 66 of determine_search_location.m.
When trying to execute line 69 "[V,D] = eig(Vr{i});" Vr{i} is NaN because cm(i,1) in the previous line is NaN. Is there an easy way to remove neurons that are listed as NaN? I'm not sure the root cause of these neurons being listed as NaN.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants