workbench/libs/mathieeesingtrans/ieeesppow.c

   1 /*
   2     Copyright © 1995-2003, The AROS Development Team. All rights reserved.
   3     $Id$
   4 */
   5
   6 #include "mathieeesingtrans_intern.h"
   7
   8 /*****************************************************************************
   9
  10     NAME */
  11
  12         AROS_LH2(float, IEEESPPow,
  13
  14 /*  SYNOPSIS */
  15         AROS_LHA(float, x, D1),
  16         AROS_LHA(float, y, D0),
  17
  18 /*  LOCATION */
  19         struct Library *, MathIeeeSingTransBase, 15, MathIeeeSingTrans)
  20
  21 /*  FUNCTION
  22       Calculate y raised to the x power (y^x)
  23
  24     INPUTS
  25
  26     RESULT
  27         IEEE single precision floating point number
  28
  29         flags:
  30         zero     : result is zero
  31         negative : result is negative
  32         overflow : result is too big
  33
  34     BUGS
  35
  36     INTERNALS
  37
  38 *****************************************************************************/
  39 {
  40     AROS_LIBFUNC_INIT
  41
  42     /*
  43         a ^ b  = e^(b * ln a )
  44         y ^ x  = e^(x * ln y )
  45     */
  46     LONG Res;
  47
  48     /* y^x is illegal if y<0 and x is not an integer-value */
  49     if (y < 0 && x != IEEESPCeil(x) ) return 0;
  50
  51     Res = IEEESPLog( y & (IEEESPMantisse_Mask + IEEESPExponent_Mask) );
  52     Res = IEEESPMul(Res, x);
  53     Res = IEEESPExp(Res);
  54
  55     /*
  56         if y < 0 and x was and even integer, the result is positive, otherwise
  57         it is negative.
  58     */
  59     if (y < 0 && TRUE == intern_IEEESPisodd(x) ) Res |= IEEESPSign_Mask;
  60
  61     if (0 == Res)
  62     {
  63         SetSR(Zero_Bit, Zero_Bit | Negative_Bit | Overflow_Bit);
  64         return 0;
  65     }
  66
  67     SetSR(0, Zero_Bit | Negative_Bit | Overflow_Bit);
  68
  69     if ( Res < 0)
  70     {
  71         SetSR(Negative_Bit, Zero_Bit | Negative_Bit | Overflow_Bit);
  72     }
  73
  74     if ( IEEESP_Pinfty == (Res & (IEEESPMantisse_Mask + IEEESPExponent_Mask)) )
  75     {
  76         /* don`t touch the Negative_Bit now!*/
  77         SetSR(Overflow_Bit, Zero_Bit | Overflow_Bit);
  78     }
  79
  80     return Res;
  81
  82     AROS_LIBFUNC_EXIT
  83 }