QFEPDE - A Simple Finite Element Program

QFEPDE is a simple finite element program, for a two dimensional spatial domain, with no time dependence, using quadratic elements.

Files you may copy include:

• qfepde.f90, the source code.
• qfepde.out, sample problem output.

The list of routines includes:

• QFEPDE is a simple finite element program
• SETINT sets up a Gauss quadrature rule.
• GEOM carries out a variety of computations related to the geometry.
• ASEM calculates the stiffness matrix and load vector.
• QBF evaluates the shape functions and their derivatives.
• AA evaluates the coefficient function A(X,Y) in the differential equation.
• UEX returns the exact solution of the differential equation.
• UEXX returns the X derivative of the exact solution.
• UEXY returns the Y derivative of the exact solution.
• RHS returns the value of the right hand side of the PDE.
• UBC returns the Dirichlet boundary condition at any point on the boundary.
• BEPP performs Gauss elimination of a banded linear system.
• EROR calculates L2 and H1 errors in the solution.
• PRISOL prints the solution as a table of 5 columns.

Return to the FORTRAN software page.