ibm.ps

fortran版本的eslack,好不容易找到的。里面有程序的索引。主要包括了常用了矩阵计算代码。

[<
007FFFFFE000007FFFFFE000007FFFFFE00000003FC0000000003FC0000000003FC0000000
003FC0000000003FC0000000003FC0000000003FC0000000003FC0000000003FC000000000
3FC0000000003FC0000000003FC0000000003FC0000000003FC0000000003FC0000000003F
C0000000003FC0000000003FC0000000003FC0000000003FC0000000003FC0000000003FC0
0000E0003FC000E0E0003FC000E0E0003FC000E0E0003FC000E0E0003FC000E0F0003FC001
E0F0003FC001E070003FC001C078003FC003C078003FC003C07E003FC007C07F803FC03FC0
7FFFFFFFFFC07FFFFFFFFFC07FFFFFFFFFC0>
	 43 40 -2 39 48] 84 @dc
[<
FFFEFFFEFFFE0FE00FE00FE00FE00FE00FE00FE00FE00FE00FE00FE00FE00FE00FE00FE00F
E00FE00FE00FE00FE00FE0FFE0FFE0FFE0000000000000000000000000000007001FC01FE0
3FE03FE03FE01FE01FC00700>
	 15 43 -3 42 19] 105 @dc
[<
007FF00003FFFE000FC01F803F0007E07E0003F0FC0001F8F80000F8F80000F8F80000F8F8
0001F87C0001F87E0007F81FFFFFF007FFFFF00FFFFFE01FFFFFC03FFFFF003FFFF8003E00
00003C0000003C0000003800000018FF80001FFFE0000FC1F8001F80FC003F80FE003F007E
007F007F007F007F007F007F007F007F007F007F007F007F003F007E103F80FE7C1F80FC7C
0FC1FE7C03FFE7FC00FF81F0>
	 30 40 -2 26 34] 103 @dc
@F10 @sf
[<
FFFFF0000F803C0007801F0007800F0007800F8007800780078007C0078007C0078007C007
8007C0078007C00780078007800F8007800F0007801E0007803C0007FFF00007803C000780
1E0007801F0007800F0007800F8007800F8007800F8007800F8007800F8007800F0007801F
0007801E000F803800FFFFE000>
	 26 31 -2 30 32] 66 @dc
[<
3F006180F0C0F0606070007000700070007000700070007000700070007000700070007000
70007000700070007000700070007000F007F0007000000000000000000000000000E001F0
01F001F000E0>
	 12 40 2 30 14] 106 @dc
[<
FFFFFF800F800F8007800180078000C0078000C00780004007800040078000400780002007
80002007800020078080000780800007808000078080000781800007FF8000078180000780
80000780800007808000078080800780008007800080078000800780018007800100078003
00078003000F800F00FFFFFF00>
	 27 31 -2 30 31] 69 @dc
[<
FFFC0FC0078007800780078007800780078007800780078007800780078007800780078007
8007800780078007800780078007800780078007800FC0FFFC>
	 14 31 -1 30 16] 73 @dc
[<
FFFC3FFC0FC007E0078007C0078007800780078007800F0007801E0007801E0007803C0007
803C00078078000780F00007C0F00007A1E0000791E0000793C000078F8000078780000783
00000781000007808000078040000780200007801000078008000780040007800200078001
00078001800FC003E0FFFC0FFC>
	 30 31 -2 30 35] 75 @dc
[<
C000C000C0006000600060003000300030001800180018000C000C000C0006000600060003
000300030001800180018000C000C000C0006000600060003000300030001800180018000C
000C000C000600060006000300030001>
	 16 45 -3 33 23] 47 @dc
[<
FFFE00000FC000000780000007800000078000000780000007800000078000000780000007
800000078080000780800007808000078080000781800007FF800007818000078080000780
80000780800007808000078000800780008007800080078000800780018007800100078003
00078003000F800F00FFFFFF00>
	 25 31 -2 30 30] 70 @dc
