split.ps

弹性波动方程正演程序

%%Page: 3 3
3 2 bop 78 -46 a Fb(Karr)n(enb)n(ach)645 b Fj(3)319 b
Fb(Splitting)19 b(the)f(wave)h(op)n(er)n(ator)78 146
y Fj(is)d(c)o(hosen:)292 220 y Fd(0)292 293 y(B)292 318
y(B)292 343 y(B)292 367 y(B)292 392 y(B)292 417 y(B)292
442 y(B)292 469 y(@)337 205 y(0)337 278 y(B)337 305 y(@)382
243 y Ff(\033)410 250 y Fa(1)382 303 y Ff(\033)410 310
y Fa(2)382 363 y Ff(\033)410 370 y Fa(3)437 205 y Fd(1)437
278 y(C)437 305 y(A)337 384 y(0)337 457 y(B)337 483 y(@)382
421 y Ff(\033)410 428 y Fa(4)382 482 y Ff(\033)410 489
y Fa(5)382 542 y Ff(\033)410 549 y Fa(6)437 384 y Fd(1)437
457 y(C)437 483 y(A)482 220 y(1)482 293 y(C)482 318 y(C)482
343 y(C)482 367 y(C)482 392 y(C)482 417 y(C)482 442 y(C)482
469 y(A)532 392 y Fj(=)584 220 y Fd(0)584 293 y(B)584
318 y(B)584 343 y(B)584 367 y(B)584 392 y(B)584 417 y(B)584
442 y(B)584 469 y(@)629 205 y(0)629 278 y(B)629 305 y(@)673
243 y Ff(c)694 250 y Fa(11)780 243 y Ff(c)801 250 y Fa(12)887
243 y Ff(c)908 250 y Fa(13)673 303 y Ff(c)694 310 y Fa(21)780
303 y Ff(c)801 310 y Fa(22)887 303 y Ff(c)908 310 y Fa(23)673
363 y Ff(c)694 370 y Fa(31)780 363 y Ff(c)801 370 y Fa(32)887
363 y Ff(c)908 370 y Fa(33)954 205 y Fd(1)954 278 y(C)954
305 y(A)1039 205 y(0)1039 278 y(B)1039 305 y(@)1084 243
y Ff(c)1105 250 y Fa(14)1191 243 y Ff(c)1212 250 y Fa(15)1298
243 y Ff(c)1319 250 y Fa(16)1084 303 y Ff(c)1105 310
y Fa(24)1191 303 y Ff(c)1212 310 y Fa(25)1298 303 y Ff(c)1319
310 y Fa(26)1084 363 y Ff(c)1105 370 y Fa(34)1191 363
y Ff(c)1212 370 y Fa(35)1298 363 y Ff(c)1319 370 y Fa(36)1365
205 y Fd(1)1365 278 y(C)1365 305 y(A)629 384 y(0)629
457 y(B)629 483 y(@)673 421 y Ff(c)694 428 y Fa(41)780
421 y Ff(c)801 428 y Fa(42)887 421 y Ff(c)908 428 y Fa(43)673
482 y Ff(c)694 489 y Fa(51)780 482 y Ff(c)801 489 y Fa(52)887
482 y Ff(c)908 489 y Fa(53)673 542 y Ff(c)694 549 y Fa(61)780
542 y Ff(c)801 549 y Fa(62)887 542 y Ff(c)908 549 y Fa(63)954
384 y Fd(1)954 457 y(C)954 483 y(A)1039 384 y(0)1039
457 y(B)1039 483 y(@)1084 421 y Ff(c)1105 428 y Fa(44)1191
421 y Ff(c)1212 428 y Fa(45)1298 421 y Ff(c)1319 428
y Fa(46)1084 482 y Ff(c)1105 489 y Fa(54)1191 482 y Ff(c)1212
489 y Fa(55)1298 482 y Ff(c)1319 489 y Fa(56)1084 542
y Ff(c)1105 549 y Fa(64)1191 542 y Ff(c)1212 549 y Fa(65)1298
542 y Ff(c)1319 549 y Fa(66)1365 384 y Fd(1)1365 457
y(C)1365 483 y(A)1409 220 y(1)1409 293 y(C)1409 318 y(C)1409
343 y(C)1409 367 y(C)1409 392 y(C)1409 417 y(C)1409 442
y(C)1409 469 y(A)1454 220 y(0)1454 293 y(B)1454 318 y(B)1454
343 y(B)1454 367 y(B)1454 392 y(B)1454 417 y(B)1454 442
y(B)1454 469 y(@)1499 205 y(0)1499 278 y(B)1499 305 y(@)1543
243 y Ff(\017)1563 250 y Fa(1)1543 303 y Ff(\017)1563
310 y Fa(2)1543 363 y Ff(\017)1563 370 y Fa(3)1591 205
y Fd(1)1591 278 y(C)1591 305 y(A)1499 384 y(0)1499 457
y(B)1499 483 y(@)1543 421 y Ff(\017)1563 428 y Fa(4)1543
482 y Ff(\017)1563 489 y Fa(5)1543 542 y Ff(\017)1563
549 y Fa(6)1591 384 y Fd(1)1591 457 y(C)1591 483 y(A)1636
220 y(1)1636 293 y(C)1636 318 y(C)1636 343 y(C)1636 367
y(C)1636 392 y(C)1636 417 y(C)1636 442 y(C)1636 469 y(A)1778
392 y Fj(\(0)p Ff(:)p Fj(3\))78 681 y(The)h(4)h(quadran)o(ts)f(of)h
(the)e(sti\013ness)i(matrix)d Ff(c)i Fj(are)g(abbreviated)g(b)o(y)1391
621 y Fd(\022)1430 651 y Ff(A)1467 658 y Fa(11)1552 651
y Ff(A)1589 658 y Fa(12)1430 712 y Ff(A)1467 719 y Fa(21)1552
712 y Ff(A)1589 719 y Fa(22)1634 621 y Fd(\023)1682 681
y Fj(and)h(their)78 763 y(lo)q(cations)j(on)f(the)g(grid)g(are)g(noted)
h(in)e(Figure)h @beginspecial @setspecial @endspecial(0.1.)33
b(Sti\013ness)20 b(v)m(alues)g(in)g(the)g(\014rst)g(diagonal)78
823 y Ff(A)115 830 y Fa(11)171 823 y Fj(quadran)o(t)f(and)h(the)e
(\014rst)h(p)q(ortion)h(of)f(the)g Ff(\033)1017 830 y
Fa(1)p Fe(::)p Fa(3)1093 823 y Fj(and)g Ff(\017)1210
830 y Fa(1)p Fe(::)p Fa(3)1286 823 y Fj(v)o(ectors)f(are)h(c)o(hosen)g
(to)g(lie)f(on)78 883 y(the)j(primary)e(grid.)36 b(The)21
b(second)g(diagonal)h(quadran)o(t)g(and)g(the)e(second)i(part)f(of)h
Ff(\033)1722 890 y Fa(4)p Fe(::)p Fa(6)1799 883 y Fj(and)78
943 y Ff(\017)98 950 y Fa(4)p Fe(::)p Fa(6)199 943 y
Fj(are)g(then)g(necessarily)e(c)o(hosen)i(to)h(lie)e(on)h(the)g
(secondary)g(grid.)39 b(The)22 b(primary)e(cub)q(e)78
1004 y(cen)o(ter)f(coincides)g(with)h(the)g(secondary)h(cub)q(e)f
(edge.)33 b(That)21 b(sc)o(heme)d(causes)j(the)f(deriv)m(ativ)o(e)78
1064 y(computation)e(to)h(ev)m(aluate)f(quan)o(tities)g(at)h(the)f
(correct)g(collo)q(cation)h(p)q(oin)o(ts.)28 b(In)19
b(anisotropic)78 1124 y(elastic)13 b(media)f(up)i(to)f(orthorhom)o(bic)
g(symmetry)l(,)d(the)k(previous)f(notion)h(w)o(orks)g(w)o(ell,)e(since)
h(the)78 1184 y(o\013-diagonal)21 b(elemen)o(ts)15 b
Ff(A)583 1191 y Fa(12)620 1184 y Ff(;)8 b(A)679 1191
y Fa(21)735 1184 y Fj(of)19 b(the)f(sti\013ness)h(matrix)e(are)i(zero)f
(and)i(th)o(us)e(do)i(not)f(ha)o(v)o(e)78 1244 y(to)c(b)q(e)g(ev)m
(aluated.)21 b(In)15 b(anisotropic)g(systems)f(b)q(ey)o(ond)h
(orthorhom)o(bic)f(those)h(elemen)o(ts)d(are)j(non)78
1304 y(zero,)f(suc)o(h)g(that)g(ev)o(en)f(in)h(non-split)g(sc)o(hemes)e
(in)o(terp)q(olation)i(to)g(collo)q(cation)g(p)q(oin)o(ts)h(has)f(to)h
(b)q(e)78 1365 y(included)i(for)i(staggered)g(grid)f(metho)q(ds.)27
b(When)18 b(splitting)g(the)g(w)o(a)o(v)o(e)g(op)q(erator,)h(one)g(has)
g(to)78 1425 y(tak)o(e)e(care)h(of)g(prop)q(er)g(collo)q(cation)g(in)o
(terp)q(olation)f(when)h(using)h(a)f(staggered)g(grid.)26
b(Th)o(us)18 b(the)78 1485 y(in)o(terp)q(olation)e(issue)h(can)g
(complicate)d(the)j(primary)e(goal)j(of)f(adapting)g(deriv)m(ativ)o(e)f
(op)q(erators)78 1545 y(to)h(the)f(quan)o(tities.)44
b(Splitting)16 b(of)g(the)h(w)o(a)o(v)o(e)e(op)q(erator)j(also)f
(a\013ects)g(the)f(n)o(umerical)e(stabilit)o(y)78 1914
y(Figure)36 b(0.1:)61 b(The)36 b(primary)e(grid)78 1962
y(and)18 b(secondary)g(grid)g(as)g(usually)f(de-)78 2011
y(\014ned)26 b(for)h(staggered)g(grid)g(applica-)78 2059
y(tions)i(for)g(elastic)f(w)o(a)o(v)o(e)g(mo)q(deling.)78
2107 y(The)34 b(primary)e(cub)q(e)i(cen)o(ter)f(coin-)78
2155 y(cides)i(with)i(the)f(secondary)g(cub)q(e)78 2203
y(edge.)70 b(The)33 b(p)q(ositions)g(of)g(stress,)78
2251 y(strain)g(and)h(sti\013ness)g(comp)q(onen)o(ts)78
2299 y(are)14 b(indicated)f(\(blac)o(k)g(triangles\),)g(to-)78
2348 y(gether)19 b(with)g(the)f(ev)m(aluation)i(p)q(oin)o(ts)78
2396 y(of)e(the)f(split)f(w)o(a)o(v)o(e)h(op)q(erator)h(sc)o(heme)78
2444 y(\(empt)o(y)c(b)q(o)o(xes\).)p 440 2408 158 2 v
440 2463 2 55 v 57 w @beginspecial @setspecial @endspecial
Fa(stagsplit)p 596 2463 V 440 2465 158 2 v 28 w Fj([NR])889
1701 y
 16056010 14208860 2433925 11117117 36377313 40916254 startTexFig
 889 1701 a
