1/* 2 * CDDL HEADER START 3 * 4 * The contents of this file are subject to the terms of the 5 * Common Development and Distribution License (the "License"). 6 * You may not use this file except in compliance with the License. 7 * 8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE 9 * or http://www.opensolaris.org/os/licensing. 10 * See the License for the specific language governing permissions 11 * and limitations under the License. 12 * 13 * When distributing Covered Code, include this CDDL HEADER in each 14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE. 15 * If applicable, add the following below this CDDL HEADER, with the 16 * fields enclosed by brackets "[]" replaced with your own identifying 17 * information: Portions Copyright [yyyy] [name of copyright owner] 18 * 19 * CDDL HEADER END 20 */ 21 22/* 23 * Copyright 2008 Sun Microsystems, Inc. All rights reserved. 24 * Use is subject to license terms. 25 */ 26 27#pragma ident "%Z%%M% %I% %E% SMI" 28 29#include "lint.h" 30#include "base_conversion.h" 31#include <sys/isa_defs.h> 32 33#define CSR_DEFAULT 0x1f80 34 35/* The following should be coded as inline expansion templates. */ 36 37/* 38 * Multiplies two normal or subnormal doubles, returns result and exceptions. 39 */ 40double 41__mul_set(double x, double y, int *pe) { 42 extern void _putmxcsr(), _getmxcsr(); 43 int csr; 44 double z; 45 46 _putmxcsr(CSR_DEFAULT); 47 z = x * y; 48 _getmxcsr(&csr); 49 if ((csr & 0x3f) == 0) { 50 *pe = 0; 51 } else { 52 /* Result may not be exact. */ 53 *pe = 1; 54 } 55 return (z); 56} 57 58/* 59 * Divides two normal or subnormal doubles x/y, returns result and exceptions. 60 */ 61double 62__div_set(double x, double y, int *pe) { 63 extern void _putmxcsr(), _getmxcsr(); 64 int csr; 65 double z; 66 67 _putmxcsr(CSR_DEFAULT); 68 z = x / y; 69 _getmxcsr(&csr); 70 if ((csr & 0x3f) == 0) { 71 *pe = 0; 72 } else { 73 *pe = 1; 74 } 75 return (z); 76} 77 78double 79__dabs(double *d) 80{ 81 /* should use hardware fabs instruction */ 82 return ((*d < 0.0) ? -*d : *d); 83} 84 85/* 86 * Returns IEEE mode/status and 87 * sets up standard environment for base conversion. 88 */ 89void 90__get_ieee_flags(__ieee_flags_type *b) { 91 extern void _getmxcsr(), _putmxcsr(); 92 93 _getmxcsr(&b->status); 94 95 /* round-to-nearest, all exceptions masked, gradual underflow */ 96 _putmxcsr(CSR_DEFAULT); 97} 98 99/* 100 * Restores previous IEEE mode/status 101 */ 102void 103__set_ieee_flags(__ieee_flags_type *b) { 104 extern void _putmxcsr(); 105 106 _putmxcsr(b->status); 107} 108