Skip to content

Support better logging for singular cases in integration#3731

Open
Shreyas-Ekanathan wants to merge 5 commits into
SciML:masterfrom
Shreyas-Ekanathan:singularity-logging
Open

Support better logging for singular cases in integration#3731
Shreyas-Ekanathan wants to merge 5 commits into
SciML:masterfrom
Shreyas-Ekanathan:singularity-logging

Conversation

@Shreyas-Ekanathan

Copy link
Copy Markdown
Contributor

Checklist

  • Appropriate tests were added
  • Any code changes were done in a way that does not break public API
  • All documentation related to code changes were updated
  • The new code follows the
    contributor guidelines, in particular the SciML Style Guide and
    COLPRAC.
  • Any new documentation only uses public API

Additional context

Allows SciMLBase to freshly calculate the Jacobian for most accurate analysis. Supports SciML/SciMLBase.jl#1379

Comment on lines +653 to +654
for i in axes(jac, 1), j in axes(jac, 2)
if !isfinite(jac[i, j]) || abs(jac[i, j]) > 1e6

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If jac is a sparse matrix, then this will be super expensive. jac[i,j] triggers a lookup on every call and we even iterate over all zeros.

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

Successfully merging this pull request may close these issues.

2 participants