Wes Lowrie

OpenMP in Pyclaw

34

I have updated step3ds.f90 in PyClaw to be able to use openmp. Some small changes were necessary to classic/solver.py to get the openmp environment variable OMP_NUM_THREADS and to properly set...

Potential problem in flux3.f90 when using f-wave solver with capacity-form differencing

2

I believe there is an issue with limiters when using an f-wave solver with capacity functions. Take a look at flux3.f90: https://github.com/clawpack/pyclaw/blob/master/src/pyclaw/classic/flux3.f90#L226-263 When a Q-wave rpn3() is called `wave`, `s`,...

Bug in pyclaw/src/pyclaw/classic/flux3.f90 with capacity function scaling when transverse_waves==21 or 22

7

At `line 418` of flux3.f90 the gadd() variable is multiplied by `dtdx1d` and `dtdz`. `dtdxtd` has been divided by the capacity function in step3.f90, while `dtdz` has not. I believe...