LocalLapl2D_LFE_Quad.m
79cd46f2bb0bmaster

LocalLapl2D_LFE_Quad.m
View Options

			function Aloc = LocalLapl2D_LFE_Quad(mapRef,x,gradshapeF,QuadRule,FHandle)

			Aloc = zeros(3,3);

			Fval = FHandle(x);

			for i =1:3
			for j = i:3
			Aloc (i,j) = sum(QuadRule.w .* Fval .* ...
			sum(gradshapeF(:, 2i-1:2i) .* (gradshapeF(:, 2j-1:2j) * mapRef.TK), 2) );
			end
			end

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