您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
xx柑橘标准化高效生态种植示范基地建设项目立项建议书.doc 51页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值，立即自动返金币，充值渠道很便利
你可能关注的文档：
··········
··········
中央现代农业生产发展资金
柑橘标准化高效生态种植示范基地建设项目
宜都市财政局
特产技术推广中心
二00九年六月
一、项目概述
1、项目名称：宜都市柑橘标准化高效生态示范基地建设项目
2、建设性质：改（扩）建
3、项目主管部门：宜都市财政局
4、项目承建单位：宜都市特产技术推广中心
5、项目负责人：裴学中　钟世泉
柑橘是世界第一大水果。目前，总产量达1.1亿吨，贸易量位居农产品第三位。2008年我国柑橘总产量达2000万吨。近几年来，宜都市的柑橘产业发展迅猛，2008年面积达28.7万亩，产量34.77万吨，柑橘系列产值达到15亿元，已成为全市29.2万农民现金收入的主要来源，柑橘面积与产量均占居湖北省第二位。利用柑橘无公害标准化综合生产技术措施实施的柑橘精品果园建设走在全省前列，柑橘精品园生产的精品蜜柑以其果形端正、皮薄光滑、色泽鲜艳、汁多化渣、风味浓甜、香气宜人、营养丰富，果实可溶性固形物达12%以上（最高达14%），受到农业部、省农业厅、华中农业大学等单位领导专家的好评，农业部、省农业厅领导多次到实地检查指导工作并给予了高度评价。
宜都市柑橘标准化建设项目区位于柑橘主产区的红花套镇、高坝洲镇和陆城镇，项目区连片，位于长江三峡出口南岸，紧邻318国道，水陆两路交通便捷。项目区形成于上世纪七十年代优质蜜柑基地的建设基础上，随着水果业的快速发展，种植规模不断扩大，已成为我省乃至全国较有名气的蜜桔生产基地。
但由于受过去“矮密早丰”栽培技术和种植技术不配套等因素的影响，普遍存在桔园密度偏大、建园质量不高、部分桔园树龄老化（有的已达40年以上，都陆续达到了柑橘栽培的经济寿命年限）等突出问题，且基地内的渠系、道路、电力配套设施老化，基地配套建设亟待升级改造。因而，建设宜都市柑橘标准化高效生态示范基地建设项目显得尤为重要，也十分必要和可行。
依据中华人民共和国财政部《关于落实好2009年财政支持现代农业生产发展政策工作的通知》（财办农〔2009〕1号）、《湖北省现代农业生产发展主导产业项目实施方案》，结合宜都市柑橘产业发展的实际，为加快柑橘标准化高效生态示范基地的建设步伐，促进我市（县）柑橘产业的可持续发展，柑橘产业生产能力和柑橘产品市场竞争力大幅提升，规模经营水平不断提高，保持宜都市经济的稳定发展和农民持续增收，编制本项目建设实施方案。
二、主导产业发展状况
（一）宜都市柑橘产业发展概述
宜都市地处武陵山余脉，位于鄂西南山区向江汉平原过渡地带,地处东径111°05′47″—111°36′06″，北纬30°05′53″—30°36′03″之间，东邻枝江，南接松滋，北与宜昌，西南与五峰、长阳毗邻。
全市版图面积1357km2，其中山区、丘陵占总面积88.4%，平原占11.6%，市内平原、丘陵、半高山呈梯级分布，西南高，东北低，海拔38-1081m。是江汉平原向鄂西山区过渡地带。属亚热带季风气候，光、热、水资源丰富，气候温和，雨量充沛，日照充足，四季分明。多年的年平均温度16.7℃，≥10℃活动积温5354.4℃，年降雨量1233.8mm，且雨热同季，无霜期273.3天，年均日照时数1641小时。地理、气候条件极适合柑橘的生产发展，是全国柑橘区划中宽皮柑橘最适宜区。
宜都地处鄂渝湘三省交界，是我国中西部地区重要的交通枢纽和物资集散中心，水陆交通发达，区位优势明显。从市政府陆城至三峡机场仅25公里，市域内有宜昌、枝城两座长江大桥，有我国九大水铁联运港之一的枝城港，贯穿南北的焦柳铁路、长江黄金水道和沪蓉高速公路构筑起立体交通网络，柑橘产品外销十分快捷、方便。
宜都市柑橘规模优势明显：截止2008年，宜都市柑橘28.7万亩（挂果面积17万亩），产量34.7万吨。2008年全市柑橘平均单产达到2040公斤/亩，全市人均占有柑橘总量达到0.88吨，2008年水果总产值4.1亿元，占种植业总产值的40%以上，柑橘系列加工产值达到15亿元。宜都柑橘曾于2001年被农业部、外经部定为首批全国园艺产品出口示范区；2002年被农业部定为全国柑橘生产优势区域，同时纳入全省九大农产品优势区域规划；2003年被省农业厅定为全省优质水果工程示范县（市）和无公害创建示范县（市）；2004年被国家林业局命名为“中国柑橘之乡”。宜都市同时还是全国柑橘标准化示范县（市）；全国柑橘富民强县项目示范县（市）；全省优质柑橘板块重点县（市）。宜都蜜柑荣获2008年度湖北省十大名牌农产品，2008年被湖北省人民政府命名为“水果大市”。
宜都市柑橘苗木繁育体系健全：宜都建有大型柑橘苗木繁育场5家，其中位于宜都市枝城镇大堰堤村的湖北省果树良种苗木繁育场，是中南地区最大的国家级果树良种苗木繁育基地。该场于2002年3月正式动工，2005年5月第一期项目建设完工，总面积1
正在加载中，请稍后...%!PS-Adobe-2.0
%%Creator: dvipsk 5.526b Copyright
Radical Eye Software
%%Title: DDall.dvi
%%Pages: 15
%%PageOrder: Ascend
%%BoundingBox: 0 0 596 842
%%EndComments
%DVIPSCommandLine: /local/TeX/bin/SunOS5.6/dvips DDall
%DVIPSParameters: dpi=600, compressed, comments removed
%DVIPSSource:
TeX output :1433
%%BeginProcSet: texc.pro
/TeXDict 250 dict def TeXDict begin /N{def}def /B{bind def}N /S{exch}N /X{S N}
B /TR{translate}N /isls false N /vsize 11 72 mul N /@rigin{isls{[0 -1 1 0 0 0]
concat}if 72 Resolution div 72 VResolution div neg scale isls{Resolution hsize
-72 div mul 0 TR}if Resolution VResolution vsize -72 div 1 add mul TR matrix
currentmatrix dup dup 4 get round 4 exch put dup dup 5 get round 5 exch put
setmatrix}N /@landscape{/isls true N}B /@manualfeed{statusdict /manualfeed
true put}B /@copies{/#copies X}B /FMat[1 0 0 -1 0 0]N /FBB[0 0 0 0]N /nn 0 N
/IE 0 N /ctr 0 N /df-tail{/nn 8 dict N nn begin /FontType 3 N /FontMatrix
fntrx N /FontBBox FBB N string /base X array /BitMaps X /BuildChar{
CharBuilder}N /Encoding IE N end dup{/foo setfont}2 array copy cvx N load 0 nn
put /ctr 0 N[}B /df{/sf 1 N /fntrx FMat N df-tail}B /dfs{div /sf X /fntrx[sf 0
0 sf neg 0 0]N df-tail}B /E{pop nn dup definefont setfont}B /ch-width{ch-data
dup length 5 sub get}B /ch-height{ch-data dup length 4 sub get}B /ch-xoff{128
ch-data dup length 3 sub get sub}B /ch-yoff{ch-data dup length 2 sub get 127
sub}B /ch-dx{ch-data dup length 1 sub get}B /ch-image{ch-data dup type
/stringtype ne{ctr get /ctr ctr 1 add N}if}B /id 0 N /rw 0 N /rc 0 N /gp 0 N
/cp 0 N /G 0 N /sf 0 N /CharBuilder{save 3 1 roll S dup /base get 2 index get
S /BitMaps get S get /ch-data X pop /ctr 0 N ch-dx 0 ch-xoff ch-yoff ch-height
sub ch-xoff ch-width add ch-yoff setcachedevice ch-width ch-height true[1 0 0
-1 -.1 ch-xoff sub ch-yoff .1 add]/id ch-image N /rw ch-width 7 add 8 idiv
string N /rc 0 N /gp 0 N /cp 0 N{rc 0 ne{rc 1 sub /rc X rw}{G}ifelse}
imagemask restore}B /G{{id gp get /gp gp 1 add N dup 18 mod S 18 idiv pl S get
exec}loop}B /adv{cp add /cp X}B /chg{rw cp id gp 4 index getinterval
putinterval dup gp add /gp X adv}B /nd{/cp 0 N rw exit}B /lsh{rw cp 2 copy get
dup 0 eq{pop 1}{dup 255 eq{pop 254}{dup dup add 255 and S 1 and or}ifelse}
ifelse put 1 adv}B /rsh{rw cp 2 copy get dup 0 eq{pop 128}{dup 255 eq{pop 127}
{dup 2 idiv S 128 and or}ifelse}ifelse put 1 adv}B /clr{rw cp 2 index string
putinterval adv}B /set{rw cp fillstr 0 4 index getinterval putinterval adv}B
/fillstr 18 string 0 1 17{2 copy 255 put pop}for N /pl[{adv 1 chg}{adv 1 chg
nd}{1 add chg}{1 add chg nd}{adv lsh}{adv lsh nd}{adv rsh}{adv rsh nd}{1 add
adv}{/rc X nd}{1 add set}{1 add clr}{adv 2 chg}{adv 2 chg nd}{pop nd}]dup{
bind pop}forall N /D{/cc X dup type /stringtype ne{]}if nn /base get cc ctr
put nn /BitMaps get S ctr S sf 1 ne{dup dup length 1 sub dup 2 index S get sf
div put}if put /ctr ctr 1 add N}B /I{cc 1 add D}B /bop{userdict /bop-hook
known{bop-hook}if /SI save N @rigin 0 0 moveto /V matrix currentmatrix dup 1
get dup mul exch 0 get dup mul add .99 lt{/QV}{/RV}ifelse load def pop pop}N
/eop{SI restore showpage userdict /eop-hook known{eop-hook}if}N /@start{
userdict /start-hook known{start-hook}if pop /VResolution X /Resolution X 1000
div /DVImag X /IE 256 array N 0 1 255{IE S 1 string dup 0 3 index put cvn put}
for 65781.76 div /vsize X 65781.76 div /hsize X}N /p{show}N /RMat[1 0 0 -1 0 0
]N /BDot 260 string N /rulex 0 N /ruley 0 N /v{/ruley X /rulex X V}B /V{}B /RV
statusdict begin /product where{pop product dup length 7 ge{0 7 getinterval
dup(Display)eq exch 0 4 getinterval(NeXT)eq or}{pop false}ifelse}{false}
ifelse end{{gsave TR -.1 -.1 TR 1 1 scale rulex ruley false RMat{BDot}
imagemask grestore}}{{gsave TR -.1 -.1 TR rulex ruley scale 1 1 false RMat{
BDot}imagemask grestore}}ifelse B /QV{gsave transform round exch round exch
itransform moveto rulex 0 rlineto 0 ruley neg rlineto rulex neg 0 rlineto fill
grestore}B /a{moveto}B /delta 0 N /tail{dup /delta X 0 rmoveto}B /M{S p delta
add tail}B /b{S p tail}B /c{-4 M}B /d{-3 M}B /e{-2 M}B /f{-1 M}B /g{0 M}B /h{
1 M}B /i{2 M}B /j{3 M}B /k{4 M}B /w{0 rmoveto}B /l{p -4 w}B /m{p -3 w}B /n{p
-2 w}B /o{p -1 w}B /q{p 1 w}B /r{p 2 w}B /s{p 3 w}B /t{p 4 w}B /x{0 S rmoveto}
B /y{3 2 roll p a}B /bos{/SS save N}B /eos{SS restore}B end
%%EndProcSet
%%BeginProcSet: special.pro
TeXDict begin /SDict 200 dict N SDict begin /@SpecialDefaults{/hs 612 N /vs
792 N /ho 0 N /vo 0 N /hsc 1 N /vsc 1 N /ang 0 N /CLIP 0 N /rwiSeen false N
/rhiSeen false N /letter{}N /note{}N /a4{}N /legal{}N}B /@scaleunit 100 N
/@hscale{@scaleunit div /hsc X}B /@vscale{@scaleunit div /vsc X}B /@hsize{/hs
X /CLIP 1 N}B /@vsize{/vs X /CLIP 1 N}B /@clip{/CLIP 2 N}B /@hoffset{/ho X}B
/@voffset{/vo X}B /@angle{/ang X}B /@rwi{10 div /rwi X /rwiSeen true N}B /@rhi
{10 div /rhi X /rhiSeen true N}B /@llx{/llx X}B /@lly{/lly X}B /@urx{/urx X}B
/@ury{/ury X}B /magscale true def end /@MacSetUp{userdict /md known{userdict
/md get type /dicttype eq{userdict begin md length 10 add md maxlength ge{/md
md dup length 20 add dict copy def}if end md begin /letter{}N /note{}N /legal{
}N /od{txpose 1 0 mtx defaultmatrix dtransform S atan/pa X newpath clippath
mark{transform{itransform moveto}}{transform{itransform lineto}}{6 -2 roll
transform 6 -2 roll transform 6 -2 roll transform{itransform 6 2 roll
itransform 6 2 roll itransform 6 2 roll curveto}}{{closepath}}pathforall
newpath counttomark array astore /gc xdf pop ct 39 0 put 10 fz 0 fs 2
F/|______Courier fnt invertflag{PaintBlack}if}N /txpose{pxs pys scale ppr
aload pop por{noflips{pop S neg S TR pop 1 -1 scale}if xflip yflip and{pop S
neg S TR 180 rotate 1 -1 scale ppr 3 get ppr 1 get neg sub neg ppr 2 get ppr 0
get neg sub neg TR}if xflip yflip not and{pop S neg S TR pop 180 rotate ppr 3
get ppr 1 get neg sub neg 0 TR}if yflip xflip not and{ppr 1 get neg ppr 0 get
neg TR}if}{noflips{TR pop pop 270 rotate 1 -1 scale}if xflip yflip and{TR pop
pop 90 rotate 1 -1 scale ppr 3 get ppr 1 get neg sub neg ppr 2 get ppr 0 get
neg sub neg TR}if xflip yflip not and{TR pop pop 90 rotate ppr 3 get ppr 1 get
neg sub neg 0 TR}if yflip xflip not and{TR pop pop 270 rotate ppr 2 get ppr 0
get neg sub neg 0 S TR}if}ifelse scaleby96{ppr aload pop 4 -1 roll add 2 div 3
1 roll add 2 div 2 copy TR .96 dup scale neg S neg S TR}if}N /cp{pop pop
showpage pm restore}N end}if}if}N /normalscale{Resolution 72 div VResolution
72 div neg scale magscale{DVImag dup scale}if 0 setgray}N /psfts{S 65781.76
div N}N /startTexFig{/psf$SavedState save N userdict maxlength dict begin
/magscale false def normalscale currentpoint TR /psf$ury psfts /psf$urx psfts
/psf$lly psfts /psf$llx psfts /psf$y psfts /psf$x psfts currentpoint /psf$cy X
/psf$cx X /psf$sx psf$x psf$urx psf$llx sub div N /psf$sy psf$y psf$ury
psf$lly sub div N psf$sx psf$sy scale psf$cx psf$sx div psf$llx sub psf$cy
psf$sy div psf$ury sub TR /showpage{}N /erasepage{}N /copypage{}N /p 3 def
@MacSetUp}N /doclip{psf$llx psf$lly psf$urx psf$ury currentpoint 6 2 roll
newpath 4 copy 4 2 roll moveto 6 -1 roll S lineto S lineto S lineto closepath
clip newpath moveto}N /endTexFig{end psf$SavedState restore}N /@beginspecial{
SDict begin /SpecialSave save N gsave normalscale currentpoint TR
@SpecialDefaults count /ocount X /dcount countdictstack N}N /@setspecial{CLIP
1 eq{newpath 0 0 moveto hs 0 rlineto 0 vs rlineto hs neg 0 rlineto closepath
clip}if ho vo TR hsc vsc scale ang rotate rwiSeen{rwi urx llx sub div rhiSeen{
rhi ury lly sub div}{dup}ifelse scale llx neg lly neg TR}{rhiSeen{rhi ury lly
sub div dup scale llx neg lly neg TR}if}ifelse CLIP 2 eq{newpath llx lly
moveto urx lly lineto urx ury lineto llx ury lineto closepath clip}if
/showpage{}N /erasepage{}N /copypage{}N newpath}N /@endspecial{count ocount
sub{pop}repeat countdictstack dcount sub{end}repeat grestore SpecialSave
restore end}N /@defspecial{SDict begin}N /@fedspecial{end}B /li{lineto}B /rl{
rlineto}B /rc{rcurveto}B /np{/SaveX currentpoint /SaveY X N 1 setlinecap
newpath}N /st{stroke SaveX SaveY moveto}N /fil{fill SaveX SaveY moveto}N
/ellipse{/endangle X /startangle X /yrad X /xrad X /savematrix matrix
currentmatrix N TR xrad yrad scale 0 0 1 startangle endangle arc savematrix
setmatrix}N end
%%EndProcSet
TeXDict begin
600 600 (DDall.dvi)
@start /Fa 42 122 df40 DI<FAFCA35C147EA314FE5C003FB612
EFCA27E6C15EC7FC5CA30F5CAF6DC8FC24
DIII<15E0EC01F0AF141FEC3FE0147FEB01FF133FDFEB3F1F
3FA47EA214FEA25CA303A25CA
130FA25CAFB512C0B612E0A36C14C01C
DI59 D<EB07FE90383FFFF838
0FFE0FC383FCC7FC5AAEEB03FEA3FF0EC
7FE0903801FFC81FFE00EB7FF890C7FCEBFCFC13
E013E113E3EBFFF85C6C5B6C5B6C90C7FC90C8FCA8EA03CFA35B120FEA07
801E D<C091390FFF83EB512F790
6C13CFFC01FE5A
141E000F92C7FC5B121F5B123F90CAFCA35A127EA412FE5AA4ED0FA2151F
5E6C143F4BC7FC6D5B001F495A0FF03FF86CB55A6C14C06C5C6C6C48C8FC
EB1FF02B D<010FB7FCD900FEC7EA3F005CAE5CA2
C7FC5CA2A487EA1B55AA25BA25EEC800F013F13
1FA24A6C5A93C8FC4990C9FCAFE160F494A7EAC7FC5BA200035D
167E007FB612FEB7FCA25E7E3 D73
D82 DI<3B1FFFE01FFFE002F014F048
5C6CB01FCCACACA
5CACAC790C7FCA25D5A007E147EA215FE12FE485CA214015D
A5AA24A5A007EEB7FFC383F83FEEBFFFC6C5B6C5BC6
48C9FC2C D97 DIIIIIIII107 DIIIII114
E /Fb 4 123 df27
D57 D63 D122
D E /Fc 1 115 df114 D E /Fd 7 119
df68 D83 D101 D108 D111 D114
/Fe 15 120 df48 D50
DI107 D114 D116
D119 D E /Ff 8 116 df53
D101 D104 D109 DI114
E /Fg 3 103 df<EC01FE0AAA2140FEB1F8F90
387FCFEF3803C07FEAE5A121E123E003C133C127C147C5A1478A2
14F85AECF0E0A21301ECE1C0780EECF807C01E1C
E /Fh 5 112 df97
DI<EB3F833801FFEFC01F
48C67E123E805A805A91C7FCA5903801FFE0AA27EA27E380FC01F3807F03F
003FE31B1A7E D108 D<EB7F803801FFE0F807C381E001E003E13
1F003C7F007CA200F814C0A7007CEB0F80AE5B6C133E6C6C
5A6C6C5A3801FFEA1A7E D E /Fi 18 120
df53 D97 DI<B512FCECFF80390FC00F
E0E7E157C81A281A2ED0F80A316C0AFE5D5D4A5A
000FEB07E0B6128002FCC7FC20 DIIIII108 DIII114 DII
118 DI E /Fj 10 100 df18 DI<ED01CFED7F00EC01FC4A5AEC07E04A5A4A
5A4AC7FC5C14FE5C1CB3B3AA1AAC8FC137E5B48
5A485AEA0FE0EA1F8F8A2127EEA1F80EA0FE0EA03F06C7E6C7E137E7F8013
1F6D7EA26D7EAB3AA6E7E6E7E6E7EEC03F86E7EEC
007FED1FCE D48
99 D E /Fk 6 121 df97 D104 D109 D115 DI120 D E /Fl 19 118 df73
D97 DIIIIIIII108 DI
114 DIII E /Fm
D<14FF381F81FC01FC133FEB07C04848
EB03E0AF8A2003F15FCA390C8FC4815FEA54815FFB3A46C15FE
A56DFCAA26C6CEB03F0A36C6CEB07ED130F6C6CEB1F
806C6CEB3F381F81F8903807FFEFC2
D E /Fn 4 64 df27 D45
D E /Fo 31 122 df19 D44 DII<EC07F0EC3FFC91B5FCF80
981F803FEB3E8EB1F5AB121FEBC1
FCEB87FF393F9E07C038F001F0D87FE07FBA290C77E5AA25AA3
1DA5A007E5C003E495AA26C495A6D48C7FC380FC0FC6CB45A00
CC8FC2 D57 D65 D67
D69 D<D9FFFC923801FFFCEFFEFC02BF150D
1303023FED19F8A26E6C070106EEC3F06E7EEF
D170CAC140F011C5BA217C06E7E38
BEC00FCF7E92C7FC5EA25EA49167E
5EFCD80FF8011E14FEBFFFFC151C150C4
7FFE3A07F00FCF1FF49C7FC013E141E5B000115
0EAF93C7FC7F6C7E6D7E14F86CEBFF806D13F815FE011F7F6D
003F7FFF0AA12784B5A
127C007E4A5A007F4AC7FC6D133ED8FBC05B39F1FC03F839F07FFFE0D8E01F
FEC8FC2 D97 D99 DII103
DII<EB03F013FF14E0A30F14C0A33F
EA513FE5BAB8
DIIIIIIIII121 D E /Fp 10 102 df48 DI64 DIII80 D83 D90 D<D91FFFF0EB01C0000301FCEB0780486DEB0F003A0F80FFC07C3A3C
001FFFFE40DA1
D E /Fq 4 83 df49 D68 D73 D82 D E
/Fr 23 120 df46
D80 D97 DIIIIIII108 DIIII114 DII119 D E /Fs 19 123 df50 D58 D65 D77 D82 D
97 DII101 D107
109 D112 D114 DIIII120 D122
D E /Ft 10 113 df0 D3 D14 D20
D48 DI54 D86 D106
D<170C171C173CF017EEE0F1C163C
EB5ADFC030E003F141ED8F3E0131C
1F05BF05DEB7C015DEB3E034A5A011F90C8FC5CEB0F8E149EEB
07DC14FC6D5A5CE327C D E /Fu 7 113
df59 D72 D76 D80 D85
D104 D112 D E /Fv 2 67 df1 D66
D E /Fw 5 94 df72 D82 D86 DI<95261FFFC01C
7FCF303FA1F047F04E0963801FFE07FF
037F01C0C7D807FCF1FFF8912607FFF8C8D801FF060FCA6C6CDDFFFCC7FC90
2603FFF0DD1FE1F90CBD80FF84BB448C8FC2601FFF0DE03FE033F13E0000F
90CDD8FFE090261FFFFEC9FCD87FF874B712E0D8FF80081F03FCCAFC00FCCE000392CBFC
00E0EC D E /Fx 3 93 df3 DI92 D E /Fy 25 123 df48 D
II66 D<B7FC16F016FC16FF2703FE00
FEE3FF0EE0FFEAAAA24C5AA24C5AEE
3FF04C5A923807FFC0B8C7FC5E16F093C8FC3 D73 D<B500F090381FFFE0A4D803FEC8EAFC5E6E1303
6C5ED6E130F013F5D6E131F011F92C7FC6E5B010F143E6E137EE
13FCFF6D5CCEF023F5B15FF6E90C8FCA26E5AA36E
5AA26E5A6E5A3 D94 D<EB3FF83801FFFE48EBFF891FE0
0FEF013F812FFA290B5FCA30180C7FCAF14786C
7E6C6C13F8CB512E0C614C090381FFE001D1A7E
E /Fz 36 121 df12
D<15C0EC0FC5C495AC49C7FC5B131E133E5B
5BAAF5BAFCA25A123EA2FCA25AB31278
AEA37E7E7FE181A
D<EC01CFEB01FF90B5FCFC7FCAC5CA513015C
AF5CA33F497EB612F0A215E01C
D58 DI68 D70 DI<017FB612E017FC17FF
EB1FE0EF0FF002FFEC07F892C7EA03FCA318FEA25B5CA418FC
F8EF0FF0A2EF1FE0AEC7F00EE01FEEE0FF891B612E094C7FCD90FF8C9FC5C
A33F5CA D82 DI87
D<14FF381F01FCFA35BEAFC
A4EC0FFF49B5FC90390FFC3F00EB7FC0FAA007F
150E90C7FC15FE48ECFC1C7C3F801C13703A1F
C0F83FE02607FFE013C0C D99 DIIIIII107 DIIIII114 DIIIIII E /FA 13 70 df0 D10 D40 DI43
54 D61 D69 D E /FB 66 122 df11 DIIII18 D20
DI<01FC14E0003FEC01FD801F814F0A200
130F16C0BED3FEAA001F495A5D49485A
4A5A003F49C7FC143E495AEB01F048485AEB0F80D87E1EC8FC13F8EAFFEC9FC
2 DIIIII30 DII58 D
III64 DIIIII71
IIIIIIIIII<49B712C05BAF80D907FCC7EA
FFAA495A91CA5A4C5A011C4A5A013C4AC7FC4B5A
CA4B5A4B5A4BC8FC4A5A4A5A4A5A4A5AA24A5A4A5A4A
5A4AC712E85C495AA5FA49C748C7FC485A4848
5CFE1F48B6FCB7FC5E3
DIIIIIIII<151C157E15FEA315FC3FECFFCE0EBF0130F
EB0E0FA215C0AA
A25CAFEA25CAB127F495AA238FE0FCFC7
FC137CEA3FF0EA0FC01F6 DIIIIII114 DIIII
120 DI E /FC 20 115 df0
D33 D50 D65 DI91 D102 DI106 DI110 D<B912FEA36C17FCA26C17F801C0C912
786C6CC6C16EC16CCCC150E171E
6D6C141C173C6D6CCD6C5CCC
C7FC5E6D6C130E161E6E131CED803C023F1338EDC0EDE0FF09138
07F1C015F15FF6E90C8FCA26E5AA4
D E /FD 3 87 df74
DI86 D E /FE 16 62
IIIIIII61 D E /FF 37 123 df14
DI64 D67 DI70 D72 D74 D78 D82 D<001FB612FCAC003C151CC5AA3C648
5AAAA449C8FCA35B381FFFFE5C2
DIIIIIIIIIIIII114 DI
E /FG 7 98 df25
D<16FE727E6D17EFFE38007FF84BEC1FFC727E72
7E023F814B7F17C05DA314FF5DA35B5D60AAA61
ACC5A4A5E4E5A4D5B013F4B90C7FC4A4A5AEF1FFC4D5A017F4AB4
5A4A010F13CFCB812FC17E06C03FCC9FC4
D73 DII78 D97 D E /FH 32 118 df40 DI<EB03FE90383FFFE090B512F848EB07FC486C7E4848EB7F
F15C0A131F007F15F0A500FF15F8BAD
133FAA26C6CEB7F806C6CEBFF006C6C485ACEBFFFD9
03FEC7FC2 DIIIIIIIII65 D68 DII
78 DI84 D97 DI100 DI103 D105
D107 DI110
E /FI 48 123 df11 DI15 D18 D20 DI26 DII30 D33 DI58 DIIII67 DI71 D74 DI
III80 D83 D97 DIII102 D104
DI107 DIIII
E /FJ 2 77 df74 D76
D E /FK 2 119 df<267FFFF0903807FFE0B516F0A26C17E0D803FEC8EA7C005F6D4A5A
5F6C4B5A4C5A6E130F94C7FC161E6C153E5E6E015E6E485AF
91C8FC151EECF03E5D011FDECFBC0EB0FFF5D92C9FC5CC6D5A
D E /FL 6 107 df0 D<13E0A438F0E1E0EAF8E3387EEFC0381F
FF00EA03F8A2EA1FFF387EEFC038F8E3E0EAF0E
D48 DI106 D E /FM 45 121 df12
D16 DI21 DIII27
65 D68 DIIIIIIII<0FB5FC496E4AD8003F
FC1A3E4A7F033F163C826F6D147C4A7E021E6EF8DA3E037FDA3C
015E836FC163FF4AC6E0
4A6D5D18F070EBF80FFC90C8FC177F725A0107ED3FFF4A6E131E199E7113BE
010F17FE91C86C5B83AF3C705AAFD801FC5FB5
00F25A5 D80
D82 DI85 DI90 D97
D<EC07FC91387FFF391FF001F39038FFC00748
EB800F48EB001F485AAF049EB0FE0003FECFC127F5BA312FF5BA4
5BAC167C6D14F8003FEC03F06C6CEB07E0000FEC3FC03A07FC01FF000001
B512FC6C6C13F0D90FFEC7FC2 DIII104
D<01FC49B4EC7FC039C003FFF848D9C03FD9F00F7F3E0F1FE0FC0FF81F83
FF3D1F0FF1F007FC7CFBC003FEFCD9FF806D5AC49
4B14CC7FCF04C7F020F5C1B004A5CA201FF021F
5C624A5C190FF4A4AECF80FAF131FF2F01E91C75B077F133C48
03FF15EB15F8F2E1FEFFFC0D801E0
023CEC00FC59 DI112
DIII117 DIII
E /FN 86 128 df11 DIII
18 DI<0FE01FC001F013E907B00FE0016013
C0A2EBE0130E000E5B485B485B485B48485A
00C05B1B D38
IIIIIIIIIIIII61 D65 DIIIIIIIIIIIIIIIIIIIIIIIIII93
DI97 DIIIIIIIIIIIIIIIIIIIIIIIIIII127 D E /FO 7 102 df48
101 D E /FP 50 122 df11 D13
D25 DIII59
I67 DI72 DIIIIII80 D83 DI87
D96 DIIIIIIIIIIIIIIII114
DIII<1FF3FCEA703E38
603C3C00C8EBF00CA200F0387F
3FE0 DI E /FQ 80 126 df
0 DIII<151E15FEEC3FC0ECFFA49
5AA25CB3AA130F5C131F495A49C7FC13FEEA03FCEA0FF0EA7FC000FEC8FCA2EA7FC0EA0F
F0EA03FCC67E137F6D7E6D7E130FAA80A26D7E6D7E806D7EEC3FC0EC0FF0EC03
FEF DIIIII<151E153E157C15F8
EC01F0EC03E0EC1FE5CA95AAAA249
C7FCAFE5B1BAFA35B121FA45B123FA548C8FCA9
12FEB3AEA20F7FA203A27F07F137EA2
7FA26D7EA26D7EA26D7EA26D7EA26D7E6D7EA0EC0FC0EC07E0
EC00FE1F DIIIII<161E167EED01FEED3FE0ED7FC0EDFF4A
5A4A5A5D140F4A5A5D143F5D147F92C7FCA25C5CB3B3B3A95AAA
AF49C8FC485A485AEA07F0EA1FEFC12FCA2B4FCEA3FC06C7EEA
07F0EA03FC6C7E6C7E6D7E133F6D7E6D7EA26D7ED7EAB3B3A38080
A41F816E7EE6E7EED7FC0ED3FE0ED0FF8ED07FE1501
ED007E161E27C DIII32 DIII40 DI46 D<177C17FCEE01F8A2EE03F0EE07E0EE0FC0A2EE1F80163FEE7F0016
7E16FE4B5A5EB5AA24B5A153F5E157F93C7FC5D5DD14075D14
0FA24A5AA24A5AA34A5AA24AC8FCACA30FA25C131FA25C133F
A35C137FAC9FCA45A5BA20FA35BA23FA75BA3
127FAD5BA212FFB3A62E DIIIIIIIIIIIIII
64 DIIIII77 D80 DIIII88 DIIII98 DIIIIIII110
DIIIIIIII<EE7F80ED0FFF157F4AB5FC140F143F5C
49B6FCC7FCC8FCEFC485AEA1FF0485A
5B485A48CAFCA25A5A5A291B DIII E /FR 45 122 df0
DIII6 DI14 D20 DI26 DI33
D48 DII<EC01C0A24A7E4A7EEC0E38EC1E3C4A7EEC380EEC
780FECF49486C7EECCECCC
011C141C013C141EFE90C91270
001C161C003C161EA27E3
D59 D62 DI65 DII69
D86 DI91 DI102 DIIIII110 D<F1C173C
FFFC5E160E161E161C163CC5C
001FC00FF85CC8FCE017E131E017F5B
6D1F1370ECC0FEFFF6D90C9FC5C13005C
147CA7B D114
D121 D E /FS 79 123 df11 DIIII17 DI<EA01F0EA
03FCEA007E7FA26D7EAD7EAD7EA26D7EA47F495A
EB039FF1F33CF314F0EA07C039
0FF14FC5A3F2
DIIIIIII30 DIIII39 D58 DIIIIIIIIIIIIIIIIIIII80 DIIIIII
DIIIIIIIIIIII
IIIIIIIIIIIIII E /FT 52 127 df0 DI3 D6 D9 DI22 D34 D40 DI43 D48 DIIIIIIIIIII61 D72 D76 D80 D84
D91 D93 DI97 D99 DIIIIII108 DIII<01FFE0E03803
F380F800D381F03A2127C12FCAE6C1307EB80
0F380FC01BFFE300A7EC3FFEA31F257E
DIIIIIII<380F838
3FFFFFC6 D E /FU 79
121 df0 DIII6
24 DIIIIIII33 D36 D39 D41 D44 D47
D49 DIIIIIIII59 DI63 D65 DIIIIIIIIIIIIIIIIIII86 DI91 DI94 D98
DIIIIIIIII110 D<F907A2F10F80A2F11F00A1A24E5AA2
4E5AA24E5AA24E5AA296C7FC60A0A24D5AA24D5AA24D5AA24D5AA24DC8FC
AF000015EEA07FE00794B5AEAE1FE00C04B5AC67E6D4A5AA26E
495A133F6E49C9FC131F6E133E130F6E5BE485A
006E485A147FED8FFCAFC15DFEC1FFEA26E5AA26E5AA26E5AA26E5A5D42547B
D114 D120 D E /FV 70 125
df11 DI14 D19
D<EA0FC0EA1FF0EA3FF8EA7FFC12FF13FEA213FFA47E7E7EEA0FEFEA000F131F131EA313
3E133C137CEA01F0A2EA03E0EA07C0EA0F80EA1F
DII<EA0FC0EA1FF0EA3FF8EA7FFC12FF13FEA213FFA47E7E7EEA0FEFEA000F
131F131EACEA01F0A2EA03E0EA07C0EA0F80EA1F
DIII49 DII
IIIIIII75 DIIIIIIIIIIII90 D97 DIIIIIIII107 DIIIIIIIIIIIIIIIII E /FW 10 58 df48 DIIIII
IIII E /FX 21 118 df[159 145 120 272 176 65 D[<BC
FCB3B3B3B3B3B3B0BC 142 121 269 87 73 D[145 140 120 267 162
84 D97 D99 D[113 144 120 270 129 II[81 144 121 271 71 II[114 143 119 270 129 I[<EC3FC0ECFFF7F497F49
A290B67EADA26D5CA26D5C6D5C6D91C8FC6D5B6D5BEC
3FC091CAFCB3A3ED1FF0017FB5FCB7FCA9EA003FB3B3B0B
144 119 271 65 I[50 143 119 270 65 108 DIIIIIII[72
132 124 258 90 II
E /FY 6 117 df73
D86 D97 D114 D116 D E /FZ 60 126 df37 D40 DI43 DIII48 DIII<
EC3FC0903801FFFF7F497F5B9039FFE07FAEA07F85B4848EB3F
C7FC485AA290C9FC5A007EFFD8FE7F7F00FDB57E
B67EFFD9EF1F07AF
6C140F16C06D131F6C7EED3F806C6C137FCB55A6C5C6C5C6D
5B011FFCC7FC2 D<121FEA3F80EA7FC0EAFFE0A5EA7FC0EA3F80
EA1F00C7FCAE121FEA3F80EA7FC0EAFFE0A5EA7FC0EA3F80EA1F000B
FE497EAEFAC7AAF8A490387E
00FCAB512FEAFEB1F80AD87FFEEBFFFC
6D5AB514FEA26C15FC497E2 DIIII<8090390FFF87C0FF90B6FC5A481383EBFC
7F485A1FAC7FC5A127EEDFC12FE5AA8913803FFF8
4A13FCA27E127E6E13F27F6C141FA27F6C6C133FA26C7E6D137F6C6C13FF
6C6C5AEBFF036C90B5FC7E6D13EF011F13CF6DEBFCC7FC2
D73 D76 DI<EB7FFF0F14FEBE0FF90C77E007E
80A300FEFB3A86C143FA07F5CA26D5B393FF007FE90B5FCA26C
5C6C5CC66C90C7FC2 D82 DIIII<387FFFFC
B512FEA414FC00FCC7FCB3B3B1B512FC14FEA46C13FC1
D<387FFFFCB512FEA47ECB1387FFFFEB5FCA46C13FC1
D<3801FFF1F6D7E15E0488081EBC01FEC03FC381FEC8127EA3
903803FFFE133F48B5FC12075A48EBC07E383FF800EA7FC090C7FC12FE5AA36C14FE1401
387F8003EBE01F6CB612FC6C15FE7E6C143F0001EBFC1F3A007FC007FC2
III<EA7FEFFFFF04A13F8
A46E13F4A5A4A5A4A5A4A5A4A5A4AC7FC14FEEBF1FC13F3EBF7FE90B5FCA2
EC9F80EC0FC001FE7FEBFCC7E16803A7FFFC0FFFCB500E113
FEA46C01C013FC27 DIIIII<D87FFEEB3FC0B53801FFF1F13FC5C6C5B39003FFFE1ED00F802
FCC5C5CA25CA391C8FCAD007FB512E0B67EA46C5C24
IIIIIIII125 D E /F0 86 128 df11 DIII<127812FCA27EA27E7E7EEA1F07EA03E
13C00B DII25
D<07F00FEF813F00F131F80030
00A6013C0EBEB48
5B485B485BB1C D38 DIII44 DIIIIIIIIIIIIIII65
DIIIIIIIIIIIIIIIIIIIIIIII91
DII97 DIII
IIIIIIIIIIIIIIIIIII127 D E /F1 9 108 df0
D E /F2 88 127 df0 D10 DII14 D19 D24 D28
DI37 DIIIIIIIIIIIIIIIIIIIII61 D64 DIIIIIIIII75
82 DIIIIII<EAFFE0A3EAE000B3B3B3A7
EAFFE0A30B DII97 DIII
IIIIIIIIIIII
126 D E /F3 22 118 df46
D67 DI73 D
77 DI97 D99 DIII104 DI108 DIIII114 DIII E /F4 2 113 df104 D112 D E /F5 30 118
df49 DI52 D
D65 D70 D73 D<DAFFEFE133C013FB6EAE0FC80
1FFFD80FFF161F6A2
7FA26D157C7F7F6D92C7FC6D7E6C13F8ECFFC015FC6CECFFC06C15F816FE6C816C826C16
E06C827E013F81EB003FFFF88182
AEA26C5D6D5D7F6D4A5A6D4A5A01FC141FD9FFF801FFE0486C
B65AD8FC1F92C7FCD8F8F
DI86 D97 D99
DIIIIII108 DIIIIIIIII E /F6 85 123 df11 DIIIIIIIIIIIIIIIII30 DIIII37
D39 D44 D58 DIIII64
DIIIIIIIIIIIIIIIIIIIIIIIIII96 DIIIIIIIIIIIIIIIIIIIIIIIIII
E /F7 91 128 df0 DI5 D11 DIII19
DIIIIIIIIIIIIIII61 D65
DIIIIIIIIIIIIIIIIIIIIIIIIII93 DII97
DIIIIIIIIIIIIIIIIIIIIIIIIIII126
DI E /F8 64 128 df[120 140 124
267 116 11 D[99 140 123 267 111 I[154 140 123 267 166 14 D[46 198 109 276 77 40 D[<E127F6C7E6C7E6C7E7F6C7E6C7E6C7E6C7E6D
7E6D7E806D7E6D7EA26D7E6D7E807F6D7FE141F816E7EA282
80A26E7FA26F7EA0F82A36F7EAC081A317E0A281
A317F0AAFCB3AA17F8AA316FFA417E0A35DA217C0
A35D0A34B5AA35E151FA25E153F5EAB5AA24A5BA25C93C7FCA2
4A5A5D140F5D141F5D143F4A5A5D14FF5DB5C495A495AAC495A49
C9FC485A485A485A485A5B485A485A48CAFC12FE5A
117 276 77 I44 DII[<F201E0F203F0F207F8
A31A0F1BF0A21A1F1BE01A3FA21BC01A7F1B80A21AFF1B00A
0F62AFAA219FF97C7FCAF61A2181F
61AFA2FCA25F60AAA3F60
AFFA295C9FC5E5FAF161F5FA67F5F16
FFA294CAFC5D5EAE151F5EA57F5EA215FF93CBFC
5CA25D1DAFA25D143F5DA4FF92CCFC5BA25C1303
5CA30F5CAFA25C137F5CA213FF91CDFCA25A5B75B
A21F5B123FA25B127F5BA212FF90CEFCA
199 117 276 100 I[<160F4C7E163F167F16FFF92B5FC3B6
FC48B7FCB8FC153F14FC14E0EBFC0048C7FCC8FCB3B3B3B3B3AC4B7FA292B57E4A80023F
ECFF80007FBA12C0A54AC 110 259 100 49 D<EA01FCEA07FF
3F13E0AB512F8A76C13F0A26C13E0AC1300EA01FCC8FCB3B3A8EA
01FCEA07FFE0AB512F8A76C13F0A26C13E0AC1300
EA01FC1 D[<1A7E1AFFA34F7FA34F7FA34F7FA24F7FA34F7FA34F7FA34F
7FA396B6FCA34EAEAE061FF
187EFAFAFAA2050F83
4E80A53F4DFE835FA
AAAACA7E4C8594BAFCA293BCFCA34B
7EA24BFE875DAA24A6CC
B57E011F02E0FC063F15F0B800E0047FB9FCA5
38F2 122 269 149 65 D[122 136 121 263 141 I[121 144 117 267 144 I[133
136 121 263 152 I[123
135 121 262 136 I[114 135 121 262 130 I[135 144 117 267 156 I[<B500F10FB591383F
FFFC91B894B7FCA5DCA6CEDC7FC021F49CB000114
E0A26E4B093BDFCA604FCCCFCB3B3B3A44A6D4E80AC07
B600F8057FECFF80BA00CFCA6C
136 121 263 149 I[60
136 122 263 72 I[83 140 120 263 102
I[139 136 121 263 155 I[109 136 121 263
166 136 120 263 183 I[<B700FC95B5903803FFFCA901ED
FE00D9000F6EDD001F14EFC785BFF3DF6D735A03CF
7F795A03C3C07F707FA807080
A17F717FA7F727FA
7F737FA7F747FA4
7F757FAAF0A27613F8
7613FCA2FFA4CFA2FFAA
A28A8AA24A6C868A4A7E8A4A6CFFFE496D48B612FEB800FC197F
1F3FA21F1F1F0FDA4 136 121 263 149 I[132 144 117 267 155 I[117 136 121 263 136 I[
133 179 117 267 155 I[<BB12FEF2FFF8F3FF801CF81CFF1DC0D0F14F0
D9007F01FCDB007F13FC021F060F13FF4CC13E00A1F7F7F76
B898B777FAC89A38CA968A3659DC9FCA267AB67
535B679AB55A5291CAFC525B525B525B0A3F13ECBFC090F13FC097F13F0080F
B512C093BACCFC1CF81CC0A21CF004F8C0313FFEF13E0757F75
7F7FAA28AABA68BA68BFA0380FA
07C0AF8BA6D727FFF725E6D5CBA
6C6F6E137E5A776C6D485AFF70F90B55AD10003
5D0C005D0D1F49C7FC9C3 140 121 263 147 I[88 144 117 267 111 I[<000FC112F8A07
804ACEC003F02F01CA150191CAFF491E1F491E0F491E
4E00A390CB1A7EA07E203FA0FC21
CC97C7FCB3B3B3B3A74E80A24EDB712E0031FBB12FEA506F8C7
121 262 144 I[134 140 121 263 149 I[<B591267FFFF894261FFFE0
EBFFF091B797B6FCA5DCCEDFC1F14C001004ACDFC6E
5B6EE775AFC80701BFEA26E711903
6E831F0F6FFF99C9FCE41E
016F1E0FFAF7097CAFCE
FC1DFE715F1CAF725E861C1F72
5E7D91CCFCCFE725D1BBCBE7735C1B
FF735CACDFCABABAAA
140 124 263 149 I[
196 140 124 263 205 I[100 136 117 263 122 90
D97 D[97 140 122 265 111 II[97
140 121 265 111 II[<EF3FF1F14C2FFFC4B9038800FFE
92391FFE001F4BB4BA5BA24A5B4A5BA24A90C7FCA24A486E13
F5A4BEC07F095C8FCAABB912F0A6C7D87FFCCAFCB3B3B3B04A7E
AF14FEC 123 267 61 II[100 138 122 265 111 I[42
134 121 261 56 I[49 175 136 261 61 I[95
138 122 265 105 I[43 138 121 265 56
IIIIIIIIIIIIIIIII[59 20 108 261 100 127
D E /F9 116 128 df0 DIII5
DIIIIIIIII16 D18 DIIIII25 D
D31 D33 DII37 DIIIIIIIII
IIIIIIIIIII61
DIIIIIIIIIIIII
IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII
IIIIIIII<D801E3E039
0FFE01C0391FFFFFF0038700FFE38E003FCB6
E /F 10 115 df40 DI
46 D68 D73 D77 D82 D103 D111 D114 D E
/Ga 72 128 df11
D<EB01F0EB03F31FEB3FF0EB7FE014C0EBFFFC5B485A485A
485A5B48C7FC127C D24 D26 D28
D31 D44 DIIIIIIII<121EA2121F13
F090B812C0A005FACB5A5F16034C
5A94C7FC485D161E5ECE5EE15074BC8FCAC157C15
A5AA7AFAF5DA414FFA65BAC6D90C9
FC143C3 D65 DIIIIIIIIIIIIIIIIIIIIIIIII95 D97 DIIIIIIIIIIIIIIIIIIIIIIIIII127 D E /Gb 9 122
D100 DI<EC3FE0EB3FFF495BA3EB007F5D143F147FA292C9FCA25CA25CA21301A2
5CA307A25CA2130FA25CED07F8011FEB3FFF92B57EFE09138C7
CFDF807FA805C137F5C5C5C13FF5CA291C7FCDA25B00
BAF5B163F000F5EAC7FC121F4C03F0201
14F84C13F0491303EEFE0A290C7EC03C016F048EE
1E177C486EB45AEE7FE04 DI112
D E /Gc 21 117 df45 D[137 169 115 291 164 67 D[150 159 118 286
159 120 286 82 73 D[<B800C00C1FB712FA3DCF702FCC7FCD9
001FA1DE301E791C8FC0F36DE303C75BA203F16DF
F06EF40F07AEACAAAA
C0AAAEA6D62A6D4F5AA2
716D4F5AAAAFCA26D60A6D60A3726D4D
5AAAAAAFCA36D5EA36D5EA373
6D4B5AAAAAAFCA46D5CA46D5C
A46D495AA46E48CAFCA375EBC01EA56D5AA375
6D5AA275EBFDE0A375EBFFC0AE7591CBFCA24A7E765A4A7E4A6C
B512CF02FCB700C00B07B612FCB900C06F484BBA12F0A3765A
A3DACF9 118 286 209 77 D[<B86CB5FC7198
B7FCC7E36EDF003F14FA91C7FC
745A03FF07F71
7F747F775808A
EFFE8AAA7814FFA28A8BA28B8BA2
8B8BA28B8BA28BA24A6C888CA24A6C888C4A7E4A6CFC8790B7
00C086B912C11F210FDAC2 159 119
286 171 I97 D99 D[111 163 120 288 126 II[74 163 122 290 70 I[<151F91B5FC0003B6FCB7FCA7EA003F1303
7F7F80A280B3B3A7F13FFF0F14FF063F15C095B77EF13F89427
07FC00037FDD1FF001007FDD3FC06E7E4DC86C7F05FE6F7F4C486F7F4C5A4D6F7F4C5A4C
486F7F4C5A94C9FC4C707F163E5EA25E757F5EA35EA35EA55EB3B3B3A24A6D4C7FA249B5
00F0F02FE043F14F8B9D8E003B91280A6DAE000F
161 121 288 126 104 D[47
155 120 282 63 I[48 161
120 288 63 108 DIIII114 DI[
71 145 124 270 89 I E /Gd 109 128 df0 DI8 D10 DIIIIII18 DII23
DI31 DI34 DI38
DIIIIIIIIIIIIIIIIIIIII61 D63
D65 DIIIIIIIIIIIIIIIII
IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII<FFFFE0F00383FFFFE38783FFC38E00FF
16 DI E /Ge 21 117 df45 D<932601FFCFC
4BB66C131E030F15EFFFFF8A01C0EB007F020790
C8381F80FEDA0FFCED07C0DA3FF8ED03E1DA7FE0ED01F34A48ED00FB0CA12
3FAF015CCCFC5A1B7E
485AAE5BAEAA25B12FFAF127F7FAEA27F
AE000F1A1CA26C7E1B3C7E6E06C6D18FCEF01E080
011FF003C06D6CEFD6CEF0F006D6D161E6D6D5E6E6C5EDA3FF84B5ADA0FFC4B
5A6EB44B5AEC1F806E01F802FEC7FC923A3FFF800FFC030F90B512F
DB003F49C8FCF DI73 D77
D99 DIII104 D
%%EndProlog
%%BeginSetup
%%Feature: *Resolution 600dpi
TeXDict begin
%%PaperSize: A4
%%EndSetup
%%Page: 251 1
251 270 bop 300 814 a FW(30)300 1212 y F8(Hierarc)-6
b(hical)63 b(Basis)i(for)h(the)300 1411 y(Con)-6 b(v)g
(ection-Di\013usion)63 b(Equation)i(on)300 1611 y(Unstructured)g
(Meshes)300 1910 y Ga(Randolph)38 b(E.)g(Bank)g(and)h(Sabine)f(Gutsc)m
(h)300 2366 y FV(1)105 b(In)m(tro)s(duction)300 2565
y F9(The)35 b(Hierarc)n(hical)f(Basis)g(Multigrid)h(Metho)r(d)h(w)n(as)
e(originally)g(dev)n(elop)r(ed)g(for)h(sequences)g(of)300
2665 y(re\014ned)d(meshes.)h(Hierarc)n(hical)e(basis)h(functions)h(can)
f(b)r(e)h(constructed)f(in)h(a)f(straigh)n(tforw)n(ard)300
2765 y(fashion)24 b(on)h(suc)n(h)f(sequences)g(of)h(nested)f(meshes.)h
(The)g(HBMG)g(iteration)f(itself)h(is)f(just)i(a)e(blo)r(c)n(k)300
2864 y(symmetric)39 b(Gau\031-Seidel)f(iteration)g(applied)h(to)g(the)g
(sti\013ness)h(matrix)e(represen)n(ted)g(in)h(the)300
2964 y(hierarc)n(hical)27 b(basis.)h(Because)f(the)i(sti\013ness)g
(matrix)f(is)h(less)f(sparse)g(than)g(when)h(the)g(standard)300
3064 y(no)r(dal)c(basis)g(functions)i(are)d(used,)i(the)g(iteration)f
(is)h(carried)e(out)i(b)n(y)f(forming)g(the)i(hierarc)n(hical)300
3163 y(basis)35 b(sti\013ness)g(matrix)g(only)f(implicitly)-7
b(.)37 b(The)e(resulting)g(algorithm)f(is)h(strongly)f(connected)300
3263 y(to)27 b(the)h(classical)e(m)n(ultigrid)h(V-cycle,)g(except)h
(that)g(only)f(a)g(subset)g(of)h(the)f(unkno)n(wns)g(on)g(eac)n(h)300
3362 y(lev)n(el)g(is)h(smo)r(othed)f(during)g(the)h(relaxation)e(steps)
i([BD)n(Y88)o(].)383 3462 y(In)e(recen)n(t)f(y)n(ears,)f(w)n(e)i(ha)n
(v)n(e)e(generalized)h(suc)n(h)g(bases)g(to)h(completely)f
(unstructured)h(meshes,)300 3562 y(not)k(just)g(those)g(arising)f(from)
g(some)g(re\014nemen)n(t)h(pro)r(cess.)f(This)h(is)g(done)f(b)n(y)h
(recognizing)e(the)300 3661 y(strong)18 b(connection)i(b)r(et)n(w)n
(een)f(the)h(Hierarc)n(hical)e(Basis)h(Multigrid)g(Metho)r(d)i(and)e
(an)g(Incomplete)300 3761 y(LU)k(factorization)e(of)i(the)g(no)r(dal)f
(sti\013ness)h(matrix.)f(This)h(can)f(b)r(est)h(b)r(e)g(understo)r(o)r
(d)f(in)h(terms)g(of)300 3861 y(the)28 b(graph)e(elimination)h(mo)r
(del)h(for)f(Gaussian)f(elimination.)h(This)h(connection)f(is)g
(explored)f(in)300 3960 y(detail)21 b(in)h([BX94)o(].)g(Once)f(this)h
(connection)f(is)g(made,)g(it)h(is)f(fairly)g(easy)f(to)i(mak)n(e)e(a)h
(sym)n(b)r(olic)g(ILU)300 4060 y(algorithm)e(for)h(unstructured)g
(meshes)g(whic)n(h)g(mimics)g(the)h(ILU)f(pro)r(cess)f(on)h(a)g
(structured)f(mesh)300 4159 y(leading)i(to)g(the)h(classical)f(hierarc)
n(hical)e(basis.)i(This)h(sym)n(b)r(olic)f(elimination)g(pro)r(cess)g
(essen)n(tially)300 4259 y(de\014nes)j(the)g F0(supp)l(orts)g
F9(of)f(the)h(hierarc)n(hical)e(basis)h(functions)h(\(or)f(the)h
(sparsit)n(y)e(structure)h(of)h(the)300 4359 y(hierarc)n(hical)i(basis)
g(sti\013ness)i(matrix\).)383 4458 y(In)39 b(the)g(classical)e(case,)h
(certain)g(linear)g(com)n(binations)f(of)i(\014ne)g(grid)f(basis)g
(functions)g(are)p 300
v 300 4608 a F3(Ninth)25
b(International)i(Confer)l(enc)l(e)f(on)g(Domain)h(De)l(c)l(omp)l
(osition)g(Metho)l(ds)300 4682 y F2(Editor)e(P)n(etter)f(E.)f
(Bj\034rstad,)h(Magne)g(S.)f(Esp)r(edal)i(and)f(Da)n(vid)h(E.)e(Key)n
(es) y(c) y F1(\015)o F2(1998)i(DDM.org)p
%%Page: 252 2
252 271 bop 300 282 a Gd(252)2305 b(BANK)25 b(&)g(GUTSCH)300
482 y F9(formed,)39 b(with)g(the)g(com)n(bination)f(co)r(e\016cien)n
(ts)g(deriv)n(ed)g(from)h(the)g(geometry)e(of)i(the)g(mesh.)300
581 y(F)-7 b(or)30 b(these)g(sp)r(ecial)g(c)n(hoices,)g(the)h(linear)e
(com)n(binations)h(simplify)h(to)f(coarse)e(grid)i(no)r(dal)g(basis)300
681 y(functions.)23 b(In)f(the)h(case)e(of)i(completely)f(unstructured)
g(meshes,)g(the)h(co)r(e\016cien)n(ts)f(for)g(the)g(linear)300
780 y(com)n(binations)28 b(can)h(also)f(b)r(e)i(sp)r(eci\014ed)f(in)h
(a)f(natural)f(w)n(a)n(y)g(using)h(the)h(geometry)e(of)h(the)g(mesh.)
300 880 y(Ho)n(w)n(ev)n(er,)23 b(since)i(there)f(is)h(no)g(coarse)e
(grid)h(as)g(suc)n(h,)g(in)i(general)d(no)i(simpli\014cation)f(of)h
(the)g(basis)300 980 y(functions)37 b(o)r(ccurs.)e(Di\013eren)n(t)i(c)n
(hoices)f(of)g(expansion)g(co)r(e\016cien)n(ts)g(t)n(ypically)g(ha)n(v)
n(e)f(no)h(e\013ect)300 1079 y(on)e(the)h(supp)r(orts)f(of)h(the)f
(basis)g(functions,)h(but)g(do)f(ha)n(v)n(e)g(a)g(profound)g(e\013ect)h
(on)f(the)h(shap)r(e)300 1179 y(of)i(the)h(basis)f(functions)g
(themselv)n(es,)g(and)g(hence)h(on)f(the)h(n)n(umerical)e(v)-5
b(alues)37 b(app)r(earing)f(in)300 1279 y(the)31 b(hierarc)n(hical)d
(basis)i(sti\013ness)h(matrix)e([BX94],)h([BX96)o(].)h(One)f(can)g(ev)n
(en)g(c)n(ho)r(ose)f(di\013eren)n(t)300 1378 y(co)r(e\016cien)n(ts)19
b(for)f(the)i(trial)f(and)g(test)g(spaces,)g(leading)f(to)h(di\013eren)
n(t)g(hierarc)n(hical)f(basis)g(functions,)300 1478 y(similar)27
b(to)g(P)n(etro)n(v-Galerkin)e(\014nite)j(elemen)n(t)g(appro)n
(ximations.)383 1577 y(In)42 b(terms)g(of)g(ILU,)h(the)f(expansion)f
(co)r(e\016cien)n(ts)h(are)f(just)i(the)f(m)n(ultipliers)g(in)h(the)f
(ILU)300 1677 y(decomp)r(osition.)29 b(This)h(leads)g(one)g(to)f
(consider)g(the)i(p)r(ossibilit)n(y)e(of)h(de\014ning)g(these)h
(expansion)300 1777 y(co)r(e\016cien)n(ts)24 b(in)h(a)g(more)f
(algebraic)f(fashion.)h(F)-7 b(or)24 b(example,)h(one)f(can)h(c)n(ho)r
(ose)e(these)i(co)r(e\016cien)n(ts)300 1876 y(to)i(eliminate)g(certain)
g(o\013)g(diagonal)e(elemen)n(ts)i(of)g(the)h(hierarc)n(hical)d(basis)h
(sti\013ness)i(matrix,)e(as)300 1976 y(is)h(done)h(in)f(the)h(case)f
(of)g(classical)f(ILU.)i(While)g(geometry)e(based)h(co)r(e\016cien)n
(ts)g(seem)g(adequate)300 2076 y(for)h(isotropic,)f(self-adjoin)n(t)h
(problems)g(\(e.g.)g FU(\000)p F9(\001)p F6(u)c F9(=)g
F6(f)9 b F9(\),)28 b(w)n(e)g(ha)n(v)n(e)f(found)i(that)g(our)e
(algebraic)300 2175 y(c)n(hoices)h(can)g(greatly)g(impro)n(v)n(e)f(the)
i(robustness)f(of)h(the)g(HBMG)g(iteration)f(for)h(other)f(t)n(yp)r(es)
h(of)300 2275 y(equations,)e(notably)g(con)n(v)n(ection)f(dominated)i
(con)n(v)n(ection-di\013usion)e(equations.)383 2374 y(In)40
b(Section)g(2,)g(w)n(e)f(review)g(the)i(connection)e(b)r(et)n(w)n(een)h
(HBMG)h(and)e(ILU)h(from)g(a)g(graph)300 2474 y(theoretical)24
b(p)r(oin)n(t)g(of)h(view.)f(In)h(Section)f(3,)h(w)n(e)f(analyze)f
(this)i(metho)r(d)g(in)g(the)g(one-dimensional)300 2574
y(case,)31 b(where)g(algebraic)e(simplicit)n(y)j(allo)n(ws)e(a)h
(fairly)g(complete)g(treatmen)n(t.)h(In)f(Section)h(4,)f(w)n(e)300
2673 y(presen)n(t)c(some)g(algorithms)f(and)i(n)n(umerical)e(results)h
(for)g(t)n(w)n(o-dimensional)f(problems.)300 2980 y FV(2)105
b(HBMG)35 b(and)g(ILU)300 3180 y F9(As)h(with)h(t)n(ypical)f(m)n
(ultigrid)g(metho)r(ds,)g(classical)f(hierarc)n(hical)f(basis)h(metho)r
(ds)i(are)e(usually)300 3279 y(de\014ned)21 b(in)h(terms)f(of)g(an)g
(underlying)f(re\014nemen)n(t)h(structure)g(for)f(a)h(sequence)f(of)h
(nested)h(meshes.)300 3379 y(In)29 b(man)n(y)g(cases,)f(this)h(is)g(no)
g(disadv)-5 b(an)n(tage,)27 b(but)j(it)f(limits)h(the)f(applicabilit)n
(y)g(of)g(the)g(metho)r(ds)300 3479 y(to)22 b(truly)h(unstructured)f
(meshes,)h(whic)n(h)f(migh)n(t)h(b)r(e)g(highly)f(non)n(uniform)g(but)h
(not)g(deriv)n(ed)f(from)300 3578 y(some)j(grid)g(re\014nemen)n(t)h
(pro)r(cess.)e(Here)i(w)n(e)f(view)g(the)h(transformation)e(of)i(the)g
(sti\013ness)g(matrix)300 3678 y F6(A)e F9(from)f(no)r(dal)f(basis)h
(represen)n(tation)f(to)h(hierarc)n(hical)e(basis)h(represen)n(tation)g
(as)h(a)f(sp)r(ecial)h(ILU)300 3778 y(decomp)r(osition.)h(This)g
(generalizes)e(the)i(construction)g(of)g(hierarc)n(hical)e(bases)h(to)h
(unstructured)300 3877 y(meshes,)k(allo)n(wing)f(HBMG)h(and)g(other)g
(hierarc)n(hical)e(basis)h(metho)r(ds)i(to)f(b)r(e)h(applied.)f(A)g
(more)300 3977 y(complete)38 b(discussion)f(of)g(this)i(p)r(oin)n(t)e
(can)h(b)r(e)g(found)g(in)g([BX94)o(])g(and)g([BX96)o(].)g(See)g([HS95)
o(],)300 4076 y([Kor96)n(],)24 b([HK94)o(],)h([KY94)o(],)g(and)f([CS94)
o(])h(for)f(related)g(algorithms.)f(See)i([Xu89)o(],)g([BPX91)o(],)g
(and)300 4176 y([Zha88)n(])j(for)f(discussions)g(of)g(non-nested)h(m)n
(ultigrid)f(algorithms.)300 4425 y F0(Gr)l(aph)k(The)l(or)l(etic)l(al)h
(Pr)l(op)l(erties)e(of)h(Hier)l(ar)l(chic)l(al)h(Bases)300
4575 y F9(W)-7 b(e)58 b(b)r(egin)f(b)n(y)g(exploring)e(the)j
(connection)f(b)r(et)n(w)n(een)g(the)g(HBMG)h(metho)r(d)f(and)g(ILU)300
4674 y(decomp)r(osition)i(in)h(terms)g(of)f(graph)g(theory)-7
b(.)59 b(W)-7 b(e)60 b(consider)f(the)h(standard)e(Gaussian)300
4774 y(elimination)44 b(and)g(the)h(classical)d(ILU)j(factorization)d
(from)i(a)g(graph)f(theoretical)g(p)r(oin)n(t)i(of)300
4873 y(view,)38 b(and)g(then)h(dev)n(elop)e(a)h(simple)g(graph)f
(elimination)h(mo)r(del)g(for)g(classical)f(hierarc)n(hical)300
4973 y(basis)43 b(metho)r(ds)h(on)f(sequences)g(of)g(nested)h(meshes.)f
(W)-7 b(e)44 b(can)f(in)n(terpret)g(this)h(mo)r(del)g(as)f(a)300
5073 y(particular)e(ILU)i(decomp)r(osition)e(and)h(generalize)f(this)i
(graph)e(elimination)h(mo)r(del)h(to)f(the)p eop
%%Page: 253 3
253 272 bop 300 282 a Gd(HIERAR)n(CHICAL)23 b(BASIS)i(F)n(OR)g(UNSTR)n
(UCTURED)f(MESHES)894 b(253)300 482 y F9(case)48 b(of)h(completely)g
(unstructured)g(meshes.)f(See)h(Rose)g([Ros72)n(])g(and)g(George)f(and)
g(Liu)300 581 y([GL81)o(])25 b(for)f(a)h(complete)f(discussion)g(of)h
(graph)f(theoretical)g(asp)r(ects)g(of)h(Gaussian)e(elimination.)300
681 y(Corresp)r(onding)44 b(to)i(a)g(sparse)f F6(n)30
b FU(\002)g F6(n)46 b F9(matrix)g F6(A)g F9(with)h(symmetric)f(pattern)
g(\(i.e.)g F6(A) y FS(ij) y FU(6)p F9(=)300
780 y(0)27 b(if)h(and)g(only)f(if)h F6(A)927 792 y FS(j)s(i)1009
780 y FU(6)p F9(=)23 b(0\),)k(let)h F6(G)p F9(\()p F6(X)r(;)14
b(E)5 b F9(\))29 b(b)r(e)f(the)g(graph)e(that)i(consists)f(of)h(a)f
(set)h(of)f F6(n)h F9(ordered)300 880 y(v)n(ertices)20
b F6(v)634 892 y FS(i)685 880 y FU(2)j F6(X)7 b F9(,)21
b(1)i FU(\024)f F6(i)h FU(\024)g F6(n)p F9(,)e(and)g(a)g(set)g(of)g
(edges)g F6(E)26 b F9(suc)n(h)21 b(that)g(the)h(edge)f(\(connecting)g
(v)n(ertices)300 980 y F6(v)340 992 y FS(i)396 980 y
F9(and)28 b F6(v)598 992 y FS(j)634 980 y F9(\))g F6(e)733
992 y FS(ij)816 980 y FU(2)c F6(E)34 b F9(if)28 b(and)g(only)g(if)h
F6(a) y FS(ij) y FU(6)p F9(=)23 b(0,)28
b F6(i)c FU(6)p F9(=)g F6(j)5 b F9(.)28 b(The)h(edges)e(in)i(the)f
(graph)g F6(G)g F9(corresp)r(ond)300 1079 y(to)k(the)g(nonzero)e(o\013)
i(diagonal)e(en)n(tries)h(of)h F6(A)p F9(.)g(If)g F6(A)g
F9(is)g(the)g(sti\013ness)g(matrix)f(for)g(the)i(space)e(of)300
1179 y(con)n(tin)n(uous)e(piecewise)g(linear)g(p)r(olynomials)g
(represen)n(ted)g(in)h(the)g(standard)f(no)r(dal)g(basis,)g(the)300
1279 y(graph)c F6(G)h F9(is)g(just)h(the)f(underlying)g(triangulation)f
(of)h(the)g(domain)g(\(with)g(minor)g(mo)r(di\014cations)300
1378 y(due)32 b(to)g(Diric)n(hlet)h(b)r(oundary)e(conditions\).)h(W)-7
b(e)32 b(de\014ne)h(for)e(a)h(v)n(ertex)f F6(v)
y FS(i) y F9(the)i(set)f(of)g(adjacen)n(t)300
1478 y(v)n(ertices)26 b F6(ad)-9 b(j)5 b F9(\()p F6(v)789
1490 y FS(i)817 1478 y F9(\))28 b(b)n(y) y F6(ad)-9
b(j)5 b F9(\()p F6(v) y FS(i) y F9(\))24
b(=)e FU(f)p F6(v) y FS(j) y FU(2)h
F6(X)7 b FU(j)p F6(e) y FS(ij) y FU(2)23
b F6(E)5 b FU(g)p F6(:)300 1727 y F9(A)48 b(clique)f
F6(C)62 b FU(\022)56 b F6(X)e F9(is)47 b(a)g(set)h(of)f(v)n(ertices)f
(whic)n(h)i(are)e(all)h(pairwise)g()g(that)h(is)300
1826 y F6(v)340 1838 y FS(i)368 1826 y F6(;)14 b(v)445
1838 y FS(j)531 1826 y FU(2)51 b F6(C)q(;)14 b(i)50 b
FU(6)p F9(=)h F6(j)k FU(\))c F6(e) y FS(ij)1299
1826 y FU(2)g F6(E)5 b F9(.)45 b(With)g(a)e(prop)r(er)h(ordering)e(of)i
(the)h(v)n(ertices,)e(a)h(clique)300 1926 y(corresp)r(onds)36
b(to)j(a)f(dense)g(submatrix)g(of)g F6(A)p F9(.)h(In)f(graph)f
(theoretic)h(terms,)g(a)g(single)g(step)h(of)300 2026
y(Gaussian)27 b(elimination)g(transforms)f F6(G)p F9(\()p
F6(X)r(;)14 b(E)5 b F9(\))29 b(to)e(a)g(new)h(graph)e
F6(G) y FR(0) y F9(\()p F6(X)
y FR(0) y F6(;)14 b(E) y FR(0)2758
2026 y F9(\))28 b(as)f(follo)n(ws:)300 2192 y(1.)41 b(Eliminate)33
b(v)n(ertex)f F6(v) y FS(i) y F9(and)h(all)g(its)g
(inciden)n(t)h(edges)e(from)h F6(G)p F9(.)g(Set)h F6(X)2655
2162 y FR(0) y F9(=)e F6(X)d FU(\000)21 b(f)p
F6(v) y FS(i) y FU(g)p F9(.)406 2291
y(Denote)28 b(the)g(resulting)f(set)g(of)h(edges)f F6(E)1683
2303 y FT(1) y FU(\022)c F6(E)5 b F9(.)300 2391
y(2.)41 b(Create)23 b(a)g(set)g F6(F)35 b F9(of)24 b(\014llin)g(edges)e
(as)h(follo)n(ws:)f(F)-7 b(or)23 b(eac)n(h)g(distinct)h(pair)e
F6(v) y FS(j) y F6(;)14 b(v)
y FS(k) y FU(2)23 b F6(ad)-9 b(j)5 b F9(\()p
F6(v) y FS(i) y F9(\))406 2491 y(in)28
b F6(G)p F9(,)g(add)f(the)h(edge)f F6(e) y FS(j)s(k)1250
2491 y F9(to)h F6(F)40 b F9(if)28 b(not)f(already)f(presen)n(t)h(in)h
F6(E) y FT(1) y F9(.)f(Set)h F6(E)2705
2461 y FR(0) y F9(=)22 b F6(E) y FT(1)2956
2491 y FU([)d F6(F)12 b F9(.)383 2657 y(Note)31 b(that)f(the)h(set)g
F6(ad)-9 b(j)5 b F9(\()p F6(v)s F9(\))31 b(in)g F6(G)g
F9(b)r(ecomes)f(a)g(clique)h(in)g F6(G) y FR(0)2333
2657 y F9(.)g(Within)h(this)e(framew)n(ork,)f(the)300
2756 y(classical)d(ILU)i(factorization)e(is)h(one)g(in)g(whic)n(h)h
F0(no)f F9(\014llin)h(edges)f(are)f(allo)n(w)n(ed,)g(i.e.)i
F6(F)35 b FU(\021)22 b(;)p F9(.)28 b(This)300 2856 y(forces)34
b(the)h(matrix)f F6(A) y FR(0) y F9(corresp)r(onding)
f(to)i(the)g(new)g(graph)e F6(G) y FR(0)2433
2856 y F9(to)i(ha)n(v)n(e)e(the)j(same)e(sparsit)n(y)300
2956 y(structure)27 b(as)g(the)h(corresp)r(onding)e(submatrix)h(of)g
F6(A)p F9(.)383 3055 y(T)-7 b(o)18 b(de\014ne)h(HBMG)g(as)f(a)h
(generalized)e(ILU)i(pro)r(cedure,)f(w)n(e)g(m)n(ust)h(\014rst)f(in)n
(tro)r(duce)h(the)g(concept)300 3155 y(of)34 b F0(vertex)i(p)l(ar)l
(ents.)f F9(W)-7 b(e)35 b(will)f(b)r(egin)h(with)g(the)f(case)g(of)g(t)
n(w)n(o)g(nested)h(meshes)f(where)f(the)i(\014ne)300
3254 y(mesh)28 b(is)g(a)g(uniform)g(re\014nemen)n(t)g(of)g(a)f(coarse)g
(mesh,)h(generated)f(b)n(y)g(pairwise)h(connecting)f(the)300
3354 y(midp)r(oin)n(ts)j(of)g(the)h(coarse)d(grid)h(edges)h(in)g(the)g
(usual)g(w)n(a)n(y)f([BD)n(Y88)o(],)i([Yse86)o(],)f([Hac85)o(].)g(Here)
300 3454 y(w)n(e)36 b(can)f(mak)n(e)h(the)g(direct)g(sum)g(decomp)r
(osition)g F6(X)43 b F9(=)37 b F6(X) y FS(c)2293
3454 y FU(\010)24 b F6(X) y FS(f) y
F9(,)37 b(where)e F6(X) y FS(c) y F9(is)h(the)g(set)h
(of)300 3553 y(coarse)26 b(grid)h(v)n(ertices)f(and)i
F6(X) y FS(f) y F9(is)f(the)h(set)g(of)g(\014ne)f
(grid)g(v)n(ertices)g(\(those)g(not)h(in)g F6(X)2964
3565 y FS(c) y F9(\).)g(F)-7 b(or)27 b(eac)n(h)300
3653 y(v)n(ertex)e F6(v)590 3665 y FS(i)641 3653 y FU(2)e
F6(X)788 3665 y FS(f)831 3653 y F9(,)j(there)g(is)f(a)h(unique)g(pair)f
(of)g(v)n(ertex)g(paren)n(ts)g F6(v) y FS(j)2386
3653 y F6(;)14 b(v) y FS(k) y FU(2)24
b F6(X) y FS(c) y F9(suc)n(h)i(that)g
F6(v) y FS(i) y F9(is)g(the)300 3753
y(midp)r(oin)n(t)i(of)g(the)g(edge)f(connecting)g F6(v)1537
3765 y FS(j) y F9(and)g F6(v) y FS(k)1870
3753 y F9(\()p F6(v) y FS(i) y F9(=)c(\()p
F6(v) y FS(j) y F9(+)18 b F6(v)2330
3765 y FS(k) y F9(\))p F6(=)p F9(2\).)383 3852
y(W)-7 b(e)35 b(no)n(w)e(view)h(HBMG)h(as)e(an)h(ILU)h(algorithm)e(in)h
(whic)n(h)g(only)g(selected)g(\014llin)h(edges)e(are)300
3952 y(allo)n(w)n(ed,)g(namely)g(those)h(connecting)g(v)n(ertex)f
(paren)n(ts.)g(In)h(this)h(algorithm,)e(w)n(e)g(sequen)n(tially)300
4051 y(eliminate)28 b(the)g(v)n(ertices)e(in)i(the)g(set)g
F6(X) y FS(f) y F9(as)f(follo)n(ws:)300
4218 y(1.)41 b(Eliminate)19 b(v)n(ertex)f F6(v)
y FS(i) y FU(2)23 b F6(X) y FS(f)1320
4218 y F9(and)c(all)g(its)h(inciden)n(t)f(edges)f(from)h
F6(G)p F9(.)h(Set)f F6(X) y FR(0) y
F9(=)j F6(X)8 b FU(\000)q(f)p F6(v) y FS(i)3100
4218 y FU(g)p F9(.)406 4317 y(Denote)28 b(the)g(resulting)f(set)g(of)h
(edges)f F6(E) y FT(1) y FU(\022)c
F6(E)5 b F9(.)300 4417 y(2.)41 b(Add)g(one)g(\014llin)g(edge)f
(connecting)h(the)g(v)n(ertex)f(paren)n(ts)f F6(v)2367
4429 y FS(j) y F6(;)14 b(v) y FS(k)2566
4417 y FU(2)45 b F6(X) y FS(c) y F9(of)c
F6(v) y FS(i) y F9(.)g(Set)406 4516
y F6(E)472 4486 y FR(0)519 4516 y F9(=)22 b F6(E)667
4528 y FT(1)723 4516 y FU([)d(f)p F6(e)878 4528 y FS(j)s(k)949
4516 y FU(g)p F9(.)300 4682 y(Note)49 b(that)g(the)g(triangulation)e
FU(T) y FS(f) y F9(is)i(the)g(graph)e(for)h(the)h
(original)e(sti\013ness)i(matrix)f F6(A)300 4782 y F9(represen)n(ted)31
b(in)g(the)i(standard)d(no)r(dal)i(basis.)f(After)h(all)f(the)h(v)n
(ertices)f(in)h F6(X) y FS(f) y F9(are)f(eliminated,)
300 4882 y(the)d(resulting)g(graph)f(is)h(just)h(the)f(triangulation)f
FU(T) y FS(c) y F9(,)h(i.e.)h(the)f(sparsit)n(y)f
(pattern)h(of)g(the)g(coarse)300 4981 y(grid)c(matrix)g(corresp)r(onds)
f(to)i(the)g(coarse)e(grid)h(triangulation.)g(F)-7 b(or)24
b(completely)h(unstructured)300 5081 y(meshes,)31 b(the)h(main)g
(problem)f(is)g(to)h(determine)f(reasonable)f(v)n(ertex)g(paren)n(ts)h
(for)g(eac)n(h)f(v)n(ertex)p eop
%%Page: 254 4
254 273 bop 300 282 a Gd(254)2305 b(BANK)25 b(&)g(GUTSCH)300
482 y F9(to)j(b)r(e)g(eliminated.)g(Once)f(this)h(is)f(done,)h(the)g
(elimination/unre\014nemen)n(t/coarsening)c(is)k(done)300
581 y(exactly)18 b(as)g(in)h(the)h(case)e(of)g(nested)h(meshes.)g(This)
g(ma)n(y)f(lead)g(to)h(graphs)e(that)i(are)f(not)h(necessarily)300
681 y(triangulations)41 b(of)h(the)h(domain,)f(but)h(t)n(ypically)f
(con)n(tain)g(p)r(olygonal)f(elemen)n(ts)h(of)h(v)-5
b(arious)300 780 y(orders.)32 b(Ev)n(en)i(if)g(the)g(graphs)f(remain)g
(triangulations,)g(they)h(will)g(generally)f(not)h(b)r(e)g(nested.)300
880 y(Nonetheless,)40 b(suc)n(h)h(a)f(sc)n(heme)g(still)h(de\014nes)g
(the)g(linear)f(com)n(binations)g(of)g(\014ne)h(grid)f(basis)300
980 y(functions)e(used)f(to)g(create)g(the)h(coarse)d(grid)i(basis)g
(functions)g(\(but)i(not)e(the)h(v)-5 b(alues)37 b(of)g(the)300
1079 y(co)r(e\016cien)n(ts)27 b(in)h(the)g(linear)f(com)n(binations\).)
383 1179 y(Algorithms)38 b(for)g(selecting)h(v)n(ertex)f(paren)n(ts)f
(are)h(curren)n(tly)g(an)g(area)g(of)h(activ)n(e)f(researc)n(h.)300
1279 y(The)31 b(sc)n(heme)g(dev)n(elop)r(ed)f(in)i([BX96)o(])f(is)g
(based)g(on)g(the)g(geometry)f(of)h(the)h(triangulation,)e(and)300
1378 y(seeks)e(to)g(coarsen)f(the)i(grid)f(in)h(a)f(fashion)g(that)h
(main)n(tains)f(the)h(shap)r(e)f(of)h(the)f(region)g(and)g(the)300
1478 y(shap)r(e)k(regularit)n(y)e(of)h(the)h(coarse)e(grid)i(elemen)n
(ts)f(to)h(the)g(exten)n(t)g(p)r(ossible.)g(F)-7 b(or)31
b(this)h(sc)n(heme,)300 1577 y(the)d(supp)r(orts)f(of)h(the)g
(resulting)f(hierarc)n(hical)e(basis)i(functions)h(gro)n(w)e(in)i(a)f
(fashion)g(analogous)300 1677 y(to)g(the)h(classical)e(case.)h(W)-7
b(e)29 b(are)f(also)f(considering)g(more)h(algebraic)f(sc)n(hemes)h
(whic)n(h)g(rely)g(only)300 1777 y(on)c(the)h(sparsit)n(y)f(structure)g
(of)g(the)h(sti\013ness)g(matrix)f(and)h(the)g(n)n(umerical)e(v)-5
b(alues)25 b(of)f(its)h(matrix)300 1876 y(elemen)n()41
b(these)g(sc)n(hemes)f(ha)n(v)n(e)g(m)n(uc)n(h)g(in)h(common)g(with)g
(more)f(classical)g(sparse)f(matrix)300 1976 y(ordering)26
b(algorithms.)300 2225 y F0(A)n(lgebr)l(aic)31 b(HBMG)f(and)h(ILU)300
2374 y F9(In)21 b(this)h(section,)f(w)n(e)f(consider)h(the)g(algebraic)
e(asp)r(ects)i(of)g(the)h(HBMG)f(metho)r(d)h(and)f(its)g(relation)300
2474 y(to)44 b(Gaussian)f(elimination.)g(Again)h(w)n(e)f(will)h
(consider)f(the)h(case)f(of)h(only)f(t)n(w)n(o)g(lev)n(els.)g(Let)300
2574 y F6(A)37 b F9(denote)g(the)h(no)r(dal)f(basis)f(sti\013ness)h
(matrix)g(for)f(the)i(\014ne)f(grid,)g(and)f(consider)h(the)g(blo)r(c)n
(k)300 2673 y(partitioning) y F6(A)23 b F9(=)1629
2787 y FQ(\022) y F6(A) y FS(cc)1954
2853 y F6(A) y FS(cf) y F6(A)
y FS(f)7 b(c) y F6(A) y FS(f)g(f)2135
2787 y FQ(\023) y F6(;)1049 b F9(\(1\))300 3135
y(where)45 b F6(A)620 3147 y FS(f)7 b(f)749 3135 y F9(corresp)r(onds)44
b(to)i(the)g(no)r(dal)g(basis)f(functions)h(of)g(the)g(\014ne)h(grid)e
(no)r(des,)h F6(A) y FS(cc)300 3235 y F9(corresp)r(onds)40
b(to)h(the)h(\(\014ne)g(grid\))f(no)r(dal)g(basis)g(functions)h(of)f
(the)h(coarse)e(grid)g(no)r(des)i(and)300 3334 y F6(A)362
3346 y FS(cf)476 3334 y F9(and)e F6(A)712 3346 y FS(f)7
b(c)825 3334 y F9(corresp)r(ond)39 b(to)i(the)g(coupling)f(b)r(et)n(w)n
(een)g(the)h(t)n(w)n(o)f(sets)g(of)g(basis)g(functions.)300
3443 y(W)-7 b(e)31 b(consider)e(transformations)g(of)h(the)h(form)f
F6(A) y FR(0) y F9(=)d F6(S)
y FS(T) y F6(A) y F9(~) y
F6(S)36 b F9(where)30 b F6(S)35 b F9(and) y(~)2771
3443 y F6(S)g F9(ha)n(v)n(e)29 b(the)i(blo)r(c)n(k)300
3542 y(structure) y F6(S)c F9(=) y
FQ(\022) y F6(I)101 b F9(0) y F6(R)83
b(I) y FQ(\023) y F6(;) y
F9(~) y F6(S)28 b F9(=) y FQ(\022)2186
3706 y F6(I)101 b F9(0) y(~) y F6(R)84
b(I) y FQ(\023) y F6(:)777 b F9(\(2\))300
3995 y(By)27 b(direct)h(calculation,)f(w)n(e)g(obtain)1061
4233 y F6(S) y FS(T) y F6(A)
y F9(~) y F6(S)h F9(=) y FQ(\022)1682
4163 y F9(^) y F6(A) y FS(cc)
y F6(A) y FS(cf) y F9(+)18 b F6(R)2328
4153 y FS(T) y F6(A) y FS(f)7 b(f)1501
4290 y F6(A) y FS(f)g(c) y F9(+)18
b F6(A) y FS(f)7 b(f) y F9(~)
y F6(R)260 b(A) y FS(f)7 b(f) y FQ(\023)3282
4233 y F9(\(3\))300 4472 y(where) y(^)
y F6(A) y FS(cc) y F9(=)23 b F6(A)1411
4669 y FS(cc) y F9(+)18 b F6(R) y FS(T)1692
4657 y F6(A) y FS(f)7 b(c) y F9(+)18
b F6(A) y FS(cf) y F9(~) y
F6(R)h F9(+)f F6(R) y FS(T) y F6(A)2407
4669 y FS(f)7 b(f) y F9(~) y F6(R)q(:)706
b F9(\(4\))300 4858 y(Di\013eren)n(t)21 b(algorithms)e(can)h(b)r(e)h(c)
n(haracterized)d(b)n(y)j(di\013eren)n(t)f(c)n(hoices)g(of)g
F6(R)h F9(and) y(~) y F6(R)p F9(.)g(F)-7
b(or)20 b(example,)300 4958 y(in)h(the)g(classical)e(blo)r(c)n(k)h
(Gaussian)f(elimination)h(w)n(e)h(ha)n(v)n(e)e F6(R)24
b F9(=)e FU(\000)p F6(A) y FR(\000)p FS(T)2426
4983 y(f)7 b(f) y F6(A) y FS(T)2592
4981 y(cf) y F9(and) y(~)
y F6(R)24 b F9(=)e FU(\000)p F6(A) y FR(\000)p
FT(1) y FS(f)7 b(f) y F6(A)
y FS(f)g(c) y F9(,)300 5081 y(and)486 5060 y(^)464
5081 y F6(A)526 5093 y FS(cc)617 5081 y F9(=)27 b F6(A)771
5093 y FS(cc)855 5081 y FU(\000)20 b F6(A) y
FS(cf) y F6(A) y FR(\000)p FT(1)1137
5106 y FS(f)7 b(f) y F6(A) y FS(f)g(c)1392
5081 y F9(is)30 b(the)g(Sc)n(h)n(ur)g(complemen)n(t.)g(In)g(this)h
(case,)e(the)i(o\013)f(diagonal)p eop
%%Page: 255 5
255 274 bop 300 282 a Gd(HIERAR)n(CHICAL)23 b(BASIS)i(F)n(OR)g(UNSTR)n
(UCTURED)f(MESHES)894 b(255)300 482 y F9(blo)r(c)n(ks)30
b(are)f(reduced)h(to)g(zero,)f(but)i(at)f(the)h(cost)f(of)g(ha)n(ving)f
(fairly)h(dense)g(matrices)g F6(R)q F9(,) y(~)3161
482 y F6(R)h F9(and)322 560 y(^)300 581 y F6(A)362 593
y FS(cc)426 581 y F9(.)c(In)h(the)g(case)e(of)h(HBMG,)h(the)g(matrices)
f F6(R)h F9(and) y(~) y F6(R)g F9(are)e(sparse)g(and)i
(con)n(tain)e(information)300 681 y(ab)r(out)39 b(c)n(hanging)f(from)i
(the)f(no)r(dal)h(to)f(the)h(hierarc)n(hical)d(basis.)i(The)g(sparsit)n
(y)g(patterns)g(of)300 780 y F6(R)d F9(and)586 759 y(~)568
780 y F6(R)g F9(are)e(the)i(same,)e(and)h(b)r(oth)h(are)e(determined)h
(b)n(y)g(the)h(v)n(ertex)e(paren)n(t)g(relationship)300
880 y(describ)r(ed)26 b(ab)r(o)n(v)n(e.)g(In)g(particular,)g(eac)n(h)f
(ro)n(w)h(of)g F6(R)h F9(and) y(~) y
F6(R)g F9(is)g(zero)e(except)h(for)h(the)f(t)n(w)n(o)g(en)n(tries)300
980 y(whic)n(h)h(corresp)r(ond)e(to)h(the)h(\(coarse)f(grid\))g(v)n
(ertex)g(paren)n(ts)f(for)i(the)g(giv)n(en)f(\014ne)h(grid)f(v)n
(ertex.)f(In)300 1079 y(the)30 b(classical)e(case,)g(where)h(eac)n(h)g
(\014ne)h(grid)e(v)n(ertex)h(is)g(the)h(midp)r(oin)n(t)g(of)f(the)h
(edge)f(connecting)300 1179 y(its)i(v)n(ertex)f(paren)n(ts,)f
F6(R)g F9(=) y(~) y F6(R)j F9(and)e(b)r(oth)h
(nonzero)f(en)n(tries)g(in)h(a)f(giv)n(en)g(ro)n(w)f(are)h(equal)g(to)h
(1)p F6(=)p F9(2.)300 1279 y(In)26 b(the)f(generalized)f(HBMG,)i(these)
f(v)-5 b(alues)25 b(are)g(replaced)f(b)n(y)h F6(\022)2349
1291 y FS(i) y F6(;) y F9(~)
y F6(\022) y FS(i) y F6(;)14 b(\027)2558
1291 y FS(i) y F9(and)30 b(~)-47 b F6(\027)2811
1291 y FS(i) y F9(.)25 b(Often)h(one)f(has)307
1356 y(~)300 1378 y F6(\022)339 1390 y FS(i)399 1378
y F9(=)31 b(1)22 b FU(\000)g F6(\022)685 1390 y FS(i)745
1378 y F9(and)38 b(~)-47 b F6(\027)953 1390 y FS(i)1013
1378 y F9(=)32 b(1)21 b FU(\000)h F6(\027) y
FS(i) y F9(.)34 b(Cho)r(osing)e F6(R)g FU(6)p
F9(=) y(~) y F6(R)i F9(corresp)r(onds)d(to)i(c)n(ho)r
(osing)e(a)i(test)g(space)300 1478 y(di\013eren)n(t)g(from)g(the)g
(trial)f(space)g(for)h(the)g(coarser)d(grids,)i(as)h(in)g(a)f(P)n(etro)
n(v-Galerkin)e(metho)r(d.)300 1577 y(Sev)n(eral)c(alternativ)n(es)g
(for)h(c)n(ho)r(osing)g F6(\022) y FS(i)1565
1577 y F6(;) y F9(~) y F6(\022)1641
1589 y FS(i) y F6(;)14 b(\027) y FS(i)1802
1577 y F9(and)32 b(~)-47 b F6(\027) y FS(i)2060
1577 y F9(are)26 b(discussed)h(in)h(Section)g(4.)300
1885 y FV(3)105 b(Analysis)35 b(of)g(a)g(One-Dimensional)f(Mo)s(del)h
(Problem)300 2084 y F0(The)c(Two-L)l(evel)g(HBMG)f(Iter)l(ation)300
2233 y F9(In)46 b(this)g(section,)g(w)n(e)g(will)g(analyze)f(the)h
(case)f(of)h(a)f(constan)n(t)h(co)r(e\016cien)n(t)f(t)n(w)n(o)g(t)n(w)n
(o-p)r(oin)n(t)300 2333 y(b)r(oundary)d(v)-5 b(alue)44
b(problem,)e(giving)h(rise)f(to)h(a)g(constan)n(t)g(co)r(e\016cien)n(t)
g(tridiagonal)e(sti\013ness)300 2433 y(matrix)26 b(when)h(discretized)f
(using)g(some)g(\014nite)h(elemen)n(t)g(appro)n(ximation)d(on)j(a)f
(uniform)g(mesh.)300 2532 y(See)j([Hac84)o(])g(and)h([BB91)n(])g(for)f
(other)f(analyses)g(of)h(m)n(ultilev)n(el)g(metho)r(ds)h(for)f
(one-dimensional)300 2632 y(mo)r(del)g(problems.)f(Let)g
F6(n)d(>)f F9(2)k(b)r(e)h(an)f(in)n(teger)g(and)g(set)g
F6(h)d F9(=)f(1)p F6(=)p F9(\(2)p F6(n)p F9(\).)k(The)g(uniform)h
(\014ne)f(mesh)300 2731 y FU(T)345 2743 y FS(h)418 2731
y F9(of)i(size)g F6(h)g F9(has)f(2)p F6(n)19 b F9(+)h(1)30
b(grid)f(p)r(oin)n(ts)h F6(x) y FS(k)
y F9(=)d F6(k)s(h;)14 b F9(0)26 b FU(\024)g F6(k)k FU(\024)d
F9(2)p F6(n)p F9(.)i(The)h(coarse)e(mesh)i(has)f F6(n)20
b F9(+)g(1)300 2831 y(grid)33 b(p)r(oin)n(ts)h F6(x)780
2843 y FT(2)p FS(k)854 2831 y F6(;)14 b F9(0)34 b FU(\024)f
F6(k)j FU(\024)d F6(n)p F9(.)h(W)-7 b(e)35 b(will)f(refer)f(to)h(the)g
(set)g(of)g(coarse)e(grid)h(p)r(oin)n(ts)h(as)f(lev)n(el)h(1)300
2931 y(v)n(ertices,)27 b(and)g F6(x)832 2943 y FT(2)p
FS(k)q FT(+1)990 2931 y F6(;)14 b F9(0)23 b FU(\024)g
F6(k)i FU(\024)e F6(n)18 b FU(\000)g F9(1,)28 b(as)f(lev)n(el)g(2)g(v)n
(ertices.)383 3030 y(Let)36 b FU(P)598 3042 y FS(h)676
3030 y F9(b)r(e)g(a)g(2)p F6(n)23 b FU(\000)g F9(1)36
b(dimensional)f(trial)g(space)g(of)h(functions)g(asso)r(ciated)e(with)i
(the)g(\014ne)300 3130 y(mesh)d FU(T)563 3142 y FS(h)639
3130 y F9(satisfying)e(the)i(b)r(oundary)f(conditions)g
F6(v)s F9(\()p F6(x) y FT(0) y F9(\))f(=)g
F6(v)s F9(\()p F6(x) y FT(2)p FS(n)
y F9(\))h(=)f(0)h(for)g(all)h F6(v)h FU(2)e(P)
y FS(h) y F9(.)g(Let)311 3217 y(^)300 3239 y
F6(\036)349 3251 y FS(k)390 3239 y F6(;)14 b F9(1)33
b FU(\024)f F6(k)k FU(\024)d F9(2)p F6(n)21 b FU(\000)h
F9(1,)34 b(denote)f(the)h(no)r(dal)f(basis)g(for)g(the)h(trial)f(space)
g FU(P) y FS(h) y F9(.)g(W)-7 b(e)34
b(assume)f(that)300 3338 y F6(suppor)r(t)p FU(f)633 3316
y F9(^)622 3338 y F6(\036)671 3350 y FS(k)712 3338 y
FU(g)f F9(=)g(\()p F6(x)962 3350 y FS(k)q FR(\000)p FT(1)1089
3338 y F6(;)14 b(x) y FS(k)q FT(+1)
y F9(\))34 b(and)f(that) y(^) y F6(\036)1766
3350 y FS(k) y F9(\()p F6(x) y FS(j)1922
3338 y F9(\))g(=)f F6(\016) y FS(k)q(j)2193
3338 y F9(.)h(Similarly)-7 b(,)33 b(w)n(e)g(de\014ne)h(the)f(2)p
F6(n)22 b FU(\000)g F9(1)300 3450 y(dimensional)27 b(test)h(space)f
FU(S) y FS(h) y F9(with)h(no)r(dal)f(basis)1900
3428 y(^) y F6( ) y FS(k)
y F6(;)14 b F9(1)22 b FU(\024)h F6(k)j FU(\024)d F9(2)p
F6(n)17 b FU(\000)h F9(1.)383 3549 y(W)-7 b(e)19 b(will)g(use)g(a)f
(discretization)g(on)h(the)g(mesh)g FU(T) y
FS(h) y F9(giv)n(en)f(b)n(y)g(a)h(bilinear)f(form)h
F6(b)p F9(\()p FU(\001)p F6(;)14 b FU(\001)p F9(\))23
b(:)g FU(P) y FS(h) y FU(\002)q(S)3240
3561 y FS(h) y FU(!)300 3649 y F9(I)-18 b(R.)30
b(The)g(details)f(of)h(the)h(bilinear)e(form)h F6(b)p
F9(\()p FU(\001)p F6(;)14 b FU(\001)p F9(\))30 b(are)f(arbitrary)e(for)
j(the)g(momen)n(t.)g(The)g(discrete)300 3749 y(system)d(of)h(equations)
f(to)g(b)r(e)h(solv)n(ed)f(is:)g(Find)h F6(u)
y FS(h) y FU(2)23 b(P) y FS(h)2130
3749 y F9(suc)n(h)k(that) y F6(b)p F9(\()p F6(u)1652
3972 y FS(h) y F6(;)14 b(v)s F9(\))23 b(=)g
F6(r)r(hs)p F9(\()p F6(v)s F9(\))300 4170 y(for)k(all)g
F6(v)g FU(2)c(S)737 4182 y FS(h)780 4170 y F9(,)28 b(where)f
F6(r)r(hs)p F9(\()p FU(\001)p F9(\))i(is)e(an)h(appropriate)e(linear)g
(functional.)383 4270 y(Supp)r(ose) y F6(b)p
F9(\() y(^) y F6(\036) y FS(k)1450
4626 y F6(;) y F9(^) y F6( )
y FS(j) y F9(\))d(=) y FQ(8)
y(>) y(>) y()
y(>) y(:) y F6(a)83 b(j)28 b F9(=)23
b F6(k)e FU(\000)d F9(1) y F6(b)87 b(j)28 b
F9(=)23 b F6(k)e F9(+)d(1) y F6(c)87 b(j)28
b F9(=)23 b F6(k) y F9(0)84 b FU(j)p F6(k)21
b FU(\000)d F6(j)5 b FU(j)23 b F6(>)g F9(1)
y(\(5\))300 4981 y(where)k F6(a)p F9(,)h F6(b)p F9(,)f(and)g
F6(c)h F9(are)f(constan)n(ts.)383 5081 y(The)35 b(resulting)f(no)r(dal)
h(basis)f(sti\013ness)h(matrix)f(in)i(natural)e(v)n(ertex)g(order)f(is)
i(the)g(constan)n(t)p eop
%%Page: 256 6
256 275 bop 300 282 a Gd(256)2305 b(BANK)25 b(&)g(GUTSCH)300
482 y F9(co)r(e\016cien)n(t)i(tridiagonal)f(matrix)1216
890 y F6(A) y FS(N)6 b(B) y F9(=)1505
599 y FQ(0) y(B) y(B) y(B)1505
894 y(B) y(B) y(@) y F6(c)118
b(b) y(a)113 b(c)144 b(b) y F9(.)1783
884 y(.) y(.) y(.) y(.)
y(.) y(.) y(.) y(.)
y F6(a)c(c)114 b(b) y(a)c(c) y FQ(1)2363
745 y(C) y(C) y(C) y(C)
y(C) y(A) y F6(:)810 b F9(\(6\))383 1299
y(No)n(w)20 b(w)n(e)f(in)n(tro)r(duce)h(the)g(generalized)f(hierarc)n
(hical)f(bases)h(for)g(the)i(spaces)e FU(P)
y FS(h) y F9(and)h FU(S) y FS(h)3116
1299 y F9(.)g(De\014ne)300 1399 y(functions)
y(~) y F6(\036) y FT(2)p FS(k)1461
1585 y F9(=)j F6(\022) y F9(^) y F6(\036)1639
1597 y FT(2)p FS(k)q FR(\000)p FT(1) y F9(+)1911
1563 y(^) y F6(\036) y FT(2)p FS(k)2042
1585 y F9(+) y(~) y F6(\022)
y F9(^) y F6(\036) y FT(2)p FS(k)q
FT(+1) y F9(\(7\))300 1788 y(for)32 b(1)f FU(\024)g
F6(k)j FU(\024)d F6(n)22 b FU(\000)f F9(1,)32 b(and)h
F6(\022)r(;) y F9(~) y F6(\022)g FU(2)f
F9(I)-18 b(R.)33 b(The)f(generalized)f(hierarc)n(hical)g(basis)h(for)g
FU(P) y FS(h) y F9(consists)300 1887
y(of)f(the)h(union)g(of)f(the)h(functions) y(~)1387
1887 y F6(\036) y FT(2)p FS(k) y F9(for)f(1)e
FU(\024)g F6(k)k FU(\024)c F6(n)21 b FU(\000)f F9(1,)31
b(and)h(the)g(basis)e(functions)i(for)f(the)300 1987
y(lev)n(el)25 b(2)g(no)r(des,)820 1965 y(^)809 1987 y
F6(\036)858 1999 y FT(2)p FS(k)q FT(+1) y F6(;)j
F9(0)22 b FU(\024)h F6(k)j FU(\024)d F6(n)14 b FU(\000)g
F9(1.)25 b(This)h(basis)f(will)g(b)r(e)h(denoted)g(b)n(y)f
F6(\036) y FS(k) y F6(;)i F9(1)c FU(\024)g
F6(k)i FU(\024)e F9(2)p F6(n)14 b FU(\000)g F9(1.)300
2087 y(The)33 b(generalized)e(hierarc)n(hical)g(basis)h(in)n(tro)r
(duces)g(a)h(natural)f(direct)g(sum)h(decomp)r(osition)g(of)300
2186 y(the)g(space)f FU(P)733 2198 y FS(h)775 2186 y
F9(.)h(If)g F6(u)e FU(2)g(P) y FS(h)
y F9(,)i(then)g(w)n(e)f(ha)n(v)n(e)f(the)i(unique)g(decomp)r(osition)f
F6(u)f F9(=)f F6(v)25 b F9(+)d F6(w)r F9(,)33 b(where)300
2286 y F6(v)e FU(2)d(V)34 b F9(=)27 b F6(span)p FU(f)859
2264 y F9(^)848 2286 y F6(\036)897 2298 y FT(2)p FS(k)q
FT(+1) y FU(g) y FS(n)p FR(\000)p FT(1)1096
2311 y FS(k)q FT(=0) y F9(and)j F6(w)h FU(2)c(W)35
b F9(=)27 b F6(span)p FU(f) y F9(~)
y F6(\036) y FT(2)p FS(k) y FU(g)2182
2250 y FS(n)p FR(\000)p FT(1) y FS(k)q FT(=1)2312
2286 y F9(.)j(The)h(generalized)e(hierarc)n(hical)300
2386 y(basis)23 b F6( )554 2398 y FS(k)595 2386 y F6(;)14
b F9(1)23 b FU(\024)f F6(k)k FU(\024)d F9(2)p F6(n)12
b FU(\000)g F9(1,)22 b(of)i FU(S) y FS(h)1417
2386 y F9(is)g(de\014ned)h(similarly)-7 b(,)23 b(but)i(with)g(constan)n
(ts)e F6(\027)29 b F9(and)g(~)-47 b F6(\027)30 b F9(instead)300
2485 y(of)e F6(\022)i F9(and)632 2463 y(~)625 2485 y
F6(\022)r F9(.)383 2585 y(Using)d(\(5\))h(and)g(\(7\),)f(w)n(e)h
(obtain)571 2983 y F6(b)p F9(\()p F6(\036)688 2995 y
FT(2)p FS(k)762 2983 y F6(;)14 b( )853 2995 y FS(j)888
2983 y F9(\))24 b(=) y FQ(8) y(>)1031
2813 y(>) y(>) y(>) y() y(>) y(>) y(>)1031
3137 y(:) y F9(^)-43 b F6(a)22 b F9(=)h F6(\022)7
b F9(~)-47 b F6(\027)6 b(c)18 b F9(+)g(\()5 b(~)-47 b
F6(\027)24 b F9(+)18 b F6(\022)r F9(\))p F6(a)863 b(j)28
b F9(=)23 b(2)p F6(k)d FU(\000)e F9(2) y F6(p)k
F9(=)h F6(\022)r(c)c F9(+)f F6(a)1164 b(j)28 b F9(=)23
b(2)p F6(k)d FU(\000)e F9(1) y(^)-43 b F6(c)23
b F9(=)f(\(1)c(+) y(~) y F6(\022)7
b F9(~)-46 b F6(\027)23 b F9(+)18 b F6(\022)r(\027)5
b F9(\))p F6(c)20 b F9(+)e(\()p F6(\022)j F9(+)h(~)-46
b F6(\027)5 b F9(\))p F6(b)18 b F9(+)g(\()p F6(\027)25
b F9(+) y(~) y F6(\022)r F9(\))p F6(a)83
b(j)28 b F9(=)23 b(2)p F6(k) y(q)j F9(=)1304
3063 y(~) y F6(\022)s(c)18 b F9(+)g F6(b)1174
b(j)28 b F9(=)23 b(2)p F6(k)d F9(+)e(1) y(^)1147
3194 y F6(b)k F9(=)h F6(\027) y F9(~)
y F6(\022)r(c)c F9(+)f(\()p F6(\027)24 b F9(+)
y(~) y F6(\022)r F9(\))p F6(b)879 b(j)28 b F9(=)23
b(2)p F6(k)d F9(+)e(2) y(\(8\))300 3377 y(for)27
b(1)c FU(\024)f F6(k)k FU(\024)d F6(n)18 b FU(\000)g
F9(1,)27 b(and) y F6(b)p F9(\()p F6(\036)1185
3670 y FT(2)p FS(k)q FT(+1) y F6(;)14 b( )1435
3670 y FS(j) y F9(\))24 b(=) y FQ(8)1612
3563 y(<) y(:) y F6(r)i F9(=)h(~)-47
b F6(\027)6 b(c)18 b F9(+)g F6(a)83 b(j)28 b F9(=)22
b(2)p F6(k) y(c)425 b(j)28 b F9(=)22 b(2)p F6(k)f
F9(+)d(1) y F6(s)23 b F9(=)f F6(\027)5 b(c)19
b F9(+)f F6(b)92 b(j)28 b F9(=)22 b(2)p F6(k)f F9(+)d(2)3282
3658 y(\(9\))300 3940 y(for)27 b(0)c FU(\024)f F6(k)k
FU(\024)d F6(n)18 b FU(\000)g F9(1.)383 4040 y(The)24
b(sti\013ness)g(matrix)f F6(A) y FS(H)t(B)1327
4040 y F9(corresp)r(onding)f(to)i(the)g(hierarc)n(hical)e(basis)h(is)h
(a)f(p)r(en)n(tadiagonal)300 4139 y(matrix)k(giv)n(en)g(b)n(y)971
4707 y F6(A) y FS(H)t(B) y F9(=)1260
4266 y FQ(0) y(B) y(B) y(B)1260
4562 y(B) y(B) y(B) y(B)1260
4761 y(B) y(B) y(B) y(B)1260
4964 y(@) y F6(c)88 b(s)113 b F9(0)
y F6(p)89 b F9(^)-44 b F6(c)115 b(q) y F9(^)1836
4430 y F6(b) y F9(0)85 b F6(r)118 b(c)142 b(s)e
F9(0) y(^)-43 b F6(a)110 b(p)143 b F9(^)-44
b F6(c)142 b(q) y F9(^) y F6(b)1630
4736 y F9(.) y(.) y(.) y(.)1842
4760 y(.) y(.) y(.) y(.)2055
4786 y(.) y(.) y(.) y(.)2350
4736 y(.) y(.) y(.) y(0)d
F6(r)145 b(c)d(s)112 b F9(0) y(^)-43 b F6(a)137
b(p)142 b F9(^)-43 b F6(c)114 b(q) y F9(0)139
b F6(r)117 b(c) y FQ(1) y(C)
y(C) y(C) y(C) y(C)
y(C) y(C) y(C) y(C)
y(C) y(C) y(A) y F6(:)p eop
%%Page: 257 7
257 276 bop 300 282 a Gd(HIERAR)n(CHICAL)23 b(BASIS)i(F)n(OR)g(UNSTR)n
(UCTURED)f(MESHES)894 b(257)383 482 y F9(No)n(w)21 b(w)n(e)f(apply)h(a)
f(simple)i(p)r(erm)n(utation)e(to)h(the)h(matrix)e F6(A)2243
451 y FS(H)t(B) y F9(,)h(in)h(whic)n(h)f(the)g(basis)f
(functions)300 581 y(asso)r(ciated)25 b(with)h(the)h(coarse)d(grid)h(p)
r(oin)n(ts)h(are)f(ordered)g(\014rst,)g(and)h(those)g(asso)r(ciated)f
(with)h(the)300 681 y(\014ne)j(grid)f(p)r(oin)n(ts)g(are)g(ordered)f
(last.)h(If)h(w)n(e)g(denote)f(the)h(relev)-5 b(an)n(t)28
b(p)r(erm)n(utation)g(matrix)g(b)n(y)h F6(P)12 b F9(,)300
780 y(the)28 b(p)r(erm)n(uted)g(matrix)f(is)g(blo)r(c)n(k)g(2)18
b FU(\002)g F9(2)28 b(of)f(the)h(form)f(giv)n(en)g(in)h(\(3\)-\(4\):)
y(\026) y F6(A) y FS(H)t(B)1363
1020 y F9(=)23 b F6(P)12 b(A) y FS(H)t(B)1694
1020 y F6(P) y FS(T) y F9(=)
y FQ(\022) y F6(A) y FS(cc) y
F6(A) y FS(cf) y F6(A) y FS(f)7
b(c) y F6(A) y FS(f)g(f) y
FQ(\023) y F6(;)300 1255 y F9(where)355 1666
y F6(A)417 1678 y FS(cc)504 1666 y F9(=)591 1374 y FQ(0)591
1520 y(B)591 1570 y(B)591 1620 y(B)591 1670 y(B)591 1720
y(B)591 1773 y(@)711 1433 y F9(^)-44 b F6(c)860 1411
y F9(^)863 1433 y F6(b)707 1543 y F9(^)h F6(a)115 b F9(^)-44
b F6(c) y F9(^) y F6(b)837 1639 y F9(.)869
1664 y(.)902 1689 y(.