https://github.com/team-pancho/HDG3D
Tip revision: af0d6cb24c63ad6fa07331f26a5733116db800b7 authored by shukaidu on 29 April 2020, 21:46:21 UTC
Merge pull request #2 from team-pancho/hdg3d_skdu
Merge pull request #2 from team-pancho/hdg3d_skdu
Tip revision: af0d6cb
checkQuadrature3d.m
function formulas=checkQuadrature3d(k,constantmass)
%formulas=checkQuadrature3d(k,constantmass)
%Imput:
% k: degree of polynomials
% constantmass: 1 constant coefficients in mass matrices
% 0 non-costant
%Output:
% formulas{1} : 3d quadrature formula for errors and var coeff
% formulas{2} : 3d quadrature formula for constant coefficients
% formulas{3} : 2d quadrature formula
% formulas{4} : 2d quadrature formula for errors
%Last update: March 14, 2013
% degrees = {3k,2k,2k,2k+2} if constant mass = 0
% {2k,2k,2k,2k+2} if constant mass = 1
% for k=0, take {2,0,0,2}
TablesQuadForm3d
TablesQuadForm
switch constantmass
case 1
switch k
case 0
formulas={tetra5,tetra1,matrix0,matrix4};
case 1
formulas={tetra3,tetra3,matrix4,matrix4};
case 2
formulas={tetra7,tetra7,matrix9,matrix9};
case 3
formulas={tetra9,tetra9,matrix11,matrix11};
end
case 0
switch k
case 0
formulas={tetra5,tetra1,matrix0,matrix4};
case 1
formulas={tetra5,tetra3,matrix4,matrix4};
case 2
formulas={tetra7,tetra7,matrix9,matrix9};
case 3
formulas={tetra9,tetra9,matrix11,matrix11};
end
end
formulas{3}(:,4)=formulas{3}(:,4)/2;
formulas{4}(:,4)=formulas{4}(:,4)/2;
return