[<
7FE7FE0700E00700E00700E00700E00700E00700E00700E00700E00700E00700E00700E007
00E00700E00700E00700E00700E00700E00700E0FFFFE00700E00700E00700E00700E00700
E00700E00700E00700E00381E001C1E000E0E0003FE0>
	 23 32 0 31 25] 13 @dc
[<
03F0001C3C00200E00400F00400780F00780F807C0F807C0F807C02007C00007C000078000
0780000F00000E00003C0003F000003800001C00000E00000F00000F00000F80380F807807
80780780780F80200F00100E000C1C0003F000>
	 18 31 -2 29 23] 51 @dc
[<
00FFE0000E00000E00000E00000E00000E00000E00000E00FFFFF0C00E00400E00200E0020
0E00100E00080E00080E00040E00020E00020E00010E00008E00008E00004E00002E00002E
00001E00000E00000E00000600000600>
	 20 30 -1 29 23] 52 @dc
[<
FFCFF80E03E00E03C00E03800E07800E07000E0E000E1E000E1C000F38000EF8000E70000E
30000E10000E08000E04000E02000E03000E03C00E0FF00E00000E00000E00000E00000E00
000E00000E00000E00000E00000E0000FE00000E0000>
	 21 32 -1 31 24] 107 @dc
[<
FFE01FFE1F8007E0070003C0020007C00100078001800F0000801F0000401E0000603C0000
207C00001078000018F8000009F0000005E0000007E0000003C00000078000000F8000000F
8000001E4000003E4000003C2000007C10000078100000F0080001F00C0001E0040003C002
0007C007000FE00FC07FF83FF8>
	 31 31 -1 30 34] 88 @dc
[<
03F0000C0C00100200200100600180C00080C000C0C000C0C000C0C001C04001C060038030
0F80103F000C7E0007FC0003F8000FF0001FC8003F06003E03007803007001806001806001
806001802001803003001006000C0C0003F000>
	 18 31 -2 29 23] 56 @dc
[<
7FE3FF00070070000700700007007000070070000700700007007000070070000700700007
00700007007000070070000700700007007000070070000700700007007000070070000700
7000FFFFFF80070070000700700007007000070070000700700007007000070070000300F0
300380F87801C0787800F06E30001F83E0>
	 29 32 0 31 27] 11 @dc
[<
7FE3FE3FF00700700700070070070007007007000700700700070070070007007007000700
70070007007007000700700700070070070007007007000700700700070070070007007007
000700700700070070070007007007000700700700FFFFFFFF000700700000070070000007
007000000700700000070070000007007000000700F00F000300F00F000380F80F0001C07C
060000F04F0400001F81F800>
	 36 32 0 31 38] 14 @dc

18 @eop0

