2.7 Assembling

For a problem described by a system of PDEs a local matrix for each element of the discretization
, *nucleus*, is constructed.
This nucleus is integrated into the *general* matrix of the system through a process called *assembling*.
In this section the assembling algorithm applied in numerical schemes of problems discussed in Chapters 3 and 4 is derived.

The Jacobi matrix needed for the Newton method in the case of a finite element discretization is,

In the further text we will call this Jacobi matrix

where , , and the discrete operators are calculated as integrals over the discretization .

The calculation of the general matrix by the inner product over the whole discretisized domain is rarely utilized in computer programs. A common approach is to construct the matrix by assembling it out of the nucleus matrix which is calculated for each element [11],

(2.42) |

with the inputs,

where and .

The basic nodal function , defined at the arbitrary point is non-zero only on the patch . Furthermore can be represented as,

at the point and at three other points of the tethraedal element .

The discrete operators obtained by testing of -th equation of the system with the basic nodal function is,

The partial derivative of (2.45) with respect to is,

Obviously in (2.46) is equal to one of the basic nodal functions at the element . Furthermore, the partial derivative is non-zero only if the stays for one nodal value of .

We now define the operator , , , which assigns a single global index to every local index of vertex belonging to the tethraedal element . The inverse function is also well-defined.

From (2.41), (2.43), and (2.46) we have,

The assembling process consists of calculating of nucleus matrix for each element of the mesh and building the inputs of the global matrix as the sum on the right side of the equation (2.47). For the implementation of the assembling algorithm into software tools the following algorithm is used,

At the end of the assembling process the general matrix contains the values given by (2.40).

H. Ceric: Numerical Techniques in Modern TCAD