Skip to content
memmett edited this page Feb 5, 2012 · 5 revisions

WENO

The focus of this group is to get WENO reconstructions working on mapped grids.

Git repo

Please use the 'hpc3-weno' branch (which is the default) of:

$ git clone [email protected]:memmett/pyclaw.git

Generating nonuniform Fortran WENO routines

The following is a short snippet of code to generate a Fortran callable WENO reconstructor to compute a 5th order accurate reconstruction (k=3) that operates on nonuniform grids. We set n=2 to say that we'll be reconstructing at two points in each cell (the left and right edges).

>>> import pyweno.symbolic
>>> import pyweno.functions
>>> f = pyweno.functions.FunctionGenerator('fortran')
>>> f.set_nonuniform(k=3, n=2)
>>> print f.generate('weno5', smoothness=True, weights=True, normalise=False)

We'll need to generate these for weno5 up to weno17 (as is done for the uniform WENO routines). To do so, we'll emulate the uniform weno.py.

Clone this wiki locally