get_area2.m

hodge decomposition for 2D vector fields.

function area2=get_area2(triangle,mesh) 
%% calculate the 2*triangle_area of a triangle 
%% triangle: the Grid Table 
%% mesh: the Node Table  
%% area2: the resulted 2*triangle_area 
x0=mesh(triangle(1),1); 
y0=mesh(triangle(1),2); 
x1=mesh(triangle(2),1); 
y1=mesh(triangle(2),2); 
x2=mesh(triangle(3),1); 
y2=mesh(triangle(3),2); 
area2=abs((x1-x0)*(y2-y0)-(x2-x0)*(y1-y0)); 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%