naginterfaces.library.mesh.dim2_​sparsity

naginterfaces.library.mesh.dim2_sparsity(nv, nnzmax, conn)

dim2_sparsity generates the sparsity pattern of a finite element matrix associated with a given mesh.

For full information please refer to the NAG Library document for d06cb

https://support.nag.com/numeric/nl/nagdoc_30/flhtml/d06/d06cbf.html

Parameters

nvint

The total number of vertices in the input mesh.

nnzmaxint

The maximum number of nonzero entries in the matrix based on the input mesh. It is the dimension of the arrays $irow$ and $icol$ as declared in the function from which dim2_sparsity is called.

connint, array-like, shape $(3,\text{nelt})$

The connectivity of the mesh between triangles and vertices. For each triangle $\text{j}$, $conn[\text{i}-1,\text{j}-1]$ gives the indices of its three vertices (in anticlockwise order), for $\text{j}=1,2,\dots ,\text{nelt}$, for $\text{i}=1,2,\dots ,3$.

Returns

nnzint

The number of nonzero entries in the matrix associated with the input mesh.

irowint, ndarray, shape $\left(nnzmax\right)$

The first $nnz$ elements contain the row and column indices of the nonzero elements supplied in the finite element matrix $A$.

icolint, ndarray, shape $\left(nnzmax\right)$

The first $nnz$ elements contain the row and column indices of the nonzero elements supplied in the finite element matrix $A$.

Raises

NagValueError

(errno $1$)

On entry, $conn(\text{I},\text{J})=⟨value⟩$, $\text{I}=⟨value⟩$, $\text{J}=⟨value⟩$ and $nv=⟨value⟩$.

Constraint: $conn(\text{I},\text{J})\ge 1$ and $conn(\text{I},\text{J})\le nv$.

(errno $1$)

On entry, $\text{nelt}=⟨value⟩$ and $nv=⟨value⟩$.

Constraint: $\text{nelt}\le 2\times nv-1$.

(errno $1$)

On entry, vertices $2$ and $3$ of the triangle $\text{K}$ have the same index $\text{I}$: $\text{K}=⟨value⟩$ and $\text{I}=⟨value⟩$.

(errno $1$)

On entry, vertices $1$ and $3$ of the triangle $\text{K}$ have the same index $\text{I}$: $\text{K}=⟨value⟩$ and $\text{I}=⟨value⟩$.

(errno $1$)

On entry, vertices $1$ and $2$ of the triangle $\text{K}$ have the same index $\text{I}$: $\text{K}=⟨value⟩$ and $\text{I}=⟨value⟩$.

(errno $1$)

On entry, $nnzmax=⟨value⟩$ and $nv=⟨value⟩$.

Constraint: $nnzmax\le {nv}^2$.

(errno $1$)

On entry, $nnzmax=⟨value⟩$, $\text{nelt}=⟨value⟩$ and $nv=⟨value⟩$.

Constraint: $nnzmax\ge (4\times \text{nelt}+nv)$.

(errno $1$)

On entry, $nv=⟨value⟩$.

Constraint: $nv\ge 3$.

(errno $2$)

A serious error has occurred in an internal call to an auxiliary function. Check the input mesh, especially the connectivity between triangles and vertices (the argument $conn$). Array dimensions should also be checked. If the problem persists, contact NAG.

Notes

dim2_sparsity generates the sparsity pattern of a finite element matrix associated with a given mesh. The sparsity pattern is returned in a coordinate storage format consistent with the sparse linear algebra functions in submodule sparse. More precisely dim2_sparsity returns the number of nonzero elements in the associated sparse matrix, and their row and column indices. This is designed to assist you in applying finite element discretization to meshes from the D06 Introduction and in solving the resulting sparse linear system using functions from submodule sparse.

The output sparsity pattern is based on the fact that finite element matrix $A$ has elements $a_{ij}$ satisfying:

$$a_{ij}\ne 0\Rightarrow i\text{and}j\text{are vertices belonging to the same triangle.}$$