[Matrix-SIG] Anybody doing/having done FE PDE solving with NumPy?

Lorenzo M. Catucci lorenzo@argon.roma2.infn.it
Mon, 23 Mar 1998 18:58:54 +0100 (CET)

So good... I managed to put two true acronyms and a false one on the same
subject line... maybe I'll do better some time... Now, I get serious: I'm
going to do some simulations, and part of my problem is to build the
static electric field inside a box, with some iregular shaped object
inside. I think that the best way to do this part, is to solve via finite
element method the Laplace equation $\nabla^2 \phi=0$ with
Dirichlet's boundary conditions on the fixed potential surfaces. Then,
I'll go to do some molecular dynamics for objects inside this field. Best
option would be to be able to describe the FEM basis directly within
python, but I still haven't managed to find a clear recipe for those
nasty functions, which need be (at least) derivable, since I'll need the
potential field gradient... Somebody could help me and tell which route
to go?