Skip to content

Conversation

@freemin7
Copy link
Owner

@freemin7 freemin7 commented Jul 7, 2020

No description provided.

ChrisRackauckas and others added 30 commits May 29, 2020 01:46
Setup solve for adjoints to deprecate concrete_solve
More explanatory error message when missing algorithm
Update License [ci skip]
Use vecvecapply in dense calculate_solution_errors
This changes the norms on dual numbers to include the dual portions. This means that the solution to the ODE under automatic differentiation is not necessarily the same as without. However, SciML/SciMLSensitivity.jl#273 nicely highlights some edge cases that show that this is somewhat necessary. In that case, the integrator has zero error, and so time steps are unaffected by tolerances which causes gradients to explode since they are not tolerance controlled. This change adds them to the summation and the length, making the solution by dual numbers equivalent to the result of forward sensitivity analysis under the default norm.

@YingboMa what do you think?

@pkofod you might want to test this with Pumas. It should give it some more robustness.
Use JuliaGPU's GitLab CI setup
Co-authored-by: Yingbo Ma <[email protected]>
Co-authored-by: Yingbo Ma <[email protected]>
Co-authored-by: Yingbo Ma <[email protected]>
ChrisRackauckas and others added 29 commits June 14, 2020 20:55
Check pointer not values when deciding whether to remake
You know how sometimes you later think that some idea was a really really bad idea? Yes, this was one of them. "It's a pair of ODEs, so it's a pair of mass matrices!". Yeah, but that breaks all implicit solver fallbacks, and even when implicit methods specialized to second order ODEs exist, mass_matrix=I will still work as the special case to catch, so it's just better...

Anyways, it didn't matter until last week when ArrayPartition suddenly became compatible with all implicit solvers, except for the fact that it would fail because the mass matrix was by default a tuple. So, let's fix this and yay, 2nd order ODEs are now better.
Relax norm test (floating point numbers are really weird)
mass matrix I for DynamicalODEFunction
Make DiffEqArrayOperator ExponentialUtilities.jl compatible
…6-17-06-10-588-1862310616

CompatHelper: bump compat for "ChainRulesCore" to "0.9"
auto-convert integer inputs
IncrementingODEProblem and IncrementingODEFunction
@freemin7 freemin7 merged commit 4a331b8 into freemin7:master Jul 7, 2020
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.

6 participants