mesh_sample.eps

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

%!PS-Adobe-3.0 EPSF-3.0
%%Title: mesh_sample.eps
%%Creator: Micrografx Graphics Engine
%%CreationDate: Tue Sep 07 05:12:00 2004
%%BoundingBox: 26 0 347 299 
%LanguageLevel: 1 
%%Pages: 1 
%%DocumentProcessColors: Black
%%DocumentSuppliedResources: font Arial
%%EndComments
%%BeginProlog
%%BeginResource: procset RangerDict 1 1
/RangerDict 150 dict def RangerDict begin /LL /languagelevel where {pop 
languagelevel} {1} ifelse def /isL2? LL 1 gt def /isL3? LL 2 gt def /bd{bind 
def}bind def /ld{load def}bd /p/closepath ld /cc/concat ld /v/curveto ld 
/gr/grestore ld /gs/gsave ld /l/lineto ld /m/moveto ld /n/newpath ld /r/rlineto 
ld /sc/scale ld /sd/setdash ld /g/setgray ld /lc/setlinecap ld /lj/setlinejoin 
ld /lw/setlinewidth ld /c/setrgbcolor ld /sp/showpage ld /s/stroke ld 
/tf/transform ld /tr/translate ld /xd{exch def}bd /MGX_Overprint false def 
/min{2 copy gt{exch} if pop}bd /arct where {pop} {/arct{arcto pop pop pop 
pop}bd} ifelse /make_ro_array { /packedarray where { pop packedarray }{array 
astore readonly } ifelse } bd /findcmykcustomcolor where { pop }{ 
/findcmykcustomcolor { 5 make_ro_array } def } ifelse /setcmykcolor where {pop 
/ck /setcmykcolor ld} {/ck {4 1 roll 3 {3 index add 1 min 1 exch sub 3 1 roll} 
repeat setrgbcolor pop} bd} ifelse /setcustomcolor where { pop } { 
/setcustomcolor { exch aload pop pop 4 { 4 index mul 4 1 roll } repeat 5 -1 
roll pop ck } def } ifelse /x {/tint xd findcmykcustomcolor MGX_Overprint 
setoverprint tint 1 exch sub setcustomcolor} def /setoverprint where 
{pop}{/setoverprint {/MGX_Overprint xd} def} ifelse /currentoverprint where 
{pop}{/currentoverprint{MGX_Overprint}def} ifelse /op{dup /MGX_Overprint xd 
setoverprint}def /_mx matrix def /smx{/_mx _mx currentmatrix def}bd /rmx{_mx 
setmatrix}bd /_ctm matrix def /sctm{/_ctm _ctm currentmatrix def /_res 72 72 
matrix defaultmatrix dtransform _ctm idtransform pop abs def}bd /rctm{_ctm 
setmatrix}bd /el{smx tr sc 1 eq {0 -1 m} if 0 0 1 -90 270 arc rmx}bd /ea{smx tr 
sc 1 eq {arcn} {arc} ifelse 1 eq {0 0 l} if rmx}bd /rt{1 eq {m} {l} ifelse 2 
copy 0 r 0 exch r neg 0 r neg 0 exch r}bd /rr{smx tr /_wy exch 3 index div def 
/_wx exch 3 index div def sc _wx 2 div 0 3 -1 roll 1 eq {m} {l} ifelse 0 0 0 
_wy 0.5 arct 0 _wy _wx _wy 0.5 arct _wx _wy _wx 0 0.5 arct _wx 0 0 0 0.5 arct 
_wx 2 div 0 l rmx}bd /clp{/_ct xd {{gs _ct 1 eq {eoclip} {clip} ifelse}stopped 
{gr currentflat 1 add setflat} {exit} ifelse} loop}bd /horz{_y1 _w 2 div add _w 
_y2 {dup _x1 exch m _x2 exch l s} for}bd /vert{_x1 _w 2 div add _w _x2 {dup _y1 
m _y2 l s} for}bd /fdiag{_x1 _w _x2 _y2 _y1 sub add{_y2 m _y2 _y1 sub neg dup r 
s} for}bd /bdiag{_x1 _y2 _y1 sub sub _w _x2{_y2 m _y2 _y1 sub dup neg r s} 
for}bd /hf{/_ht xd clp rctm /_w _res 12 div def pathbbox _w div cvi 1 add _w 
mul /_y2 xd _w div cvi 1 add _w mul /_x2 xd _w div cvi _w mul /_y1 xd _w div 
cvi _w mul /_x1 xd 1 lw [] 0 sd n _ht 0 eq {horz} if _ht 1 eq {vert} if _ht 2 
eq {fdiag} if _ht 3 eq {bdiag} if _ht 4 eq {horz vert} if _ht 5 eq {fdiag 
bdiag} if gr}bd /bf{/_im xd clp rctm /_w _res 12 div round def pathbbox _w div 
cvi 1 add _w mul /_y2 xd _w div cvi 1 add _w mul /_x2 xd _w div cvi _w mul /_y1 
xd _w div cvi _w mul /_x1 xd /_wx _x2 _x1 sub def /_wy _y2 _y1 sub def /_swx 
_wx _w div cvi 8 mul def /_swy _wy _w div cvi 8 mul def /_ic 0 def _x1 _y1 tr 
_wx _wy sc _swx _swy false [_swx 0 0 _swy 0 0] /_swx _swx 8 div def {_im _ic 
_swx div cvi 8 mod 1 getinterval /_ic _ic 1 add def} imagemask gr}bd /accbds 
{dup _uy gt {dup /_uy xd} if dup _ly lt {/_ly xd} {pop} ifelse  dup _ux gt {dup 
/_ux xd} if dup _lx lt {/_lx xd} {pop} ifelse} bd /pathbounds {/_uy -65535 def 
/_ux -65535 def /_ly 65535 def /_lx 65535 def {accbds} {accbds} {} {} 
pathforall} bd /gradrgb { _c 0 get _c 3 get 1 index sub 2 index mul add _c 1 
get _c 4 get 1 index sub 3 index mul add _c 2 get _c 5 get 1 index sub 4 index 
mul add setrgbcolor pop} def /gradcmyk { _c 0 get _c 4 get 1 index sub 2 index 
mul add _c 1 get _c 5 get 1 index sub 3 index mul add _c 2 get _c 6 get 1 index 
sub 4 index mul add _c 3 get _c 7 get 1 index sub 5 index mul add setcmykcolor 
pop} def /lincomp {_s0 sub abs _yp 50 ge {_s0} {_ns _s0 sub} ifelse div} bd 
/lingrad {/_h _uy _ly sub def /_w _ux _lx sub def /_ns _h abs _rc div .5 add 
cvi dup 2 lt {pop 2} if dup 128 gt {pop 128} if 1 sub def /_s0 _ns _yp mul 100 
div .5 add cvi dup 0 lt {pop 0} if def 0 1 _ns {dup lincomp gradcolor _ns 1 add 
div _h mul _lx exch _ly add _w _h 2 mul _ns 1 add div rectfill} for} bd 
/getmaxr {dup mul exch dup mul add sqrt dup _r gt {/_r xd} {pop} ifelse} bd 
/getmaxd {abs exch abs 2 copy lt {exch} if pop dup _r gt {/_r xd} {pop} ifelse} 
bd/squradgrad {_ux _lx sub _xp mul 100 div _lx add _uy _ly sub _yp mul 100 div 
_ly add translate _a neg rotate /_r 0 def {_smx} {_smx} {} {} pathforall 
newpath /_ns _r abs _rc div .5 add cvi dup 2 lt {pop 2} if dup 128 gt {pop 128} 
if def _ns -1 0 {dup _ns div gradcolor .5 add _r mul _ns div _srf} for } bd 
/congrad {_ux _lx sub _xp mul 100 div _lx add _uy _ly sub _yp mul 100 div _ly 
add translate _a rotate /_r 0 def {getmaxr} {getmaxr} {} {} pathforall newpath 
/_ns 60 def 0 1 _ns {dup _ns div gradcolor 3 mul 1.5 sub dup 3 add 0 0 moveto 0 
0 _r 4 index 4 index arc fill 0 0 moveto 0 0 _r 360 5 index sub 360 5 index sub 
arcn fill pop pop} for } bd /gradbase {/_a xd /_or xd /_yp xd /_xp xd /_ct xd 
/gradcolor exch 1 eq {{gradcmyk}} {{gradrgb}} ifelse bd /_c xd /_rc _res 16 div 
dup 18 lt {pop 18} if def _ct clp flattenpath} bd /grad2linear {gradbase _a _or 
add neg rotate pathbounds newpath lingrad gr} bd /grad2square {gradbase _or neg 
rotate pathbounds /_smx {getmaxd} bd /_srf {dup neg dup 2 index 2 mul dup 
rectfill pop} bd squradgrad gr} bd /grad2radial {gradbase _or neg rotate 
pathbounds /_smx {getmaxr} bd /_srf {0 0 2 index 0 360 arc fill pop} bd 
squradgrad gr} bd /gradconical {gradbase pathbounds congrad gr} bd /_lnclrfunc 
{/_c0 xd /_c1 xd 10 dict begin /Decode [_c0 _c1] def /FunctionType 0 def 
/Domain [0 1] def /Range [0 1] def /Size [2] def /BitsPerSample 8 def 
/DataSource <ff00> def currentdict end } bd /_sqclrfunc {/_c0 xd /_c1 xd 10 
dict begin /FunctionType 0 def /Domain [0 1 0 1] def /Range [0 1] def /Decode 
[_c0 _c1] def /Size [17 17] def /BitsPerSample 4 def /DataSource 
<0000000000000000002222222222222220 0244444444444442002466666666666420 
0246888888888642002468999999986420 024689bbbbb986420024689bdddb986420 
024689bdfdb986420024689bdddb986420 024689bbbbb98642002468999999986420 
0246888888888642002466666666666420 0244444444444442002222222222222220 
000000000000000000 > def currentdict end } bd /grad3rgb {_c 0 get _c 3 get 
_clrfunc _c 1 get _c 4 get _clrfunc _c 2 get _c 5 get _clrfunc} bd /grad3cmyk 
{_c 0 get _c 4 get _clrfunc _c 1 get _c 5 get _clrfunc _c 2 get _c 6 get 
_clrfunc _c 3 get _c 7 get _clrfunc} bd /grad3base {/_a xd /_or xd 100 div /_yp 
xd 100 div /_xp xd /_ct xd dup /_isc xd /grad3color exch 1 eq {{grad3cmyk}} 
{{grad3rgb}} ifelse bd /_c xd _ct clp flattenpath} bd /lin3pgrad {/_y0 xd /_y1 
xd 8 dict begin /ShadingType 2 def /Domain [0 1] def /ColorSpace [_isc 1 eq 
{/DeviceCMYK} {/DeviceRGB} ifelse] def /BBox [_lx _ly _ux _uy] def /Coords [_lx 
_y0 _lx _y1] def /Function [grad3color] def currentdict end shfill} bd 
/lin3grad {_yp 0 le {_ly _uy lin3pgrad} {_yp 1 ge {_uy _ly lin3pgrad} {/_h _uy 
_ly sub def /_s0 _h _yp .5 ge {_yp} {1 _yp sub} ifelse mul def /_ym _ly _h _yp 
mul add def _ym _s0 sub _ym lin3pgrad _ym _s0 add _ym lin3pgrad } ifelse } 
ifelse } bd /grad3linear {/_clrfunc {_lnclrfunc} bd grad3base _a _or add neg 
rotate pathbounds newpath lin3grad gr} bd /rad3pgrad {8 dict begin /ShadingType 
3 def /Domain [0 1] def /ColorSpace [_isc 1 eq {/DeviceCMYK} {/DeviceRGB} 
ifelse] def /BBox [_r neg dup _r dup] def /Coords [0 0 0 0 0 _r] def /Function 
[grad3color] def currentdict end shfill} bd /rad3grad {_ux _lx sub _xp mul _lx 
add _uy _ly sub _yp mul _ly add translate _a neg rotate /_r 0 def {getmaxr} 
{getmaxr} {} {} pathforall newpath rad3pgrad} bd /grad3radial {/_clrfunc 
{_lnclrfunc} bd grad3base _or neg rotate pathbounds rad3grad gr} bd/sq3pgrad {8 
dict begin /ShadingType 1 def /Domain [0 1 0 1] def /ColorSpace [_isc 1 eq 
{/DeviceCMYK} {/DeviceRGB} ifelse] def /BBox [_r neg dup _r dup] def /Matrix 
[_r dup add 0 0 _r dup add _r neg dup] def /Function [grad3color] def 
currentdict end shfill} bd /sq3grad {_ux _lx sub _xp mul _lx add _uy _ly sub 
_yp mul _ly add translate _a neg rotate /_r 0 def {getmaxd} {getmaxd} {} {} 
pathforall newpath sq3pgrad} bd /grad3square {/_clrfunc {_sqclrfunc} bd 
grad3base _or neg rotate pathbounds sq3grad gr} bd /gradlinear isL3? 
{{grad3linear}} {{grad2linear}} ifelse bd /gradradial isL3? {{grad3radial}} 
{{grad2radial}} ifelse bd /gradsquare isL3? {{grad3square}} {{grad2square}} 
ifelse bd userdict /RangerImage 40 dict dup begin put /set_rdfile {pop /_rdfile 
currentfile def} bd /readimg{_rdfile exch readstring pop}bd /readimgx{_rdfile 
exch readhexstring pop}bd /i{gs cc /_imgstr exch string def [4 index 0 0 6 
index 0 0] exch dup set_rdfile 1 eq {{_imgstr readimg}} {{_imgstr readimgx}} 
ifelse image gr}bd /im{gs cc /_imgstr exch string def 1 eq {true} {false} 
ifelse 0 set_rdfile [3 index 0 0 5 index 0 0] {_imgstr readimgx} imagemask 
gr}bd /_image systemdict /image get def /_setgray systemdict /setgray get def 
/_currentgray systemdict /currentgray get def /_settransfer systemdict 
/settransfer get def /_currenttransfer systemdict /currenttransfer get def 
/updatetransfer {_currenttransfer band 1 ge band 4 le and {{1 exch sub /dummy 
exec} dup 3 _currenttransfer put _settransfer} if band 6 eq{negative {{pop 
0}}{{pop 1}}ifelse _settransfer} if } def/restoretransfer {_settransfer} def 
/clr? /deviceinfo where {pop deviceinfo /Colors known {deviceinfo /Colors get 
exec 1 gt}{false}ifelse }{ /statusdict where { pop statusdict /processcolors 
known { statusdict /processcolors get exec 1 gt }{ false }ifelse }{ false 
}ifelse }ifelse def /blank 0 _currenttransfer exec 1 _currenttransfer exec eq 
def /negative blank{ 0 _currenttransfer exec 0.5 lt }{ 0 _currenttransfer exec 
1 _currenttransfer exec gt }ifelse def /inverted? negative def /l2? 
/languagelevel where{pop languagelevel 2 ge}{false} ifelse def /4= { 4 index eq 
8 1 roll 4 index eq 8 1 roll 4 index eq 8 1 roll 4 index eq 8 1 roll pop pop 
pop pop and and and }def clr? {/band 0 def}{/band 5 def}ifelse /setcmykcolor 
where { pop 1 0 0 0 setcmykcolor _currentgray 1 exch sub 0 1 0 0 setcmykcolor 
_currentgray 1 exch sub 0 0 1 0 setcmykcolor _currentgray 1 exch sub 0 0 0 1 
setcmykcolor _currentgray 1 exch sub 4 {4 copy} repeat 1 0 0 0 4= {/band 1 
store} if 0 1 0 0 4= {/band 2 store} if 0 0 1 0 4= {/band 3 store} if 0 0 0 1 
4= {/band 4 store} if 0 0 0 0 4= {/band 6 store} if } if blank {/band 6 store} 
if l2? {band 0 eq {/DeviceCMYK}{/DeviceGray} ifelse setcolorspace } if 
/sepimage l2? { /image ld }{ { begin Width Height BitsPerComponent ImageMatrix 
Decode length 2 eq {/DataSource load image} if Decode length 6 eq {DataSource 0 
get DataSource 1 get DataSource 2 get true 3 colorimage }if Decode length 8 eq 
{ DataSource 0 get DataSource 1 get DataSource 2 get DataSource 3 get true 4 
colorimage }if end }def }ifelse /_image2 l2? { /_image ld }{ { begin Width 
Height BitsPerComponent ImageMatrix /DataSource load _image end } def }ifelse 
/doimage { band 0 eq band 5 eq or { sepimage } if band 1 ge band 4 le and { 
_image2 } if band 6 eq { _image2 }if } def /ic { gsave cc /scanlen xd /ncomp xd 
/IsBin xd /bpcomp xd /rows xd /cols xd 0 set_rdfile /picstr1 cols string def 
/picstr2 cols string def /picstr3 cols string def /picstr4 cols string def 
/picstr5 cols string def l2? {band 0 eq {/DeviceCMYK}{/DeviceGray} ifelse 
setcolorspace} if 12 dict begin /ImageType 1 def /Width cols def /Height rows 
def /ImageMatrix [cols 0 0 rows 0 0] def /BitsPerComponent bpcomp def band 0 eq 
{ /Decode [0 1 0 1 0 1 0 1] def } { /Decode [0 1] def } ifelse band 0 eq { 
/MultipleDataSources true def /DataSource [ { picstr1 readimgx } { picstr2 
readimgx } { picstr3 readimgx } { picstr4 readimgx picstr5 readimgx pop } ] def 
} if band 1 eq { /DataSource { picstr1 readimgx picstr2 readimgx pop picstr3 
readimgx pop picstr4 readimgx pop picstr5 readimgx pop } def } if band 2 eq { 
/DataSource { picstr1 readimgx pop picstr2 readimgx picstr3 readimgx pop 
picstr4 readimgx pop picstr5 readimgx pop } def } if band 3 eq { /DataSource { 
picstr1 readimgx pop picstr2 readimgx pop picstr3 readimgx picstr4 readimgx pop 
picstr5 readimgx pop } def } if band 4 eq { /DataSource { picstr1 readimgx pop 
picstr2 readimgx pop picstr3 readimgx pop picstr4 readimgx picstr5 readimgx pop 
} def } if band 5 ge { /DataSource { picstr1 readimgx pop picstr2 readimgx pop 
picstr3 readimgx pop picstr4 readimgx pop picstr5 readimgx } def } if