%%BeginDocument: ././../split/Fig/stagsplit.ps.save

/arrowhead {
0 begin
transform originalCTM itransform
/taily exch def
/tailx exch def
transform originalCTM itransform
/tipy exch def
/tipx exch def
/dy tipy taily sub def
/dx tipx tailx sub def
/angle dx 0 ne dy 0 ne or { dy dx atan } { 90 } ifelse def
gsave
originalCTM setmatrix
tipx tipy translate
angle rotate
newpath
arrowHeight neg arrowWidth 2 div moveto
0 0 lineto
arrowHeight neg arrowWidth 2 div neg lineto
patternNone not {
originalCTM setmatrix
/padtip arrowHeight 2 exp 0.25 arrowWidth 2 exp mul add sqrt brushWidth mul
arrowWidth div def
/padtail brushWidth 2 div def
tipx tipy translate
angle rotate
padtip 0 translate
arrowHeight padtip add padtail add arrowHeight div dup scale
arrowheadpath
ifill
} if
brushNone not {
originalCTM setmatrix
tipx tipy translate
angle rotate
arrowheadpath
istroke
} if
grestore
end
} dup 0 9 dict put def

/arrowheadpath {
newpath
arrowHeight neg arrowWidth 2 div moveto
0 0 lineto
arrowHeight neg arrowWidth 2 div neg lineto
} def