0 0 18 @bop1 59 192 a @F10 @sf(The)10 b(b)q(est)h(examples)f(of)f(the)i
(latter)e(tec)o(hnique)h(are)g @F5 @sf(elmbak)f @F10 @sf(and)i @F5 @sf
(combak)p @F10 @sf(.)17 b(With)9 b(only)h(v)o(ery)g(minor)-1711 70 y(mo)q
(di\014cations)h(and)h(the)f(application)g(of)g(a)g(few)g(v)o(ector)g
(directiv)o(es,)h(the)f(co)q(de)h(ran)g(appro)o(ximately)-1712 71 y(6)h
(and)h(7)f(\(resp)q(ectiv)o(ely\))g(times)f(faster)h(for)g @F9 @sf(N)k 
@F10 @sf(=)c(200.)19 b(Notice)13 b(that)f(the)i(imp)q(ortan)o(t)e
(subroutines)-1710 70 y @F5 @sf(hqr)j @F10 @sf(and)g @F5 @sf(hqr2)f @F10 @sf
(p)q(erformed)h(70%)f(and)h(90%,)f(resp)q(ectiv)o(ely)l(,)h(b)q(etter)g
(at)f(dimension)h(100)f(with)g(the)-1710 71 y(PVS)20 b(v)o(ersions)e(than)h
(the)g(with)g(NA)l(TS.)g(Alternativ)o(ely)l(,)g(PVS)g @F5 @sf(orthes)f 
@F10 @sf(p)q(erformed)i(only)f(20%)-1711 70 y(b)q(etter)13 b(at)g
(dimension)g(100)f(and)i @F5 @sf(ortran)e @F10 @sf(p)q(erformed)h(no)h(b)q
(etter)f(than)g(its)f(NA)l(TS)i(v)o(ersion.)k(Since)-1709 71 y(impro)o(v)o
(emen)o(ts)c(for)h(small)f(dimension)h(matrices)g(w)o(ere)g(realized)g
(with)g(P)l(ASC)h(double)f(precision)-1710 70 y(v)o(ersions)g(of)g @F5 @sf
(eltran)p @F10 @sf(,)f @F5 @sf(htribk)p @F10 @sf(,)g(and)i @F5 @sf(ortran)f 
@F10 @sf(but)g(not)g(with)g(the)h(PVS)f(v)o(ersions,)g(w)o(e)g(are)g(led)h
(to)-1711 71 y(b)q(eliev)o(e)i(that)g(the)g(signi\014can)o(t)g(stride-one)g
(reorien)o(tati)o(on)f(that)g(w)o(as)h(emplo)o(y)o(ed)g(for)f(the)i(P)l
(ASC)-1710 71 y(co)q(des)d(are)f(the)g(correct)g(approac)o(h)g(to)f(emplo)o
(y)h(ev)o(en)g(for)g(single)f(precision.)-1324 162 y @F1 @sf(7)69 b
(Correctness)22 b(T)-6 b(esting)-698 115 y @F10 @sf(Both)15 b(P)l(ASC)h
(and)f(PVS)h(co)q(des)g(w)o(ere)f(sub)s(jected)h(to)e(tests)h(to)g(ensure)h
(their)e(correctness.)21 b(These)-1710 71 y(tests)16 b(w)o(ere)g(created)g
(b)o(y)g(mo)q(difying)f(the)h(thirteen)g(NA)l(TS-supplied)h(driv)o(ers.)
23 b(These)16 b(programs)-1711 70 y(originall)o(y)f(read)h(a)f(large)g
(set)h(of)g(test)f(matrices,)g(called)g(a)h(sequence)h(of)f(EISP)l(A)o(CK)h
(mo)q(dules)f(to)-1711 71 y(pro)q(duce)11 b(eigen)o(v)m(alues)g(and/or)f
(eigen)o(v)o(ectors,)g(compared)g(results)g(across)g(alternativ)o(e)f
(paths,)i(and)-1710 70 y(computed)i(residuals.)19 b(Eac)o(h)12 b(test)h
(pro)q(duces)g(a)g(\014gure)g(of)f(merit)g(that)g(essen)o(tially)f(is)i
(the)g(residual)-1711 71 y(norm)k(scaled)g(b)o(y)g(the)g(the)h(pro)q(duct)f
(of)g(ten,)g(the)h(mac)o(hine)f(precision,)g(the)g(matrix)f(dimension,)-1711 
70 y(and)j(the)h(matrix)d(norm.)32 b(These)19 b(tests)g(ha)o(v)o(e)g(b)q
(een)h(emplo)o(y)o(ed)f(since)g(1972)f(to)h(determine)g(the)-1710 71 y
(correctness)c(of)g(v)o(ersions)f(of)h(EISP)l(A)o(CK)h(on)f(v)m(arious)g
(mac)o(hines.)-1067 71 y(The)d(mo)q(di\014cations)f(emplo)o(y)o(ed)g(allo)o
(w)o(ed)g(for)g(a)h(set)g(of)f(random)h(p)q(erturbations)f(to)g(b)q(e)i
(applied)-1710 70 y(to)j(eac)o(h)h(of)g(the)g(input)g(matrices.)24 b(F)l
(or)16 b(certain)h(situations)e(\(e.)i(g.)f(symmetric)g(matrices,)g
(tridi-)-1711 71 y(agonal)f(matrices\))f(the)i(p)q(erturbations)f(had)h
(to)f(re\015ect)h(the)g(c)o(haracter)f(of)g(the)h(original)o 15 x(matrix.)
-1712 70 y(A)f(large)g(and)g(v)m(aried)g(set)g(of)g(tests)f(could)i(b)q(e)f
(p)q(erformed)h(with)e(this)h(pro)q(cedure.)-1343 71 y(Although)h(the)h
(co)q(des)g(w)o(ere)g(targeted)f(for)g(the)h(3090-VF)e(they)i(w)o(ere)g
(also)e(run)i(\(emplo)o(ying)-1711 70 y(sim)o(ulatio)o(n\))f(on)i(Sun)g(3)f
(and)h(Sun)g(4)f(w)o(orkstat)o(ions)f(and)i(on)f(a)g(V)-5 b(AX)18 b
(11/780.)o 26 x(T)l(esting)e(on)i(the)-1710 71 y(alternativ)o(e)c
(equipmen)o(t)h(with)g(di\013eren)o(t)g(\015oating)g(p)q(oin)o(t)g
(precisions)g(and)h(ranges)f(as)g(w)o(ell)g(as)g(on)-1710 71 y(the)g
(3090-VF)e(w)o(as)h(in)o(tended)h(to)e(rev)o(eal)h(subtle)h(errors.)j(In)e
(fact,)d(they)i(did)f(since)h(sev)o(eral)f(of)g(the)-1710 70 y(original)k
(P)l(ASC)i(co)q(des)g(\(since)f(reco)q(ded\))h(sho)o(w)o(ed)g(di\016culties)e
(with)h(o)o(v)o(er\015o)o(ws)g(that)g(had)g(not)-1710 71 y(previously)c(b)q
(een)h(unco)o(v)o(ered.)-465 70 y(The)k(\014nal)f(v)o(ersions)g(of)h(all)f
(P)l(ASC)h(and)g(PVS)g(routines)f(p)q(erform)h(within)f(the)h(acceptable)
-1710 71 y(tolerances)15 b(on)g(all)f(tests.)388 319 y(18)-937 b
18 @eop1

