|
The above results are obtained using a computer program that was developed
for the Maxwell equations and the eigenvalue
problem in the frequency domain.
It is written in C++, in the framework of the
AggieFEM finite element
library which supports complex geometries, local refinement, multigrid
preconditioning and OpenGL visualization.
The code is based on the solvers for the magnetostatic and electrostatic
problems.
It works on triangular, tetrahedral and hexahedral meshes.
It provides an eigenvalue solver, which allows for computations of blocks
of eigenvalues, and a solver for the full time-harmonic system.
A parallel version of this code was implemented based on MPI and CASC's
hypre preconditioning library.
It uses matrices build in parcsr parallel format,
Lobpcg
as a parallel eigensolver and
BoomerAMG, or
ParaSails as a
preconditioner.
In these examples we used piecewise constant for each component of the
solution space and piecewise linears+face bubbles for each component of the
test spaces.
|
|
