laspack2.ps

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

1179 y(LambdaMax)e(=)j(GetMaxEigenval\(&A,)c(SSORPrecond,)h(1.0\);)141
1229 y(LambdaMin)g(=)j(GetMinEigenval\(&A,)c(SSORPrecond,)h(1.0\);)141
1278 y(printf\("with)g(SSOR)i(preconditioning:)84 b(k)21
b(=)h(\04512.5e\\n",)e(LambdaMax)g(/)h(LambdaMin\);)141
1378 y(LambdaMax)e(=)j(GetMaxEigenval\(&A,)c(ILUPrecond,)i(1.0\);)141
1428 y(LambdaMin)f(=)j(GetMinEigenval\(&A,)c(ILUPrecond,)i(1.0\);)141
1478 y(printf\("with)f(ILU)i(preconditioning:)106 b(k)21
b(=)h(\04512.5e\\n",)e(LambdaMax)g(/)h(LambdaMin\);)141
1577 y(Q_Destr\(&A\);)-9 1710 y Fi(SEE)16 b(ALSO)141
1760 y Fj(qmatrix\(3LAS\),)11 b(precond\(3LAS\))-9 2594
y(6)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: 7 7
7 6 bop -9 -260 a Fj(ERRHANDL\(3LAS\))249 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)250 b(ERRHANDL\(3LAS\))-9 -67 y
Fi(NAME)141 -17 y Fh(LASError)p Fj(,)11 b Fh(LASBreak)p
Fj(,)h Fh(LASResult)p Fj(,)f Fh(WriteLASErrDescr)g Fj({)j(error)g
(handling)f(routines)-9 82 y Fi(SYNOPSIS)141 143 y Fh(#include)20
b(<laspack/errhand)o(l.h>)141 242 y(typedef)g(enum)h({)228
292 y(LASOK,)228 342 y(LASMemAllocErr,)228 392 y(LASLValErr,)228
441 y(LASDimErr,)228 491 y(LASRangeErr,)228 541 y(LASSymStorErr,)228
591 y(LASMatrCombErr,)228 641 y(LASMulInvErr,)228 690
y(LASElNotSortedEr)o(r,)228 740 y(LASZeroInDiagErr)o(,)228
790 y(LASZeroPivotErr,)228 840 y(LASILUStructErr,)228
890 y(LASBreakdownErr,)228 940 y(LASUserBreak)141 989
y(})g(LASErrType;)141 1089 y(void)g(LASError\(LASErr)o(IdTyp)o(e)e
(ErrId,)i(char)f(*ProcName,)g(char)h(*Object1Name,)446
1139 y(char)g(*Object2Name,)e(char)i(*Object3Name\);)141
1189 y(void)g(LASBreak\(void\);)141 1238 y(LASErrIdType)e
(LASResult\(void\);)141 1288 y(void)i(WriteLASErrDesc)o(r\(FIL)o(E)e
(*File\);)-9 2045 y Fi(DESCRIPTION)141 2095 y Fh(LASPack)e
Fj(incorp)q(orates)k(a)d(system)i(for)e(error)j(handling.)32
b(In)20 b(most)e(of)g(the)i(routines)g(of)f(the)h(library)m(,)f(the)141
2145 y(consistency)h(of)f(passed)h(parameters)f(as)h(w)o(ell)e(as)h
(the)h(app)q(earance)h(of)d(some)h(exceptional)g(situations)g(are)141
2195 y(c)o(hec)o(k)o(ed.)141 2294 y(The)14 b(pro)q(cedure)i
Fh(LASError)c Fj(is)i(in)o(tended)h(for)f(in)o(ternal)f(usage)i(in)f
Fh(LASPack)n Fj(.)19 b(In)14 b(case)h(of)e(errors,)i(it)f(transmits)141
2344 y(information)g(ab)q(out)j(their)g(cause)i(to)e(the)g(error)h
(handling)e(system)h(whic)o(h)g(prev)o(en)o(ts)i(further)f(execution)g
(of)141 2394 y Fh(LASPack)e Fj(routines.)29 b(As)18 b(parameter)f
Fh(ErrId)p Fj(,)g(one)h(of)e(the)j(follo)o(wing)14 b(sym)o(b)q(olic)i
(v)n(alues)h(dep)q(ending)h(on)g(the)141 2444 y(reason)c(of)f(the)i
(exceptional)f(situation)f(is)h(used:)-9 2594 y(V)m(ersion:)k(1.12.3)
410 b(Last)14 b(c)o(hange:)19 b(Jan)o(uary)14 b(17,)e(1996)661
b(7)p eop
%%Page: 8 8
8 7 bop -9 -260 a Fj(ERRHANDL\(3LAS\))249 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)250 b(ERRHANDL\(3LAS\))219 -61 y
Fh(LASMemAllocErr)112 b Fj(Not)14 b(enough)g(memory)d(is)j(a)o(v)n
(ailable.)219 -11 y Fh(LASLValErr)200 b Fj(L-v)n(alue)13
b(parameter)g(is)h(exp)q(ected.)219 39 y Fh(LASDimErr)222
b Fj(Ob)r(jects)15 b(ha)o(v)o(e)f(incompatible)e(dimensions.)219
89 y Fh(LASRangeErr)178 b Fj(Indices)15 b(are)f(out)g(of)f(range.)219
139 y Fh(LASSymStorErr)134 b Fj(Some)13 b(elemen)o(ts)g(are)i(stored)f
(in)g(the)g(lo)o(w)o(er)g(triangular)f(part)639 188 y(of)g(a)h
(symmetric)e(matrix.)219 238 y Fh(LASMatrCombErr)112
b Fj(Matrices)15 b(cannot)f(b)q(e)g(com)o(bined.)219
288 y Fh(LASMulInvErr)156 b Fj(In)o(v)o(erse)15 b(m)o(ultiplicatio)o(n)
c(cannot)j(b)q(e)h(carried)f(out.)219 338 y Fh(LASElNotSortedE)o(rr)47
b Fj(Matrix)13 b(elemen)o(ts)h(are)g(not)g(sorted.)219
388 y Fh(LASZeroInDiagEr)o(r)69 b Fj(Zero)15 b(elemen)o(ts)e(in)h(the)g
(matrix)e(diagonal)g(are)j(not)e(allo)o(w)o(ed.)219 438
y Fh(LASZeroPivotErr)90 b Fj(F)m(actorization)13 b(pro)q(duces)i(zero)g
(piv)o(ot)e(elemen)o(ts.)219 487 y Fh(LASILUStructErr)90
b Fj(Matrix)13 b(has)h(a)g(structure)i(whic)o(h)e(is)g(not)f(allo)o(w)o
(ed)g(for)g(ILU)h(factorization.)219 537 y Fh(LASBreakdownErr)90
b Fj(Iterativ)o(e)14 b(solv)o(er)g(fails.)141 652 y(In)h(addition,)f
(the)h(name)f(of)h(the)g(routine)h(where)g(the)g(error)g(has)f(o)q
(ccurred,)i Fh(ProcName)p Fj(,)c(and)i(names)g(of)f(con-)141
702 y(cerning)g(ob)r(jects,)g Fh(Object1Name)p Fj(,)d
Fh(Object2Name)g Fj(and)j Fh(Object3Name)d Fj(\(e.g.)i(names)g(of)g(v)o
(ector)h(or)g(matrices\),)141 752 y(should)f(b)q(e)i(passed,)f(to)q(o.)
k(The)c(un)o(used)h(parameters)f(should)g(b)q(e)g(set)h(to)f
Fh(NULL)p Fj(.)141 851 y(F)m(rom)j(an)h(application,)h(the)g(pro)q
(cessing)h(of)e Fh(LASPack)g Fj(routines)h(ma)o(y)e(b)o(y)i(terminated)
f(b)o(y)h(the)g(pro)q(cedure)141 901 y Fh(LASBreak)12
b Fj(whic)o(h)h(set)i(the)g(error)g(status)f(to)g Fh(LASUserBreak)p
Fj(.)141 1001 y(The)j(curren)o(t)i(error)f(status)g(can)g(b)q(e)g(got)f
(b)o(y)g(the)h(pro)q(cedure)h Fh(LASResult)p Fj(.)27
b(If)17 b(all)f Fh(LASPack)f Fj(routines)j(ha)o(v)o(e)141
1051 y(b)q(een)d(runs)g(faultless,)f(the)h(usual)f(v)n(alue)g
Fh(LASOK)f Fj(b)o(y)h(whic)o(h)h(the)g(library)e(is)i(initialized)e(is)
h(returned.)21 b(In)15 b(other)141 1100 y(cases,)f(the)h(pro)q(cedure)h
Fh(LASResult)c Fj(returns)j(one)f(of)f(the)i(sym)o(b)q(olic)d(v)n
(alues)i(listed)f(ab)q(o)o(v)o(e.)141 1200 y(V)m(erbal)22
b(description)h(the)g(the)g(error)h(o)q(ccasion)f(can)g(b)q(e)g(put)g
(out)f(in)g(the)h(\014le)g Fh(File)f Fj(b)o(y)g(the)h(pro)q(cedure)141
1250 y Fh(WriteLASErrDesc)o(r)p Fj(.)-9 1353 y Fi(FILES)141
1394 y(errhandl.h)530 b Fj(header)14 b(\014le)141 1444
y Fi(errhandl.c)536 b Fj(source)15 b(\014le)-9 1547 y
Fi(EXAMPLES)141 1597 y Fj(The)f(co)q(de)h(fragmen)o(t)141
1697 y Fh(Vector)20 b(a,)h(b;)141 1796 y(V_Constr\(&a,)e("a",)i(6,)g
(Normal,)g(True\);)141 1846 y(V_Constr\(&b,)e("b",)i(7,)g(Normal,)g
(True\);)141 1946 y(V_SetAllCmp\(&a,)d(0.0\);)141 2045
y(Asgn_VV\(&b,)h(&a\);)141 2145 y(if)i(\(LASResult\(\))e(!=)j(LASOK\))e
({)228 2195 y(printf\("LASPack)e(error:)j("\);)228 2245
y(WriteLASErrDescr)o(\(stdi)o(n\);)141 2294 y(})141 2394
y(V_Destr\(&a\);)141 2444 y(V_Destr\(&b\);)-9 2594 y
Fj(8)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: 9 9
9 8 bop -9 -260 a Fj(ERRHANDL\(3LAS\))249 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)250 b(ERRHANDL\(3LAS\))141 -67 y(pro)q(duces)15
b(the)g(follo)o(wing)c(error)k(message:)141 66 y Fh(LASPack)20
b(error:)g(in)i(Asgn_VV)e(for)h(a,)h(b:)141 116 y(Objects)e(have)h
(incompatible)e(dimensions.)-9 2594 y Fj(V)m(ersion:)f(1.12.3)410
b(Last)14 b(c)o(hange:)19 b(Jan)o(uary)14 b(17,)e(1996)661
b(9)p eop
%%Page: 10 10
10 9 bop -9 -260 a Fj(F)-5 b(A)o(CTOR\(3LAS\))315 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)315 b(F)-5 b(A)o(CTOR\(3LAS\))-9
-67 y Fi(NAME)141 -17 y Fh(ILUFactor)12 b Fj({)h(incomplete)g
(factorization)g(of)g(quadratic)h(matrices)-9 82 y Fi(SYNOPSIS)141
143 y Fh(#include)20 b(<laspack/factor.)o(h>)141 242
y(QMatrix)g(*ILUFactor\(QMatri)o(x)f(*Q\);)-9 357 y Fi(DESCRIPTION)141
407 y Fj(A)o(t)c(the)g(\014rst)h(call)e(of)g(the)i(pro)q(cedure)g
Fh(ILUFactor)d Fj(for)i(a)g(giv)o(en)f(matrix)f Fb(Q)p
Fj(,the)i(incomplete)f(factorization)g(of)141 456 y(the)g(matrix)e(is)i
(carried)g(out.)k(F)m(or)c(non-symmetric)e(matrices,)h(the)h(ILU)g
(factorization)f(with)740 547 y Fb(Q)f Fj(=)f(\()p Fb(D)g
Fj(+)f Fb(L)p Fj(\))i Fb(D)1022 530 y Fd(\000)p Fc(1)1078
547 y Fj(\()p Fb(D)f Fj(+)f Fb(U)5 b Fj(\))k(+)h Fb(R)141
638 y Fj(is)20 b(p)q(erformed,)g(where)i Fb(D)q Fj(,)f
Fb(U)5 b Fj(,)21 b(and)f Fb(L)h Fj(are)f(certain)h(diagonal,)e(upp)q
(er,)j(and)e(lo)o(w)o(er)g(triangular)f(matrices,)141
687 y(resp)q(ectiv)o(ely)m(.)j(The)15 b(remainder)g(matrix)e
Fb(R)i Fj(con)o(tains)g(\014ll)f(elemen)o(ts,)g(whic)o(h)h(ha)o(v)o(e)g
(b)q(een)h(ignored)f(during)g(the)141 737 y(factorization)10
b(pro)q(cess.)19 b(F)m(or)11 b(symmetric)f(matrices,)h(the)h
(incomplete)e(Cholesky)i(factorization)e(with)i Fb(L)f
Fj(=)h Fb(U)1886 722 y Fe(T)141 787 y Fj(is)j(applied.)22
b(In)15 b(order)h(to)f(b)q(e)h(able)g(to)f(p)q(erform)f(the)i
(incomplete)f(factorization)f(for)h(singular)g(matrices)g(to)q(o,)141
837 y(regularization)e(b)o(y)g(means)g(of)h(increasing)g(of)f(diagonal)
f(en)o(tries)j(is)f(applied.)141 936 y(In)e(the)g(curren)o(t)i(implem)o
(en)o(tation,)9 b(matrices)j Fb(L)g Fj(and)g Fb(U)k Fj(are)d(laid)e(do)
o(wn)g(with)h(the)h(same)e(p)q(osition)g(of)h(non-zero)141
986 y(elemen)o(ts)h(as)h(in)g(matrix)e Fb(Q)p Fj(.)141
1086 y(Matrices)j(are)g(stored)g Fb(D)q Fj(,)f Fb(U)5
b Fj(,)14 b(and)g Fb(L)h Fj(in)f(connection)h(with)f
Fb(Q)p Fj(as)g(matrix)e Fb(D)f Fj(+)f Fb(L)g Fj(+)g Fb(U)19
b Fj(of)13 b(the)i(t)o(yp)q(e)g Fh(QMatrix)p Fj(.)141
1136 y(This)j(one)g(is)g(also)g(the)h(one)f(returned)i(b)o(y)f(the)f
(pro)q(cedure)j Fh(ILUFactor)p Fj(.)29 b(Ho)o(w)18 b(to)g(extract)h
(the)g(particular)141 1186 y(matrices)13 b(from)f(it,)h(is)h(sho)o(wn)g
(in)f(the)i(example)d(b)q(elo)o(w.)-9 1285 y Fi(REFERENCES)141
1335 y Fj(The)i(incomplete)f(factorization)g(is)h(comprehensiv)o(ely)f
(describ)q(ed)j(and)e(analyzed)f(e.g.)g(in:)244 1435
y(W.)19 b(Hac)o(kbusc)o(h:)30 b(Iterativ)o(e)21 b(Solution)d(of)h
(Large)h(Sparse)g(Systems)g(of)f(Equations,)h(Springer-)244
1484 y(V)m(erlag,)13 b(Berlin,)h(1994.)-9 1617 y Fi(FILES)141
1658 y(factor.h)588 b Fj(header)14 b(\014le)141 1708
y Fi(factor.c)594 b Fj(source)15 b(\014le)-9 1808 y Fi(EXAMPLES)141
1858 y Fj(The)g(follo)o(wing)d(example)i(sho)o(ws)h(the)h(usage)f(of)g
(the)g(pro)q(cedure)i Fh(ILUFactor)c Fj(in)i(the)h(implem)o(en)o
(tation)c(of)i(an)141 1907 y(ILU)f(preconditioner.)19
b(This)14 b(has)g(to)g(solv)o(e)g(the)g(system)g(of)f(equations)953
1998 y Fb(W)g(y)g Fj(=)f Fb(c)141 2089 y Fj(whic)o(h)h(arises)i(during)
e(the)i(solution)e(of)g(preconditioned)i(systems)866
2179 y Fb(W)911 2162 y Fd(\000)p Fc(1)955 2179 y Fb(A)7
b(x)12 b Fj(=)g Fb(W)1118 2162 y Fd(\000)p Fc(1)1169
2179 y Fb(b)141 2270 y Fj(at)h(whic)o(h)770 2320 y Fb(W)k
Fj(=)12 b(\()p Fb(D)f Fj(+)e Fb(L)p Fj(\))j Fb(D)1063
2303 y Fd(\000)p Fc(1)1120 2320 y Fj(\()p Fb(D)f Fj(+)e
Fb(U)c Fj(\))p Fb(:)141 2444 y Fj(The)14 b(corresp)q(onding)h
Fh(LASPack)d Fj(routine)i(could)g(b)q(e)g(built)g(as)g(follo)o(ws:)-9
2594 y(10)643 b(Last)14 b(c)o(hange:)f(Jan)o(uary)h(17,)f(1996)419
b(V)m(ersion:)13 b(1.12.3)p eop
%%Page: 11 11
11 10 bop -9 -260 a Fj(F)-5 b(A)o(CTOR\(3LAS\))315 b(LASP)m(A)o(CK)14
b(REFERENCE)g(MANUAL)315 b(F)-5 b(A)o(CTOR\(3LAS\))141
-67 y Fh(Vector)20 b(*ILUPrecond\(QMatri)o(x)f(*A,)i(Vector)g(*y,)g
(Vector)f(*c,)i(double)e(Omega\))141 -17 y({)228 33 y(Q_Lock\(A\);)228
83 y(V_Lock\(y\);)228 133 y(V_Lock\(c\);)228 232 y(Asgn_VV\(y,)f
(MulInv_QV\(Add_QQ\(Dia)o(g_Q\()o(ILUFa)o(ctor\()o(A\)\),)f
(Upper_Q\(ILUFactor\(A)o(\)\)\),)315 282 y(Mul_QV\(Diag_Q\(ILU)o(Facto)
o(r\(A\)\))o(,)315 332 y(MulInv_QV\(Add_QQ\()o(Diag_)o(Q\(ILU)o(Fact)o
(or\(A\))o(\),)h(Lower_Q\(ILUFactor)o(\(A\)\))o(\),)g(c\)\)\)\);)228
431 y(Q_Unlock\(A\);)228 481 y(V_Unlock\(y\);)228 531
y(V_Unlock\(c\);)228 631 y(return\(y\);)141 681 y(})-9
813 y Fi(SEE)d(ALSO)141 863 y Fj(qmatrix\(3LAS\),)11
b(op)q(erats\(3LAS\),)k(errhandl\(3LAS\))-9 963 y Fi(BUGS)141
1013 y Fj(In)h(the)g(curren)o(t)i(implemen)o(tatio)o(n,)c(it)h(is)h
(assumed)g(that)g(the)h(matrix)d Fb(Q)i Fj(has)g(a)g(symmetric)e
(structure)k(with)141 1062 y(regard)i(to)f(non-zero)i(elemen)o(ts.)35
b(Because)21 b(during)f(discretization)g(of)f(di\013eren)o(tial)g
(equations)h(ev)o(en)g(suc)o(h)141 1112 y(matrices)12
b(arise,)g(this)h(restriction)g(is)g(not)f(gra)o(v)o(e)h(for)f(man)o(y)
e(applications.)17 b(F)m(or)12 b(matrices)g(whic)o(h)h(did)f(not)g(ha)o
(v)o(e)141 1162 y(satisfy)h(the)i(ab)q(o)o(v)o(e)e(condition,)g(the)h
(error)h Fh(LASILUStructErr)c Fj(is)j(raised.)-9 2594
y(V)m(ersion:)k(1.12.3)410 b(Last)14 b(c)o(hange:)19
b(Jan)o(uary)14 b(17,)e(1996)641 b(11)p eop
%%Page: 12 12
12 11 bop -9 -260 a Fj(ITERSOL)-5 b(V)q(\(3LAS\))278
b(LASP)m