17 @bop0
@F10 @sf
[<
03E0000C3800100E00200600400700400380E00380F003C0F003C07003C00003C00003C000
03C0000380000380100780100700180E00161C0011F0001000001000001000001000001000
001000001FE0001FF8001FFC001FFE00180300>
	 18 31 -2 29 23] 53 @dc
@F5 @sf
[<
FFFFF0FFFFF0FFFFF07800703C00701E00700F007007800003C00001F00000F800003C0000
1E00000F007007807003C07001E07FFFF07FFFF03FFFF0>
	 20 20 -1 19 24] 122 @dc
@F10 @sf
[<
FFF8FFF80F800F8007800F0007800F0007800F0007800F0007800F0007800F0007800F0007
800F0007800F0007800F0007800F0007800F0007800F0007800F0007FFFF0007800F000780
0F0007800F0007800F0007800F0007800F0007800F0007800F0007800F0007800F0007800F
0007800F000F800F80FFF8FFF8>
	 29 31 -2 30 34] 72 @dc
@F5 @sf
[<
7FFCFFFE7FFC03800380038003800380038003800380038003800380038003800380038003
804380FB80FF803F800F800780038003800180>
	 15 28 -5 27 24] 49 @dc
@F5 @sf
[<
00700000F80000F80001DC0001DC0001DC00038E00038E00038E00038E0007070007070007
07000E03800E03800E03801E03C07F8FF0FF8FF87F8FF0>
	 21 20 -1 19 24] 118 @dc
/@F3 @newfont
@F3 @sf
[<
00FF0003FFC00783E00F00F01E00F83E007C3C007C7C007E7C007E7C007EFC007EFC007EFC
007EFC007EFE007CFE007CFF00F8FD80F0FCFFE0FC7FC0FC08007C00007C00007E00F03E01
F83E01F81F01F80F01F807C0F003F07000FFE0001F80>
	 23 32 -2 31 28] 54 @dc
[<
387CFEFEFE7C38>
	 7 7 -4 6 16] 46 @dc
[<
FFFFF0FFFFF07FFFF03FFFF01FFFF00FFFF00E003807001803801801C01800E01800700000
3800003E00001F00000F800007C00007E00003F00003F00003F83803F87C01F8FE01F8FE03
F8FE03F8FE03F07C07F07007E03C1FC00FFF0003FC00>
	 21 32 -3 31 28] 50 @dc
