Dr. Dobb's is part of the Informa Tech Division of Informa PLC

This site is operated by a business or businesses owned by Informa PLC and all copyright resides with them. Informa PLC's registered office is 5 Howick Place, London SW1P 1WG. Registered in England and Wales. Number 8860726.


Channels ▼
RSS

Floating-point Summation


September 1996/Floating-point Summation/Listing 8

Listing 8: Tests performance of summation methods


int main(void)
{
  float flt_arr[ARR_SIZE];
  double dbl_arr[ARR_SIZE];
  double val;
  double x;
  int j, k;
  long i;

  /* seed random number generator from time */
  srand((unsigned int)time(0));

  /* initialize arrays for test 1: uniform random 0.0 - 1.0 */

  printf("\n%s     Simple      Extended"
         "    Sorted      Pairwise    Kahan\n", "Unif 0-1   ");
  for (j = 0; j < OUTER_ITER; j++)
  {
    for (i = 0; i < ARR_SIZE; i++)
      flt_arr[i] = dbl_arr[i] = 
          (  rand()*(double)RAND_MAX*RAND_MAX*(double)RAND_MAX
           + rand()*(double)RAND_MAX*RAND_MAX
           + rand()*(double)RAND_MAX
           + rand())
          /
          (  RAND_MAX*(double)RAND_MAX*RAND_MAX*(double)RAND_MAX
           + RAND_MAX*(double)RAND_MAX*RAND_MAX
           + RAND_MAX*(double)RAND_MAX
           + RAND_MAX);

    loop(flt_arr, dbl_arr, "Unif 0-1   ");
  }
  print_and_reset_minmaxerr();

  /* reinitialize for test 2: uniform random -1.0 - 1.0 */
  printf("\n%s     Simple      Extended"
         "    Sorted      Pairwise    Kahan\n", "Unif -1-1  ");
  for (j = 0; j < OUTER_ITER; j++)
  {
    for (i = 0; i < ARR_SIZE; i++)
      flt_arr[i] = dbl_arr[i] =
          (  rand()*(double)RAND_MAX*RAND_MAX*(double)RAND_MAX
           + rand()*(double)RAND_MAX*RAND_MAX
           + rand()*(double)RAND_MAX
           + rand()) * 2.0
          /
          (  RAND_MAX*(double)RAND_MAX*RAND_MAX*(double)RAND_MAX
           + RAND_MAX*(double)RAND_MAX*RAND_MAX
           + RAND_MAX*(double)RAND_MAX
           + RAND_MAX)
          - 1.0;

    loop(flt_arr, dbl_arr, "Unif -1-1  ");
  }
  print_and_reset_minmaxerr();

  /* reinitialize for test 3: more or less Gaussian */
  /*
   * Simple but not very accurate way to make a gaussian: just
   * sum lots (here 12) uniform random variables.
   */
  printf("\n%s     Simple      Extended"
         "    Sorted      Pairwise    Kahan\n", "Gaussian   ");
  for (j = 0; j < OUTER_ITER; j++)
  {
    for (i = 0; i < ARR_SIZE; i++)
    {
      val = 0.0;
      for (k = 0; k < 12; k++)
        val += rand();
      flt_arr[i] = dbl_arr[i] = val / (double)RAND_MAX - 6.0;
    }

    loop(flt_arr, dbl_arr, "Gaussian   ");
  }
  print_and_reset_minmaxerr();
  /* reinitialize for test 4:
   * numerical integration to calculate pi */
  /*
   * The area of a circle is pi * radius^2.  When the radius is
   * 1.0 the area is pi.  So we approximate pi as 4 times the
   * area of one quarter of the unit circle, and get the area
   * of that quadrant using numerical integration.
   */
  printf("\n%s     Simple      Extended"
         "    Sorted      Pairwise    Kahan\n", "Pi         ");
  for (j = 0; j < OUTER_ITER; j++)
  {
    for (i = 0; i < ARR_SIZE; i++)
    {
      x = (1 + 2 * i) / (double)(2 * ARR_SIZE);
      flt_arr[i] = dbl_arr[i] =
          sqrt(1.0 - x * x) * (4.0 / ARR_SIZE);

    loop(flt_arr, dbl_arr, "Pi         ");
  }
  print_and_reset_minmaxerr();

  /*
   * Omitted from listing: code to print times to nearest tenth
   * of a second
   * ...
   */

  return EXIT_SUCCESS;
}
/* End of File */

Related Reading


More Insights






Currently we allow the following HTML tags in comments:

Single tags

These tags can be used alone and don't need an ending tag.

<br> Defines a single line break

<hr> Defines a horizontal line

Matching tags

These require an ending tag - e.g. <i>italic text</i>

<a> Defines an anchor

<b> Defines bold text

<big> Defines big text

<blockquote> Defines a long quotation

<caption> Defines a table caption

<cite> Defines a citation

<code> Defines computer code text

<em> Defines emphasized text

<fieldset> Defines a border around elements in a form

<h1> This is heading 1

<h2> This is heading 2

<h3> This is heading 3

<h4> This is heading 4

<h5> This is heading 5

<h6> This is heading 6

<i> Defines italic text

<p> Defines a paragraph

<pre> Defines preformatted text

<q> Defines a short quotation

<samp> Defines sample computer code text

<small> Defines small text

<span> Defines a section in a document

<s> Defines strikethrough text

<strike> Defines strikethrough text

<strong> Defines strong text

<sub> Defines subscripted text

<sup> Defines superscripted text

<u> Defines underlined text

Dr. Dobb's encourages readers to engage in spirited, healthy debate, including taking us to task. However, Dr. Dobb's moderates all comments posted to our site, and reserves the right to modify or remove any content that it determines to be derogatory, offensive, inflammatory, vulgar, irrelevant/off-topic, racist or obvious marketing or spam. Dr. Dobb's further reserves the right to disable the profile of any commenter participating in said activities.

 
Disqus Tips To upload an avatar photo, first complete your Disqus profile. | View the list of supported HTML tags you can use to style comments. | Please read our commenting policy.