/leftarrow {
0 begin
y exch get /taily exch def
x exch get /tailx exch def
y exch get /tipy exch def
x exch get /tipx exch def
brushLeftArrow { tipx tipy tailx taily arrowhead } if
end
} dup 0 4 dict put def

/rightarrow {
0 begin
y exch get /tipy exch def
x exch get /tipx exch def
y exch get /taily exch def
x exch get /tailx exch def
brushRightArrow { tipx tipy tailx taily arrowhead } if
end
} dup 0 4 dict put def


/arrowHeight 10 def
/arrowWidth 5 def

/IdrawDict 52 dict def
IdrawDict begin

/reencodeISO {
dup dup findfont dup length dict begin
{ 1 index /FID ne { def }{ pop pop } ifelse } forall
/Encoding ISOLatin1Encoding def
currentdict end definefont
} def

/ISOLatin1Encoding [
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
/space/exclam/quotedbl/numbersign/dollar/percent/ampersand/quoteright
/parenleft/parenright/asterisk/plus/comma/minus/period/slash
/zero/one/two/three/four/five/six/seven/eight/nine/colon/semicolon
/less/equal/greater/question/at/A/B/C/D/E/F/G/H/I/J/K/L/M/N
/O/P/Q/R/S/T/U/V/W/X/Y/Z/bracketleft/backslash/bracketright
/asciicircum/underscore/quoteleft/a/b/c/d/e/f/g/h/i/j/k/l/m
/n/o/p/q/r/s/t/u/v/w/x/y/z/braceleft/bar/braceright/asciitilde
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
/.notdef/dotlessi/grave/acute/circumflex/tilde/macron/breve
/dotaccent/dieresis/.notdef/ring/cedilla/.notdef/hungarumlaut
/ogonek/caron/space/exclamdown/cent/sterling/currency/yen/brokenbar
/section/dieresis/copyright/ordfeminine/guillemotleft/logicalnot
/hyphen/registered/macron/degree/plusminus/twosuperior/threesuperior
/acute/mu/paragraph/periodcentered/cedilla/onesuperior/ordmasculine
/guillemotright/onequarter/onehalf/threequarters/questiondown
/Agrave/Aacute/Acircumflex/Atilde/Adieresis/Aring/AE/Ccedilla
/Egrave/Eacute/Ecircumflex/Edieresis/Igrave/Iacute/Icircumflex
/Idieresis/Eth/Ntilde/Ograve/Oacute/Ocircumflex/Otilde/Odieresis
/multiply/Oslash/Ugrave/Uacute/Ucircumflex/Udieresis/Yacute
/Thorn/germandbls/agrave/aacute/acircumflex/atilde/adieresis
/aring/ae/ccedilla/egrave/eacute/ecircumflex/edieresis/igrave
/iacute/icircumflex/idieresis/eth/ntilde/ograve/oacute/ocircumflex
/otilde/odieresis/divide/oslash/ugrave/uacute/ucircumflex/udieresis
/yacute/thorn/ydieresis
] def
/Helvetica-Bold reencodeISO def
/Symbol dup findfont def

