MULTIGRID and ADAPTIVE MESHING

This page is to figure out what we have to do to fipy to get multigrid and adaptinve meshing working.

Questions that need answering

• Should we pull in an algebraic multigrid solver before

embarking on a geometric multigrid?

• Is it useful to implement geometric multigrid without

having adaptive griding?

• Have we noticed a trend in code use in the phase field

community?

Review of packages

The first step has to be to figure out which code base to start from. Are they Python wrapped? How would they fit in? Would it require updating Variable, Mesh and Terms in a big way? Can the package work mostly inside our paradigm?

PyTrilinos

​http://software.sandia.gov/trilinos/packages/pytrilinos/

Trilinos comes with the ML package that Roman wrapped at one point and so is probably worth looking at.

MG net

​http://www.mgnet.org/

Paramesh

​http://www.physics.drexel.edu/~olson/paramesh-doc/Users_manual/amr.html ​http://www.physics.drexel.edu/~olson/paramesh-doc/Users_manual/multigrid.html ​http://sourceforge.net/projects/paramesh/

The multigrid support is still new. Distributed under the NASA license, like NIST probably. Licence is completly open.

Chombo

​http://seesar.lbl.gov/ANAG/chombo/

Not Python.

PyTrilinos

​http://software.sandia.gov/trilinos/packages/pytrilinos/index.html

Has mutigrid and python wrappers

Pyro

​http://www.ucolick.org/~zingale/pyro/

Gerris

​http://gfs.sourceforge.net/

Has multigrid.

PETSc

​http://www-unix.mcs.anl.gov/petsc/petsc-as/

Don't see any multigrid solvers in Petc.

Other Resources

Need a link here StoreHere?