adptivesupg.m

来自「C++编译指导」· M 代码 · 共 41 行

function [mesh,eta, res, osc]=ADPTIVESUPG(mesh, eps, a1, a2, b, f, g_D)
   mesh=SUPG(mesh, eps, a1, a2, b, f, g_D);
   [eta, res, osc]= estimate(mesh,eps,a1,a2,b,f,0);
for i=1:7
   mesh = bisection(mesh,eta,res,osc,0.5,0.5);
   mesh=SUPG(mesh, eps, a1, a2, b, f, g_D);
   [eta, res, osc]= estimate(mesh,eps,a1,a2,b,f,0);
end   























% subplot(1,2,1)
% hold off
% trisurf(mesh.elem, mesh.node(:,1), mesh.node(:,2), u', ...
%         'FaceColor', 'interp', 'EdgeColor', 'interp');
% view(40,15);
% title('Solution to the Poisson equation', 'FontSize', 14)

%--------------------------------------------------------------------------
% end of function CDFM
%--------------------------------------------------------------------------