I just pushed python wrappers + a simple schroedinger solver for our
C++ hermes. See the README for instructions how to turn them on.
Example:
$ examples/schroedinger/schroedinger
-------------------------------------------
This is Hermes1D - a free ODE solver
based on the hp-FEM and Newton's method,
developed by the hp-FEM group at UNR
and distributed under the BSD license.
For more details visit http://hpfem.org/.
-------------------------------------------
Importing hermes1d
Python initialized
inverting
solving
sorting the eigenvalues
eigenvalues:
-0.124999914632
-0.0556812382878
-0.0351854495991
-0.0111430096155
Output written to solution.gp.
Done.
The main program is in C++, but it uses numpy for solving the
eigenproblem and matplotlib for plotting (a plot will automatically
arise on the screen, see the attachement).
Ondrej