bijiao.m

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

function u=bijiao
     u=sparse(5,2);
     for i=1:5 
        mesh=getmesh([0 0;1 0;1 1;0 1],[1 2 3;3 4 1],[1 2;2 3;3 4;4 1],[]);
        mesh=uniformrefine(mesh,i);
        
        A=MCDFM(mesh,0.0005,@a1,@a2,@b,@f,@g_D);
        mesh=getmesh([0 0;1 0;1 1;0 1],[1 2 3;3 4 1],[1 2;2 3;3 4;4 1],[]);
        
        mesh=uniformrefine(mesh,i);
        B=SUPG(mesh,0.0005,@a1,@a2,@b,@f,@g_D);
        mesh=getmesh([0 0;1 0;1 1;0 1],[1 2 3;3 4 1],[1 2;2 3;3 4;4 1],[]);
        mesh=uniformrefine(mesh,i);
        
        A=WMCDFM(mesh,0.0005,@a1,@a2,@b,@f,@g_D);
        u(i,1)=condest(A);
        u(i,2)=condest(B);
     end