PyBem

BEM code for solving the heat equation.

More precisely, we solve the Dirichlet problem in a smooth domain Omega with boundary Gamma and a time interval I = (0,4)

Find u: Omega x I --> R satisfying:

(partial_t-Delta)u = 0, in Omega x I

u = 0, at Omega x {t=0}

u|_{Gamma x I} = g, in Gamma x I,

Prerequisites

Running nosetests in the root directory will run all the unit tests.

If you want to see the library in action, you can run:

PYTHONPATH=.:test python -i test/functional.py

This will run a longer calculation and keep a plot of the result open.

Name		Name	Last commit message	Last commit date
Latest commit History 16 Commits
test		test
.gitignore		.gitignore
L2proj.py		L2proj.py
README.md		README.md
assembleMass.py		assembleMass.py
assembleStiffness.py		assembleStiffness.py
assembleVector.py		assembleVector.py
basis.py		basis.py
calcSol.py		calcSol.py
integrate.py		integrate.py
integrate1d.py		integrate1d.py
memoize.py		memoize.py
plotHeatMap.py		plotHeatMap.py
projection.py		projection.py
solve.py		solve.py
timekernel.py		timekernel.py
util.py		util.py