COMPILATION LISTING OF SEGMENT dexerfc_ Compiled by: Multics PL/I Compiler, Release 28d, of September 14, 1983 Compiled at: Honeywell LCPD Phoenix, System M Compiled on: 10/03/83 1400.9 mst Mon Options: optimize map 1 /* ****************************************************** 2* * * 3* * * 4* * Copyright (c) 1972 by Massachusetts Institute of * 5* * Technology and Honeywell Information Systems, Inc. * 6* * * 7* * * 8* ****************************************************** */ 9 10 dexerfc_: procedure (number) returns (float binary (63)); 11 12 /* compute the special error function of a double-precision floating-point number */ 13 declare (number, f, n, hn, dh, dl, dm, ph, rh, rl, rm, sh, sl, sm, th, tl, tm) float binary (63), 14 (abs, exp, erfc) builtin; 15 dcl code_ ext entry (fixed bin); 16 n = number + 0.0e0; 17 if n < -9.30630096e0 then err: do; 18 call code_ (66); 19 return (170141182.e30); 20 end err; 21 f = abs (n); 22 if f < 2.5e0 then th = erfc (n) * exp (f*f); 23 else large: do; 24 rm, th = 0.5e0 / f; 25 if f >= 1.e11 then go to done; 26 ph = f; 27 rl = 0.0e0; 28 hn = 0.5e0; 29 th = -1.e10; 30 loop: dm = hn / ph; 31 ph = dm + f; 32 rh = (rl * dm + rm * f) / ph; 33 dl = rl - rm; 34 dh = rh - rm; 35 dm = dh + dl; 36 if dm = 0.0e0 then go to dvc1; 37 sh = (dh/dm) * dl + rm; 38 if hn < 1.25e0 then go to step; 39 dl = sl - sm; 40 dh = sh - sm; 41 dm = dh + dl; 42 if dm = 0.0e0 then go to dvc2; 43 th = (dh/dm) * dl + sm; 44 if th = tm then go to done; 45 if th = tl then go to done; 46 step: hn = hn + 0.5e0; 47 rl = rm; 48 sl = sm; 49 tl = tm; 50 rm = rh; 51 sm = sh; 52 tm = th; 53 go to loop; 54 dvc1: sh = rh; 55 dvc2: th = sh; 56 done: th = 1.128379167095512574e0 * th; 57 if n < 0.0e0 then th = 2.e0 * exp (f*f) - th; 58 end large; 59 return (th); 60 end dexerfc_; SOURCE FILES USED IN THIS COMPILATION. LINE NUMBER DATE MODIFIED NAME PATHNAME 0 10/03/83 1004.9 dexerfc_.pl1 >spec>on>pl128d>dexerfc_.pl1 NAMES DECLARED IN THIS COMPILATION. IDENTIFIER OFFSET LOC STORAGE CLASS DATA TYPE ATTRIBUTES AND REFERENCES (* indicates a set context) NAMES DECLARED BY DECLARE STATEMENT. abs builtin function dcl 13 ref 21 code_ 000012 constant entry external dcl 15 ref 18 dh 000106 automatic float bin(63) dcl 13 set ref 34* 35 37 40* 41 43 dl 000110 automatic float bin(63) dcl 13 set ref 33* 35 37 39* 41 43 dm 000112 automatic float bin(63) dcl 13 set ref 30* 31 32 35* 36 37 41* 42 43 erfc builtin function dcl 13 ref 22 exp builtin function dcl 13 ref 22 57 f 000100 automatic float bin(63) dcl 13 set ref 21* 22 22 22 24 25 26 31 32 57 57 hn 000104 automatic float bin(63) dcl 13 set ref 28* 30 38 46* 46 n 000102 automatic float bin(63) dcl 13 set ref 16* 17 21 22 57 number parameter float bin(63) dcl 13 ref 10 16 ph 000114 automatic float bin(63) dcl 13 set ref 26* 30 31* 32 rh 000116 automatic float bin(63) dcl 13 set ref 32* 34 50 54 rl 000120 automatic float bin(63) dcl 13 set ref 27* 32 33 47* rm 000122 automatic float bin(63) dcl 13 set ref 24* 32 33 34 37 47 50* sh 000124 automatic float bin(63) dcl 13 set ref 37* 40 51 54* 55 sl 000126 automatic float bin(63) dcl 13 set ref 39 48* sm 000130 automatic float bin(63) dcl 13 set ref 39 40 43 48 51* th 000132 automatic float bin(63) dcl 13 set ref 22* 24* 29* 43* 44 45 52 55* 56* 56 57* 57 59 tl 000134 automatic float bin(63) dcl 13 set ref 45 49* tm 000136 automatic float bin(63) dcl 13 set ref 44 49 52* NAMES DECLARED BY EXPLICIT CONTEXT. dexerfc_ 000022 constant entry external dcl 10 done 000220 constant label dcl 56 ref 25 44 45 dvc1 000214 constant label dcl 54 set ref 36 dvc2 000216 constant label dcl 55 ref 42 err 000035 constant label dcl 17 large 000077 constant label dcl 23 loop 000116 constant label dcl 30 ref 53 step 000174 constant label dcl 46 ref 38 THERE WERE NO NAMES DECLARED BY CONTEXT OR IMPLICATION. STORAGE REQUIREMENTS FOR THIS PROGRAM. Object Text Link Symbol Defs Static Start 0 0 300 314 240 310 Length 456 240 14 126 40 0 BLOCK NAME STACK SIZE TYPE WHY NONQUICK/WHO SHARES STACK FRAME dexerfc_ 142 external procedure is an external procedure. STORAGE FOR AUTOMATIC VARIABLES. STACK FRAME LOC IDENTIFIER BLOCK NAME dexerfc_ 000100 f dexerfc_ 000102 n dexerfc_ 000104 hn dexerfc_ 000106 dh dexerfc_ 000110 dl dexerfc_ 000112 dm dexerfc_ 000114 ph dexerfc_ 000116 rh dexerfc_ 000120 rl dexerfc_ 000122 rm dexerfc_ 000124 sh dexerfc_ 000126 sl dexerfc_ 000130 sm dexerfc_ 000132 th dexerfc_ 000134 tl dexerfc_ 000136 tm dexerfc_ THE FOLLOWING EXTERNAL OPERATORS ARE USED BY THIS PROGRAM. call_ext_out return ext_entry dexp THE FOLLOWING EXTERNAL ENTRIES ARE CALLED BY THIS PROGRAM. code_ derf_$derfc_ NO EXTERNAL VARIABLES ARE USED BY THIS PROGRAM. LINE LOC LINE LOC LINE LOC LINE LOC LINE LOC LINE LOC LINE LOC 10 000016 16 000027 17 000033 18 000035 19 000045 21 000051 22 000055 24 000077 25 000103 26 000106 27 000110 28 000112 29 000114 30 000116 31 000121 32 000123 33 000133 34 000136 35 000141 36 000143 37 000144 38 000150 39 000153 40 000156 41 000161 42 000163 43 000164 44 000170 45 000172 46 000174 47 000177 48 000201 49 000203 50 000205 51 000207 52 000211 53 000213 54 000214 55 000216 56 000220 57 000223 59 000234 ----------------------------------------------------------- Historical Background This edition of the Multics software materials and documentation is provided and donated to Massachusetts Institute of Technology by Group BULL including BULL HN Information Systems Inc. as a contribution to computer science knowledge. This donation is made also to give evidence of the common contributions of Massachusetts Institute of Technology, Bell Laboratories, General Electric, Honeywell Information Systems Inc., Honeywell BULL Inc., Groupe BULL and BULL HN Information Systems Inc. to the development of this operating system. Multics development was initiated by Massachusetts Institute of Technology Project MAC (1963-1970), renamed the MIT Laboratory for Computer Science and Artificial Intelligence in the mid 1970s, under the leadership of Professor Fernando Jose Corbato. Users consider that Multics provided the best software architecture for managing computer hardware properly and for executing programs. Many subsequent operating systems incorporated Multics principles. Multics was distributed in 1975 to 2000 by Group Bull in Europe , and in the U.S. by Bull HN Information Systems Inc., as successor in interest by change in name only to Honeywell Bull Inc. and Honeywell Information Systems Inc. . ----------------------------------------------------------- Permission to use, copy, modify, and distribute these programs and their documentation for any purpose and without fee is hereby granted,provided that the below copyright notice and historical background appear in all copies and that both the copyright notice and historical background and this permission notice appear in supporting documentation, and that the names of MIT, HIS, BULL or BULL HN not be used in advertising or publicity pertaining to distribution of the programs without specific prior written permission. Copyright 1972 by Massachusetts Institute of Technology and Honeywell Information Systems Inc. Copyright 2006 by BULL HN Information Systems Inc. Copyright 2006 by Bull SAS All Rights Reserved