Today I added preliminary support for 3D bilinear and linear forms using tetrahedral P1 elements. Unfortunately, the visualization is based on MATLAB using MATLAB.jl Julia package (https://github.com/lindahua/MATLAB.jl) and the mesh generation is also based on a MATLAB script package iso2mesh (http://iso2mesh.sourceforge.net/cgi-bin/index.cgi) which is again based on tetgen (http://wias-berlin.de/software/tetgen/).
Thus, I haven't yet figured out how I should package the code properly for distribution. Meanwhile, here is a screenshot of solving the Poisson equation
$$\Delta u = 1$$ on a unit cube $[0,1]^3$ with zero Dirichlet boundary conditions and the corresponding mesh:
 |
| Cross section at $z=0.5$. |
 |
| The mesh. Color indicates simply the $z$-coordinate. |
