Reorganize and commit all the modified TACLeBench code and run scripts

1 files changed, 423 insertions, 0 deletions

diff --git a/all_pairs/source/audiobeam/audiobeamlibm.c b/all_pairs/source/audiobeam/audiobeamlibm.c
new file mode 100644
index 0000000..259f40b
--- /dev/null
+++ b/all_pairs/source/audiobeam/audiobeamlibm.c
@@ -0,0 +1,423 @@
+/*
+
+  This program is part of the TACLeBench benchmark suite.
+  Version V 2.0
+
+  Name: audiobeamlibm.c
+
+  Author: Ian Lance Taylor and J.T. Conklin
+
+  Function: IEEE754 software library routines.
+
+  Source: Sun Microsystems and Cygnus
+
+  Original name: Unknown
+
+  Changes: No major functional changes.
+
+  License: See the terms below.
+
+*/
+
+/*
+   ====================================================
+   Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+
+   Developed at SunPro, a Sun Microsystems, Inc. business.
+   Permission to use, copy, modify, and distribute this
+   software is freely granted, provided that this notice
+   is preserved.
+   ====================================================
+*/
+
+
+#include "audiobeamlibm.h"
+#include "audiobeamlibmath.h"
+
+
+static const int audiobeam_npio2_hw[] = {
+  0x3fc90f00, 0x40490f00, 0x4096cb00, 0x40c90f00, 0x40fb5300, 0x4116cb00,
+  0x412fed00, 0x41490f00, 0x41623100, 0x417b5300, 0x418a3a00, 0x4196cb00,
+  0x41a35c00, 0x41afed00, 0x41bc7e00, 0x41c90f00, 0x41d5a000, 0x41e23100,
+  0x41eec200, 0x41fb5300, 0x4203f200, 0x420a3a00, 0x42108300, 0x4216cb00,
+  0x421d1400, 0x42235c00, 0x4229a500, 0x422fed00, 0x42363600, 0x423c7e00,
+  0x4242c700, 0x42490f00
+};
+
+static const float
+audiobeam_invpio2 =  6.3661980629e-01f, /* 0x3f22f984 */
+audiobeam_pio2_1  =  1.5707855225e+00f, /* 0x3fc90f80 */
+audiobeam_pio2_1t =  1.0804334124e-05f, /* 0x37354443 */
+audiobeam_pio2_2  =  1.0804273188e-05f, /* 0x37354400 */
+audiobeam_pio2_2t =  6.0770999344e-11f, /* 0x2e85a308 */
+audiobeam_pio2_3  =  6.0770943833e-11f, /* 0x2e85a300 */
+audiobeam_pio2_3t =  6.1232342629e-17f; /* 0x248d3132 */
+
+static const float
+audiobeam_C1  =  4.1666667908e-02f, /* 0x3d2aaaab */
+audiobeam_C2  = -1.3888889225e-03f, /* 0xbab60b61 */
+audiobeam_C3  =  2.4801587642e-05f, /* 0x37d00d01 */
+audiobeam_C4  = -2.7557314297e-07f, /* 0xb493f27c */
+audiobeam_C5  =  2.0875723372e-09f, /* 0x310f74f6 */
+audiobeam_C6  = -1.1359647598e-11f; /* 0xad47d74e */
+
+static const float
+audiobeam_S1  = -1.6666667163e-01f, /* 0xbe2aaaab */
+audiobeam_S2  =  8.3333337680e-03f, /* 0x3c088889 */
+audiobeam_S3  = -1.9841270114e-04f, /* 0xb9500d01 */
+audiobeam_S4  =  2.7557314297e-06f, /* 0x3638ef1b */
+audiobeam_S5  = -2.5050759689e-08f, /* 0xb2d72f34 */
+audiobeam_S6  =  1.5896910177e-10f; /* 0x2f2ec9d3 */
+
+static const float
+audiobeam_two25   =  3.355443200e+07f,  /* 0x4c000000 */
+audiobeam_twom25  =  2.9802322388e-08f; /* 0x33000000 */
+
+
+int audiobeam___ieee754_rem_pio2f( float x, float *y )
+{
+  float z, w, t, r, fn;
+  int i, j, n = 0, ix, hx;
+
+  AUDIOBEAM_GET_FLOAT_WORD( hx, x );
+  ix = hx & 0x7fffffff;
+  if ( ix <= 0x3f490fd8 ) {
+    y[0] = x;
+    y[1] = 0;
+    return 0;
+  }
+  if ( ix < 0x4016cbe4 ) {
+    if ( hx > 0 ) {
+      z = x - audiobeam_pio2_1;
+      if ( ( ix & 0xfffffff0 ) != 0x3fc90fd0 ) {
+        y[0] = z - audiobeam_pio2_1t;
+        y[1] = ( z - y[0] ) - audiobeam_pio2_1t;
+      } else {
+        z -= audiobeam_pio2_2;
+        y[0] = z - audiobeam_pio2_2t;
+        y[1] = ( z - y[0] ) - audiobeam_pio2_2t;
+      }
+      return 1;
+    } else {
+      z = x + audiobeam_pio2_1;
+      if ( ( ix & 0xfffffff0 ) != 0x3fc90fd0 ) {
+        y[0] = z + audiobeam_pio2_1t;
+        y[1] = ( z - y[0] ) + audiobeam_pio2_1t;
+      } else {
+        z += audiobeam_pio2_2;
+        y[0] = z + audiobeam_pio2_2t;
+        y[1] = ( z - y[0] ) + audiobeam_pio2_2t;
+      }
+      return -1;
+    }
+  }
+  if ( ix <= 0x43490f80 ) {
+    t  = audiobeam_fabsf( x );
+    n  = ( int ) ( t * audiobeam_invpio2 + audiobeam_half );
+    fn = ( float )n;
+    r  = t - fn * audiobeam_pio2_1;
+    w  = fn * audiobeam_pio2_1t;
+    if ( n < 32 && ( int )( ix & 0xffffff00 ) != audiobeam_npio2_hw[n - 1] )
+      y[0] = r - w;
+    else {
+      unsigned int high;
+      j  = ix >> 23;
+      y[0] = r - w;
+      AUDIOBEAM_GET_FLOAT_WORD( high, y[0] );
+      i = j - ( ( high >> 23 ) & 0xff );
+      if ( i > 8 ) {
+        t  = r;
+        w  = fn * audiobeam_pio2_2;
+        r  = t - w;
+        w  = fn * audiobeam_pio2_2t - ( ( t - r ) - w );
+        y[0] = r - w;
+        AUDIOBEAM_GET_FLOAT_WORD( high, y[0] );
+        i = j - ( ( high >> 23 ) & 0xff );
+        if ( i > 25 )  {
+          t  = r;
+          w  = fn * audiobeam_pio2_3;
+          r  = t - w;
+          w  = fn * audiobeam_pio2_3t - ( ( t - r ) - w );
+          y[0] = r - w;
+        }
+      }
+    }
+    y[1] = ( r - y[0] ) - w;
+    if ( hx < 0 )  {
+      y[0] = -y[0];
+      y[1] = -y[1];
+      return -n;
+    } else   return n;
+  }
+  if ( ix >= 0x7f800000 ) {
+    y[0] = y[1] = x - x;
+    return 0;
+  }
+
+  return n;
+}
+
+
+float audiobeam___kernel_cosf( float x, float y )
+{
+  float a, hz, z, r, qx;
+  int ix;
+  AUDIOBEAM_GET_FLOAT_WORD( ix, x );
+  ix &= 0x7fffffff;
+  if ( ix < 0x32000000 ) {
+    if ( ( ( int )x ) == 0 ) return audiobeam_one;
+  }
+  z  = x * x;
+  r  = z * ( audiobeam_C1 + z * ( audiobeam_C2 + z * ( audiobeam_C3 + z *
+                                  ( audiobeam_C4 + z *
+                                    ( audiobeam_C5 + z * audiobeam_C6 ) ) ) ) );
+  if ( ix < 0x3e99999a )
+    return audiobeam_one - ( ( float )0.5f * z - ( z * r - x * y ) );
+  else {
+    if ( ix > 0x3f480000 )
+      qx = ( float )0.28125f;
+    else
+      AUDIOBEAM_SET_FLOAT_WORD( qx, ix - 0x01000000 );
+    hz = ( float )0.5f * z - qx;
+    a  = audiobeam_one - qx;
+    return a - ( hz - ( z * r - x * y ) );
+  }
+}
+
+
+float audiobeam___kernel_sinf( float x, float y, int iy )
+{
+  float z, r, v;
+  int ix;
+  AUDIOBEAM_GET_FLOAT_WORD( ix, x );
+  ix &= 0x7fffffff;