72 	ldcomplex ans;  in casinl()  local
88 			LD_IM(ans) = x;  in casinl()
92 					LD_RE(ans) = pi_4 + pi_4_l;  in casinl()
94 					LD_RE(ans) = y + y;  in casinl()
96 				LD_RE(ans) = pi_2 + pi_2_l;  in casinl()
105 				LD_IM(ans) = y + y;  in casinl()
106 				LD_RE(ans) = x + x;  in casinl()
111 				LD_IM(ans) = LD_RE(ans) = x + y;  in casinl()
115 			LD_RE(ans) = -LD_RE(ans);  in casinl()
117 			LD_IM(ans) = -LD_IM(ans);  in casinl()
118 		return (ans);  in casinl()
127 			LD_IM(ans) = y;  in casinl()
128 			LD_RE(ans) = zero;  in casinl()
130 			LD_IM(ans) = x + y;  in casinl()
132 				LD_RE(ans) = x;  in casinl()
134 				LD_RE(ans) = y;  in casinl()
137 			LD_RE(ans) = -LD_RE(ans);  in casinl()
139 			LD_IM(ans) = -LD_IM(ans);  in casinl()
140 		return (ans);  in casinl()
145 			LD_RE(ans) = asinl(x);  in casinl()
146 			LD_IM(ans) = zero;  in casinl()
148 			LD_RE(ans) = pi_2 + pi_2_l;  in casinl()
150 				LD_IM(ans) = ln2 + logl(x);  in casinl()
152 				LD_IM(ans) = logl(x + sqrtl((x - one) * (x +  in casinl()
156 				LD_IM(ans) = log1pl(xm1 + sqrtl(xm1 * (x +  in casinl()
162 			LD_RE(ans) = asinl(x);  in casinl()
163 			LD_IM(ans) = y / sqrtl((one + x) * (one - x));  in casinl()
165 			LD_RE(ans) = pi_2 + pi_2_l;  in casinl()
167 				LD_IM(ans) = ln2 + logl(x);  in casinl()
169 				LD_IM(ans) = logl(x + sqrtl((x - one) * (x +  in casinl()
172 				LD_IM(ans) = log1pl((x - one) + sqrtl((x -  in casinl()
177 		LD_RE(ans) = pi_2 - (t - pi_2_l);  in casinl()
178 		LD_IM(ans) = t;  in casinl()
180 		LD_RE(ans) = x / y;	/* need to fix underflow cases */  in casinl()
181 		LD_IM(ans) = ln2 + logl(y);  in casinl()
185 		LD_RE(ans) = atanl(t);  in casinl()
186 		LD_IM(ans) = ln2 + logl(y) + half * log1pl(t * t);  in casinl()
190 		LD_RE(ans) = x / A;	/* may underflow */  in casinl()
192 			LD_IM(ans) = logl(y + A);  in casinl()
194 			LD_IM(ans) = half * log1pl((y + y) * (y + A));  in casinl()
204 			LD_RE(ans) = asinl(B);  in casinl()
208 				LD_RE(ans) = atanl(x / sqrtl(half * Apx * (y2 /  in casinl()
211 				LD_RE(ans) = atanl(x / (y * sqrtl(half * (Apx /  in casinl()
220 			LD_IM(ans) = log1pl(Am1 + sqrtl(Am1 * (A + one)));  in casinl()
222 			LD_IM(ans) = logl(A + sqrtl(A * A - one));  in casinl()
227 		LD_RE(ans) = -LD_RE(ans);  in casinl()
229 		LD_IM(ans) = -LD_IM(ans);  in casinl()
231 	return (ans);  in casinl()