[<
FFFF8000FFFF800007F0000007F0000007F0000007F0000007F0000007F0000007F0000007
F0000007F0000007F0000007F0000007F0000007F0000007FFFF0007FFFFE007F007F007F0
01FC07F000FC07F0007E07F0007E07F0007F07F0007F07F0007F07F0007F07F0007F07F000
7E07F0007E07F000FC07F001FC07F007F0FFFFFFE0FFFFFF00>
	 32 34 -2 33 38] 80 @dc
[<
0000380000000038000000007C000000007C00000000FE00000000FE00000000FE00000001
FF00000001FF00000003FD80000003F980000007F9C0000007F0C0000007F0C000000FF060
00000FE06000001FE03000001FC03000003FC03800003F801800003F801800007F800C0000
7F000C0000FF00060000FE00060001FE00070001FC00030001FC00030003F800018003F800
018007F80000C007F00000C0FFFF800FFEFFFF800FFE>
	 39 34 -1 33 42] 86 @dc
[<
80FF80C7FFE0FF00F8FC003CF0003CE0001EE0001EC0001FC0001FC0001F00003F00003F00
007F0003FF003FFE03FFFE0FFFFC1FFFF83FFFF07FFFC07FFF00FFE000FF0000FC0000FC00
0CF8000CF8000C78001C78001C7C003C3C007C1F03FC07FF8C01FC04>
	 24 34 -3 33 31] 83 @dc
[<
000E0000000E0000001F0000001F0000003F8000003F8000007FC000007CC000007CC00000
F8600000F8600001F0300001F0300003E0180003E0180007E01C0007C00C000FC00E000F80
06000F800600FFE01FE0FFE01FE0>
	 27 22 -1 21 30] 118 @dc
[<
C7F000EFFC00FC0E00F00700E00380E00380C00380C00780000F8001FF000FFF003FFE007F
F8007FE000FC0000F00300E00300E00300600700781F003FFF000FF300>
	 17 22 -2 21 22] 115 @dc
[<
FFF00000C0FFF00001C006000003C006000003C006000007C00600000FC00600001FC00600
003FC00600007FC00600007FC0060000FFC0060001FEC0060003FCC0060007F8C0060007F0
C006000FF0C006001FE0C006003FC0C006007F80C00600FF00C00600FE00C00601FE00C006
03FC00C00607F800C0060FF000C0060FE000C0061FE000C0063FC000C0067F8000C006FF00
00C007FE0000C007FC0000C0FFFC001FFEFFF8001FFE>
	 39 34 -2 33 44] 78 @dc
[<
FFE00FFFF8FFE00FFFF80600007F000600007F00030000FE00030000FE00038001FE000180
01FC00018001FC0000FFFFF80000FFFFF80000E007F800006007F00000600FF00000300FE0
0000300FE00000381FE00000181FC00000183FC000000C3F8000000C3F8000000E7F800000
067F00000006FF00000003FE00000003FE00000003FE00000001FC00000001FC00000000F8
00000000F800000000F8000000007000000000700000>
	 37 34 -2 33 42] 65 @dc
[<
03FFFFF80003FFFFF8000003F800000003F800000003F800000003F800000003F800000003
F800000003F800000003F800000003F800000003F800000003F800000003F800000003F800
000003F800000003F800000003F800000003F800000003F800000003F800000003F80000C0
03F800C0C003F800C0C003F800C0C003F800C0E003F801C0E003F801C06003F801807003F8
03807803F807807E03F80F807FFFFFFF807FFFFFFF80>
	 34 34 -2 33 39] 84 @dc
[<
00180030006000C001C00380070007000E001E001C003C003C003C007800780078007800F8
00F000F000F000F000F000F000F000F000F000F000F000F80078007800780078003C003C00
3C001C001E000E0007000700038001C000C0006000300018>
	 13 49 -5 36 22] 40 @dc
