# How to Improve Convergence for the MLFMM /FEM

The hybrid multilevel fast multipole method (MLFMM) / finite element method (FEM) is an iterative solution method and under certain conditions, the iterative solution may fail to converge. Several model or solution settings are presented that could improve the model's convergence behaviour.

ERROR 4673: Iterative solution of the system of linear equations failed, maybe try another pre-conditioner (solution settings).

WARNING 830: Maximum number of iterations reached without convergence, using in the following the solution with the smallest residuum.

- Adjust the mesh.
- Change the default box size.
- Use the CFIE for metallic structures.
- Use double precision.
- Change the FEM to use first order basis functions.

## Adjusting the Mesh

Slight adjustments to the mesh size (smaller or larger elements) could lead to improved convergence. If a model is discretised too finely or too coarsely, convergence could be negatively affected.

## Changing the Default MLFMM Box Size

The MLFMM uses a boxing algorithm that encloses the entire computational space in a single box at the highest level, dividing this box in three dimensions into a maximum of eight child boxes and repeating the process iteratively until the side length of each child box is approximately a quarter wavelength at the finest level. Using a different box size at the finest level can sometimes facilitate convergence, although memory consumption could increase if the box size is increased.

The default box size is 0.23. A lower value decreases memory consumption while a higher value increases the memory consumption.

## Using the CFIE on MLFMM Metallic Surfaces

- Use the default preconditioner with the CFIE.
- CFIE can only be applied to surfaces bounding closed PEC structures.
- A mixture of CFIE and EFIE surfaces can be used.
- Sharp corners on CFIE surfaces can lead to inaccurate results.

Sharp corners should be meshed finer if CFIE is applied. If there is uncertainty, the EFIE should rather be applied around sharp corners. A rule of thumb is to apply the EFIE up to a few meshed triangles away from the sharp corners.

## Using Double Precision

The Solver uses single precision by default- a single byte is used to store a complex number.

- Double precision requires twice the memory compared to single precision.
- Double precision does not improve convergence for the stabilised MLFMM.

## Changing the FEM to Use First Order Basis Functions

The FEM uses higher order (order two) basis functions by default. For large volumes, the higher order results in a much smaller number of tetrahedra in the mesh.