Assembling Mass Weight Matrix by $P_1$ -Lagrange finite elements using "OptV0" version (see report). More...

Functions
function M =	MassWAssemblingP1OptV0 (nq, nme, me, areas, Tw)
	Assembling Mass Weight Matrix by -Lagrange finite elements using "OptV0" version (see report).

Detailed Description

Assembling Mass Weight Matrix by $P_1$ -Lagrange finite elements using "OptV0" version (see report).

Function Documentation

function M = MassWAssemblingP1OptV0	(	nq,
		nme,
		me,
		areas,
		Tw
	)

Assembling Mass Weight Matrix by $P_1$ -Lagrange finite elements using "OptV0" version (see report).

The Mass Weight Matrix is given by

$\MasseF{w}_{i,j}=\int_\DOMH w(\q)\FoncBase_i(\q) \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);
    w=@(x,y) cos(x+y);
    Tw=w(Th.q(1,:),Th.q(2,:));
    Mw=MassWAssemblingP1OptV0(Th.nq,Th.nme,Th.me,Th.areas,Tw);

Parameters

nq	total number of nodes of the mesh, also noted $\nq$ ,
nme	number of triangles, also noted $\nme$ ,
me	$3\times\nme$ `int32` array, $\me(\jl,k)$ index of storage, in the array , of the $\jl$ -th vertex of the triangle of index , $\jl\in\{1,2,3\}$ and $k\in\{1,\hdots,\nme\}.$ Also noted $\me$ .
areas	$1\times\nme$ array, areas(k) is the area of triangle k.
Tw	$\nme$ `double` array, $Tw(i)=w(\q^i),$ $\forall i\in\ENS{1}{\nq}$ .

Return values

M	$\nq\times\nq$ sparse matrix

Definition at line 17 of file MassWAssemblingP1OptV0.m.