Open64 (mfef90, whirl2f, and IR tools): osprey1.0/libF77/pow

Go to the documentation of this file.
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.1 of the GNU Lesser General Public License 
00007   as 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 Lesser General Public 
00021   License along with this program; if not, write the Free Software 
00022   Foundation, Inc., 59 Temple Place - Suite 330, Boston MA 02111-1307, 
00023   USA.
00024 
00025   Contact information:  Silicon Graphics, Inc., 1600 Amphitheatre Pky,
00026   Mountain View, CA 94043, or:
00027 
00028   http://www.sgi.com
00029 
00030   For further information regarding this notice, see:
00031 
00032   http://oss.sgi.com/projects/GenInfo/NoticeExplan
00033 
00034 */
00035 
00036 
00037 /* $Header: /m_home/m_utkej/Argonne/cvs2svn/cvs/Open64/osprey1.0/libF77/pow_di.c,v 1.1.1.1 2002-05-22 20:09:13 dsystem Exp $ */
00038 #include <math.h>
00039 #include "cmplrs/host.h"
00040 #include "moremath.h"
00041 
00042 double_t __powdi(double_t x, int32 n)
00043 {
00044 double_t pow;
00045 
00046 pow = 1;
00047 
00048 if(n != 0)
00049         {
00050         if(n < 0)
00051                 {
00052                 if(x == 0)
00053                         {
00054                         return(pow);
00055                         }
00056                 n = -n;
00057                 x = 1/x;
00058                 }
00059         for( ; ; )
00060                 {
00061                 if(n & 01)
00062                         pow *= x;
00063                 if(n >>= 1)
00064                         x *= x;
00065                 else
00066                         break;
00067                 }
00068         }
00069 return(pow);
00070 }
00071 
00072 double_t __powdl(double_t x, int64 n)
00073 {
00074 double_t pow;
00075 
00076 pow = 1;
00077 
00078 if(n != 0)
00079         {
00080         if(n < 0)
00081                 {
00082                 if(x == 0)
00083                         {
00084                         return(pow);
00085                         }
00086                 n = -n;
00087                 x = 1/x;
00088                 }
00089         for( ; ; )
00090                 {
00091                 if(n & 01)
00092                         pow *= x;
00093                 if(n >>= 1)
00094                         x *= x;
00095                 else
00096                         break;
00097                 }
00098         }
00099 return(pow);
00100 }
00101 
00102 /* By-reference versions for backward compatibility. */
00103 
00104 double_t pow_di(double_t *ap, int32 *bp)
00105 {
00106 double_t pow;
00107 pow=__powdi(*ap,*bp);
00108 return pow;
00109 }
00110 
00111 double_t pow_dl(double_t *ap, int64 *bp)
00112 {
00113 double_t pow;
00114 pow=__powdl(*ap,*bp);
00115 return pow;
00116 }