1 /* ****************************************************** 2 * * 3 * * 4 * Copyright (c) 1972 by Massachusetts Institute of * 5 * Technology and Honeywell Information Systems, Inc. * 6 * * 7 * * 8 ****************************************************** */ 9 10 dcxp2_: proc (base, exponent) returns (complex float bin (63)); 11 12 dcl (base, exponent, a, b) complex float bin (63); 13 14 dcl code_ ext entry (fixed bin (17)); 15 16 a = base; 17 b = exponent; 18 19 if a = 0.0e0 20 then do; 21 if real (b)>0.0e0 22 & imag (b) = 0.0e0 23 then goto ret; 24 25 call code_ (57); 26 goto ret; 27 end; 28 29 if b = 0.0e0 then return (1.0e0); 30 31 a = exp (log (a)*b); 32 33 ret: 34 return (a); 35 36 end dcxp2_;