laspack2.ps

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

(usage)g(of)f(temp)q(orary)g(v)n(ariables)g(allo)o(ws)f(to)i(c)o(hain)f
(these)i(routines)f(up)f(so)h(that)g(n)o(umerical)141
1796 y(algorithms)c(could)j(b)q(e)g(transcrib)q(ed)i(directly)e(in)f
Fh(LASPack)o Fj(.)141 1896 y(The)h(mo)q(dule)e Fi(F)-5
b(A)o(CTOR)14 b Fj(consists)h(of)e(pro)q(cedures)j(for)e(incomplete)f
(factorization)g(of)g(quadratic)h(matrices.)141 1996
y(In)22 b(top-lev)o(el)g(mo)q(dules,)i(sev)o(eral)f(iterativ)o(e)g
(solv)o(ers)g(are)g(implemen)o(ted:)33 b(the)24 b(classical)e(iterativ)
o(e,)i(semi-)141 2045 y(iterativ)o(e,)11 b(and)g(conjugate)g(gradien)o
(t)g(metho)q(ds,)f(selected)j(pro)q(cedures)h(for)c(non-symmetric)f
(systems)i(\(mo)q(dule)141 2095 y Fi(ITERSOL)-5 b(V)p
Fj(\),)12 b(and)i(m)o(ultilev)o(el)d(solv)o(ers)k(\(mo)q(dule)d
Fi(MLSOL)-5 b(V)p Fj(\).)141 2195 y(In)13 b(addition,)f(the)i(most)e
(common)f(t)o(yp)q(es)k(of)d(preconditioners)j(are)f(a)o(v)n(ailable)d
(in)i(the)h(mo)q(dule)f Fi(PRECOND)p Fj(.)141 2245 y(The)j(mo)q(dule)f
Fi(EIGENV)-5 b(AL)16 b Fj(con)o(tains)h(pro)q(cedures)h(for)e(the)h
(estimation)d(of)i(extremal)f(eigen)o(v)n(alues)h(whic)o(h)141
2294 y(are)e(required)h(e.g.)e(b)o(y)g(the)i(Cheb)o(yshev)g(metho)q(d.)
141 2394 y(Whereas)d(the)g(ab)q(o)o(v)o(e)f(mo)q(dules)g(de\014ne)h
(basic)g(ob)r(jects)g(of)f(linear)g(algebra)g(and)g(the)i(corresp)q
(onding)f(op)q(erations)141 2444 y(concerning)i(them,)f(the)h
(remaining)e(mo)q(dules)h Fi(R)l(TC)h Fj(and)g Fi(ERRHANDL)g
Fj(are)g(of)g(a)f(global)f(kind.)-9 2594 y(2)661 b(Last)14
b(c)o(hange:)19 b(Jan)o(uary)14 b(17,)e(1996)412 b(V)m(ersion:)18
b(1.12.3)p eop
%%Page: 3 3
3 2 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\))141 -67 y(Con)o(v)o
(ergence)10 b(of)f(iterativ)o(e)h(solution)e(is)i(con)o(trolled)f(b)o
(y)g(routines)h(of)f(the)h(mo)q(dule)e Fi(R)l(TC)p Fj(.)i(In)f(order)h
(to)g(terminate)141 -17 y(the)k(solution)f(pro)q(cess,)i(a)f(residual)g
(criterion)g(is)g(applied.)141 83 y(Routines)h(in)f(the)i(mo)q(dule)e
Fi(ERRHANDL)h Fj(supp)q(ort)h(the)g(error)g(handling.)21
b(They)16 b(sup)q(ervise)h(some)d(excep-)141 133 y(tional)e(states)j
(whic)o(h)e(ma)o(y)f(happ)q(en)i(during)f(generating)h(or)f(solving)f
(systems)i(of)f(equations)h(and)f(allo)o(w)f(their)141
182 y(treatmen)o(t)h(in)h(the)g(application)f(co)q(de.)141
282 y(F)m(rom)c(a)j(global)e(\(ob)r(ject-orien)o(ted\))j(p)q(oin)o(t)e
(of)h(view,)f Fh(LASPack)f Fj(consists)j(of)e(de\014nition)h(of)f(the)h
(ob)r(jects)h Fh(Vector)p Fj(,)141 332 y Fh(Matrix)p
Fj(,)20 b(and)g Fh(QMatrix)f Fj(and)h(of)g(a)g(set)h(of)f(more)f(or)h
(less)h(complicated)e(op)q(erations)i(on)f(them.)37 b(F)m(or)20
b(the)141 382 y(extensibilit)o(y)12 b(of)g(the)i(library)m(,)e(there)i
(are)g(the)f(strict)h(hiding)e(of)h(in)o(ternal)f(represen)o(tation)j
(of)d(these)j(ob)r(jects)f(as)141 431 y(w)o(ell)c(as)h(the)h(comm)o
(unicatio)o(n)c(b)q(et)o(w)o(een)13 b(the)f(hierarc)o(hically)e(built)g
(mo)q(dules)g(b)o(y)h(simple)f(structured)j(in)o(terfaces)141
481 y(in)h(the)h(form)d(of)i(pro)q(cedure)i(calls.)k(This)14
b(mak)o(es)f(it)h(p)q(ossible)h(to)f(extend)i(or)e(exc)o(hange)h
Fh(LASPack)e Fj(mo)q(dules)h(in)141 531 y(a)f(simple)g(manner.)141
631 y Fh(LASPack)f Fj(is)i(written)h(in)e(ANSI)h(C)g(and)g(th)o(us)h(p)
q(ortable)f(to)g(most)f(computer)g(platforms.)k(It)d(w)o(as)g
(successfully)141 681 y(tested)k(on)f(Sun)g(Sparc5)g(\(SunOS)h(4.1.3,)d
(gcc)j(2.4.5\),)d(HP)j(9000/735)d(\(HP-UX)i(9.05\),)f(IBM)i(RS/6000)d
(550)141 730 y(\(AIX)g(3.2.5\),)f(DEC)h(3000/800)e(M)i(\(OSF/1)g
(3.0\),)f(SGI)h(IRIS)g(Indigo)f(\(IRIX)h(5.3\),)f(and)h(PC)h(486)e
(\(MS-DOS)141 780 y(6.2,)e(BC)i(2.0/MSC)f(6.0;)f(Lin)o(ux)i(1.1.54,)d
(gcc)j(2.5.8/lib)q(c)e(4.5.26\).)-9 881 y Fi(REFERENCES)244
976 y Fj([1])g(R.)f(Barrett,)j(M.)e(Berry)m(,)h(T.)e(F.)h(Chan,)g(J.)h
(Demmel,)c(J.)j(Donato,)g(J.)g(Dongarra,)f(V.)h(Eijkhout,)244
1026 y(R.)f(P)o(ozo,)g(Ch.)g(Romine,)f(H.)h(v)n(an)g(der)h(V)m(orst:)35
b(T)m(emplates)10 b(for)h(the)h(Solution)f(of)g(Linear)g(Systems:)244
1076 y(Building)j(Blo)q(c)o(ks)h(for)g(Iterativ)o(e)g(Solv)o(ers,)g
(SIAM,)g(Philadelphia,)e(1994.)20 b(\(A)15 b(HTML)g(v)o(ersion)h(is)244
1126 y(also)e(a)o(v)n(ailable\).)244 1191 y([2])f(G.)g(Golub,)g(C.)h(v)
n(an)f(Loan:)18 b(Matrix)c(Computations,)e(second)j(edition,)e(The)i
(Johns)g(Hopkins)244 1241 y(Univ)o(ersit)o(y)f(Press,)h(Baltimore,)d
(1989.)244 1307 y([3])i(W.)f(Hac)o(kbusc)o(h:)20 b(Iterativ)o(e)15
b(Solution)e(of)h(Large)h(Sparse)g(Systems)f(of)g(Equations,)g
(Springer-)244 1356 y(V)m(erlag,)f(Berlin,)h(1994.)244
1422 y([4])20 b(W.)h(Hac)o(kbusc)o(h:)33 b(Multi-Grid)20
b(Metho)q(ds)i(and)f(Applications,)h(Springer-V)m(erlag,)h(Berlin,)244
1472 y(1985.)244 1537 y([5])13 b(S.)g(F.)h(McCormic)o(k:)j(Multigrid)12
b(Metho)q(ds,)j(SIAM,)e(Philadelphia,)f(1987.)-9 1666
y Fi(A)o(UTHOR)141 1716 y Fj(The)i(library)f Fh(LASPack)f
Fj(w)o(as)i(dev)o(elop)q(ed)h(b)o(y)e(T)m(om\023)-21
b(a)n(\024)i(s)12 b(Sk)n(alic)o(k)q(\023)-21 b(y.)-9
1816 y Fi(COPYRIGHT)141 1866 y Fj(Cop)o(yrigh)o(t)13
b(\(C\))h(1992-1996)d(T)m(om\023)-21 b(a)n(\024)i(s)12
b(Sk)n(alic)o(k)q(\023)-21 b(y.)16 b(All)d(righ)o(ts)h(reserv)o(ed.)375
1965 y Fi(ANY)j(USE)e(OF)h(THIS)g(CODE)h(CONSTITUTES)f(A)o(CCEPT)l
(ANCE)479 2015 y(OF)g(THE)h(TERMS)f(OF)g(THE)g(COPYRIGHT)g(NOTICE)141
2144 y Fj(See)e(the)h(\014le)f Fh(laspack/copyrght)o(.h)d
Fj(for)i(details.)-9 2245 y Fi(FILES)203 2394 y Fg(\017)20
b Fj(installed:)244 2444 y Fi(/usr/lo)q(cal/incl)o(ud)o(e/l)o(aspac)o
(k/*.h)103 b Fh(LASPack)13 b Fj(header)i(\014les)-9 2594
y(V)m(ersion:)j(1.12.3)410 b(Last)14 b(c)o(hange:)19
b(Jan)o(uary)14 b(17,)e(1996)661 b(3)p eop
%%Page: 4 4
4 3 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\))244 -64 y
Fi(/usr/lo)q(cal/lib)o(/l)o(ib)o(l)o(aspac)o(k.a)198
b Fh(LASPack)13 b Fj(library)203 19 y Fg(\017)20 b Fj(in)14
b(the)g(distribution:)244 69 y Fi(laspac)o(k/*.h)511
b Fh(LASPack)13 b Fj(header)i(\014les)244 118 y Fi(laspac)o(k/*.c)517
b Fh(LASPack)13 b Fj(sources)244 168 y Fi(laspac)o(k/examples/*)335
b Fj(some)13 b(examples)g(and)h(test)h(programs)244 218
y Fi(laspac)o(k/do)q(c/*.ps)395 b Fj(reference)17 b(man)o(ual)11
b(incl.)i(these)i(man)o(ual)d(pages)-9 351 y Fi(LIST)j(OF)h(MANUAL)g(P)
l(A)o(GES)269 436 y Fj(eigen)o(v)n(al\(3LAS\))144 b(estimation)12
b(of)h(extremal)g(eigen)o(v)n(alues)269 486 y(errhandl\(3LAS\))136
b(error)15 b(handling)d(routines)269 535 y(factor\(3LAS\))183
b(incomplete)13 b(factorization)g(of)g(quadratic)h(matrices)269
585 y(itersolv\(3LAS\))155 b(classical)14 b(iterativ)o(e,)f
(semi-iterativ)o(e,)f(CG,)h(and)h(CG-lik)o(e)e(solv)o(ers)269
635 y(matrix\(3LAS\))165 b(t)o(yp)q(e)14 b Fh(Matrix)f
Fj(for)g(general)i(rectangular)f(sparse)h(matrices)269
685 y(mlsolv\(3LAS\))169 b(m)o(ultilev)o(el)11 b(solv)o(ers)269
735 y(op)q(erats\(3LAS\))156 b(basic)14 b(op)q(erations)g(of)f(linear)h
(algebra)269 785 y(precond\(3LAS\))146 b(pre-de\014ned)16
b(preconditioners)269 834 y(qmatrix\(3LAS\))143 b(t)o(yp)q(e)14
b Fh(QMatrix)f Fj(for)g(quadratic)h(sparse)h(matrices)269
884 y(rtc\(3LAS\))238 b(residual)14 b(termination)e(con)o(trol)i(of)f
(iterativ)o(e)h(solv)o(ers)269 934 y(v)o(ector\(3LAS\))178
b(t)o(yp)q(e)14 b Fh(Vector)-9 2594 y Fj(4)661 b(Last)14
b(c)o(hange:)19 b(Jan)o(uary)14 b(17,)e(1996)412 b(V)m(ersion:)18
b(1.12.3)p eop
%%Page: 5 5
5 4 bop -9 -260 a Fj(EIGENV)-5 b(AL\(3LAS\))272 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)271 b(EIGENV)-5 b(AL\(3LAS\))-9
-67 y Fi(NAME)141 -17 y Fh(SetEigenvalAccu)o(racy)p Fj(,)10
b Fh(GetMaxEigenval)p Fj(,)g Fh(GetMinEigenval)h Fj({)j(estimation)e
(of)h(extremal)g(eigen)o(v)n(alues)-9 82 y Fi(SYNOPSIS)141
143 y Fh(#include)20 b(<laspack/eigenva)o(l.h>)141 242
y(void)h(SetEigenvalAccu)o(racy\()o(doubl)o(e)e(Eps\);)141
292 y(double)h(GetMinEigenval\(QMa)o(trix)e(*Q,)k(PrecondProcType)c
(PrecondProc,)620 342 y(double)j(OmegaPrecond\);)141
392 y(double)f(GetMaxEigenval\(QMa)o(trix)e(*Q,)k(PrecondProcType)c
(PrecondProc,)620 441 y(double)j(OmegaPrecond\);)-9 568
y Fi(DESCRIPTION)141 618 y Fj(The)10 b(pro)q(cedures)j
Fh(GetMinEigenval)8 b Fj(and)i Fh(GetMaxEigenval)d Fj(estimate)j
(extremal)g(eigen)o(v)n(alues)g(of)g(the)h(matrix)141
668 y Fh(Q)p Fj(.)g(In)h Fh(LASPack)18 b Fj(b)q(oth)13
b(eigen)o(v)n(alues)f(are)g(required)i(in)o(ternally)d(for)h(the)h
(Cheb)o(yshev)g(metho)q(d.)k(Otherwise)d(they)141 718
y(could)e(b)q(e)h(applied)e(in)h(an)g(application,)f(e.g.)h(for)g(the)h
(estimation)d(of)i(the)h(condition)f(n)o(um)o(b)q(er)f(of)h(a)g(matrix)
f(\(see)141 768 y(example)h(b)q(elo)o(w\).)141 867 y(The)i(algorithm)e
(used)j(is)f(based)h(on)f(the)h(Lanczos)g(metho)q(d)f(for)g(symmetric)e
(matrices.)19 b(Eigen)o(v)n(alues)14 b(of)f(non-)141
917 y(symmetric)h(matrices)i(are)g(estimated)g(b)o(y)g(means)g(of)f
Fh(Q)1038 902 y Fe(T)1064 917 y Fh(Q)p Fj(.)25 b(F)m(or)16
b(singular)f(matrices,)h(if)g(the)h(n)o(ull)e(space)i(has)141
967 y(b)q(een)e(sp)q(eci\014ed,)g(the)f(zero)h(eigen)o(v)n(alue)e(is)h
(not)g(considered.)141 1067 y(If)g(an)o(y)g(pro)q(cedure)i(for)e
(preconditioning)g(is)g(giv)o(en)g(b)o(y)h(the)g(parameter)f
Fh(PrecondProc)e Fj(\(i.e.)i(if)f Fh(PrecondProc)141
1116 y(!=)21 b(NULL)p Fj(\),)16 b(eigen)o(v)n(alues)h(are)h(computed)e
(for)h(the)h(preconditioned)g(matrix.)26 b(The)18 b(pro)q(cedure)h
(used)f(should)141 1166 y(ha)o(v)o(e)13 b(the)i(protot)o(yp)q(e)184
1266 y Fh(Vector)21 b(*\(*PrecondProcTy)o(pe\)\(Q)o(Matri)o(x)e(*,)i
(Vector)g(*,)g(Vector)g(*,)g(double\))141 1404 y Fj(The)14
b(parameter)f Fh(OmegaPrecond)f Fj(will)g(b)q(e)j(passed)g(as)f(a)f
(relaxation)g(parameter)h(to)f(the)i(preconditioner.)141
1504 y(The)d(accuracy)i(for)e(the)h(stopping)f(criterion)h(can)g(b)q(e)
g(set)g(b)o(y)g Fh(SetEigenvalAccur)o(acy)p Fj(.)h(The)f(default)f(is)h
(10)1857 1489 y Fd(\000)p Fc(4)1901 1504 y Fj(.)141 1604
y(Eac)o(h)h(of)g(the)g(pro)q(cedures)j Fh(GetMinEigenval)11
b Fj(and)j Fh(GetMaxEigenval)d Fj(causes)16 b(the)e(estimation)f(of)g
(b)q(oth,)h(mi-)141 1654 y(nim)o(um)f(and)k(maxim)n(um)c(eigen)o(v)n
(alues.)27 b(They)18 b(are)f(stored)h(together)g(with)f(the)g
(parameters)g Fh(PrecondProc)141 1703 y Fj(and)12 b Fh(OmegaPrecond)e
Fj(in)i(connection)i(with)e(the)h(matrix)e Fh(Q)p Fj(,)h(and)g(are)h
(returned)i(at)d(the)h(next)g(call)f(of)g(the)h(ab)q(o)o(v)o(e)141
1753 y(routines)h(if)g(the)h(same)e(parameters)i(are)g(used.)20
b(Therefore,)15 b(the)g(extremal)e(eigen)o(v)n(alues)i(need)g(not)f(b)q
(e)h(stored)141 1803 y(additionally)m(,)10 b(but)k(could)g(b)q(e)g
(queried)h(on)f(demand)e(without)i(extra)g(computational)e(costs.)-9
1909 y Fi(REFERENCES)141 1959 y Fj(The)i(implemen)o(tatio)o(n)d(of)i
(the)i(Lanczos)g(metho)q(d)e(follo)o(ws)f(algorithms)g(describ)q(ed)j
(in:)244 2058 y(G.)j(Golub,)g(C.)h(v)n(an)f(Loan:)27
b(Matrix)18 b(Computations,)g(second)h(edition,)g(The)g(Johns)h
(Hopkins)244 2108 y(Univ)o(ersit)o(y)14 b(Press,)h(Baltimore,)d(1989.)
-9 2247 y Fi(FILES)141 2288 y(eigen)o(v)m(al.h)543 b
Fj(header)14 b(\014le)141 2338 y Fi(eigen)o(v)m(al.c)549
b Fj(source)15 b(\014le)-9 2444 y Fi(EXAMPLES)-9 2594
y Fj(V)m(ersion:)j(1.12.3)410 b(Last)14 b(c)o(hange:)19
b(Jan)o(uary)14 b(17,)e(1996)661 b(5)p eop
%%Page: 6 6
6 5 bop -9 -260 a Fj(EIGENV)-5 b(AL\(3LAS\))272 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)271 b(EIGENV)-5 b(AL\(3LAS\))141
-67 y(In)15 b(the)h(follo)o(wing)c(co)q(de)17 b(fragmen)o(t,)c
(condition)i(n)o(um)o(b)q(ers)g(of)g(a)g(matrix)e(and)i(their)h
(preconditioned)g(v)n(arian)o(ts)141 -17 y(are)e(estimated.)141
83 y Fh(QMatrix)20 b(A;)141 133 y(double)g(LambdaMax,)g(LambdaMin;)141
182 y(size_t)g(Dim;)141 332 y(Dim)h(=)g(...)g(;)141 382
y(Q_Constr\(&A,)e("A",)i(Dim,)g(True,)g(Rowws,)f(Normal,)g(True\);)141
481 y(/*)h(generation)f(of)h(the)g(matrix)g(A)g(*/)141
581 y(...)141 681 y(printf\("Conditi)o(on)e(numbers)h(of)i(the)f
(matrix)f(\045s:\\n\\n",)g(Q_GetName\(&A\)\);)141 780
y(LambdaMax)f(=)j(GetMaxEigenval\(&A,)c(NULL,)j(1.0\);)141
830 y(LambdaMin)e(=)j(GetMinEigenval\(&A,)c(NULL,)j(1.0\);)141
880 y(printf\("without)d(preconditioning:)128 b(k)21
b(=)h(\04512.5e\\n",)e(LambdaMax)g(/)h(LambdaMin\);)141
979 y(LambdaMax)e(=)j(GetMaxEigenval\(&A,)c(JacobiPrecond,)h(1.0\);)141
1029 y(LambdaMin)g(=)j(GetMinEigenval\(&A,)c(JacobiPrecond,)h(1.0\);)
141 1079 y(printf\("with)g(Jacobi)h(preconditioning:)41
b(k)21 b(=)h(\04512.5e\\n",)e(LambdaMax)g(/)h(LambdaMin\);)141