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

Possible incorrectness in VDP example Jacobian ordering / docs #6

Open
3 tasks
skyreflectedinmirrors opened this issue May 2, 2018 · 3 comments
Open
3 tasks

Comments

@skyreflectedinmirrors
Copy link
Collaborator

CC: @kyleniemeyer @andrewalferman

So Andrew's work has show that the Jacobian ordering in the VDP problem is possibly incorrect (and at the very least the docs are confusing / misleading)

Todo:

  • Verify expected Jacobian order from pyJac is what is described in accelerInt docs
  • Verify that VDP example uses correct Jacobian order
  • Check that "Column major" references in docs correspond to expected VDP order

I think we actually ran into this issue with @stonecp a few months ago?

see the fixed oregonator Jacobian here

@stonecp
Copy link

stonecp commented May 3, 2018

I did have trouble using the Jacobian. The matrix index ordering is confusing as it seems to depend upon the target platform. It only stores K-1 species (using mass conservation) but that is not consistent with the RHS function itself.

@kyleniemeyer
Copy link
Member

Bump?

@skyreflectedinmirrors
Copy link
Collaborator Author

@kyleniemeyer So I had updated the pyJac docs previously, and I believe that the documentation branch of accelerInt has fixed this problem, but I believe I was waiting on merging Andrew's Oregonator example before I merged this?

I'll try to take a look at it in the next few days

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

3 participants