ElemMassWMatP1.m

Go to the documentation of this file.

function AElem=ElemMassWMatP1(area,w)
% function AElem=ElemMassWMatP1(area,w)
%   Computation of the element weight mass matrix for a triangle with area 'area'.
%   
% Parameters:
%  area: triangle area
%  w: weight function values at the triangle vertices 
%  
% Return values:
%  AElem: `3\times 3` matrix
%
  AElem=(area/30)*[ 3*w(1)+w(2)+w(3), w(1)+w(2)+w(3)/2, w(1)+w(2)/2+w(3); ...
                    w(1)+w(2)+w(3)/2, w(1)+3*w(2)+w(3), w(1)/2+w(2)+w(3); ...
                    w(1)+w(2)/2+w(3), w(1)/2+w(2)+w(3), w(1)+w(2)+3*w(3)];