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_mul(v, z, w) sets *v = *z * *w computed by 31*7c478bd9Sstevel@tonic-gate * the textbook formula without regard to exceptions or special cases. 32*7c478bd9Sstevel@tonic-gate * 33*7c478bd9Sstevel@tonic-gate * On SPARC V9, _Q_cplx_lr_mul(z, w) returns *z * *w. 34*7c478bd9Sstevel@tonic-gate * 35*7c478bd9Sstevel@tonic-gate * This code is intended to be used only when CX_LIMITED_RANGE is ON; 36*7c478bd9Sstevel@tonic-gate * otherwise use _Q_cplx_mul. 37*7c478bd9Sstevel@tonic-gate */ 38*7c478bd9Sstevel@tonic-gate 39*7c478bd9Sstevel@tonic-gate #if !defined(sparc) && !defined(__sparc) 40*7c478bd9Sstevel@tonic-gate #error This code is for SPARC only 41*7c478bd9Sstevel@tonic-gate #endif 42*7c478bd9Sstevel@tonic-gate 43*7c478bd9Sstevel@tonic-gate #ifdef __sparcv9 44*7c478bd9Sstevel@tonic-gate long double _Complex 45*7c478bd9Sstevel@tonic-gate _Q_cplx_lr_mul(const long double _Complex *z, const long double _Complex *w) 46*7c478bd9Sstevel@tonic-gate { 47*7c478bd9Sstevel@tonic-gate long double _Complex v; 48*7c478bd9Sstevel@tonic-gate #else 49*7c478bd9Sstevel@tonic-gate void 50*7c478bd9Sstevel@tonic-gate _Q_cplx_lr_mul(long double _Complex *v, const long double _Complex *z, 51*7c478bd9Sstevel@tonic-gate const long double _Complex *w) 52*7c478bd9Sstevel@tonic-gate { 53*7c478bd9Sstevel@tonic-gate #endif 54*7c478bd9Sstevel@tonic-gate long double a, b, c, d; 55*7c478bd9Sstevel@tonic-gate 56*7c478bd9Sstevel@tonic-gate a = ((long double *)z)[0]; 57*7c478bd9Sstevel@tonic-gate b = ((long double *)z)[1]; 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 #ifdef __sparcv9 62*7c478bd9Sstevel@tonic-gate ((long double *)&v)[0] = (a * c - b * d); 63*7c478bd9Sstevel@tonic-gate ((long double *)&v)[1] = (a * d + b * c); 64*7c478bd9Sstevel@tonic-gate return (v); 65*7c478bd9Sstevel@tonic-gate #else 66*7c478bd9Sstevel@tonic-gate ((long double *)v)[0] = (a * c - b * d); 67*7c478bd9Sstevel@tonic-gate ((long double *)v)[1] = (a * d + b * c); 68*7c478bd9Sstevel@tonic-gate #endif 69*7c478bd9Sstevel@tonic-gate } 70