LocalMass2D_LFE_Quad.m
79cd46f2bb0bmaster

LocalMass2D_LFE_Quad.m
View Options

			function Mloc = LocalMass2D_LFE_Quad(mapRef,x,shapeF,QuadRule,FHandle)

			Mloc = zeros(3,3);

			Fval = FHandle(x);

			for i =1:3
			for j = i:3
			Mloc (i,j) = mapRef.detJK * sum(QuadRule.w .* Fval .* shapeF(:,i) .* shapeF(:,j));
			end
			end

			% symetry
			Mloc (3 ,1) = Mloc (1 ,3);
			Mloc (3 ,2) = Mloc (2 ,3);
			Mloc (2 ,1) = Mloc (1 ,2);
			end