


INPUT function md=cielomesh(md) mesh solution


0001 function md=cielomesh(md) 0002 %INPUT function md=cielomesh(md) 0003 %mesh solution 0004 0005 %timing 0006 t1=clock; 0007 0008 analysis='mesh'; 0009 0010 %configure, and initialise data model and parameter defaults: 0011 SetParameterDefaults; 0012 0013 %first build an initial triangulation, using the domain outline and qhull. 0014 [md.elements,md.x,md.y,md.z]=MeshTriangulate(md); 0015 md.numberofgrids=length(md.x); 0016 md.numberofelements=length(md.elements); 0017 0018 %timing 0019 t2=clock; 0020 disp(sprintf('\n%s\n',[' solution lasted ' num2str(etime(t2,t1)) ' seconds']));