quadrilateral
FENICS Solution of Poisson Equation

quadrilateral, a FENICS code which uses the finite element method to solve a version of the Poisson equation over the unit square, adapted from the FENICS tutorial by Langtangen and Logg.

Note that I have installed FENICS using Docker, and so to run this script I issue the commands:

1. cd $HOME/fenicsproject/quadrilateral
2. fenicsproject run
3. python3 quadrilateral.py
4. exit

Licensing:

The GNU LGPL license.

Source Code:

• quadrilateral.py a script which sets up and solves the problem.
• quadrilateral.txt the output file.
• quadrilateral000000.vtu a VTK plot of the solution.
• quadrilateral.pvd an index file for the plots.
• quadrilateral_mesh.png a plot of the mesh, created by running ParaView on quadrilateral.pvd, and selecting "Surface with Edges". The point of this plot is to notice that the elements are, indeed, quadrilateral.