From ae38@iamk4524.mathematik.uni-karlsruhe.de Wed Jan 5 17:27:28 1994 Received: from iamk4524.mathematik.uni-karlsruhe.de by dkuug.dk with SMTP id AA16704 (5.65c8/IDA-1.4.4j for ); Wed, 5 Jan 1994 16:26:01 +0100 Received: from localhost (ae38@localhost) by iamk4524.mathematik.uni-karlsruhe.de (8.6.4/8.6.4) id QAA25136 for SC22WG5@dkuug.dk; Wed, 5 Jan 1994 16:27:29 +0100 From: Wolfgang Walter Message-Id: <199401051527.QAA25136@iamk4524.mathematik.uni-karlsruhe.de> Subject: (Postscript version) Some Comments on a Note by Mary Payne To: SC22WG5@dkuug.dk (Fortran 90) Date: Wed, 5 Jan 1994 16:27:28 +0100 (MET) X-Mailer: ELM [version 2.4 PL21] Content-Type: text Content-Length: 64826 X-Charset: ASCII X-Char-Esc: 29 %!PS-Adobe-2.0 %%Creator: dvips, version 5.37 (C) 1986-90 Radical Eye Software %%Title: comments.dvi %%Pages: 3 1 %%BoundingBox: 0 0 612 842 %%EndComments %%BeginProcSet: tex.pro /TeXDict 200 dict def TeXDict begin /N /def load def /B{bind def}N /S /exch load def /X{S N}B /isls false N /vsize 10 N /@rigin{isls{[0 1 -1 0 0 0]concat} if 72 Resolution div 72 VResolution div neg scale Resolution VResolution vsize neg mul translate}B /@letter{/vsize 10 N}B /@landscape{/isls true N /vsize -1 N}B /@a4{/vsize 10.6929133858 N}B /@legal{/vsize 13 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 /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 /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 /E{pop nn dup definefont setfont}B /ch-image{ch-data dup type /stringtype ne{ctr get /ctr ctr 1 add N}if}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 /ctr 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 add]{ch-image}imagemask restore}B /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}B /eop{clear SI restore showpage userdict /eop-hook known{eop-hook}if}B /@start{userdict /start-hook known{ start-hook}if /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}B /p /show load N /RMat[1 0 0 -1 -.1 -.1]N /BDot 8 string N /v{/ruley X /rulex X V}B /V{gsave translate rulex ruley false RMat{BDot}imagemask 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 /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 /w{0 rmoveto}B /x{0 S rmoveto}B /y{3 2 roll p a}B /bos{/SS save N}B /eos{clear SS restore}B end %%EndProcSet TeXDict begin 1000 300 300 @start /Fa 16 122 df<000100020004000800100020006000 C00180018003000300060006000C000C001C001800180038003000300070007000700060006000 60006000E000E000E000E000E000E000600060006000600060003000300030001000180008000C 0004000200010010327BA413>40 D<0080004000200030001000180008000C000C000C00060006 0006000600060007000700070007000700070006000600060006000E000E000E000C000C001C00 1800180038003000300060006000C000C001800180030006000400080010002000400080001032 7FA413>I<03FC000606000F03000F03800601800001C0000380000380007F8003E3800F03801C 0380380700780700F00708F00708F00F08F00F08F017107867A01F83C015157D9418>97 D<00FF000381C00603C00C03C01C0180380000780000700000F00000F00000F00000F00000F000 00E00000F00000F000807001007001003806001C180007E00012157C9416>99 D<00FE000383800701C00C00E01C00E03800E07800E07000E0FFFFE0F00000F00000F00000F000 00E00000E00000F000407000803000801803000E0C0003F00013157D9416>101 D<00780003F80000700000700000700000700000700000700000E00000E00000E00000E00000E0 0000E00001C3F001CC1801D00C01E00E01E00E01C00E03C01C03801C03801C03801C03801C0380 1C0700380700380700380700380700380700380E00700F0078FFE7FF18237FA21B>104 D<007000F001F000F000E00000000000000000000000000000000001C00FC001C001C001C001C0 0380038003800380038003800700070007000700070007000E000F00FFE00C227FA10E>I<0003 00000780000F80000F80000700000000000000000000000000000000000000000000000000000E 0000FE00001E00000E00000E00000E00001C00001C00001C00001C00001C00001C000038000038 0000380000380000380000380000700000700000700000700000700000700000E00000E00060E0 00F1C000F18000E300003E0000112C83A10F>I<01C1F807E01FC60C183001D80E603801E00780 1C01E007801C01C007001C03C00F003803800E003803800E003803800E003803800E003803800E 003807001C007007001C007007001C007007001C007007001C007007001C00700E003800E00F00 3C00F0FFE3FF8FFE27157F942A>109 D<01C3F01FCC1801D00C01E00E01E00E01C00E03C01C03 801C03801C03801C03801C03801C0700380700380700380700380700380700380E00700F0078FF E7FF18157F941B>I<007E000383800600C00C00E01C0070380070780078700078F00078F00078 F00078F00078E000F0E000F0E000E0F001E07001C07003803807001C1C0007F00015157D9418> I<00E1F8000FE60C0000E8060000F0070000E0038000E0038001C003C001C003C001C003C001C0 03C001C003C001C003C003800780038007800380070003800F0003801E0003801C000740380007 30E000070F80000700000007000000070000000E0000000E0000000E0000000E0000000E000000 1E000000FFC000001A1F80941B>I<01C7C01FC8E001D1E001E1E001E0C001C00003C000038000 0380000380000380000380000700000700000700000700000700000700000E00000F0000FFF000 13157F9413>114 D<01F906070C0318031801180118021C001FE00FF807FC007E000E40064006 40066006600CE008D83087C010157E9413>I<008000800080018001000300030007000F001F00 FFF80E000E000E000E000E001C001C001C001C001C001C0038103810381038103810382038201C 4007800D1F7C9E13>I<0FFC0FE001E0078000E0030000E0020000E0040000E004000070080000 700800007010000070300000782000003840000038400000388000003C8000001D0000001D0000 001E0000001C0000000C0000000800000008000000100000001000000020000000400000F04000 00F0800000F1000000C20000003C0000001B1F80941A>121 D E /Fb 1 15 df<0F801FC0306060306030C018C018C018C0186030603030601FC00F800D0E7E8E12>14 D E /Fc 12 104 df0 D<70F8F8F87005057C8E0E>I<00C000 00C00000C00000C00000C000C0C0C0F0C3C038C7000EDC0003F00000C00003F0000EDC0038C700 F0C3C0C0C0C000C00000C00000C00000C00000C00012157D9619>3 D<03F0000FFC001C0E0030 0300600180600180C000C0C000C0C000C0C000C0C000C0C000C0C000C0C000C060018060018030 03001C0E000FFC0003F00012147D9519>14 D<0000000C0000003C000000F0000003C000000F00 00003C000000F0000007C000001F00000078000001E00000078000001E00000078000000E00000 00780000001E0000000780000001E0000000780000001F00000007C0000000F00000003C000000 0F00000003C0000000F00000003C0000000C000000000000000000000000000000000000000000 00000000000000000000007FFFFFF8FFFFFFFC1E277C9F27>20 D<000000006000000000003000 000000003000000000001800000000001800000000000C00000000000600000000000380FFFFFF FFFFE0FFFFFFFFFFC0000000000380000000000600000000000C00000000001800000000001800 0000000030000000000030000000000060002B127D9432>33 D<00000060000000000060000000 00003000000000001800000000001800000000000C0000000000060000000000030000FFFFFFFF 8000FFFFFFFFC000000000007000000000001C00000000000F800000000003E000000000078000 0000001E00000000003800000000006000FFFFFFFFC000FFFFFFFF800000000007000000000006 00000000000C00000000001800000000003000000000003000000000006000000000006000002B 1C7D9932>41 D<001FFF007FFF01E0000380000600000C00001800003000003000006000006000 00600000C00000C00000FFFFFFFFFFFFC00000C000006000006000006000003000003000001800 000C000006000003800001E000007FFF001FFF181E7C9A21>50 D<000030000000003000000000 78000000007800000000CC00000000CC0000000186000000018600000003030000000601800000 06018000000C00C000000C00C00000180060000018006000003000300000300030000060001800 006000180000C0000C0000C0000C00018000060001800006000300000300030000030006000001 800C000000C00C000000C018000000601800000060300000003030000000306000000018600000 0018FFFFFFFFFC7FFFFFFFFC26247DA32D>52 D<7FFFFFFFFCFFFFFFFFFC600000001860000000 1830000000303000000030180000006018000000600C000000C00C000000C00600000180030000 030003000003000180000600018000060000C0000C0000C0000C00006000180000600018000030 003000003000300000180060000018006000000C00C000000C00C0000006018000000601800000 03030000000186000000018600000000CC00000000CC0000000078000000007800000000300000 000030000026247D982D>I<000F0038007000E001C001C001C001C001C001C001C001C001C001 C001C001C001C001C001C001C001C0038007001E00F0001E000700038001C001C001C001C001C0 01C001C001C001C001C001C001C001C001C001C001C001C000E000700038000F10317CA419> 102 DI E /Fd 6 99 df<70F8FCFC7404040404080810102040060F7C840E>59 D<000080000180000180 000300000300000300000600000600000600000C00000C00000C00001800001800001800003000 00300000300000600000600000600000C00000C00000C000018000018000018000018000030000 0300000300000600000600000600000C00000C00000C0000180000180000180000300000300000 300000600000600000600000C00000C00000C0000011317DA418>61 D<007FC00001FF0007C000 03E00007C00005E00007C00005E00009E0000BC00009E0000BC00009E00013C00009E00023C000 11E00027800011E00047800011E00047800011E00087800021E0010F000020F0010F000020F002 0F000020F0040F000040F0041E000040F0081E000040F0081E000040F0101E000080F0203C0000 8078203C00008078403C00008078803C0001007880780001007900780001007900780001007A00 780002007C00F00002007C00F00002003800F00006003800F0000F003001F000FFE0203FFF0030 227EA12F>77 D<007FC003FF0007C000780007C000600005E000200009E000400009E000400008 F000400008F000400010F800800010780080001078008000103C008000203C010000203E010000 201E010000201E010000400F020000400F020000400F020000400782000080078400008007C400 008003C400008003C400010001E800010001E800010001F800010000F800020000F00002000070 00020000700006000070000F00002000FFE000200028227EA127>I<00786001C4E00302E00601 C00E01C01C01C03C01C0380380780380780380780380F00700F00700F00700F00708F00E10700E 10701E1030262018C6200F01C015157E941A>97 D<03C0003F8000038000038000038000070000 0700000700000700000E00000E00000E00000E00001C00001C78001D8E001E07003C0700380380 380380380780700780700780700780700780E00F00E00F00E00F00E01E00E01C00601C00603800 30700030C0000F000011237DA215>I E /Fe 75 124 df<001F83E000706E3000C07C780180F8 780380F07807007000070070000700700007007000070070000700700007007000070070000700 7000FFFFFFC0070070000700700007007000070070000700700007007000070070000700700007 007000070070000700700007007000070070000700700007007000070070000700700007007000 070078007FE3FF801D2380A21C>11 D<001FC0000070200000C010000180380003807800070078 000700300007000000070000000700000007000000070000000700000007000000FFFFF8000700 780007003800070038000700380007003800070038000700380007003800070038000700380007 00380007003800070038000700380007003800070038000700380007003800070038007FE1FF80 192380A21B>I<001FD8000070380000C078000180780003807800070038000700380007003800 070038000700380007003800070038000700380007003800FFFFF8000700380007003800070038 000700380007003800070038000700380007003800070038000700380007003800070038000700 38000700380007003800070038000700380007003800070038007FF3FF80192380A21B>I<000F C07F00007031C08000E00B004001801E00E003803E01E007003C01E007001C00C007001C000007 001C000007001C000007001C000007001C000007001C000007001C0000FFFFFFFFE007001C01E0 07001C00E007001C00E007001C00E007001C00E007001C00E007001C00E007001C00E007001C00 E007001C00E007001C00E007001C00E007001C00E007001C00E007001C00E007001C00E007001C 00E007001C00E007001C00E07FF1FFCFFE272380A229>I<7038F87CFC7EFC7E743A0402040204 0204020804080410081008201040200F0F7EA218>34 D<70F8FCFC740404040408081010204006 0F7CA20E>39 D<00200040008001000300060004000C000C001800180030003000300070006000 60006000E000E000E000E000E000E000E000E000E000E000E000E000E000E00060006000600070 00300030003000180018000C000C0004000600030001000080004000200B327CA413>I<800040 002000100018000C000400060006000300030001800180018001C000C000C000C000E000E000E0 00E000E000E000E000E000E000E000E000E000E000E000C000C000C001C0018001800180030003 000600060004000C00180010002000400080000B327DA413>I<00018000000180000001800000 018000000180000001800000018000000180000001800000018000000180000001800000018000 000180000001800000018000FFFFFFFEFFFFFFFE00018000000180000001800000018000000180 000001800000018000000180000001800000018000000180000001800000018000000180000001 8000000180001F227D9C26>43 D<70F8FCFC7404040404080810102040060F7C840E>II<70F8F8F87005057C840E>I<0000800001800001800003000003000003000006 00000600000600000C00000C00000C000018000018000018000030000030000030000060000060 0000600000C00000C00000C0000180000180000180000180000300000300000300000600000600 000600000C00000C00000C00001800001800001800003000003000003000006000006000006000 00C00000C00000C0000011317DA418>I<01F000071C000C06001803003803803803807001C070 01C07001C07001C0F001E0F001E0F001E0F001E0F001E0F001E0F001E0F001E0F001E0F001E0F0 01E0F001E0F001E0F001E07001C07001C07001C07803C03803803803801C07000C0600071C0001 F00013227EA018>I<008003800F80F38003800380038003800380038003800380038003800380 038003800380038003800380038003800380038003800380038003800380038007C0FFFE0F217C A018>I<03F0000C1C001007002007804003C04003C08003E0F003E0F801E0F801E0F801E02003 E00003E00003C00003C0000780000700000E00001C0000180000300000600000C0000180000100 000200200400200800201800603000403FFFC07FFFC0FFFFC013217EA018>I<03F8000C1E0010 07002007804007C07807C07803C07807C03807C0000780000780000700000F00000E0000380003 F000001C00000F000007800007800003C00003C00003E02003E07003E0F803E0F803E0F003C040 03C0400780200780100F000C1C0003F00013227EA018>I<000200000600000E00000E00001E00 001E00002E00004E00004E00008E00008E00010E00020E00020E00040E00040E00080E00100E00 100E00200E00200E00400E00800E00FFFFF8000E00000E00000E00000E00000E00000E00000E00 001F0001FFF015217FA018>I<1000801E07001FFF001FFE001FF80013E0001000001000001000 0010000010000010000010F800130E001407001803801003800001C00001C00001E00001E00001 E00001E07001E0F001E0F001E0E001C08001C04003C04003802007001006000C1C0003F0001322 7EA018>I<4000006000007FFFE07FFFC07FFFC0400080C0010080010080020080020000040000 080000080000100000300000200000600000600000600000E00000C00000C00001C00001C00001 C00001C00003C00003C00003C00003C00003C00003C00003C00003C00001800013237DA118>55 D<01F800060E000803001001802001802000C06000C06000C06000C07000C07801803E01003F02 001FC4000FF80003F80003FC00067F00083F80100F803007C06001C06000E0C000E0C00060C000 60C00060C000606000406000C03000801803000E0E0003F00013227EA018>I<01F000060C000C 0600180700380380700380700380F001C0F001C0F001C0F001E0F001E0F001E0F001E0F001E070 01E07003E03803E01805E00C05E00619E003E1E00001C00001C00001C000038000038030030078 0700780600700C002018001030000FC00013227EA018>I<70F8F8F87000000000000000000000 0070F8F8F87005157C940E>I61 D<07E01838201C400E80 0FF00FF00FF00F000F000E001C00380030006000C000C000800080018001000100010001000100 010000000000000000000000038007C007C007C0038010237DA217>63 D<000180000001800000 0180000003C0000003C0000003C0000005E0000005E000000DF0000008F0000008F0000010F800 001078000010780000203C0000203C0000203C0000401E0000401E0000401E0000800F0000800F 0000FFFF000100078001000780030007C0020003C0020003C0040003E0040001E0040001E00C00 00F00C0000F03E0001F8FF800FFF20237EA225>65 DI<0007E0100038183000E0063001C00170038000F0070000F00E0000701E00 00701C0000303C0000303C0000307C0000107800001078000010F8000000F8000000F8000000F8 000000F8000000F8000000F8000000F800000078000000780000107C0000103C0000103C000010 1C0000201E0000200E000040070000400380008001C0010000E0020000381C000007E0001C247D A223>IIII<0007F008003C0C1800E0021801C001B8038000F8070000780F0000381E0000381E0000 183C0000183C0000187C0000087800000878000008F8000000F8000000F8000000F8000000F800 0000F8000000F8000000F8001FFF780000F8780000787C0000783C0000783C0000781E0000781E 0000780F00007807000078038000B801C000B800E00318003C0C080007F00020247DA226>III75 DIII< 000FE00000783C0000E00E0003C00780078003C00F0001E00E0000E01E0000F03C0000783C0000 787C00007C7C00007C7800003C7800003CF800003EF800003EF800003EF800003EF800003EF800 003EF800003EF800003EF800003E7800003C7C00007C7C00007C3C0000783E0000F81E0000F00F 0001E00F0001E0078003C003C0078000E00E0000783C00000FE0001F247DA226>II82 D<03F0200C0C601802603001E07000E0600060E00060E00060E00020E00020 E00020F00000F000007800007F00003FF0001FFE000FFF0003FF80003FC00007E00001E00000F0 0000F0000070800070800070800070800070C00060C00060E000C0F000C0C80180C6070081FC00 14247DA21B>I<7FFFFFF87807807860078018400780084007800840078008C007800C80078004 800780048007800480078004000780000007800000078000000780000007800000078000000780 000007800000078000000780000007800000078000000780000007800000078000000780000007 800000078000000780000007800000078000000FC00003FFFF001E227EA123>IIII<7FF807FF0007 E001F80003C000E00003E000C00001E000800000F001000000F80300000078020000007C040000 003E0C0000001E080000001F100000000FB000000007A000000007C000000003E000000001E000 000001F000000003F80000000278000000047C0000000C3E000000081E000000101F000000300F 80000020078000004007C00000C003E000008001E000010001F000030000F000070000F8001F80 01FC00FFE007FFC022227FA125>I<0804100820102010402040208040804080408040B85CFC7E FC7E7C3E381C0F0F7AA218>92 D<0FE0001838003C0C003C0E0018070000070000070000070000 FF0007C7001E07003C0700780700700700F00708F00708F00708F00F087817083C23900FC1E015 157E9418>97 D<0E0000FE00001E00000E00000E00000E00000E00000E00000E00000E00000E00 000E00000E00000E00000E1F000E61C00E80600F00300E00380E003C0E001C0E001E0E001E0E00 1E0E001E0E001E0E001E0E001E0E001C0E003C0E00380F00700C80600C41C0083F0017237FA21B >I<01FE000703000C07801C0780380300780000700000F00000F00000F00000F00000F00000F0 0000F000007000007800403800401C00800C010007060001F80012157E9416>I<0000E0000FE0 0001E00000E00000E00000E00000E00000E00000E00000E00000E00000E00000E00000E001F8E0 0704E00C02E01C01E03800E07800E07000E0F000E0F000E0F000E0F000E0F000E0F000E0F000E0 7000E07800E03800E01801E00C02E0070CF001F0FE17237EA21B>I<01FC000707000C03801C01 C03801C07801E07000E0F000E0FFFFE0F00000F00000F00000F00000F000007000007800203800 201C00400E008007030000FC0013157F9416>I<003C00C6018F038F030F070007000700070007 000700070007000700FFF807000700070007000700070007000700070007000700070007000700 070007000700070007807FF8102380A20F>I<00007001F198071E180E0E181C07001C07003C07 803C07803C07803C07801C07001C07000E0E000F1C0019F0001000001000001800001800001FFE 000FFFC00FFFE03800F0600030400018C00018C00018C000186000306000303800E00E038003FE 0015217F9518>I<0E0000FE00001E00000E00000E00000E00000E00000E00000E00000E00000E 00000E00000E00000E00000E1F800E60C00E80E00F00700F00700E00700E00700E00700E00700E 00700E00700E00700E00700E00700E00700E00700E00700E00700E00700E0070FFE7FF18237FA2 1B>I<1C001E003E001E001C00000000000000000000000000000000000E00FE001E000E000E00 0E000E000E000E000E000E000E000E000E000E000E000E000E000E000E00FFC00A227FA10E>I< 01C003E003E003E001C00000000000000000000000000000000001E00FE001E000E000E000E000 E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E0 60E0F0C0F18061803E000B2C82A10F>I<0E0000FE00001E00000E00000E00000E00000E00000E 00000E00000E00000E00000E00000E00000E00000E03FC0E01F00E01C00E01800E02000E04000E 08000E10000E38000EF8000F1C000E1E000E0E000E07000E07800E03C00E01C00E01E00E00F00E 00F8FFE3FE17237FA21A>I<0E00FE001E000E000E000E000E000E000E000E000E000E000E000E 000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E000E00 0E00FFE00B237FA20E>I<0E1FC07F00FE60E183801E807201C00F003C00E00F003C00E00E0038 00E00E003800E00E003800E00E003800E00E003800E00E003800E00E003800E00E003800E00E00 3800E00E003800E00E003800E00E003800E00E003800E00E003800E00E003800E0FFE3FF8FFE27 157F942A>I<0E1F80FE60C01E80E00F00700F00700E00700E00700E00700E00700E00700E0070 0E00700E00700E00700E00700E00700E00700E00700E00700E0070FFE7FF18157F941B>I<01FC 000707000C01801800C03800E0700070700070F00078F00078F00078F00078F00078F00078F000 787000707800F03800E01C01C00E038007070001FC0015157F9418>I<0E1F00FE61C00E80600F 00700E00380E003C0E001C0E001E0E001E0E001E0E001E0E001E0E001E0E001E0E003C0E003C0E 00380F00700E80E00E41C00E3F000E00000E00000E00000E00000E00000E00000E00000E00000E 0000FFE000171F7F941B>I<01F8200704600E02601C01603801E07800E07800E0F000E0F000E0 F000E0F000E0F000E0F000E0F000E07000E07800E03801E01C01E00C02E0070CE001F0E00000E0 0000E00000E00000E00000E00000E00000E00000E00000E0000FFE171F7E941A>I<0E3CFE461E 8F0F0F0F060F000E000E000E000E000E000E000E000E000E000E000E000E000E000F00FFF01015 7F9413>I<0F8830786018C018C008C008E008F0007F803FE00FF001F8003C801C800C800CC00C C008E018D0308FC00E157E9413>I<02000200020002000600060006000E001E003E00FFF80E00 0E000E000E000E000E000E000E000E000E000E000E040E040E040E040E040E040708030801F00E 1F7F9E13>I<0E0070FE07F01E00F00E00700E00700E00700E00700E00700E00700E00700E0070 0E00700E00700E00700E00700E00700E00F00E00F006017003827800FC7F18157F941B>IIIII<3FFFC0380380300780200700600E00401C00403C0040 380000700000E00001E00001C0000380400700400F00400E00C01C0080380080780180700780FF FF8012157F9416>II E /Ff 47 124 df<00180030006000C001C0038007 0007000E001E001C003C003C003C007800780078007800F800F000F000F000F000F000F000F000 F000F000F000F000F80078007800780078003C003C003C001C001E000E0007000700038001C000 C00060003000180D317BA416>40 DI<387CFEFFFF7F3B03030706060C1C18702008117C8610>44 D<387CFEFEFE7C3807077C86 10>46 D<0000180000380000380000700000700000E00000E00000E00001C00001C00003800003 80000380000700000700000700000E00000E00001C00001C00001C000038000038000070000070 0000700000E00000E00001C00001C00001C0000380000380000700000700000700000E00000E00 000E00001C00001C0000380000380000380000700000700000E00000E00000C0000015317DA41C >I<00180000780001F800FFF800FFF80001F80001F80001F80001F80001F80001F80001F80001 F80001F80001F80001F80001F80001F80001F80001F80001F80001F80001F80001F80001F80001 F80001F80001F80001F80001F8007FFFE07FFFE013207C9F1C>49 D<03FC000FFF003C1FC07007 E07C07F0FE03F0FE03F8FE03F8FE01F87C01F83803F80003F80003F00003F00007E00007C0000F 80001F00003E0000380000700000E01801C0180380180700180E00380FFFF01FFFF03FFFF07FFF F0FFFFF0FFFFF015207D9F1C>I<00FE0007FFC00F07E01E03F03F03F03F81F83F81F83F81F81F 03F81F03F00003F00003E00007C0001F8001FE0001FF000007C00001F00001F80000FC0000FC3C 00FE7E00FEFF00FEFF00FEFF00FEFF00FC7E01FC7801F81E07F00FFFC001FE0017207E9F1C>I< 1000201E01E01FFFC01FFF801FFF001FFE001FF8001BC00018000018000018000018000019FC00 1FFF001E0FC01807E01803E00003F00003F00003F80003F83803F87C03F8FE03F8FE03F8FC03F0 FC03F07007E03007C01C1F800FFF0003F80015207D9F1C>53 D<01FE0007FF800F83E01E01F03E 00F07C00F87C0078FC007CFC007CFC007CFC007EFC007EFC007EFC007E7C00FE7C00FE3E01FE1E 037E0FFE7E07FC7E00207E00007C00007C1E007C3F00F83F00F83F00F03F01E01E03C01C0F800F FE0003F80017207E9F1C>57 D<000070000000007000000000F800000000F800000000F8000000 01FC00000001FC00000003FE00000003FE00000003FE00000006FF000000067F0000000E7F8000 000C3F8000000C3F800000183FC00000181FC00000381FE00000300FE00000300FE00000600FF0 00006007F00000E007F80000FFFFF80000FFFFF800018001FC00018001FC00038001FE00030000 FE00030000FE000600007F000600007F00FFE00FFFF8FFE00FFFF825227EA12A>65 D<0003FE0080001FFF818000FF01E38001F8003F8003E0001F8007C0000F800F800007801F8000 07803F000003803F000003807F000001807E000001807E00000180FE00000000FE00000000FE00 000000FE00000000FE00000000FE00000000FE00000000FE000000007E000000007E000001807F 000001803F000001803F000003801F800003000F8000030007C000060003F0000C0001F8003800 00FF00F000001FFFC0000003FE000021227DA128>67 DIII<0003FE0040001FFFC0C0007F00F1C001F8003FC003 F0000FC007C00007C00FC00003C01F800003C03F000001C03F000001C07F000000C07E000000C0 7E000000C0FE00000000FE00000000FE00000000FE00000000FE00000000FE00000000FE000000 00FE000FFFFC7E000FFFFC7F00001FC07F00001FC03F00001FC03F00001FC01F80001FC00FC000 1FC007E0001FC003F0001FC001FC003FC0007F80E7C0001FFFC3C00003FF00C026227DA12C>I< FFFFE0FFFFE003F80003F80003F80003F80003F80003F80003F80003F80003F80003F80003F800 03F80003F80003F80003F80003F80003F80003F80003F80003F80003F80003F80003F80003F800 03F80003F80003F80003F80003F80003F800FFFFE0FFFFE013227FA115>73 D<03FFFF03FFFF0007F00007F00007F00007F00007F00007F00007F00007F00007F00007F00007 F00007F00007F00007F00007F00007F00007F00007F00007F00007F00007F00007F01807F07E07 F0FF07F0FF07F0FF07F0FF07E07E0FE07C1FC01FFF0007F80018227FA11D>II77 DI<0007FC0000003FFF800000FC07E00003F0 01F80007E000FC000FC0007E001F80003F001F80003F003F00001F803F00001F807F00001FC07E 00000FC07E00000FC0FE00000FE0FE00000FE0FE00000FE0FE00000FE0FE00000FE0FE00000FE0 FE00000FE0FE00000FE0FE00000FE07E00000FC07F00001FC07F00001FC03F00001F803F80003F 801F80003F000FC0007E0007E000FC0003F001F80000FC07E000003FFF80000007FC000023227D A12A>II<01FC0407FF8C1F03FC 3C007C7C003C78001C78001CF8000CF8000CFC000CFC0000FF0000FFE0007FFF007FFFC03FFFF0 1FFFF80FFFFC03FFFE003FFE0003FF00007F00003F00003FC0001FC0001FC0001FE0001EE0001E F0003CFC003CFF00F8C7FFE080FF8018227DA11F>83 D<7FFFFFFF807FFFFFFF807E03F80F8078 03F807807003F803806003F80180E003F801C0E003F801C0C003F800C0C003F800C0C003F800C0 C003F800C00003F800000003F800000003F800000003F800000003F800000003F800000003F800 000003F800000003F800000003F800000003F800000003F800000003F800000003F800000003F8 00000003F800000003F800000003F800000003F800000003F8000003FFFFF80003FFFFF8002222 7EA127>IIII<07FC001FFF803F07C03F03E03F01E03F01F01E01 F00001F00001F0003FF003FDF01FC1F03F01F07E01F0FC01F0FC01F0FC01F0FC01F07E02F07E0C F81FF87F07E03F18167E951B>97 DI<00FF8007FFE00F83F01F03F03E03F07E03F07C01E07C0000FC0000FC0000FC0000 FC0000FC0000FC00007C00007E00007E00003E00301F00600FC0E007FF8000FE0014167E9519> I<0001FE000001FE0000003E0000003E0000003E0000003E0000003E0000003E0000003E000000 3E0000003E0000003E0000003E0001FC3E0007FFBE000F81FE001F007E003E003E007E003E007C 003E00FC003E00FC003E00FC003E00FC003E00FC003E00FC003E00FC003E00FC003E007C003E00 7C003E003E007E001E00FE000F83BE0007FF3FC001FC3FC01A237EA21F>I<00FE0007FF800F87 C01E01E03E01F07C00F07C00F8FC00F8FC00F8FFFFF8FFFFF8FC0000FC0000FC00007C00007C00 007E00003E00181F00300FC07003FFC000FF0015167E951A>I<003F8000FFC001E3E003C7E007 C7E00F87E00F83C00F80000F80000F80000F80000F80000F8000FFFC00FFFC000F80000F80000F 80000F80000F80000F80000F80000F80000F80000F80000F80000F80000F80000F80000F80000F 80000F80000F80007FF8007FF80013237FA211>I<03FC1E0FFF7F1F0F8F3E07CF3C03C07C03E0 7C03E07C03E07C03E07C03E03C03C03E07C01F0F801FFF0013FC003000003000003800003FFF80 1FFFF00FFFF81FFFFC3800FC70003EF0001EF0001EF0001EF0001E78003C7C007C3F01F80FFFE0 01FF0018217E951C>I I<1C003F007F007F007F003F001C000000000000000000000000000000FF00FF001F001F001F00 1F001F001F001F001F001F001F001F001F001F001F001F001F001F001F00FFE0FFE00B247EA310 >I 108 DII<00FE0007FFC00F83E01E00F03E00F87C007C7C007C7C007CFC007EFC007EFC007EFC007E FC007EFC007EFC007E7C007C7C007C3E00F81F01F00F83E007FFC000FE0017167E951C>I114 D<0FF3003FFF00781F00600700E00300E00300F00300FC00007FE0007FF8003FFE000FFF0001FF 00000F80C00780C00380E00380E00380F00700FC0E00EFFC00C7F00011167E9516>I<01800001 80000180000180000380000380000780000780000F80003F8000FFFF00FFFF000F80000F80000F 80000F80000F80000F80000F80000F80000F80000F80000F80000F81800F81800F81800F81800F 81800F830007C30003FE0000F80011207F9F16>II121 D123 D E /Fg 31 124 df<1C007F007F00FF80FFC0FFC07FC07FC01CC000C000C00180018001800300 030006000C001800300020000A157B8813>44 D<000E00001E00007E0007FE00FFFE00FFFE00F8 FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000 FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000FE0000 FE0000FE0000FE0000FE007FFFFE7FFFFE7FFFFE17277BA622>49 D<00000E0000001E0000003E 0000007E000000FE000000FE000001FE000003FE0000077E00000E7E00000E7E00001C7E000038 7E0000707E0000E07E0000E07E0001C07E0003807E0007007E000E007E000E007E001C007E0038 007E0070007E00E0007E00FFFFFFF8FFFFFFF8FFFFFFF80000FE000000FE000000FE000000FE00 0000FE000000FE000000FE000000FE00007FFFF8007FFFF8007FFFF81D277EA622>52 D<0C0003000F803F000FFFFE000FFFFC000FFFF8000FFFF0000FFFE0000FFFC0000FFE00000E00 00000E0000000E0000000E0000000E0000000E0000000E7FC0000FFFF8000F80FC000E003E000C 003F0000001F8000001FC000001FC000001FE000001FE018001FE07C001FE0FE001FE0FE001FE0 FE001FE0FE001FC0FC001FC078003F8078003F803C007F001F01FE000FFFF80003FFF00000FF80 001B277DA622>I<380000003E0000003FFFFFF03FFFFFF03FFFFFF07FFFFFE07FFFFFC07FFFFF 807FFFFF0070000E0070000E0070001C00E0003800E0007000E000E0000000E0000001C0000003 80000007800000078000000F0000000F0000001F0000001F0000003F0000003E0000003E000000 7E0000007E0000007E0000007E000000FE000000FE000000FE000000FE000000FE000000FE0000 00FE000000FE0000007C0000003800001C297CA822>55 D<000003800000000007C00000000007 C0000000000FE0000000000FE0000000000FE0000000001FF0000000001FF0000000003FF80000 00003FF8000000003FF80000000073FC0000000073FC00000000F3FE00000000E1FE00000000E1 FE00000001C0FF00000001C0FF00000003C0FF80000003807F80000007807FC0000007003FC000 0007003FC000000E003FE000000E001FE000001E001FF000001C000FF000001FFFFFF000003FFF FFF800003FFFFFF80000780007FC0000700003FC0000700003FC0000E00001FE0000E00001FE00 01E00001FF0001C00000FF0001C00000FF00FFFE001FFFFEFFFE001FFFFEFFFE001FFFFE2F297E A834>65 D<00003FF001800003FFFE0380000FFFFF8780003FF007DF8000FF8001FF8001FE0000 7F8003FC00003F8007F000001F800FF000000F801FE0000007801FE0000007803FC0000007803F C0000003807FC0000003807F80000003807F8000000000FF8000000000FF8000000000FF800000 0000FF8000000000FF8000000000FF8000000000FF8000000000FF8000000000FF80000000007F 80000000007F80000000007FC0000003803FC0000003803FC0000003801FE0000003801FE00000 07000FF00000070007F000000E0003FC00001E0001FE00003C0000FF8000F800003FF007E00000 0FFFFFC0000003FFFF000000003FF8000029297CA832>67 D69 DI73 D75 DIII<007F806003FFF0E007FFF9E00F807FE01F001FE03E0007E07C0003E07C0001E0FC0001E0FC 0001E0FC0000E0FE0000E0FE0000E0FF000000FFC000007FFE00007FFFE0003FFFFC001FFFFE00 0FFFFF8007FFFFC003FFFFE000FFFFE00007FFF000007FF000000FF8000007F8000003F8600001 F8E00001F8E00001F8E00001F8F00001F0F00001F0F80003F0FC0003E0FF0007C0FFE01F80F3FF FF00E0FFFE00C01FF0001D297CA826>83 D<01FF800007FFF0000F81F8001FC07E001FC07E001F C03F000F803F8007003F8000003F8000003F8000003F80000FFF8000FFFF8007FC3F800FE03F80 3F803F803F003F807F003F80FE003F80FE003F80FE003F80FE003F807E007F807F00DF803F839F FC0FFF0FFC01FC03FC1E1B7E9A21>97 D<001FF80000FFFE0003F01F0007E03F800FC03F801F80 3F803F801F007F800E007F0000007F000000FF000000FF000000FF000000FF000000FF000000FF 000000FF0000007F0000007F0000007F8000003F8001C01F8001C00FC0038007E0070003F01E00 00FFFC00001FE0001A1B7E9A1F>99 D<00003FF80000003FF80000003FF800000003F800000003 F800000003F800000003F800000003F800000003F800000003F800000003F800000003F8000000 03F800000003F800000003F800001FE3F80000FFFBF80003F03FF80007E00FF8000FC007F8001F 8003F8003F8003F8007F0003F8007F0003F8007F0003F800FF0003F800FF0003F800FF0003F800 FF0003F800FF0003F800FF0003F800FF0003F8007F0003F8007F0003F8007F0003F8003F8003F8 001F8003F8000F8007F80007C00FF80003F03BFF8000FFF3FF80003FC3FF80212A7EA926>I<00 3FE00001FFF80003F07E0007C01F000F801F801F800F803F800FC07F000FC07F0007C07F0007E0 FF0007E0FF0007E0FFFFFFE0FFFFFFE0FF000000FF000000FF0000007F0000007F0000007F0000 003F8000E01F8000E00FC001C007E0038003F81F0000FFFE00001FF0001B1B7E9A20>I104 D<07000FC01FE03FE03FE03FE01FE00FC00700 0000000000000000000000000000FFE0FFE0FFE00FE00FE00FE00FE00FE00FE00FE00FE00FE00F E00FE00FE00FE00FE00FE00FE00FE00FE00FE00FE00FE0FFFEFFFEFFFE0F2B7DAA14>I107 DIII<003FE00001FFFC0003F07E000FC01F801F800FC03F 800FE03F0007E07F0007F07F0007F07F0007F0FF0007F8FF0007F8FF0007F8FF0007F8FF0007F8 FF0007F8FF0007F8FF0007F87F0007F07F0007F03F800FE03F800FE01F800FC00FC01F8007F07F 0001FFFC00003FE0001D1B7E9A22>I114 D<03FE300FFFF01E03F03800 F0700070F00070F00070F80070FC0000FFE0007FFE007FFF803FFFE01FFFF007FFF800FFF80003 FC0000FC60007CE0003CF0003CF00038F80038FC0070FF01E0F7FFC0C1FF00161B7E9A1B>I<00 700000700000700000700000F00000F00000F00001F00003F00003F00007F0001FFFF0FFFFF0FF FFF007F00007F00007F00007F00007F00007F00007F00007F00007F00007F00007F00007F00007 F00007F03807F03807F03807F03807F03807F03803F03803F87001F86000FFC0001F8015267FA5 1B>II123 D E end %%EndProlog %%BeginSetup %%Feature: *Resolution 300 TeXDict begin @a4 %%EndSetup %%Page: 1 1 bop 514 75 a Fg(Some)23 b(Commen)n(ts)e(on)i(a)h(Note)e(on)59 178 y(F)-6 b(ortran,)25 b(LIA{1,)f(IEEE)g(754)g(and)f(Kulisc)n(h{Mirank)n(er) f(Arithmetic)294 303 y Ff(\(F)-5 b(ortran)19 b(do)r(cumen)n(t)e(ISO/IEC)i (JTC1/SC22/W)n(G5)11 b({)e(N953,)559 385 y(b)n(y)19 b(Mary)g(P)n(a)n(yne,)f (August)i(31,)e(1993\))451 534 y(b)n(y)h(Ulric)n(h)f(Kulisc)n(h)h(and)g(W)-5 b(olfgang)19 b(V.)f(W)-5 b(alter)742 616 y(Decem)n(b)r(er)16 b(31,)i(1993)0 768 y Fe(The)23 b(note)h(b)o(y)f(Mary)g(P)o(a)o(yne)g (addresses)h(an)g(imp)q(ortan)o(t)e(sub)s(ject)h(whic)o(h)g(is)g(cen)o(tral)f (to)i(the)f(whole)0 828 y(of)c(n)o(umerical)e(computation:)26 b(the)19 b(prop)q(erties)g(of)h(computer)e(arithmetic)e(and)k(their)f(in)o (terpla)o(y)e(with)0 888 y(programming)12 b(languages.)21 b(The)14 b(note)g(is)f(highly)g(relev)m(an)o(t)f(to)i(con)o(temp)q(orary)e(dev)o (elopmen)o(ts)f(of)j(arith-)0 948 y(metic)e(standards)j(and)g(hardw)o(are.)21 b(W)l(e)14 b(w)o(ould)g(therefore)f(lik)o(e)f(to)j(commen)o(t)c(on)j(it,)g (in)f(particular)h(with)0 1009 y(resp)q(ect)i(to)h(what)f(is)h(there)e (called)g(Kulisc)o(h{Mirank)o(er)g(\(K-M\))h(arithmetic.)0 1094 y(K-M)i(arithmetic)d(has)k(b)q(een)f(around)h(for)f(ab)q(out)i(20)e(y)o (ears.)26 b(It)18 b(has)h(often)f(b)q(een)g(misin)o(terprete)o(d)d(and)0 1154 y(misundersto)q(o)q(d.)28 b(Most)19 b(misunderstandings)f(stem)f(from)g (the)i(fact)f(that)h(K-M)f(arithmetic)e(is)j(based)0 1214 y(on)14 b(a)h(clear)e(mathematical)e(theory)j(whereas)g(the)g(usual)g(discussion)g (of)h(computer)d(arithmetic)g(in)h(more)0 1274 y(collo)q(quial)h(terms)e (generally)i(fails)g(to)h(meet)d(a)j(rigorous)g(mathematical)c(de\014nition.) 20 b(The)15 b(misin)o(terpre-)0 1335 y(tations)i(are)f(mostly)f(of)h(a)h(p)q (ersonal)g(nature.)22 b(W)l(e)16 b(will)f(not)h(commen)o(t)d(on)k(them.)0 1420 y(K-M)e(arithmetic)e(w)o(as)j(essen)o(tially)e(dev)o(elop)q(ed)h(to)h (enhance)f(the)g(mathematical)e(p)q(o)o(w)o(er)i(of)h(the)g(digital)0 1480 y(computer.)23 b(T)l(raditionally)l(,)16 b(the)h(computer)f(is)h(only)h (capable)f(of)g(computing)g(appro)o(ximate)e(solutions)0 1540 y(to)21 b(a)g(problem.)32 b(A)20 b(lot)g(of)g(mathematics)e(ma)o(y)h(b)q(e)h (required)g(to)g(dev)o(elop)g(an)g(algorithm)g(for)g(suc)o(h)h(a)0 1600 y(task)h(and)g(to)f(deriv)o(e)f(error)h(estimates)f(for)h(the)g (computed)f(results.)37 b(Ho)o(w)o(ev)o(er,)20 b(this)h(in)g(itself)f(do)q (es)0 1660 y(not)c(imply)c(that)k(a)f(strict)g(mathematical)d(answ)o(er)j (can)g(b)q(e)h(deduced)e(b)q(ecause)h(classical)g(error)g(analysis)0 1721 y(do)q(es)i(not)f(tak)o(e)g(in)o(to)g(accoun)o(t)g(all)f(computational)g (errors.)22 b(Th)o(us,)16 b(ev)o(en)f(though)i(appro)o(ximations)e(are)0 1781 y(often)i(satisfactory)l(,)f(they)g(ma)o(y)f(nev)o(ertheless)g(b)q(e)h (inaccurate)g(or)h(completely)d(wrong)j(in)f(a)h(signi\014can)o(t)0 1841 y(n)o(um)o(b)q(er)c(of)h(cases.)21 b(With)14 b(traditional)h (\015oating-p)q(oin)o(t)g(arithmetic,)d(it)i(is)g(usually)g(v)o(ery)f (di\016cult)g(to)i(tell)0 1901 y(when)h(something)g(has)h(gone)f(wrong)i (simply)c(b)o(y)i(lo)q(oking)g(at)h(the)f(n)o(umerical)d(results.)0 1986 y(In)h(con)o(trast,)g(K-M)g(arithmetic)d(aims)i(to)i(do)f(rigorous)h (mathematics)d(with)h(the)h(computer.)19 b(F)l(or)14 b(exam-)0 2047 y(ple,)i(rigorous,)i(highly)f(accurate)g(b)q(ounds)i(can)e(b)q(e)g (computed)f(for)h(the)g(solution)h(of)f(systems)f(of)i(equa-)0 2107 y(tions)e(or)g(initial,)e(b)q(oundary)l(,)i(or)g(eigen)o(v)m(alue)f (problems)f(of)i(ODEs)h(or)f(in)o(tegral)f(equations.)21 b(A)15 b(pro)q(of)i(of)0 2167 y(the)h(existence)e(and)j(uniqueness)f(of)g(the)g (solution)h(within)e(the)h(computed)f(b)q(ounds)j(is)e(automatically)0 2227 y(p)q(erformed)f(b)o(y)i(the)f(computer.)27 b(This)19 b(is)f(often)h(referred)e(to)i(as)h(\\automatic)e(result)g(v)o (eri\014cation")g(or)0 2287 y(\\self-v)m(alidation)13 b(of)h(computational)f (results".)20 b(Ev)o(en)13 b(when)h(used)f(to)h(compute)e(appro)o(ximations)h (only)l(,)0 2347 y(K-M)k(arithmetic)e(generally)i(deliv)o(ers)e(results)i(of) h(higher)f(accuracy)g(than)h(traditional)g(\015oating-p)q(oin)o(t)0 2408 y(arithmetic)g(due)h(to)i(its)f(strict)f(mathematical)e(requiremen)o (ts.)29 b(More)20 b(imp)q(ortan)o(tly)l(,)f(ho)o(w)o(ev)o(er,)g(K-M)0 2468 y(arithmetic)14 b(enables)i(the)g(transition)g(from)f(Numerical)f (Mathematics)g(to)j(Mathematical)d(Numerics.)0 2553 y(What)k(is)e(K-M)h (arithmetic?)22 b(Instead)17 b(of)g(reducing)g(all)f(calculations)h(to)h(the) e(four)i(elemen)o(tary)c(arith-)0 2613 y(metic)k(op)q(erations)k(for)g (\015oating-p)q(oin)o(t)g(n)o(um)o(b)q(ers,)e(K-M)h(arithmetic)d(requires)i (t)o(w)o(elv)o(e)e(fundamen)o(tal)0 2673 y(data)f(t)o(yp)q(es)g(\(or)g (mathematical)c(spaces\))k(in)f(a)h(computing)e(en)o(vironmen)o(t.)20 b(These)c(are)h(the)f(four)h(basic)0 2734 y(data)j(t)o(yp)q(es)e(real,)h (complex,)e(in)o(terv)m(al)h(and)h(complex)e(in)o(terv)m(al)h(and)h(the)g(v)o (ectors)f(and)i(matrices)d(o)o(v)o(er)p eop %%Page: 2 2 bop 0 75 a Fe(these)19 b(t)o(yp)q(es.)32 b(All)18 b(computer)g(op)q(erations) j(for)f(these)g(t)o(yp)q(es)f(should)h(b)q(e)g(de\014ned)g(b)o(y)f(the)g (concept)h(of)0 135 y(\\semimorphism)o(":)0 218 y(If)c Fd(M)22 b Fe(is)17 b(an)o(y)g(one)g(of)g(these)f(t)o(w)o(elv)o(e)f(spaces)i(and)g Fd(N)22 b Fe(its)17 b(computer-represen)o(table)d(subset,)j(and)g(if)g Fc(\016)f Fe(is)0 278 y(an)o(y)h(arithmetic)d(op)q(eration)k(in)e Fd(M)5 b Fe(,)17 b(then)g(the)g(corresp)q(onding)g(computer)f(op)q(eration)p 1641 252 29 2 v 1641 277 2 25 v 1647 272 a Fb(\016)p 1668 277 V 1641 279 29 2 v 1693 278 a Fe(in)g Fd(N)23 b Fe(should)0 338 y(b)q(e)16 b(de\014ned)g(b)o(y)16 442 y(\(R)o(G\))159 b Fd(a)p 319 416 V 319 440 2 25 v 325 435 a Fb(\016)p 345 440 V 319 442 29 2 v 353 442 a Fd(b)14 b Fe(:=)p 460 416 V 460 440 2 25 v 486 440 V 460 442 29 2 v 55 w(\()p Fd(a)c Fc(\016)h Fd(b)p Fe(\))317 b(for)16 b(all)g Fd(a;)8 b(b)13 b Fc(2)h Fd(N)323 b Fa(\(semimorphism)o(\))0 541 y Fe(where)p 152 516 V 152 540 2 25 v 178 540 V 152 542 29 2 v 82 w(:)20 b Fd(M)27 b Fc(!)20 b Fd(N)26 b Fe(is)20 b(a)g(mapping)g(from)f Fd(M)26 b Fe(on)o(to)21 b Fd(N)k Fe(\(called)19 b(a)i(rounding\))g(with)f(the)g(follo)o(wing)0 602 y(prop)q(erties:)16 705 y(\(R1\))p 293 679 V 293 704 2 25 v 319 704 V 293 706 29 2 v 214 w(\()p Fd(a)p Fe(\))13 b(=)h Fd(a)460 b Fe(for)16 b(all)g Fd(a)d Fc(2)h Fd(N)458 b Fa(\(pro)s(jection\))16 783 y Fe(\(R2\))172 b Fd(a)13 b Fc(\024)h Fd(b)30 b Fc(\))p 515 757 V 515 781 2 25 v 542 781 V 515 783 29 2 v 71 w Fe(\()p Fd(a)p Fe(\))14 b Fc(\024)p 687 757 V 687 781 2 25 v 713 781 V 687 783 29 2 v 55 w Fe(\()p Fd(b)p Fe(\))162 b(for)16 b Fd(a;)8 b(b)13 b Fc(2)h Fd(M)411 b Fa(\(monotonicit)o(y\))16 860 y Fe(\(R3\))p 293 834 V 293 859 2 25 v 319 859 V 293 861 29 2 v 214 w(\()p Fc(\000)p Fd(a)p Fe(\))13 b(=)g Fc(\000)p 541 834 V 541 859 2 25 v 567 859 V 541 861 29 2 v 42 w Fe(\()p Fd(a)p Fe(\))303 b(for)16 b(all)g Fd(a)d Fc(2)h Fd(M)371 b Fa(\(an)o(tisymme)o(try\))0 960 y Fe(This)18 b(means)e(that)i(ev)o(ery)e (computer)g(op)q(eration)i(should)g(b)q(e)g(p)q(erformed)e(in)h(suc)o(h)g(a)h (w)o(a)o(y)f(that)h(it)f(pro-)0 1020 y(duces)i(the)g(same)e(result)i(as)h(if) e(the)h(mathematicall)o(y)d(correct)i(op)q(eration)i(w)o(ere)e(\014rst)i(p)q (erformed)d(and)0 1080 y(the)g(exact)g(result)g(then)g(rounded)h(in)o(to)f (the)g(computer-represen)o(table)d(subset.)25 b(This)18 b(rule)e(is)h (optimal)0 1140 y(in)e(the)g(sense)g(that)h(there)e(is)h(no)h(b)q(etter)f (computer-represen)o(table)e(appro)o(ximation)h(to)i(the)f(true)g(result)0 1201 y(\(with)h(resp)q(ect)g(to)h(the)f(prescrib)q(ed)f(rounding\).)0 1283 y(A)o(t)20 b(the)g(end)h(of)g(the)f(in)o(tro)q(duction)h(of)g(the)f(b)q (o)q(ok)i(on)f(Computer)f(Arithmetic)e(b)o(y)i(Kulisc)o(h)f(and)i(Mi-)0 1343 y(rank)o(er,)c(it)h(is)g(already)f(stated)h(that)h(this)f(concept)f(of)h (semimorphism)c(should)k(and)h(\\can)f(b)q(e)g(used)g(as)0 1403 y(an)g(axiomatic)f(de\014nition)g(of)i(computer)d(arithmetic)g(in)h(the) h(con)o(text)f(of)h(programming)f(languages".)0 1463 y(A)d(careful)g (analysis)h(of)g(the)f(requiremen)o(ts)d(is)k(giv)o(en)e(in)i(the)f(b)q(o)q (ok.)22 b(All)13 b(semimorphic)e(computer)i(op)q(er-)0 1524 y(ations)i(are)g(of)f(1)h(ulp)g(accuracy)f(\(or)h(ev)o(en)e(1/2)j(ulp)e (accuracy)g(in)g(the)g(case)h(of)g(rounding)g(to)g(nearest\).)20 b(It)0 1584 y(is)14 b(ultimately)d(sho)o(wn)k(that)g(semimorphic)c(op)q (erations)k(can)g(b)q(e)f(pro)o(vided)g(on)h(a)f(computer)f(if,)h(on)h(a)f (lo)o(w)0 1644 y(lev)o(el)f(\(preferrably)h(in)h(hardw)o(are\),)g(15)h (fundamen)o(tal)e(op)q(erations)i(are)g(a)o(v)m(ailable:)k(the)15 b(\014v)o(e)f(op)q(erations)0 1704 y(+)p Fd(;)8 b Fc(\000)p Fd(;)16 b Fc(\003)8 b Fd(;)16 b(=)8 b(;)17 b Fc(\001)8 b Fe(,)21 b(eac)o(h)e(one)h(with)g(the)g(three)g(roundings)p 1085 1678 V 1085 1703 2 25 v 1112 1703 V 1085 1705 29 2 v 62 w Fd(;)8 b Fc(5)p Fd(;)g Fc(4)p Fe(.)31 b(Here)19 b Fc(\001)h Fe(means)f(the)h(dot)g (pro)q(duct)0 1764 y(of)f(t)o(w)o(o)g(v)o(ectors,)p 342 1739 V 342 1764 2 25 v 369 1764 V 342 1766 29 2 v 79 w(is)g(a)g(monotone,)g(an)o (tisymmetri)o(c)d(rounding)k(\(e.)e(g.)h(to)g(nearest\))g(and)h Fc(5)e Fe(and)i Fc(4)0 1825 y Fe(are)c(the)g(monotone)g(do)o(wn)o(w)o(ardly)g (and)h(up)o(w)o(ardly)e(directed)g(roundings.)23 b(All)14 b(15)j(op)q (erations)p 1788 1799 V 1788 1824 2 25 v 1794 1819 a Fb(\016)p 1815 1824 V 1788 1826 29 2 v 1823 1825 a Fd(;)12 b Fc(5)1864 1819 y Fb(\016)1887 1825 y Fd(;)g Fc(4)1927 1819 y Fb(\016)0 1885 y Fe(with)k Fc(\016)e(2)g(f)p Fe(+)p Fd(;)8 b Fc(\000)p Fd(;)15 b Fc(\003)8 b Fd(;)17 b(=)8 b(;)17 b Fc(\001)8 b(g)16 b Fe(are)g(de\014ned)g(b)o(y)g(\(R)o(G\).)0 1967 y(The)e(IEEE)g(754)h (standard)h(prescrib)q(es)d(12)i(of)f(these)g(15)h(op)q(erations)g(\(for)f Fc(\016)g(2)g(f)p Fe(+)p Fd(;)8 b Fc(\000)p Fd(;)15 b Fc(\003)8 b Fd(;)17 b(=)8 b Fc(g)p Fe(\),)14 b(but)g(fails)0 2027 y(to)j(de\014ne)f (the)h(dot)g(pro)q(duct)g(op)q(erations,)h(whic)o(h)e(are)h(the)f(most)g (crucial)g(op)q(erations)h(in)g(the)f(common)0 2087 y(v)o(ector)f(spaces.)22 b(IEEE)16 b(754)i(also)f(prescrib)q(es)e(sp)q(eci\014c)h(data)h(formats.)0 2170 y(The)j(LIA{1)h(standard)g(is)f(more)f(general)h(and)h(do)q(es)g(not)g (sp)q(ecify)e(an)o(y)h(particular)g(data)h(formats)f(or)0 2230 y(rounding)25 b(mo)q(des.)43 b(It)24 b(only)f(requires)g(that)i(the)e(four)i (elemen)o(tary)c(\015oating-p)q(oin)o(t)k(op)q(erations)g(b)q(e)0 2290 y(accurate)d(to)g(1)h(ulp)f(\(unit)g(in)g(the)f(last)i(place\).)38 b(As)22 b(a)h(consequence,)f(in)o(terv)m(al)f(arithmetic)e(cannot)0 2350 y(b)q(e)h(implem)o(en)o(te)o(d)e(in)h(a)i(straigh)o(tforw)o(ard)f (manner)f(\(in)h(con)o(trast)g(to)h(the)f(situation)g(with)g(IEEE)g(754)0 2411 y(arithmetic\).)f(Also,)d(addition)g(and)h(subtraction)g(are)g(not)g (required)e(to)h(satisfy)h(prop)q(ert)o(y)f(\(R)o(G\).)g(This)0 2471 y(unfortunate)d(allo)o(w)o(ance)f(w)o(as)h(apparen)o(tly)g(made)e(to)i (accomo)q(date)f(/370)j(arithmetic,)10 b(whic)o(h)i(is)g(of)h(1)g(ulp)0 2531 y(accuracy)l(,)i(but)i(whose)g(addition)f(and)h(subtraction)g(are)f(not) h(semim)o(orphic.)0 2613 y(Mary)12 b(P)o(a)o(yne's)f(note)i(states)f(that)h (\\lik)o(e)e(IEEE)h(754,)i(K-M)e(arithmetic)d(can)k(conform)e(to)h(LIA{1)q (".)20 b(Since)0 2673 y(K-M)g(arithmetic)e(sp)q(eci\014es)h(the)h (mathematical)d(requiremen)o(ts,)h(the)i(opp)q(osite)g(question)g(should)h(b) q(e)0 2734 y(ask)o(ed:)g(do)16 b(IEEE)f(754)i(and)f(LIA{1)g(conform)e(to)i (K-M)f(arithmetic.)k(The)c(answ)o(er)h(is)f(no)h(in)f(b)q(oth)i(cases)963 2858 y(2)p eop %%Page: 3 3 bop 0 75 a Fe(since)19 b(neither)f(standard)i(pro)o(vides)f(the)g(dot)h(pro)q (duct)g(op)q(erations)p 1303 49 29 2 v 1303 74 2 25 v 1311 73 a Fc(\001)p 1329 74 V 1303 76 29 2 v 1338 75 a Fd(;)12 b Fc(5)-29 b(\001)7 b Fd(;)12 b Fc(4)-29 b(\001)16 b Fe(.)30 b(F)l(urthermore,)18 b(LIA{1)0 135 y(do)q(es)i(not)g(require)f(semim)o (orphic)e(addition)i(and)i(subtraction)f(and)g(do)q(es)g(not)g(pro)o(vide)f (the)g(directed)0 195 y(roundings)e Fc(5)f Fe(and)h Fc(4)8 b Fe(.)0 276 y(Tw)o(o)13 b(serious)g(disadv)m(an)o(tages)g(of)g(most)f (existing)g(IEEE)g(pro)q(cessors)i(are)f(that)g(full)e(double-length)i(pro)q (d-)0 337 y(ucts)18 b(cannot)h(b)q(e)f(obtained)g(from)f(the)h(outside)g (\(ev)o(en)f(though)j(they)d(are)h(computed)f(in)o(ternally\))f(and)0 397 y(that)i(the)f(op)q(erations)h(with)f(directed)f(roundings)23 b Fc(5)975 391 y Fb(\016)1015 397 y Fe(and)f Fc(4)1130 391 y Fb(\016)1170 397 y Fe(where)17 b Fc(\016)f(2)f(f)p Fe(+)p Fd(;)8 b Fc(\000)p Fd(;)16 b Fc(\003)8 b Fd(;)16 b(=)8 b Fc(g)18 b Fe(are)g(not)f(pro-)0 457 y(vided)i(as)h(single)g(instructions)f(with)h (the)g(rounding)g(mo)q(de)f(incorp)q(orated.)33 b(The)20 b(former)e (restriction)0 517 y(complicates)f(the)i(implem)o(e)o(n)o(tation)d(of)k (accurate)f(complex)d(arithmetic)h(and)i(of)h(accurate)e(dot)i(pro)q(d-)0 577 y(ucts)i(unnecessarily)e(whereas)i(the)g(latter)f(slo)o(ws)h(do)o(wn)g (in)o(terv)m(al)e(arithmetic)g(considerably)l(.)36 b(Man)o(y)0 638 y(pro)q(cessors)15 b(need)e(as)i(man)o(y)d(mac)o(hine)f(cycles)i(to)h (set)f(the)h(rounding)g(mo)q(de)f(as)h(to)g(p)q(erform)f(the)h(rounded)0 698 y(arithmetic)21 b(op)q(eration)j(itself.)41 b(Since)23 b(in)g(in)o(terv)m(al)f(calculations,)i(the)g(rounding)g(mo)q(de)e(has)i(to)g (b)q(e)0 758 y(switc)o(hed)15 b(frequen)o(tly)l(,)f(this)i(in)o(tro)q(duces)g (a)h(sev)o(ere)e(and)i(unnecessary)f(p)q(enalt)o(y)l(.)0 839 y(A)40 b(\\Prop)q(osal)i(for)e(Accurate)f(Floating-P)o(oin)o(t)h(V)l(ector)f (Arithmetic")e(w)o(as)k(drafted)f(b)o(y)f(the)0 899 y(\\IMA)o(CS/GAMM)15 b(W)l(orking)h(Group)g(on)h(Enhanced)f(Computer)f(Arithmeti)o(c")e(and)k (extensiv)o(ely)c(dis-)0 959 y(cussed)j(during)g(the)g(past)h(t)o(w)o(o)f(y)o (ears.)k(In)c(its)g(\014nal)g(form,)f(the)g(prop)q(osal)j(w)o(as)e (unanimously)f(appro)o(v)o(ed)0 1019 y(and)20 b(adopted)h(b)o(y)e(GAMM)f(and) j(b)o(y)e(IMA)o(CS)f(in)h(1993)j(and)e(is)f(no)o(w)h(an)g(o\016cial)f(do)q (cumen)o(t)f(of)i(these)0 1080 y(mathematical)15 b(so)q(cieties.)27 b(In)18 b(essence,)f(it)h(requires)f(K-M)h(arithmetic)e(and)j(giv)o(es)e (some)h(hin)o(ts)g(on)g(its)0 1140 y(realization.)0 1221 y(K-M)g(arithmetic)e (has)j(b)q(een)f(implem)o(en)o(te)o(d)e(in)h(hardw)o(are,)i(in)f(\014rm)o(w)o (are,)f(and)i(in)e(soft)o(w)o(are)i(on)f(man)o(y)0 1281 y(di\013eren)o(t)23 b(platforms)h(and)h(has)g(b)q(een)f(a)o(v)m(ailable)g(since)f(the)h(early)g (1980's.)47 b(The)24 b(accum)o(ulation)f(of)0 1341 y(pro)q(ducts)13 b(in)e(a)i(long)f(\014xed-p)q(oin)o(t)g(register)f(allo)o(ws)h(an)h(easy)f (realization)f(of)h(the)g(optimal)e(dot)j(pro)q(duct)f(as)0 1401 y(de\014ned)j(b)o(y)f(\(R)o(G\).)h(Besides)f(b)q(eing)h(highly)f (accurate,)h(\014xed-p)q(oin)o(t)g(accum)o(ulation)e(in)i(hardw)o(are)h (turns)0 1462 y(out)h(to)f(b)q(e)h(at)f(least)g(as)h(fast)g(as)g(a)g (customary)e(implem)o(en)o(tati)o(on)f(using)j(\015oating-p)q(oin)o(t)g(op)q (erations.)0 1543 y(Numerical)f(libraries)i(and)i(programming)d(languages)k (whic)o(h)d(supp)q(ort)i(K-M)f(arithmetic)e(ha)o(v)o(e)h(b)q(een)0 1603 y(a)o(v)m(ailable)12 b(for)i(a)f(n)o(um)o(b)q(er)f(of)h(y)o(ears)g(\(e.) f(g.)h(A)o(CRITH,)f(A)o(CRITH{XSC,)f(P)l(ASCAL{XSC,)i(C{XSC,)g(and)0 1663 y(F)o(OR)l(TRAN{XSC\).)21 b(Besides)h(the)g(elemen)o(tary)e(op)q (erations)k(discussed)f(ab)q(o)o(v)o(e,)h(m)o(ultiple-prec)o(ision)0 1723 y(\(\\long"\))15 b(arithmetic)d(for)i(real)f(and)i(complex)d(n)o(um)o(b) q(ers)g(and)j(in)o(terv)m(als)e(of)h(v)m(arying)g(length)g(is)g(also)g(pro-)0 1783 y(vided.)20 b(It)15 b(can)h(easily)e(b)q(e)i(implem)o(en)n(ted)c(via)k (the)f(optimal)f(dot)i(pro)q(duct.)21 b(Using)15 b(these)g(programming)0 1843 y(systems,)d(problem-solving)g(routines)i(with)f(automatic)g(result)g(v) o(eri\014cation)f(ha)o(v)o(e)h(b)q(een)g(dev)o(elop)q(ed)g(for)0 1904 y(a)k(large)f(v)m(ariet)o(y)f(of)i(standard)g(problems)e(of)i(n)o (umerical)c(analysis)k(as)g(w)o(ell)d(as)j(for)g(man)o(y)e(applications.)0 1985 y(The)k(mathematical)d(insigh)o(t)i(in)o(to)h(the)g(fundamen)o(tal)e (principles)h(of)h(computer)e(arithmetic)g(and)i(the)0 2045 y(tec)o(hnology)g(to)g(realize)f(them)f(in)i(hardw)o(are)h(ha)o(v)o(e)e(b)q (een)h(a)o(v)m(ailable)g(for)g(man)o(y)f(y)o(ears)h(and)g(should)h(b)q(e)0 2105 y(common)14 b(kno)o(wledge)j(b)o(y)f(no)o(w.)22 b(Not)17 b(to)g(incorp)q(orate)g(them)e(in)o(to)h(arithmetic)e(standards)k(has)g(pro)o (v)o(en)0 2165 y(to)h(b)q(e)h(a)f(great)g(hindrance)g(to)g(further)g (progress)h(to)o(w)o(ards)g(making)d(n)o(umerical)g(computations)h(more)0 2225 y(reliable)d(and)i(to)o(w)o(ards)f(making)f(computers)g(do)i(real)f (mathematics.)0 2306 y(Man)o(y)g(man)o(ufacturers)f(b)q(eliev)o(e)g(that)i (all)f(that)h(is)f(necessary)g(is)h(to)f(impleme)o(n)o(t)d(computer)j (arithmetic)0 2367 y(according)f(to)g(one)f(of)h(the)f(existing)g(standards.) 22 b(Ev)o(en)14 b(though)i(this)e(is)h(a)f(safe)h(approac)o(h)h(for)e(the)h (man-)0 2427 y(ufacturer,)e(it)g(do)q(es)h(not)g(serv)o(e)e(the)h(user)h(as)g (w)o(ell)e(as)h(it)g(could.)20 b(W)l(e)13 b(hop)q(e)h(that)g(forthcoming)f (arithmetic)0 2487 y(standards)k(will)e(tak)o(e)g(the)h(necessary)g(steps)g (to)g(turn)g(computers)e(in)o(to)i(rigorous)h(mathematical)c(to)q(ols.)963 2858 y(3)p eop %%Trailer end userdict /end-hook known{end-hook}if %%EOF