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

Polynomial reconstruction in SharpClaw #492

Open
wants to merge 16 commits into
base: master
Choose a base branch
from

Conversation

sanromd
Copy link
Contributor

@sanromd sanromd commented Dec 28, 2014

This PR adds a routine for polynomial based reconstruction in SharpClaw. In detail,

  • 4th and 6th order Polynomial reconstrcution
  • 8th and 10th order reconstruction are also included but have not been fully tested yet
  • renames "weno_order" to "reconstruction_order"
  • deallocate calls for WENO5 and Poly_comp
  • sets number of ghost cells for either Weno or Poly in solver
  • documentation

@ketch
Copy link
Member

ketch commented Dec 30, 2014

Just a reminder for anybody testing this -- you'll need to rebuild:

cd clawpack
pip install -e .

@@ -9,7 +9,7 @@ module ClawParams
integer :: num_dim, num_waves, index_capa

! Method-related parameters:
integer :: char_decomp, lim_type, weno_order
integer :: char_decomp, lim_type, reconstruction_order
Copy link
Member

Choose a reason for hiding this comment

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

Just a note: renaming weno_order to reconstruction_order is backward-incompatible, but it is a small enough change that I think it's okay for this to go into 5.3.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

@ketch, if needed for 5.4 I can add backward-compatibility...

Copy link
Member

Choose a reason for hiding this comment

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

Maybe that's a good idea; we can just emit a warning saying that the old attribute will be deprecated in future versions.

@ketch
Copy link
Member

ketch commented Dec 30, 2014

I just merged in the fix for the doctests (from master) so this is passing now.


Order of the WENO reconstruction. From 1st to 17th order (PyWENO)
Order of the reconstruction. From 1st to 17th order (PyWENO), and from
4 to 10 (even) in Poly
Copy link
Member

Choose a reason for hiding this comment

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

This should say

From 5th to 17th order (odd) for WENO reconstruction and from 4th to 10th order (even) for piecewise polynomial reconstruction.

Copy link
Member

Choose a reason for hiding this comment

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

Ah, I just noticed that it's wrong in master...

@ketch
Copy link
Member

ketch commented Dec 30, 2014

This looks great and will be very useful. It would of course be nice to have:

  • A test that uses the new reconstruction
  • At least some comments at the top of poly.f90 documenting the functions there.

@coveralls
Copy link

Coverage Status

Coverage remained the same when pulling 5f6ee25 on sanromd:dev-poly into e20a787 on clawpack:master.

@sanromd
Copy link
Contributor Author

sanromd commented Dec 30, 2014

@ketch working on

  • acoustics_1d test
  • comments on poly
  • documentation on solver

@mandli
Copy link
Member

mandli commented Oct 24, 2015

What's the status of this PR?

@sanromd
Copy link
Contributor Author

sanromd commented Oct 24, 2015

I made some modifications to it in my fork and rebased it... I haven't got
time to finalize testing. I will update it soon.

On Saturday, October 24, 2015, Kyle Mandli notifications@github.com wrote:

What's the status of this PR?


Reply to this email directly or view it on GitHub
#492 (comment).

Damián

ubi dubium ibi libertas

@mandli
Copy link
Member

mandli commented Oct 24, 2015

Thanks for the update. Would you like this included in the 5.3.1 release or would this be better for the 5.4.0 release?

@ketch
Copy link
Member

ketch commented Oct 25, 2015

If @sanromd can add a test or two soon, it would be great to get this into 5.4.0. It probably shouldn't go into 5.3.1 since it's a major feature and 5.4 is coming very soon anyway.

@coveralls
Copy link

Coverage Status

Changes Unknown when pulling 458933d on sanromd:dev-poly into * on clawpack:master*.

@coveralls
Copy link

coveralls commented Oct 5, 2016

Coverage Status

Coverage decreased (-0.01%) to 69.671% when pulling c0095d2 on sanromd:dev-poly into 314fd93 on clawpack:master.

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.

None yet

4 participants