FETK  
The Finite Element ToolKit  
Developed by the MCP Research Group at the UCSD Center for Computational Mathematics 


Overview PMG is an algebraic multilevel code written in several languages (FORTRAN, C, C++, and CC++). The original FORTRAN/C/C++ versions were written for sequential computers, whereas the more recent CC++ port is a distributed memory parallel implementation. All of the versions of PMG are designed to numerically approximate the solutions of linear and nonlinear elliptic partial differential equations in threedimensional (logically) bricklike domains in an extremely efficient and robust way. To accomplish this task, PMG employs a nonuniform Cartesian mesh of the user's choice, boxmethods or finite element methods for discretization, algebraic (or geometric) construction of lowerdimensional subspace problems, dampedinexact Global Newton methods for nonlinearities, and a jacobian multilevel iteration based on the algebraic (or geometric) subspace hierarchy. The code is applicable to threedimensional nonlinear Poissonlike scalar equations, allowing for nonlinearities that depend on the scalar unknown (but not on derivatives); nonlinearities of this type occur in the PoissonBoltzmann equation as well as in other applications, such as the driftdifussion equations in semiconductor modeling. The class of problems for which the code is applicable can be extended with suitable modifications to the code (the methods need not be modified). For strongly elliptic equations with selfadjoint differential operators and smooth coefficients, the solution method is provably optimal (meaning O(N), where N mesh points are used on the finest userspecified mesh). For strongly elliptic selfadjoint operators with simple coefficient discontinuities, the solution method is provably O(N log N). For strongly elliptic selfadjoint operators with arbitrary discontinuities, the solution method is provably convergent, and emperically appears to retain the O(N log N) behavior. At the heart of the code is the algebraic construction of the lowerdimensional subspace problems, which is performed by a Maplegenerated code. This code was generated by Maple using a symbolic stencil calculus. Usage To use the code as provided, the user defines the problem coefficients, the mesh points lying on the finest nonuniform Cartesian mesh along with the basic domain geometry, as well as the boundary conditions and types, in a very simple and welldefined way. All necessary workspace is passed in by the user; the code will indicate exactly how much workspace is required for a particular problem configuration. The nonlinearity must be defined in a subroutine supplied by the user which is to be called by the code, and a second subroutine must also be supplied which provides the derivative of the nonlinear term for use in forming the Jacobian required by the Newton solver. Using the solver involves a single subroutine call with a simple and understandable argument list. Discretization The code is restricted to ``logically'' nonuniform Cartesian threedimensional meshes due to the multilevel solver employed, and the boxmethod discretization routine we provide requires a ``physically'' nonuniform Cartesian mesh as well. This discretization routine allows arbitrary Dirichlet or Neumann boundary conditions at any boundary point on a threedimensional rectangular box, allowing for specification of arbitrary contact placements and geometries, as may be necessary in semiconductor modeling. Although we provide a boxmethod discretization with the code, the multilevel method can form the coarse level equations algebraically, the only requirement on the fine level discretization being that the equations be representable as stencils, and that the nonlinear term be diagonal. This is always true for either box or finite element (with mass lumping) discretization of the above class of equations on a logically nonuniform Cartesian mesh. Algebraic and geometric construction of coarse problems The linear multilevel method at the core of the code is designed specifically for discontinuous coefficients as occur in material interface problems. Techniques employed include coefficient averaging methods or algebraric (Galerkin) methods, combined with operatorinduced prolongation procedures which attempt to enforce flux conservation at box boundaries when a grid function is interpolated from a coarse to a fine mesh. The averaging methods are usually sufficient for many problems (such as the PoissonBoltzmann equation), but particularly complex discontinuities (i.e., complex molecular surfaces) seem to require the Galerkin approach for robustness; this is computationally more costly and requires more memory than the averaging approach (sevendiagonal matrices produced by the boxmethod on the fine mesh expand to twentysevendiagonal matrices on coarse meshes), but it always converges rapidly. Galerkin coefficient expressions are extremely complex in three dimensions, and were obtained with the help of the MAPLE symbolic manipulation package. Multilevel linear solver The user is given many choices in configuring the multilevel method, including choice of smoother, number of levels employed, iteration strategy, choice of prolongation operator, coarse problem formulation, as well as iterative (CG) or direct (banded LINPACK) solve on the coarsest mesh. A conjugate gradient method is also employed to accelerate the multilevel methods, and while this approach is mathematically guaranteed to improve the convergence rate of the multilevel method with a corresponding increase in computational labor, it is also the case that in practice the CPU time to solve the problem is usually reduced using the acceleration (i.e., the acceleration more than pays for itself). Nonlinear solver Nonlinear problems are handled with a robust (essentially globally convergent) inexactNewton solver coupled with the fast linear multilevel method. Configuration options include a choice of several stopping criteria, degree of inexactness in the approximate Newton solves, and damping strategy, in addition to all possible configuration options for the multilevel Jacobian system solver. A nonlinear operatorbased prolongation procedure has been developed for the nonlinear case as well. Also included in the code are several implementations of more traditional nonlinear methods, including nonlinear SOR (SORouter iteration with onedimensional Newtonsolves), nonlinear CG (the FletherReeves nonlinear extension to the linear HestenesSteifel algorithm), and a ``true'' nonlinear multilevel method (the ``FAS'' or ``nonlinear multigrid'' algorithm). PoissonBoltzmann equationspecific code An interface routine is provided for using the code specifically as a nonlinear PoissonBoltzmann equation solver. To handle the severe numerical problems occurring with nonlinearities of exponentialtype present in the PoissonBoltzmann equation, we developed argumentcapping functions which avoid nonvectorizable statements. Calls to the standard intrinsic functions are replaced by these modified functions, and overflows are successfully avoided during early transient iterations without loosing the execution efficiency of the intrinsic functions. The C/C++/CC++ version and parallel computing The C/C++/CC++ version of PMG runs sequentially or in parallel on a wide range of hardware. Parallelism is obtained on workstation clusters as well as supercomputers such as the SP2 through the CC++ (Compositional C++) extensions. The (sequential) FORTRAN 77 version There is also a sequential FORTRAN 77 version of PMG (which we refer to as MG) with much of the core numerical capabilities as the C/C++/CC++ version. The F77 version installs itself in single or double precision on the specified machine, inserting the appropriate compiler directives automatically, and runs without modification on most current systems, including the Cray C90, Cray YMP, Convex C3, Convex C240, Convex C220, IBM RS/6000, Sun 4, Sun 3, and most other systems with a standard FORTRAN 77 compiler. The X Interface in C An Xinterface has been written (in C) for running the numerical code. The interface is completely disconnected from the numerical code in the sense that it provides only input files which are read by the numerical code, executes the numerical code, and then reads and processes the output files produced by the numerical code. The Xinterface runs as is on any machine that has a standard C compiler and the standard X11 distribution Release 4 or newer. The interface uses only the Athena Widget set (along with Xlib and the Xt intrinsics), so that the interface is also quite portable. Obtaining PMG The C++ version of PMG was developed by Michael Holst at Caltech, based on his initial FORTRAN version developed at the University of Illinois. The CC++ extension of PMG was developed by John Garnett and Carl Kesselman at Caltech. PMG is copylefted and distributed under the GNU General Public License (GPL). (A copy of the GPL should be included with your copy of PMG.) The PMG source can be downloaded from the FETK Download Page. 