1 /* ******************************************************
 2    *                                                    *
 3    *                                                    *
 4    * Copyright (c) 1972 by Massachusetts Institute of   *
 5    * Technology and Honeywell Information Systems, Inc. *
 6    *                                                    *
 7    *                                                    *
 8    ****************************************************** */
 9 
10 dcasin_: proc (number) returns (complex float bin (63));
11 
12 dcl (number, a, b, c) complex float bin (63);
13 dcl       (imag, log, real, sqrt) builtin;
14 
15           real (a) = -imag (number);
16           imag (a) = real (number);
17 
18           b = 1.0e0;
19 
20 trig:
21           c = -1.0e0i;
22 
23 ret:
24           return (log (sqrt (a*a+b)+a)*c);
25 
26 dcacos_: entry (number) returns (complex float bin (63));
27 
28           a = number;
29           b = -1.0e0;
30 
31           goto trig;
32 
33 dcasinh_: entry (number) returns (complex float bin (63));
34 
35           b = 1.0e0;
36 
37 hyper:
38           a = number;
39 
40           c = 1.0e0;
41 
42           goto ret;
43 
44 dcacosh_: entry (number) returns (complex float bin (63));
45 
46           b = -1.0e0;
47           goto hyper;
48 
49      end dcasin_;