[<
FFE0FFE01F001F001F001F001F001F001F001F001F001F001F001F001F001F001F001F001F
001F00FF00FF0000000000000000000000000000001C003F007F007F007F003F001C00>
	 11 36 -2 35 16] 105 @dc
[<
FFE1FFC0FFE1FFC01F003E001F003E001F003E001F003E001F003E001F003E001F003E001F
003E001F003E001F003E001F003E001F003E001F003E001F003E001F803E001F803E001F40
3C001F307C00FF1FF800FF07E000>
	 26 22 -2 21 31] 110 @dc
[<
01FF000FFFE03F01F87C007C78003CF0001EF0001EF0001EF0001E70003E3800FC1FFFFC0F
FFF81FFFF03FFF8038000030000030000013FC001FFF001F0F803E07C03C03C07C03E07C03
E07C03E07C03E07C03E03C03C03E07CF1F0F8F0FFF7F03FC1E>
	 24 33 -2 21 28] 103 @dc
[<
FFE0FFE01F001F001F001F001F001F001F001F001F001F001F001F001F001F001F001F001F
001F001F001F001F001F001F001F001F001F001F001F001F001F001F00FF00FF00>
	 11 35 -2 34 16] 108 @dc
[<
00FF0003FFC00FC0701F00303E00187E00007C00007C0000FC0000FC0000FC0000FFFFF8FF
FFF8FC00F8FC00F87C00F87C00F03E01F01E01E00F87C007FF8000FE00>
	 21 22 -2 21 26] 101 @dc
[<
FFE00000FFE000001F0000001F0000001F0000001F0000001F0000001F0000001F0000001F
0000001F0FC0001F3FF8001FE0FC001FC03E001F803F001F001F801F001F801F000FC01F00
0FC01F000FC01F000FC01F000FC01F000FC01F000FC01F000FC01F001F801F001F801F001F
001F803E001FF07C00FF3FF800FF0FE000>
	 26 32 -2 21 31] 112 @dc
[<
FFF000FFF0001F00001F00001F00001F00001F00001F00001F00001F00001F00001F00001F
00001F00001F00001F83C01E87E01E87E01EC7E01E67E0FE3FC0FE1F00>
	 19 22 -2 21 23] 114 @dc
[<
00FE0007FF800FC0E01F00603E00307E00007E00007C0000FC0000FC0000FC0000FC0000FC
0000FC00007C00007C01E07E03F03E03F01F03F00F83F007FFE000FF80>
	 20 22 -2 21 25] 99 @dc
[<
00FE0007FFC00F83E01F01F03E00F87C007C7C007CFC007EFC007EFC007EFC007EFC007EFC
007EFC007E7C007C7C007C7C007C3E00F81E00F00F83E007FFC000FE00>
	 23 22 -2 21 28] 111 @dc
[<
C0006000300018001C000E0007000700038003C001C001E001E001E000F000F000F000F000
F80078007800780078007800780078007800780078007800F800F000F000F000F001E001E0
01E001C003C00380070007000E001C00180030006000C000>
	 13 49 -3 36 22] 41 @dc
@F10 @sf
[<
FEFEC0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0C0
C0C0C0C0C0C0FEFE>
	 7 45 -4 33 13] 91 @dc
/@F8 @newfont
@F8 @sf
[<
FFFFFFC0FFFFFFC0>
	 26 2 -4 12 35] 0 @dc
@F10 @sf
[<
FEFE0606060606060606060606060606060606060606060606060606060606060606060606
060606060606FEFE>
	 7 45 -1 33 13] 93 @dc
[<
001F800000F0F00001C0380007801E000F000F000E0007001E0007803C0003C03C0003C07C
0003E07C0003E0780001E0F80001F0F80001F0F80001F0F80001F0F80001F0F80001F0F800
01F0F80001F0F80001F0780001E0780001E07C0003E03C0003C03C0003C01E0007800E0007
000F000F0007801E0001C0380000F0F000001F8000>
	 28 33 -3 31 35] 79 @dc
[<
000003E0FFFC0F100FC01E0807803E0407807E0407807C0407807C0007807C0007807C0007
8078000780780007807800078070000780F0000780E0000781C00007FF800