xref: /illumos-gate/usr/src/lib/libc/sparc/fp/_Q_cplx_lr_div_rx.c (revision 7c478bd95313f5f23a4c958a745db2134aa03244)
1*7c478bd9Sstevel@tonic-gate /*
2*7c478bd9Sstevel@tonic-gate  * CDDL HEADER START
3*7c478bd9Sstevel@tonic-gate  *
4*7c478bd9Sstevel@tonic-gate  * The contents of this file are subject to the terms of the
5*7c478bd9Sstevel@tonic-gate  * Common Development and Distribution License, Version 1.0 only
6*7c478bd9Sstevel@tonic-gate  * (the "License").  You may not use this file except in compliance
7*7c478bd9Sstevel@tonic-gate  * with the License.
8*7c478bd9Sstevel@tonic-gate  *
9*7c478bd9Sstevel@tonic-gate  * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
10*7c478bd9Sstevel@tonic-gate  * or http://www.opensolaris.org/os/licensing.
11*7c478bd9Sstevel@tonic-gate  * See the License for the specific language governing permissions
12*7c478bd9Sstevel@tonic-gate  * and limitations under the License.
13*7c478bd9Sstevel@tonic-gate  *
14*7c478bd9Sstevel@tonic-gate  * When distributing Covered Code, include this CDDL HEADER in each
15*7c478bd9Sstevel@tonic-gate  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
16*7c478bd9Sstevel@tonic-gate  * If applicable, add the following below this CDDL HEADER, with the
17*7c478bd9Sstevel@tonic-gate  * fields enclosed by brackets "[]" replaced with your own identifying
18*7c478bd9Sstevel@tonic-gate  * information: Portions Copyright [yyyy] [name of copyright owner]
19*7c478bd9Sstevel@tonic-gate  *
20*7c478bd9Sstevel@tonic-gate  * CDDL HEADER END
21*7c478bd9Sstevel@tonic-gate  */
22*7c478bd9Sstevel@tonic-gate /*
23*7c478bd9Sstevel@tonic-gate  * Copyright 2003 Sun Microsystems, Inc.  All rights reserved.
24*7c478bd9Sstevel@tonic-gate  * Use is subject to license terms.
25*7c478bd9Sstevel@tonic-gate  */
26*7c478bd9Sstevel@tonic-gate 
27*7c478bd9Sstevel@tonic-gate #pragma ident	"%Z%%M%	%I%	%E% SMI"
28*7c478bd9Sstevel@tonic-gate 
29*7c478bd9Sstevel@tonic-gate /*
30*7c478bd9Sstevel@tonic-gate  * On SPARC V8, _Q_cplx_lr_div_rx(v, a, w) sets *v = *a / *w computed
31*7c478bd9Sstevel@tonic-gate  * by the textbook formula without regard to exceptions or special
32*7c478bd9Sstevel@tonic-gate  * cases.
33*7c478bd9Sstevel@tonic-gate  *
34*7c478bd9Sstevel@tonic-gate  * On SPARC V9, _Q_cplx_lr_div_rx(a, w) returns *a / *w.
35*7c478bd9Sstevel@tonic-gate  *
36*7c478bd9Sstevel@tonic-gate  * This code is intended to be used only when CX_LIMITED_RANGE is ON;
37*7c478bd9Sstevel@tonic-gate  * otherwise use _Q_cplx_div_rx.
38*7c478bd9Sstevel@tonic-gate  */
39*7c478bd9Sstevel@tonic-gate 
40*7c478bd9Sstevel@tonic-gate #if !defined(sparc) && !defined(__sparc)
41*7c478bd9Sstevel@tonic-gate #error This code is for SPARC only
42*7c478bd9Sstevel@tonic-gate #endif
43*7c478bd9Sstevel@tonic-gate 
44*7c478bd9Sstevel@tonic-gate #ifdef __sparcv9
45*7c478bd9Sstevel@tonic-gate long double _Complex
46*7c478bd9Sstevel@tonic-gate _Q_cplx_lr_div_rx(const long double *pa, const long double _Complex *w)
47*7c478bd9Sstevel@tonic-gate {
48*7c478bd9Sstevel@tonic-gate 	long double _Complex	v;
49*7c478bd9Sstevel@tonic-gate #else
50*7c478bd9Sstevel@tonic-gate void
51*7c478bd9Sstevel@tonic-gate _Q_cplx_lr_div_rx(long double _Complex *v, const long double *pa,
52*7c478bd9Sstevel@tonic-gate 	const long double _Complex *w)
53*7c478bd9Sstevel@tonic-gate {
54*7c478bd9Sstevel@tonic-gate #endif
55*7c478bd9Sstevel@tonic-gate 	long double	a, c, d;
56*7c478bd9Sstevel@tonic-gate 
57*7c478bd9Sstevel@tonic-gate 	a = *pa;
58*7c478bd9Sstevel@tonic-gate 	c = ((long double *)w)[0];
59*7c478bd9Sstevel@tonic-gate 	d = ((long double *)w)[1];
60*7c478bd9Sstevel@tonic-gate 
61*7c478bd9Sstevel@tonic-gate 	a /= (c * c + d * d);
62*7c478bd9Sstevel@tonic-gate 
63*7c478bd9Sstevel@tonic-gate #ifdef __sparcv9
64*7c478bd9Sstevel@tonic-gate 	((long double *)&v)[0] = a * c;
65*7c478bd9Sstevel@tonic-gate 	((long double *)&v)[1] = a * -d;
66*7c478bd9Sstevel@tonic-gate 	return (v);
67*7c478bd9Sstevel@tonic-gate #else
68*7c478bd9Sstevel@tonic-gate 	((long double *)v)[0] = a * c;
69*7c478bd9Sstevel@tonic-gate 	((long double *)v)[1] = a * -d;
70*7c478bd9Sstevel@tonic-gate #endif
71*7c478bd9Sstevel@tonic-gate }
72