/none null def
/numGraphicParameters 17 def
/stringLimit 65535 def

/Begin {
save
numGraphicParameters dict begin
} def

/End {
end
restore
} def

/SetB {
dup type /nulltype eq {
pop
false /brushRightArrow idef
false /brushLeftArrow idef
true /brushNone idef
} {
/brushDashOffset idef
/brushDashArray idef
0 ne /brushRightArrow idef
0 ne /brushLeftArrow idef
/brushWidth idef
false /brushNone idef
} ifelse
} def

/SetCFg {
/fgblue idef
/fggreen idef
/fgred idef
} def

/SetCBg {
/bgblue idef
/bggreen idef
/bgred idef
} def

/SetF {
/printSize idef
/printFont idef
} def

/SetP {
dup type /nulltype eq {
pop true /patternNone idef
} {
dup -1 eq {
/patternGrayLevel idef
/patternString idef
} {
/patternGrayLevel idef
} ifelse
false /patternNone idef
} ifelse
} def

/BSpl {
0 begin
storexyn
newpath
n 1 gt {
0 0 0 0 0 0 1 1 true subspline
n 2 gt {
0 0 0 0 1 1 2 2 false subspline
1 1 n 3 sub {
/i exch def
i 1 sub dup i dup i 1 add dup i 2 add dup false subspline
} for
n 3 sub dup n 2 sub dup n 1 sub dup 2 copy false subspline
} if
n 2 sub dup n 1 sub dup 2 copy 2 copy false subspline
patternNone not brushLeftArrow not brushRightArrow not and and { ifill } if
brushNone not { istroke } if
0 0 1 1 leftarrow
n 2 sub dup n 1 sub dup rightarrow
} if
end
} dup 0 4 dict put def

/Circ {
newpath
0 360 arc
patternNone not { ifill } if
brushNone not { istroke } if
} def

/CBSpl {
0 begin
dup 2 gt {
storexyn
newpath
n 1 sub dup 0 0 1 1 2 2 true subspline