test4.di

来自「算断裂的」· DI 代码 · 共 118 行

17:24:00.2 27-Aug-1999
-- >> %% Test 4: mesh a cylinder with two small spheres removed (like bubbles)

17:24:00.2 27-Aug-1999
-- >> 

17:24:00.2 27-Aug-1999
-- >> cyl1 = gmmake_cyl(4);

17:24:00.5 27-Aug-1999
-- >> global QMG_DATA

17:24:00.5 27-Aug-1999
-- >> s = gm_read(fullfile(QMG_DATA, 'sphere3.geo'));

17:24:01.1 27-Aug-1999
-- >> bubble1 = gmapply([0.01, 0, 0, -0.02; 0, 0.01, 0, 0; 0, 0, 0.01, 0], s);

17:24:01.7 27-Aug-1999
-- >> bubble2 = gmapply([0.01, 0, 0, 0.02; 0, 0.01, 0, 0; 0, 0, 0.01, 0], s);

17:24:01.7 27-Aug-1999
-- >> cyl2 = gmcavity(cyl1, bubble1);

17:24:01.9 27-Aug-1999
-- >> cyl3 = gmcavity(cyl2, bubble2);

17:24:02.1 27-Aug-1999
-- >> show = 0;

17:24:02.1 27-Aug-1999
-- >> global interactive

17:24:02.2 27-Aug-1999
-- >> if length(interactive) > 0
  show = 1;
end

17:24:02.2 27-Aug-1999
-- >> 

17:24:02.2 27-Aug-1999
-- >> mesh = gmmeshgen(cyl3, 'curve', '(const 1.2)', 'show', show);

17:26:49.8 27-Aug-1999
-- >> asp = gmchecktri(cyl3, mesh);
Maximum aspect ratio =        49.0823 achieved in
simplex #2922 of topological entity cylinder (3:0) which has vertices 179 152 209 739 
Maximum global side length =  1.03923
Minimum global altitude =     0.000131217
Number of nodes = 1179 number of elements = 6054

17:26:51.4 27-Aug-1999
-- >> 

17:26:51.4 27-Aug-1999
-- >> global aspprod

17:26:51.4 27-Aug-1999
-- >> global meshsizesum

17:26:51.4 27-Aug-1999
-- >> 

17:26:51.4 27-Aug-1999
-- >> [scrap, numvtx] = size(mesh{4});

17:26:51.4 27-Aug-1999
-- >> 

17:26:51.4 27-Aug-1999
-- >> if length(aspprod) > 0
    aspprod  = aspprod * asp;
    meshsizesum = meshsizesum + numvtx;
end

17:26:51.4 27-Aug-1999
-- >> 

17:26:51.4 27-Aug-1999
-- >> 

17:26:51.4 27-Aug-1999
-- >> 

17:26:51.4 27-Aug-1999
-- >> % ------------------------------------------------------------------

17:26:51.4 27-Aug-1999
-- >> % Copyright (c) 1999 by Cornell University.  All rights reserved

17:26:51.4 27-Aug-1999
-- >> % See the accompanying file 'Copyright' for authorship information,

17:26:51.4 27-Aug-1999
-- >> % the terms of the license governing this software, and disclaimers

17:26:51.4 27-Aug-1999
-- >> % concerning this software.

17:26:51.5 27-Aug-1999
-- >> % ------------------------------------------------------------------

17:26:51.5 27-Aug-1999
-- >> % This file is part of the QMG software.  

17:26:51.5 27-Aug-1999
-- >> % Version 2.0 of QMG, release date RELDATE.

17:26:51.5 27-Aug-1999
-- >> % ------------------------------------------------------------------

17:26:51.5 27-Aug-1999
-- >> 

17:26:51.5 27-Aug-1999
-- >>