code.ps

有限内存BFGS约束优化算法Fortran程序及其说明

%!PS-Adobe-2.0
%%Creator: dvipsk 5.58f Copyright 1986, 1994 Radical Eye Software
%%Title: tomspaper.dvi
%%Pages: 12
%%PageOrder: Ascend
%%BoundingBox: 0 0 612 792
%%DocumentPaperSizes: Letter
%%EndComments
%DVIPSCommandLine: dvips tomspaper
%DVIPSParameters: dpi=600, compressed, comments removed
%DVIPSSource:  TeX output 1996.10.08:0746
%%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 /hsize 8.5 72
mul N /landplus90{false}def /@rigin{isls{[0 landplus90{1 -1}{-1 1}
ifelse 0 0 0]concat}if 72 Resolution div 72 VResolution div neg scale
isls{landplus90{VResolution 72 div vsize mul 0 exch}{Resolution -72 div
hsize mul 0}ifelse TR}if Resolution VResolution vsize -72 div 1 add mul
TR[matrix currentmatrix{dup dup round sub abs 0.00001 lt{round}if}
forall round exch round exch]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 sub]/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 userdict
/eop-hook known{eop-hook}if showpage}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 newpath 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
TeXDict begin 40258431 52099146 1000 600 600 (tomspaper.dvi)
@start /Fa 7 85 df<EC01FE4A7EA34A7FA391380FDFC0159F021F7F158FA291383F0F
F0A21507027E7FA2150302FC7FA201016D7E14F8A201036D7E5CA201076E7E5C010F6E7E
A25C011F6E7E5CA2013F6E7E91C7FC49811607137E90B77EA34882A2D803F8C7EA01FF82
5B00071780177F5B000FEE3FC05B001F17E0171F5B003F17F0170F90C9FC48EE07F8A200
FE17FC17035A363F7DBE3D>65 D<EDFFF0020FEBFF80023F14F091B612FC5B1307499038
803FF890391FFC0003D93FF01300D97FC01438494814184890C9FC485A5B12075B485AA2
485AA2485AA3485AA448CAFCAD6C7EA46C7EA36C7EA26C7EA26C7E7F12037F6C6C15046C
6D140C6D6C143CD93FF0147CD91FFCEB01FE90390FFF801F6D90B5FC010115FC6D15F002
3F14C0020F1400020013F02F437AC03A>67 D<B8FCA690CAFCB3A390B612FCA690CAFCB3
A590B71280A7293F78BE36>69 D<B4FCB3B3B3A390B612F0A6243F78BE31>76
D<D8FFF015FEA27FA27FA27F12FE137FA280133F80131F80130F80130780130380130180
1300808081143F81141F81140F81140781140381140181140081811680153F16C0151F16
E0150F16F0150716F8150316FC1501A2ED00FE16FF167FA2163FA2161FA22F3F78BE40>
78 D<4AB47E020F13F0027F13FE91B6FC010315C04981011F010013F8D93FF8EB1FFCD9
7FE0EB07FE4A130349486D7E4890C813804848ED7FC049153F4848ED1FE04848ED0FF0A2
4848ED07F8A2491503003F17FCA2491501007F17FEA390CAFC4817FFAC6D5D007F17FEA4
6D1503003F17FCA26D1507001F17F86D150F000F17F06D151F6C6CED3FE0A26C6CED7FC0
6C6CEDFF806C6D4913006E5BD97FF0EB0FFE6D6C495A6DB4EBFFF8010790B512E06D5D01
0092C7FC6E5B020F13F00201138038437BC043>I<B912FEA6C890C8FCB3B3B3A2373E7D
BD3E>84 D E /Fb 2 110 df<16F8ED03FEED0F8792381F0F80ED3E3F167F157CA215FC
1700161C4A48C7FCA414035DA414075DA20107B512F0A39026000FE0C7FC5DA4141F5DA4
143F92C8FCA45C147EA514FE5CA413015CA4495AA45C1307A25C121E123F387F8F80A200
FF90C9FC131E12FEEA7C3CEA7878EA1FF0EA07C0294C7CBA29>102
D<D803E0017F14FE3D07F801FFE003FFC03D0E3C0781F00F03E03D1C3E1E00F83C01F026
383F38D9FC707F00304914E04A90387DC000007049EB7F8000604991C7FCA200E090C700
FE1301485A017E5CA200000201140301FE5F495CA203031407000160495C180F03075D12
03494A011F13601980030F023F13E00007F000C0495C1901031F023E1380000F1803494A
150061033F150E001FEF1E1C4991C7EA0FF80007C7000EEC03E043267EA449>109
D E /Fc 1 4 df<1338A50060130C00F8133E00FC137E00FE13FE383FBBF83807FFC000
011300EA007C48B4FC000713C0383FBBF838FE38FE00FC137E00F8133E0060130C000013
00A517197B9A22>3 D E /Fd 6 56 df<B612F0A600FCC8FCB3B3B3B3B3B3B3B01C9466
8137>50 D<B612F0A6C71203B3B3B3B3B3B3B3B01C94808137>I<12FCB3B3B3B3B3B3B3
B0B612F0A61C94668237>I<EC03F0B3B3B3B3B3B3B3B0B6FCA61C94808237>I<12FCB3B3
B00634668037>I<12FCB3B3B006346A8037>I E /Fe 2 50 df<B812C0A32A037A9137>
0 D<D93F80EC07F0D9FFF0EC3FFC000301FC91B5FC4801FF903901F80780D80F80903A80
07C001C03D1E003FC00F8000E0486D6C48C71270003890260FF03E143048902607F83C14
386E6C48141800606D6C5A00E06D6C48141C48027F150C5E153F6F7E6F7E82151F6C6F14
1C0060DA3DFE1418DB78FF143800704A6C7E003001016D6C1370003849486C6C13F06C90
3B07C00FF001E06C903B0F8007FC07C02807807E0003B512806CB44801001400C601F0EC
3FFCD93F80EC07F03E1F7C9D47>49 D E /Ff 5 108 df<147C49B4FC903803C7809038
0783C090381F03E0EB1E01133E017C13F013F8A2EA01F0120313E01207A2EA0FC01403A2
EA1F80A21407003F14E0130090B5FCA2397F000FC0127EA2141F1580127C00FC14005CA2
147EA248137C14FC00785B495AA2387C03E0383C07C0495A001C90C7FCEA1E3EEA0FF8EA
03E01C307DAE21>18 D<1670A216F01501A24B7EA21507150DA2151915391531ED61FC15
6015C0EC0180A2EC03005C14064A7F167E5C5CA25C14E05C4948137F91B6FC5B0106C712
3FA25B131C1318491580161F5B5B120112031207000FED3FC0D8FFF8903807FFFEA22F2F
7DAE35>65 D<000FB8FCA23B1FC003F8003F0100151F001C4A130E123C00380107140612
3000704A130EA20060010F140C12E0485CA2141FC715005DA2143FA292C8FCA25CA2147E
A214FEA25CA21301A25CA21303A25CA21307A25C130F131F001FB512F0A2302D7FAC29>
84 D<0107B612F8A2903A0FFC0007F002E0EB0FE00280EB1FC049C71380011E143F011C
EC7F004914FE4B5A0130495A0170495A0160495A4B5A4B5A90C748C7FCA215FE4A5A4A5A
4A5A4A5A4A5A4A5A4AC8FC14FE5C130149481306495A4948130E4948130C495A49C7121C
01FE141848481438485A5E484814F048481301484813034848495A48C7127FB7FC5E2D2D
7CAC30>90 D<131FEA03FFA2EA003FA2133EA2137EA2137CA213FCA25BA2120115F89038
F003FCEC0F0E0003EB1C1EEC387EEBE07014E03807E1C09038E3803849C7FC13CEEA0FDC
13F8A2EBFF80381F9FE0EB83F0EB01F81300481404150C123EA2007E141C1518007CEBF0
38ECF83000FC1470EC78E048EB3FC00070EB0F801F2F7DAD25>107
D E /Fg 6 108 df<007FB812FEBAFCA26C17FE3804799847>0 D<EB03C0A2805CA600F0
140F00FC143F00FE147F00FF14FF393FC3C3FC390FE187F03903F18FC03900FDBF00EB3F
FCEB0FF0EB03C0EB0FF0EB3FFCEBFDBF3903F18FC0390FE187F0393FC3C3FC39FF03C0FF
00FE147F00FC143F00F0140F00001400A6805CA220277AA92D>3
D<180E183F18FFEF03FEEF0FF8EF3FE0EFFF80933803FE00EE0FF8EE3FE0EEFF80DB03FE
C7FCED0FF8ED7FE0913801FF80DA07FEC8FCEC1FF8EC7FC04948C9FCEB07FCEB1FF0EB7F
C04848CAFCEA07FCEA1FF0EA7FC048CBFCA2EA7FC0EA1FF0EA07FCEA01FF38007FC0EB1F
F0EB07FCEB01FF9038007FC0EC1FF0EC07FE913801FF809138007FE0ED1FF8ED03FE9238
00FF80EE3FE0EE0FF8EE03FE933800FF80EF3FE0EF0FF8EF03FEEF00FF183F180E1800AE
007FB812FEBAFCA26C17FE384879B947>20 D<127012FCB4FCEA7FC0EA1FF0EA07FCEA01
FF38007FC0EB1FF0EB07FCEB01FF9038007FC0EC1FF0EC07FE913801FF809138007FE0ED
1FF8ED03FE923800FF80EE3FE0EE0FF8EE03FE933800FF80EF3FE0EF0FF8EF03FEEF00FF
A2EF03FEEF0FF8EF3FE0EFFF80933803FE00EE0FF8EE3FE0EEFF80DB03FEC7FCED0FF8ED
7FE0913801FF80DA07FEC8FCEC1FF8EC7FC04948C9FCEB07FCEB1FF0EB7FC04848CAFCEA
07FCEA1FF0EA7FC048CBFC12FC1270CCFCAE007FB812FEBAFCA26C17FE384879B947>I<
126012F0B3B3B3B3B11260045B76C319>106 D<0060131800F0133CB3B3B3B3B0006013
18165A75C32D>I E /Fh 23 122 df<15FCEC03FF91380F87C091383E03E0EC7C0102F8
13F01301903903F000F8495A010F14FC5C495A133F91C7FC4914FE13FEA212015B120349
13011207A25B000F15FC1503121F5BA21507003F15F890B6FCA33A7FC0000FF05BA2151F
16E048C7FCA2ED3FC0A2481580157F1600A215FEA24A5AA24A5A007E5C14075D4A5A003E
5C141F4AC7FC6C137E5C380F81F03807C3E03801FF80D8007EC8FC27417DBF2B>18
D<121EEA7F80A2EAFFC0A4EA7F80A2EA1E000A0A798919>58 D<121EEA7F8012FF13C0A2
13E0A3127FEA1E601200A413E013C0A312011380120313005A120E5A1218123812300B1C
798919>I<17075F84171FA2173F177FA217FFA25E5EA24C6C7EA2EE0E3F161E161C1638
A21670A216E0ED01C084ED0380171FED07005D150E5DA25D157815705D844A5A170F4A5A
4AC7FC92B6FC5CA2021CC7120F143C14384A81A24A140713015C495AA249C8FC5B130E13
1E4982137C13FED807FFED1FFEB500F00107B512FCA219F83E417DC044>65
D<49B712F818FF19C0D9000190C7EA3FF0F00FF84BEC03FCF000FE197F0203EE3F805DF1
1FC0A20207EE0FE05D1AF0A2020F16075DA21AF8141F5DA2190F143F5DA21AF0147F4B15
1FA302FF17E092C9123FA21AC049177F5C1A8019FF010318005C4E5A61010716034A5E4E
5A180F010F4C5A4A5E4E5A4EC7FC011F16FE4A4A5AEF07F8013FED0FE0EF3FC04A49B4C8
FC017FEC0FFCB812F017C004FCC9FC453E7DBD4B>68 D<49B612C05BA2D90001EB800093
C7FC5DA314035DA314075DA3140F5DA3141F5DA3143F5DA3147F5DA314FF92C8FCA35B5C
A313035CA313075CA3130F5CA3131F5CA2133FA25CEBFFE0B612E0A32A3E7DBD28>73
D<49B612F0A3D900010180C7FC93C8FC5DA314035DA314075DA3140F5DA3141F5DA3143F
5DA3147F5DA314FF92C9FCA35B5C180C181E0103161C5C183C183813074A1578187018F0
130F4AEC01E0A21703011FED07C04A140F171F013FED3F8017FF4A1303017F021F1300B9
FCA25F373E7DBD3E>76 D<49B56C93383FFFF05113E098B5FCD90001F1E000704B5B03DF
933803BF80A2F2077F1403039F040E90C7FC1A1CDB8FE05E02075F030F4C5AA21AE1020F
EE01C1020E606F6CEC03811A83021EEE0703021C040E5BA2F11C07023C16380238606F6C
1470F1E00F14780270DB01C05BA2953803801F02F0ED07004A6C6C5E180E4E133F130102
C04B5C601A7F01036D6C5B4A95C8FC4D5A4D485B130791C749C75A170E047F1401495D01
0E4B5CA24D1303131E011C4B5C5F013C023F1407017C5D01FE92C75BD803FF4D7EB500FC
013E011FB512F8163C041C5E5C3E7DBD58>I<49B77E18F818FFD90001D900017F943800
3FE04BEC0FF0727E727E14034B6E7EA30207825DA3020F4B5A5DA24E5A141F4B4A5A614E
5A023F4B5A4B4A5A06FEC7FCEF03FC027FEC0FF04BEBFF8092B500FCC8FC5F9139FF8001
FE92C7EA7F80EF1FC084496F7E4A1407A28413035CA2170F13075C60171F130F5CA3011F
033F5B4AEE038018E0013F17071A004A021F5B496C160EB600E090380FF01E05075B716C
5ACBEAFFE0F03F8041407DBD45>82 D<DB07FC131C92383FFF8092B5EAE038913A03F803
F078913A0FE000F8F8DA1F80133D4AC7EA1FF0027E140F5C494814074A15E01303494814
03A2010F16C05CA3011F1680A38094C7FC808014FE90380FFFC015FC6DEBFFC016F86D14
FE6D806D81023F800207801400030F7F1500163F707E160F1607A2160312075A5F120EA2
001E15075FA24C5A123E003F4B5AA26D4AC7FC007F157E6D5C6D495AD87DF0495AD8F8FC
EB0FE090393F803F8027F01FFFFEC8FCD8E00713F839C0007FC036427BBF38>I<B6017F