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> -- replaced with ieeefp for a portable interface */ 00047 00048 #if !defined(__alpha) /* FIXME: need a replacement for alpha-Tru64 */ 00049 # include <ieeefp.h> 00050 #endif 00051 #include "quad.h" 00052 00053 /* complex *32 division */ 00054 00055 /* ==================================================================== 00056 * 00057 * FunctionName cq_div 00058 * 00059 * Description computes complex quotient of args 00060 * 00061 * ==================================================================== 00062 */ 00063 00064 qcomplex __cqdiv(long double xqreal, long double xqimag, long double yqreal, long double yqimag) 00065 { 00066 qcomplex result; 00067 00068 #if defined(__alpha) /* FIXME: need a replacement for alpha-Tru64 */ 00069 00070 result.qreal = 0; 00071 result.qimag = 0; 00072 00073 #else 00074 00075 long double tmp; 00076 unsigned int m, n; 00077 fp_except oldmask; 00078 fp_except oldsticky; 00079 00080 if ((xqreal != xqreal) || (xqimag != xqimag) || 00081 (yqreal != yqreal) || (yqimag != yqimag)) { 00082 result.qreal = __libm_qnan_ld; 00083 result.qimag = __libm_qnan_ld; 00084 return result; 00085 } 00086 00087 if ((yqreal == 0.0L) && (yqimag == 0.0L)) { 00088 result.qreal = xqreal/__libm_zero_ld; 00089 result.qimag = xqimag/__libm_zero_ld; 00090 return result; 00091 } 00092 00093 if (yqreal == 0.0L) { 00094 result.qreal = xqimag/yqimag; 00095 result.qimag = -(xqreal/yqimag); 00096 return result; 00097 } 00098 00099 if (yqimag == 0.0L) { 00100 result.qreal = xqreal/yqreal; 00101 result.qimag = xqimag/yqreal; 00102 return result; 00103 } 00104 00105 if (__qabs(yqreal) <= __qabs(yqimag)) { 00106 /* turn off traps on underflow while computing 00107 yqreal/yqimag 00108 */ 00109 00110 /* record fpu exception handling state */ 00111 oldmask = fpgetmask(); 00112 oldsticky = fpgetsticky(); 00113 00114 /* suppress underflow exceptions */ 00115 fpsetmask(oldmask & ~FP_X_UFL); 00116 00117 tmp = yqreal/yqimag; 00118 00119 /* restore underflow handling to its previous state */ 00120 fpsetmask(oldmask); 00121 fpsetsticky(oldsticky); 00122 00123 result.qreal = (xqimag + xqreal*tmp)/(yqimag + yqreal*tmp); 00124 result.qimag = (-xqreal + xqimag*tmp)/(yqimag + yqreal*tmp); 00125 return result; 00126 } 00127 00128 /* turn off traps on underflow while computing 00129 yqimag/yqreal 00130 */ 00131 /* record fpu exception handling state */ 00132 oldmask = fpgetmask(); 00133 oldsticky = fpgetsticky(); 00134 00135 /* suppress underflow exceptions */ 00136 fpsetmask(oldmask & ~FP_X_UFL); 00137 00138 tmp = yqimag/yqreal; 00139 00140 /* restore underflow handling to its previous state */ 00141 fpsetmask(oldmask); 00142 fpsetsticky(oldsticky); 00143 00144 result.qreal = (xqreal + xqimag*tmp)/(yqreal + yqimag*tmp); 00145 result.qimag = (xqimag - xqreal*tmp)/(yqreal + yqimag*tmp); 00146 00147 #endif /* __alpha */ 00148 00149 return result; 00150 } 00151 00152 void __cq_div(qcomplex *result, qcomplex *x, qcomplex *y) 00153 { 00154 *result = __cqdiv(x->qreal, x->qimag, y->qreal, y->qimag); 00155 }