Use the integral() routine for L2 errors
Implement func2quad
Refactor the specific functions out
Refactor poisson3d_assembly.f90
Refactor a bit
Make poisson3d a test
Test exact numbers
Add another test
Export box dimensions
Add another test
Precalculate local element matrix