Skip to content

Conversation

@brownbaerchen
Copy link
Contributor

@brownbaerchen brownbaerchen commented Jan 19, 2023

Simple script that solves the non-linear equations implicitly. I also added a scipy reference solution analogously to the van der Pol problem. I used this reference solution to verify the implementation of the Jacobian. The solution also has the familiar butterfly shape.
I tried to write everything such that the changes going to v6 will be minimal, that is to say I already put in defaults instead of essential keys.

I added a hook for computing the global error at the end of the run. Because the step sizes get redistributed between the last post_step hook and the post_run hook, this is a bit more messy than anticipated, but it is also tested.

@lisawim, now is your turn to point out the many places where I can improve ;)

@brownbaerchen brownbaerchen requested a review from lisawim January 19, 2023 17:59
@codecov
Copy link

codecov bot commented Jan 19, 2023

Codecov Report

Base: 68.50% // Head: 68.66% // Increases project coverage by +0.15% 🎉

Coverage data is based on head (5dfabbf) compared to base (c2c066a).
Patch coverage: 91.72% of modified lines in pull request are covered.

Additional details and impacted files
@@            Coverage Diff             @@
##           master     #261      +/-   ##
==========================================
+ Coverage   68.50%   68.66%   +0.15%     
==========================================
  Files         226      229       +3     
  Lines       19510    19643     +133     
==========================================
+ Hits        13365    13487     +122     
- Misses       6145     6156      +11     
Impacted Files Coverage Δ
pySDC/projects/Resilience/Lorenz.py 84.50% <84.50%> (ø)
pySDC/implementations/hooks/log_errors.py 100.00% <100.00%> (ø)
pySDC/implementations/problem_classes/Lorenz.py 100.00% <100.00%> (ø)

Help us with your feedback. Take ten seconds to tell us how you rate us. Have a feature suggestion? Share it here.

☔ View full report at Codecov.
📢 Do you have feedback about the report comment? Let us know in this issue.

@lisawim
Copy link
Contributor

lisawim commented Jan 20, 2023

Looks really good to me! As usual, it looks all very clean coded. But there is one thing I noticed: Due to consistency it makes sense to rename the file Lorentz.py to Lorenz.py ;) That's all!

@pancetta pancetta merged commit 8b987cd into Parallel-in-Time:master Jan 20, 2023
@pancetta
Copy link
Member

Merged it. You could consider testing the MPI part or removing it (if you don't need it anyway).

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.

3 participants