Next: Package FL.MULTIGRID, Previous: Package FL.NAVIER-STOKES-FE, Up: Reference manual
The FL.ITERATION package includes the
definition for the abstract classes <solver>,
<iterative-solver>, as well as the generic functions
iterate and solve which constitutes the interface for
linear and non-linear solving. Both functions work on a blackboard which
is passed together with the iteration used as argument.
Several instances of iterative solvers are implemented, e.g. Gauss-Seidel,
SOR, ILU (in linit.lisp) and CG (in krylow.lisp). A larger
block of code is contained in a separate package FL.MULTIGRID and
contains the multigrid iteration. From this class, an algebraic multigrid
iteration is derived in amg.lisp and a geometric multigrid iteration
in a separate file geomg.lisp and package FL.GEOMG.
Observe time during an iteration. This should be used as element in the observe list of an iteration.
Preconditioned Bi-CGStab iteration
Superclasses: <LINEAR-ITERATION>
Direct slots:
- PRECONDITIONER: NIL
Subspace correction scheme generated by collecting overlapping blocks of unknowns.
Superclasses: <LINEAR-ITERATION>
Direct slots:
- INNER-ITERATION: Iteration which is used to solve for each block.
- ORDERING: NIL
- STORE-P: T if diagonal can be stored.
Preconditioned conjugate gradient iteration
Superclasses: <LINEAR-ITERATION>
Direct slots:
- PRECONDITIONER: NIL
PSC with custom BLOCK-SETUP function slot.
Superclasses: <SETUP-BLOCKS-MIXIN> <PSC>
SSC with custom BLOCK-SETUP function slot.
Superclasses: <SETUP-BLOCKS-MIXIN> <SSC>
The <function> class is an abstract class for a general function. This function will usually accept vector arguments, the dimensions of domain and image are fixed when defining the function. If the function is differentiable, the gradient matrix can be obtained by evaluating the gradient slot.
Direct slots:
- DOMAIN-DIMENSION: NIL
- IMAGE-DIMENSION: NIL
Incomplete LU iteration. omega is the modification parameter, eta is the diagonal enhancement.
Superclasses: <LINEAR-ITERATION>
Direct slots:
- OMEGA: NIL
- ETA: NIL
- ORDERING: NIL
Direct slots:
- OBSERVE: The initform depends on the subclass.
- OUTPUT: A boolean indicating if output is to be done.
- SUCCESS-IF: A form specifying a success criterion.
- FAILURE-IF: A form specifying a failure criterion.
- START-TIME: Start time of the iteration.
- SUCCESS-IF-FN: Compiled success-if form.
- FAILURE-IF-FN: Compiled failure-if form.
Base class of all iterative solvers and solution strategies.
Superclasses: <ITERATION> <SOLVER>
An <iterator> object contains functions doing iteration work or flags indicating which work has or has not to be done for calling that iterator. It is generated by the generic function make-iterator.
Direct slots:
- MATRIX: The matrix on which the iterator is working.
- INITIALIZE: NIL or an initialization function which is called with the matrix as argument.
- ITERATE: A function of the arguments solution, rhs, and residual which performs an update step.
- RESIDUAL-BEFORE: T if the residual has to be up-to-date before the iteration step.
- RESIDUAL-AFTER: T if the residual is updated through the iteration step.
The <linear-iteration> class. Linear iterations are e.g. <gauss-seidel> or <multigrid>.
Direct slots:
- DAMP: NIL
Class for linear iterative solvers.
Superclasses: <DISCRETE-ITERATIVE-SOLVER>
Direct slots:
- ITERATION: The inner iteration.
A linear iteration interface for the LU exact solver.
Superclasses: <LINEAR-ITERATION>
Direct slots:
- STORE-P: Store decomposition for multiple applications.
One step of this iteration performs nr-steps of the base iteration.
Superclasses: <LINEAR-ITERATION>
Direct slots:
- BASE: Base iteration.
- NR-STEPS: Number of steps with which the base iteration is performed.
If failure occurs, object of this class try an alternative iteration. Usually this will be a direct decomposition.
Superclasses: <LINEAR-SOLVER>
Direct slots:
- FALLBACK: The fallback iteration.
If you happen to have a problem-adapted solver given as a function, you may use this base class.
Superclasses: <ITERATIVE-SOLVER>
Direct slots:
- SOLVER-FUNCTION: NIL
Successive subspace correction scheme.
Superclasses: <BLOCK-ITERATION>
Direct slots:
- OMEGA: NIL
Generic evaluation of functions on an argument. Numbers and arrays are treated as constants. Special evaluation is defined for multivariate polynomials on vectors and for <function> objects.
Generic evaluation of gradients of differentiable functions.
Iterates on the data in the blackboard according to the iteration iter.
Constructs an iterator object given a linear iteration and a matrix.
Returns an iterator which does several steps of the given iterator.
Setup routine for determining the blocking of unknowns. Returns a list of blocks where each block is a vector of keys. May return a second value which is a list of pair. Each pair is of the form start-index/end-index and can be used to filter out different fe components.