Open64 (mfef90, whirl2f, and IR tools)
TAG: version-openad; SVN changeset: 916
|
00001 /* 00002 00003 Copyright (C) 2000, 2001 Silicon Graphics, Inc. All Rights Reserved. 00004 00005 This program is free software; you can redistribute it and/or modify it 00006 under the terms of version 2 of the GNU General Public License as 00007 published by the Free Software Foundation. 00008 00009 This program is distributed in the hope that it would be useful, but 00010 WITHOUT ANY WARRANTY; without even the implied warranty of 00011 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. 00012 00013 Further, this software is distributed without any warranty that it is 00014 free of the rightful claim of any third person regarding infringement 00015 or the like. Any license provided herein, whether implied or 00016 otherwise, applies only to this software file. Patent licenses, if 00017 any, provided herein do not apply to combinations of this program with 00018 other software, or any other product whatsoever. 00019 00020 You should have received a copy of the GNU General Public License along 00021 with this program; if not, write the Free Software Foundation, Inc., 59 00022 Temple Place - Suite 330, Boston MA 02111-1307, USA. 00023 00024 Contact information: Silicon Graphics, Inc., 1600 Amphitheatre Pky, 00025 Mountain View, CA 94043, or: 00026 00027 http://www.sgi.com 00028 00029 For further information regarding this notice, see: 00030 00031 http://oss.sgi.com/projects/GenInfo/NoticeExplan 00032 00033 */ 00034 00035 00036 /* ==================================================================== 00037 * ==================================================================== 00038 * 00039 * 00040 * Description: source code for complex quad division 00041 * 00042 * ==================================================================== 00043 * ==================================================================== 00044 */ 00045 00046 #include <sys/fpu.h> 00047 #include "quad.h" 00048 00049 /* complex *32 division */ 00050 00051 /* ==================================================================== 00052 * 00053 * FunctionName cq_div 00054 * 00055 * Description computes complex quotient of args 00056 * 00057 * ==================================================================== 00058 */ 00059 00060 qcomplex __cqdiv(long double xqreal, long double xqimag, long double yqreal, long double yqimag) 00061 { 00062 long double tmp; 00063 unsigned int m, n; 00064 qcomplex result; 00065 00066 if ((xqreal != xqreal) || (xqimag != xqimag) || 00067 (yqreal != yqreal) || (yqimag != yqimag)) { 00068 result.qreal = __libm_qnan_ld; 00069 result.qimag = __libm_qnan_ld; 00070 return result; 00071 } 00072 00073 if ((yqreal == 0.0L) && (yqimag == 0.0L)) { 00074 result.qreal = xqreal/__libm_zero_ld; 00075 result.qimag = xqimag/__libm_zero_ld; 00076 return result; 00077 } 00078 00079 if (yqreal == 0.0L) { 00080 result.qreal = xqimag/yqimag; 00081 result.qimag = -(xqreal/yqimag); 00082 return result; 00083 } 00084 00085 if (yqimag == 0.0L) { 00086 result.qreal = xqreal/yqreal; 00087 result.qimag = xqimag/yqreal; 00088 return result; 00089 } 00090 00091 if (__qabs(yqreal) <= __qabs(yqimag)) { 00092 /* turn off traps on underflow while computing 00093 yqreal/yqimag 00094 */ 00095 00096 m = get_fpc_csr(); 00097 n = (m & 0xfffffeff); 00098 (void)set_fpc_csr(n); 00099 tmp = yqreal/yqimag; 00100 (void)set_fpc_csr(m); 00101 result.qreal = (xqimag + xqreal*tmp)/(yqimag + yqreal*tmp); 00102 result.qimag = (-xqreal + xqimag*tmp)/(yqimag + yqreal*tmp); 00103 return result; 00104 } 00105 00106 /* turn off traps on underflow while computing 00107 yqimag/yqreal 00108 */ 00109 00110 m = get_fpc_csr(); 00111 n = (m & 0xfffffeff); 00112 (void)set_fpc_csr(n); 00113 00114 tmp = yqimag/yqreal; 00115 (void)set_fpc_csr(m); 00116 00117 result.qreal = (xqreal + xqimag*tmp)/(yqreal + yqimag*tmp); 00118 result.qimag = (xqimag - xqreal*tmp)/(yqreal + yqimag*tmp); 00119 00120 return result; 00121 } 00122 00123 void __cq_div(qcomplex *result, qcomplex *x, qcomplex *y) 00124 { 00125 *result = __cqdiv(x->qreal, x->qimag, y->qreal, y->qimag); 00126 }