laspack2.ps

一个用来实现偏微分方程中网格的计算库

003C140C1404A214061402A213011400A35B13FF13037FA3140113001402A31406A2140C
143CB512FC181C7E9B1C>I<B512F8380F007814181408A2140C1404A213011400A35B13
FF13037FA490C7FCA8EAFFF8161C7E9B1B>I<90381F8080EBE061380180193807000700
0E13035A14015A00781300A2127000F01400A6ECFFF0EC0F80007013071278A212387EA2
7E6C130B380180113800E06090381F80001C1E7E9C21>I<39FFF0FFF0390F000F00AC90
B5FCEB000FAD39FFF0FFF01C1C7F9B1F>I<EAFFF0EA0F00B3A8EAFFF00C1C7F9B0F>I<38
07FF8038007C00133CB3127012F8A21338EA7078EA4070EA30E0EA0F80111D7F9B15>I<
39FFF01FE0390F000780EC060014045C5C5C5C5C49C7FC13021306130FEB17801327EB43
C0EB81E013016D7E1478A280143E141E80158015C039FFF03FF01C1C7F9B20>I<EAFFF8
EA0F8090C7FCB01408A31418A2141014301470EB01F0B5FC151C7E9B1A>I<B46CEBFF80
000FECF800A2390BC00178A33809E002A23808F004A3EB7808A3EB3C10A3EB1E20A3EB0F
40A2EB0780A3EB0300121C3AFF8307FF80211C7E9B26>I<B4EB1FF0390F8003809038C0
0100120BEA09E0EA08F0A21378133C133E131E130F14811307EB03C1EB01E114F1130014
79147D143D141FA2808080121CB46C7E1C1C7F9B1F>I<EB3F80EBE0E03803803848487E
000E7F487F003C148000381303007814C0A20070130100F014E0A8007014C000781303A2
00381480003C1307001C14006C130E6C5B6C6C5A3800E0E0EB3F801B1E7E9C20>I<B512
80380F00E01478143C141C141EA5141C143C147814E0EBFF8090C7FCACEAFFF0171C7E9B
1C>I<EB3F80EBE0E03803803848487E000E7F487F003C148000381303007814C0A20070
130100F014E0A8007014C000781303A200381480383C0E07D81C111300380E208E000713
5C3803A0783900F0E020133FEB0060EC3060EC38E0EC3FC0A2EC1F80EC0F001B257E9C20
>I<B5FC380F01E0EB007880141C141EA4141C143C5CEB01E001FFC7FCEB01E0EB007014
7880A61510A2141CA239FFF00E20C7EA03C01C1D7E9B1F>I<3807E080EA1C19EA300513
03EA600112E01300A36C13007E127CEA7FC0EA3FF8EA1FFEEA07FFC61380130FEB07C013
0313011280A300C01380A238E00300EAD002EACC0CEA83F8121E7E9C17>I<007FB512C0
38700F010060130000401440A200C014201280A300001400B1497E3803FFFC1B1C7F9B1E
>I<39FFF01FF0390F000380EC0100B3A26C1302138000035BEA01C03800E018EB7060EB
0F801C1D7F9B1F>I<39FFE00FF0391F0003C0EC01806C1400A238078002A213C000035B
A2EBE00C00011308A26C6C5AA213F8EB7820A26D5AA36D5AA2131F6DC7FCA21306A31C1D
7F9B1F>I<3AFFE1FFC0FF3A1F003E003C001E013C13186C6D1310A32607801F1320A33A
03C0278040A33A01E043C080A33A00F081E100A39038F900F3017913F2A2017E137E013E
137CA2013C133C011C1338A20118131801081310281D7F9B2B>I<39FFF07FC0390FC01E
003807800CEBC00800035B6C6C5A13F000005BEB7880137C013DC7FC133E7F7F80A2EB13
C0EB23E01321EB40F0497E14783801007C00027F141E0006131F001F148039FF807FF01C
1C7F9B1F>I<387FFFF0EA7C01007013E0386003C0A238400780130F1400131E12005B13
7C13785BA2485A1203EBC010EA0780A2EA0F00481330001E13205A14604813E0EAF803B5
FC141C7E9B19>90 D<12FEA212C0B3B312FEA207297C9E0C>I<EA0808EA1010EA2020EA
4040A2EA8080A3EAB0B0EAF8F8EA7878EA30300D0C7A9C15>I<12FEA21206B3B312FEA2
0729809E0C>I<EA1FC0EA3070EA78387F12301200A2EA01FCEA0F1C12381270126000E0
1340A3EA603C38304E80381F870012127E9115>97 D<12FC121CAA137CEA1D87381E0180
381C00C014E014601470A6146014E014C0381E018038190700EA10FC141D7F9C17>I<EA
03F8EA0C0CEA181E1230EA700CEA600012E0A61260EA70021230EA1804EA0C18EA03E00F
127F9112>I<EB1F801303AAEA03F3EA0E0BEA1807EA30031270126012E0A61260127012
30EA1807EA0E1B3803E3F0141D7F9C17>I<EA07E0EA0C30EA1818EA300CEA700EEA6006
12E0EAFFFEEAE000A41260EA70021230EA1804EA0C18EA03E00F127F9112>I<13F8EA01
8CEA071E1206EA0E0C1300A6EAFFE0EA0E00B0EA7FE00F1D809C0D>I<EB03803807C4C0
EA1C783838388038301800EA701CA4EA3018EA38386C5AEA27C00020C7FCA21230EA3FF8
6CB4FC1480EA2003386001C0EAC000A33860018038300300EA1C0EEA07F8121C7F9215>
I<12FC121CAA137C1387EA1D03001E1380121CAD38FF9FF0141D7F9C17>I<1218123CA2
1218C7FCA712FC121CB0EAFF80091D7F9C0C>I<13C0EA01E0A2EA00C01300A7EA07E012
00B3A21260EAF0C012F1EA6180EA3E000B25839C0D>I<12FC121CAAEB0FE0EB0780EB06
005B13105B5B13E0121DEA1E70EA1C781338133C131C7F130F148038FF9FE0131D7F9C16
>I<12FC121CB3A9EAFF80091D7F9C0C>I<39FC7E07E0391C838838391D019018001EEBE0
1C001C13C0AD3AFF8FF8FF8021127F9124>I<EAFC7CEA1C87EA1D03001E1380121CAD38
FF9FF014127F9117>I<EA03F0EA0E1CEA1806487E00701380EA600100E013C0A6006013
80EA700300301300EA1806EA0E1CEA03F012127F9115>I<EAFC7CEA1D87381E0180001C
13C0EB00E0A21470A614E0A2EB01C0001E1380381D0700EA1CFC90C7FCA7B47E141A7F91
17>I<3803E080EA0E19EA1805EA3807EA7003A212E0A61270A2EA38071218EA0E1BEA03
E3EA0003A7EB1FF0141A7F9116>I<EAFCE0EA1D38EA1E78A2EA1C301300ACEAFFC00D12
7F9110>I<EA1F90EA2070EA4030EAC010A212E0EAF800EA7F80EA3FE0EA0FF0EA00F8EA
8038131812C0A2EAE010EAD060EA8FC00D127F9110>I<1204A4120CA2121C123CEAFFE0
EA1C00A91310A5120CEA0E20EA03C00C1A7F9910>I<38FC1F80EA1C03AD1307120CEA0E
1B3803E3F014127F9117>I<38FF07E0383C0380381C0100A2EA0E02A2EA0F06EA0704A2
EA0388A213C8EA01D0A2EA00E0A3134013127F9116>I<39FF3FC7E0393C0703C0001CEB
01801500130B000E1382A21311000713C4A213203803A0E8A2EBC06800011370A2EB8030
000013201B127F911E>I<38FF0FE0381E0700EA1C06EA0E046C5AEA039013B0EA01E012
007F12011338EA021C1204EA0C0E487E003C138038FE1FF014127F9116>I<38FF07E038
3C0380381C0100A2EA0E02A2EA0F06EA0704A2EA0388A213C8EA01D0A2EA00E0A31340A2
5BA212F000F1C7FC12F312661238131A7F9116>I<EA7FFCEA70381260EA407013F013E0
EA41C012031380EA0700EA0F04120E121CEA3C0CEA380812701338EAFFF80E127F9112>
I<B512F01401808B15>I E end
%%EndProlog
%%BeginSetup
%%Feature: *Resolution 300dpi
TeXDict begin
%%PaperSize: Letter

