Assembly of the Stiffness Matrix using $P_1$ -Lagrange finite elements. More...

Functions
function R =	StiffAssemblingP1OptV2 (nq, nme, q, me, areas)
	Assembly of the Stiffness Matrix using -Lagrange finite elements.

Detailed Description

Assembly of the Stiffness Matrix using $P_1$ -Lagrange finite elements.

Function Documentation

function R = StiffAssemblingP1OptV2	(	nq,
		nme,
		q,
		me,
		areas
	)

Assembly of the Stiffness Matrix using $P_1$ -Lagrange finite elements.

The Stiffness Matrix $\Stiff$ is given by

$\Stiff_{i,j}=\int_\DOMH \DOT{\GRAD\FoncBase_i(\q)}{\GRAD\FoncBase_j(\q)}d\q,\ \forall (i,j)\in{\ENS{1}{\nq}}^2$

where $\FoncBase_i$ are $P_1$ -Lagrange basis functions.

Example

    Th=SquareMesh(10);
    R=StiffAssemblingP1OptV2(Th.nq,Th.nme,Th.me,Th.q,Th.areas);

Parameters

nq	total number of nodes of the mesh, also denoted by $\nq$ .
nme	total number of triangles, also denoted by $\nme$ .
q	Array of vertices coordinates, $2\times\nq$ array. ${\q}(\il,j)$ is the $\il$ -th coordinate of the -th vertex, $\il\in\{1,2\}$ and $j\in\ENS{1}{\nq}$
me	Connectivity array, $3\times\nme$ array. $\me(\jl,k)$ is the storage index of the $\jl$ -th vertex of the -th triangle in the array $\q$ of vertices coordinates, $\jl\in\{1,2,3\}$ and $k\in{\ENS{1}{\nme}}$ .
areas	Array of areas, $1\times\nme$ array. areas(k) is the area of the -th triangle.

Return values

R	Stiffness Matrix, $\nq\times\nq$ sparse matrix.

Definition at line 17 of file StiffAssemblingP1OptV2.m.