OptFEM2D  0.1
Matlab optimized FEM2D
 All Files Functions Groups Pages
ComputeArea.m
Go to the documentation of this file.
1 %****f* FEM/ComputeArea
2 % NAME
4 % USE
5 % area=ComputeArea(q,me)
6 % DESCRIPTION
7 % Computation of areas of triangles from the mesh [q,me]
8 % INPUTS
9 %
10 % OUTPUT
11 %
12 % RETURN VALUE
13 %
14 % EXAMPLE
15 %
16 % CREATION DATE
17 % January 14th 2004
18 % MODIFICATION HISTORY
19 % None
20 % NOTES
21 % Unoptimized
22 % AUTHOR
23 % Cuvelier F., cuvelier@math.univ-paris13.fr
24 % SOURCE
25 %
26 function area=ComputeArea(q,me)
27 nq=size(q,2);
28 nme=size(me,2);
29 area=zeros(1,nme);
30 for k=1:nme
31  area(k).5*abs(det([ q(1,me(1,k)) q(2,me(1,k)) 1 ; ...
32  q(1,me(2,k)) q(2,me(2,k)) 1 ; ...
33  q(1,me(3,k)) q(2,me(3,k)) 1 ]));
34 end
35 %*****