%%EndSetup
%%Page: 1 1
1 0 bop -9 -260 a Fj(INTR)o(O\(3LAS\))352 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)352 b(INTR)o(O\(3LAS\))-9 -67 y
Fi(NAME)141 -17 y Fj(in)o(tro)13 b({)h(in)o(tro)q(duction)f(to)h
Fh(LASPack)e Fj(\(P)o(ac)o(k)n(age)i(for)g(Linear)f(Algebra)h(with)g
(Sparse)h(Matrices\))-9 82 y Fi(DESCRIPTION)141 132 y
Fh(LASPack)h Fj(is)h(a)g(pac)o(k)n(age)g(for)h(solving)e(large)h
(sparse)i(systems)f(of)f(linear)g(equations)g(lik)o(e)g(those)h(whic)o
(h)g(arise)141 182 y(from)12 b(discretization)i(of)f(partial)g
(di\013eren)o(tial)h(equations.)141 282 y(Main)f(features:)203
381 y Fg(\017)20 b Fj(The)g(primary)d(aim)g(of)i Fh(LASPack)e
Fj(is)i(the)h(implem)o(en)o(tation)c(of)j(e\016cien)o(t)g(iterativ)o(e)
g(metho)q(ds)g(for)g(the)244 431 y(solution)d(of)g(systems)g(of)g
(linear)g(equations.)25 b(All)16 b(routines)h(and)f(data)g(structures)j
(are)e(optimized)e(for)244 481 y(e\013ectiv)o(e)g(usage)f(of)f
(resources)j(esp)q(ecially)e(with)f(regard)h(to)f(large)h(sparse)g
(matrices.)k(The)c(pac)o(k)n(age)f(can)244 531 y(b)q(e)k(accessed)i
(from)14 b(an)i(application)f(through)h(a)g(straigh)o(tforw)o(ard)g(in)
o(terface)g(de\014ned)i(in)e(the)g(form)f(of)244 581
y(pro)q(cedure)h(calls.)203 667 y Fg(\017)k Fj(Beside)c(the)e
(obligatory)f(Jacobi,)g(successiv)o(e)j(o)o(v)o(er-relaxation,)d(Cheb)o
(yshev,)h(and)g(conjugate)g(gradien)o(t)244 716 y(solv)o(ers,)h
Fh(LASPack)e Fj(con)o(tains)h(selected)j(state-of-the-art)e(algorithms)
d(\(cf.)i(Barrett)i(et)f(al.)k([1],)13 b(McCor-)244 766
y(mic)o(k[5]\))e(whic)o(h)j(are)h(commonl)o(y)c(used)k(for)e(large)h
(sparse)h(systems:)291 852 y Fi({)21 b Fj(CG-lik)o(e)13
b(metho)q(ds)h(for)h(non-symmetric)e(systems:)20 b(CGN,)14
b(GMRES,)f(BiCG,)h(QMR,)g(CGS,)g(and)336 902 y(BiCGStab,)291
970 y Fi({)21 b Fj(m)o(ultilev)o(el)6 b(metho)q(ds)j(suc)o(h)i(as)e
(the)h(m)o(ultigrid)d(and)i(conjugate)g(gradien)o(t)h(metho)q(ds)f
(preconditioned)336 1020 y(b)o(y)k(m)o(ultigrid)e(and)j(BPX)h
(preconditioners.)244 1105 y(All)g(the)h(ab)q(o)o(v)o(e)f(solv)o(ers)g
(are)h(applicable)f(not)g(only)f(to)i(the)f(p)q(ositiv)o(e)g
(de\014nite)h(or)g(non-symmetric)d(ma-)244 1155 y(trices,)h(but)e(are)h
(also)f(adopted)g(for)g(singular)g(systems)h(\(e.g.)e(arising)h(from)f
(discretization)h(of)g(Neumann)244 1205 y(b)q(oundary)i(v)n(alue)f
(problems\).)203 1291 y Fg(\017)20 b Fj(The)d(implem)o(en)o(tation)c
(is)j(based)g(on)g(an)g(ob)r(ject-orien)o(ted)h(approac)o(h)f
(\(although)f(C)h(as)g(programmi)o(ng)244 1341 y(language)g(is)g
(used\).)27 b(V)m(ectors)18 b(and)e(matrices)g(are)h(de\014ned)g(as)g
(new)g(data)f(t)o(yp)q(es)h(in)f(connection)h(with)244
1391 y(the)d(corresp)q(onding)g(supp)q(orting)f(routines.)19
b(The)13 b(basic)h(op)q(erations)f(are)h(implem)o(en)o(ted)d(in)i(suc)o
(h)h(a)f(w)o(a)o(y)244 1440 y(that)h(they)h(allo)o(w)d(the)i
(programming)d(of)i(linear)g(algebra)h(algorithms)d(in)j(a)f(natural)h
(w)o(a)o(y)m(.)203 1526 y Fg(\017)20 b Fh(LASPack)11
b Fj(is)h(extensible)h(in)f(a)f(simple)g(manner.)16 b(An)d(access)h(to)
e(the)g(in)o(ternal)g(represen)o(tation)i(of)d(v)o(ectors)244
1576 y(and)18 b(matrices)f(is)h(not)f(necessary)j(and)e(is)f(a)o(v)o
(oided,)h(as)g(required)g(of)g(the)g(ob)r(ject-orien)o(ted)h(program-)
244 1626 y(ming.)j(This)16 b(allo)o(ws)e(an)i(impro)o(v)o(emen)o(t)d
(of)i(algorithms)f(or)h(a)h(mo)q(di\014cation)d(of)j(data)f(structures)
j(with)244 1676 y(no)c(adjustmen)o(t)f(of)g(application)g(programs)f
(using)i(the)h(pac)o(k)n(age.)203 1762 y Fg(\017)20 b
Fh(LASPack)13 b Fj(is)g(written)i(in)e(ANSI)h(C)g(and)g(is)g(th)o(us)g
(largely)f(p)q(ortable.)141 1896 y(The)h(structure)i(of)d
Fh(LASPack)g Fj(is)g(depicted)i(in)f(\014gure)g(1.)k(The)c(library)g
(consists)g(of)g(sev)o(eral)g(mo)q(dules)f(of)g(whic)o(h)141
1946 y(some)g(ma)o(y)e(b)q(e)k(group)q(ed)f(together)h(forming)d
(larger)i(units.)141 2045 y(The)h(mo)q(dules)g Fi(VECTOR)p
Fj(,)g Fi(MA)l(TRIX)p Fj(,)h(and)f Fi(QMA)l(TRIX)h Fj(as)f(w)o(ell)g
(as)g Fi(OPERA)l(TS)p Fj(,)f(and)h Fi(F)-5 b(A)o(CTOR)141
2095 y Fj(are)15 b(the)g(main)d(comp)q(onen)o(t)i(of)g(the)h(library)m
(.)k(The)c(basic)g(ob)r(jects)g(of)f(linear)g(algebra,)g(v)o(ectors)i
(and)e(matrices,)141 2145 y(are)f(here)h(implemen)o(ted)d(as)i(data)f
(t)o(yp)q(es.)19 b(Their)13 b(de\014nition)g(consists)h(not)e(only)h
(of)f(data)g(structures,)k(but)d(also)141 2195 y(of)f(the)h(corresp)q
(onding)h(managemen)o(t)c(routines,)j(e.g.)f(for)h(allo)q(cation)e(and)
i(release)h(of)e(v)n(ariables)g(of)g(the)h(ab)q(o)o(v)o(e)141
2245 y(t)o(yp)q(es,)18 b(or)f(storage)h(and)f(querying)g(of)g(v)o
(ector)h(comp)q(onen)o(ts)f(and)g(matrix)f(elemen)o(ts.)28
b(F)m(or)17 b(matrix)e(storage,)141 2294 y(the)i(compressed)h(ro)o(w)e
(format)f(or)i(the)g(compressed)h(column)d(format)g(can)i(b)q(e)g
(used.)28 b(Tw)o(o)16 b(di\013eren)o(t)i(t)o(yp)q(es)141
2344 y(of)c(matrices)g(are)h(implem)o(en)o(ted)e(in)h
Fh(LASPack)o Fj(.)20 b(The)15 b(t)o(yp)q(e)g Fh(Matrix)e
Fj(is)i(a)f(data)h(t)o(yp)q(e)g(for)f(general)h(rectangular)141
2394 y(sparse)f(matrices.)k(V)m(ariables)12 b(of)h(suc)o(h)h(a)f(kind)g
(in)g Fh(LASPack)f Fj(are)i(in)o(tended)g(esp)q(ecially)g(for)f(the)h
(description)g(of)141 2444 y(restriction)g(and)f(prolongation)e(op)q
(erators)j(within)f(m)o(ultigrid)d(algorithms.)16 b(The)d(t)o(yp)q(e)h
Fh(QMatrix)d Fj(is)i(designed)-9 2594 y(V)m(ersion:)18
b(1.12.3)410 b(Last)14 b(c)o(hange:)19 b(Jan)o(uary)14
b(17,)e(1996)661 b(1)p eop
%%Page: 2 2
2 1 bop -9 -260 a Fj(INTR)o(O\(3LAS\))352 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)352 b(INTR)o(O\(3LAS\))p 420 78
1213 2 v 420 172 2 95 v 932 143 a Ff(MLSOLV)p 1631 172
V 420 174 1213 2 v 420 238 V 420 333 2 95 v 901 304 a(ITERSOLV)p
1631 333 V 420 335 1213 2 v 979 432 a(RTC)p 420 559 911
2 v 420 654 2 95 v 765 625 a(PRECOND)p 1329 654 V 420
656 911 2 v 1327 559 306 2 v 1327 654 2 95 v 372 w(EIGENVAL)p
1631 654 V 1327 656 306 2 v 420 720 760 2 v 420 815 2
95 v 690 786 a(OPERATS)p 1178 815 V 420 817 760 2 v 1176
720 457 2 v 1176 815 2 95 v 403 w(FACTOR)p 1631 815 V
1176 817 457 2 v 420 815 306 2 v 420 909 2 95 v 479 880
a(VECTOR)p 724 909 V 420 911 306 2 v 722 815 V 722 909
2 95 v 116 w(MATRIX)p 1026 909 V 722 911 306 2 v 1025
815 609 2 v 1025 909 2 95 v 1219 876 a(QMATRIX)p 1631
909 V 1025 911 609 2 v 901 1008 a(ERRHANDL)p 418 77 1217
4 v 418 171 4 95 v 1631 171 V 418 174 1217 4 v 418 237
V 418 332 4 95 v 1631 332 V 418 335 1217 4 v 352 11 1349
4 v 352 492 4 482 v 1697 492 V 352 495 1349 4 v 418 559
1217 4 v 418 653 4 95 v 1631 653 V 418 656 1217 4 v 418
719 V 418 908 4 189 v 1631 908 V 418 911 1217 4 v 286
-56 1481 4 v 286 1069 4 1125 v 1763 1069 V 286 1072 1481
4 v 743 1203 a Fj(Figure)14 b(1:)k(Structure)e(of)d Fh(LASPack)141
1398 y Fj(for)e(quadratic)g(sparse)i(matrices)e(of)g(systems)g(of)g
(linear)g(equations.)18 b(With)10 b(regard)i(to)g(the)g(t)o(yp)q(e)g
Fh(Matrix)p Fj(,)e(some)141 1448 y(additional)k(prop)q(erties)k(are)f
(a)o(v)n(ailable)d(for)i(this)g(data)g(t)o(yp)q(e)h(suc)o(h)g(as)g
(symmetry)m(,)d(in)o(v)o(ertabilit)o(y)g(etc.)27 b(These)141
1497 y(are)14 b(also)f(tak)o(en)h(in)o(to)f(accoun)o(t)i(for)e(storage)
i(and)e(in)h(n)o(umerical)e(algorithms.)141 1597 y(The)17
b(basic)h(op)q(erations)f(of)g(linear)g(algebra)f(are)i(implemen)o(ted)
d(in)i(the)h(mo)q(dule)e Fi(OPERA)l(TS)p Fj(.)g(It)h(con)o(tains)141
1647 y(pro)q(cedures)j(for)d(e.g.)g(assignmen)o(t)g(and)g(addition)g
(of)g(v)o(ectors,)i(m)o(ultiplication)c(of)i(v)o(ectors)i(b)o(y)e
(scalars)h(and)141 1697 y(matrices,)e(or)h(generation)f(of)h(transp)q
(osed,)h(diagonal,)d(and)h(triangular)g(matrices.)26
b(F)m(or)16 b(complex)g(algebraic)141 1747 y(expressions,)f(the)g