17c478bd9Sstevel@tonic-gate /*
27c478bd9Sstevel@tonic-gate * CDDL HEADER START
37c478bd9Sstevel@tonic-gate *
47c478bd9Sstevel@tonic-gate * The contents of this file are subject to the terms of the
57c478bd9Sstevel@tonic-gate * Common Development and Distribution License, Version 1.0 only
67c478bd9Sstevel@tonic-gate * (the "License"). You may not use this file except in compliance
77c478bd9Sstevel@tonic-gate * with the License.
87c478bd9Sstevel@tonic-gate *
97c478bd9Sstevel@tonic-gate * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
107c478bd9Sstevel@tonic-gate * or http://www.opensolaris.org/os/licensing.
117c478bd9Sstevel@tonic-gate * See the License for the specific language governing permissions
127c478bd9Sstevel@tonic-gate * and limitations under the License.
137c478bd9Sstevel@tonic-gate *
147c478bd9Sstevel@tonic-gate * When distributing Covered Code, include this CDDL HEADER in each
157c478bd9Sstevel@tonic-gate * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
167c478bd9Sstevel@tonic-gate * If applicable, add the following below this CDDL HEADER, with the
177c478bd9Sstevel@tonic-gate * fields enclosed by brackets "[]" replaced with your own identifying
187c478bd9Sstevel@tonic-gate * information: Portions Copyright [yyyy] [name of copyright owner]
197c478bd9Sstevel@tonic-gate *
207c478bd9Sstevel@tonic-gate * CDDL HEADER END
217c478bd9Sstevel@tonic-gate */
227c478bd9Sstevel@tonic-gate /*
237c478bd9Sstevel@tonic-gate * Copyright 2003 Sun Microsystems, Inc. All rights reserved.
247c478bd9Sstevel@tonic-gate * Use is subject to license terms.
257c478bd9Sstevel@tonic-gate */
267c478bd9Sstevel@tonic-gate
277c478bd9Sstevel@tonic-gate /*
287c478bd9Sstevel@tonic-gate * On SPARC V8, _Q_cplx_lr_mul(v, z, w) sets *v = *z * *w computed by
297c478bd9Sstevel@tonic-gate * the textbook formula without regard to exceptions or special cases.
307c478bd9Sstevel@tonic-gate *
317c478bd9Sstevel@tonic-gate * On SPARC V9, _Q_cplx_lr_mul(z, w) returns *z * *w.
327c478bd9Sstevel@tonic-gate *
337c478bd9Sstevel@tonic-gate * This code is intended to be used only when CX_LIMITED_RANGE is ON;
347c478bd9Sstevel@tonic-gate * otherwise use _Q_cplx_mul.
357c478bd9Sstevel@tonic-gate */
367c478bd9Sstevel@tonic-gate
377c478bd9Sstevel@tonic-gate #if !defined(sparc) && !defined(__sparc)
387c478bd9Sstevel@tonic-gate #error This code is for SPARC only
397c478bd9Sstevel@tonic-gate #endif
407c478bd9Sstevel@tonic-gate
417c478bd9Sstevel@tonic-gate #ifdef __sparcv9
427c478bd9Sstevel@tonic-gate long double _Complex
_Q_cplx_lr_mul(const long double _Complex * z,const long double _Complex * w)437c478bd9Sstevel@tonic-gate _Q_cplx_lr_mul(const long double _Complex *z, const long double _Complex *w)
447c478bd9Sstevel@tonic-gate {
45*1d9fea2aSToomas Soome long double _Complex v = 0;
467c478bd9Sstevel@tonic-gate #else
477c478bd9Sstevel@tonic-gate void
487c478bd9Sstevel@tonic-gate _Q_cplx_lr_mul(long double _Complex *v, const long double _Complex *z,
497c478bd9Sstevel@tonic-gate const long double _Complex *w)
507c478bd9Sstevel@tonic-gate {
517c478bd9Sstevel@tonic-gate #endif
527c478bd9Sstevel@tonic-gate long double a, b, c, d;
537c478bd9Sstevel@tonic-gate
547c478bd9Sstevel@tonic-gate a = ((long double *)z)[0];
557c478bd9Sstevel@tonic-gate b = ((long double *)z)[1];
567c478bd9Sstevel@tonic-gate c = ((long double *)w)[0];
577c478bd9Sstevel@tonic-gate d = ((long double *)w)[1];
587c478bd9Sstevel@tonic-gate
597c478bd9Sstevel@tonic-gate #ifdef __sparcv9
607c478bd9Sstevel@tonic-gate ((long double *)&v)[0] = (a * c - b * d);
617c478bd9Sstevel@tonic-gate ((long double *)&v)[1] = (a * d + b * c);
627c478bd9Sstevel@tonic-gate return (v);
637c478bd9Sstevel@tonic-gate #else
647c478bd9Sstevel@tonic-gate ((long double *)v)[0] = (a * c - b * d);
657c478bd9Sstevel@tonic-gate ((long double *)v)[1] = (a * d + b * c);
667c478bd9Sstevel@tonic-gate #endif
677c478bd9Sstevel@tonic-gate }
68