17:19:56.2 27-Aug-1999-- >> %% QMG test 117:19:56.3 27-Aug-1999-- >> %% This is an approximation to the mesh on the cover of Claes 17:19:56.3 27-Aug-1999-- >> %% Johnson's book ``Numerical solutions of partial differential17:19:56.3 27-Aug-1999-- >> %% equations by the finite element method'' , Cambridge U. Press 1987.17:19:56.3 27-Aug-1999-- >> %% It has greater refinement near the origin.17:19:56.3 27-Aug-1999-- >> global interactive17:19:56.3 27-Aug-1999-- >> wedgeangle = pi / 6.0;17:19:56.3 27-Aug-1999-- >> %% Make a circular approximation with 4 Bezier curves from pi/6 to 2pi.17:19:56.3 27-Aug-1999-- >> [verts,codes] = gm_circ_approx(wedgeangle, 2*pi, 4);17:19:56.4 27-Aug-1999-- >> newverts = [0,0;verts];17:19:56.4 27-Aug-1999-- >> newcodes = [0;codes];17:19:56.4 27-Aug-1999-- >> domain = gm_cpoly(newverts, newcodes);17:19:56.5 27-Aug-1999-- >> % The size function is of the form h0*r^p+f, where p is determined17:19:56.5 27-Aug-1999-- >> % by (4.31) in Johnson, and f is determined so that the H1 error on17:19:56.5 27-Aug-1999-- >> % the innermost triangles matches the H1 error everywhere else.17:19:56.5 27-Aug-1999-- >> h0 =  0.2;17:19:56.5 27-Aug-1999-- >> gamma = pi / (2*pi - wedgeangle);17:19:56.5 27-Aug-1999-- >> f  = h0 ^ (1.0/gamma);17:19:56.5 27-Aug-1999-- >> p = 1.0 - gamma;17:19:56.5 27-Aug-1999-- >> sizectl = sprintf('(formula (%d * (%%0^2+%%1^2)^(%d/2.0) + %d))', h0, p, f);17:19:56.5 27-Aug-1999-- >> show =  0;17:19:56.5 27-Aug-1999-- >> if length(interactive) > 0  disp(sprintf('size control string = %s', sizectl));  show = 1;end17:19:56.5 27-Aug-1999-- >> mesh = gmmeshgen(domain, 'size', sizectl, 'show', show);17:20:01.4 27-Aug-1999-- >> [scrap, numvtx] = size(mesh{4});17:20:01.4 27-Aug-1999-- >> asp  = gmchecktri(domain,mesh);Maximum aspect ratio =        10.2397 achieved insimplex #359 of topological entity mregion (2:0) which has vertices 13 11 230 Maximum global side length =  0.200806Minimum global altitude =     0.00727931Number of nodes = 259 number of elements = 44717:20:02.1 27-Aug-1999-- >> global aspprod17:20:02.1 27-Aug-1999-- >> global meshsizesum17:20:02.1 27-Aug-1999-- >> if length(aspprod) > 0  aspprod = aspprod * asp;  meshsizesum = meshsizesum + numvtx;end17:20:02.1 27-Aug-1999-- >> if length(interactive) > 0  gmviz(mesh)end17:20:02.1 27-Aug-1999-- >> % ------------------------------------------------------------------17:20:02.1 27-Aug-1999-- >> % Copyright (c) 1999 by Cornell University.  All rights reserved17:20:02.1 27-Aug-1999-- >> % See the accompanying file 'Copyright' for authorship information,17:20:02.1 27-Aug-1999-- >> % the terms of the license governing this software, and disclaimers17:20:02.2 27-Aug-1999-- >> % concerning this software.17:20:02.2 27-Aug-1999-- >> % ------------------------------------------------------------------17:20:02.2 27-Aug-1999-- >> % This file is part of the QMG software.  17:20:02.2 27-Aug-1999-- >> % Version 2.0 of QMG, release date RELDATE.17:20:02.2 27-Aug-1999-- >> % ------------------------------------------------------------------17:20:02.2 27-Aug-1999-- >>