Function xsum_large_add_sqnorm

xs/sources/xsum.c:1645–1801

/* ADD SQUARED NORM OF VECTOR OF FLOATING-POINT NUMBERS TO LARGE ACCUMULATOR. */

void xsum_large_add_sqnorm (xsum_large_accumulator *restrict lacc,
                            const xsum_flt *restrict vec,
                            xsum_length n)
{
  if (xsum_debug) c_printf("\nLARGE ADD_SQNORM OF %ld VALUES\n",(long)n);

# if OPT_LARGE_SQNORM
  {
    xsum_lcount count;
    xsum_expint ix;
    xsum_uint uintv;
    double fltv;

    while (n > 3)
    {
      fltv = *vec * *vec;
      COPY64 (uintv, fltv);
      vec += 1;

      ix = uintv >> XSUM_MANTISSA_BITS;

      count = lacc->count[ix] - 1;

      if (count < 0)
      { xsum_large_add_value_inf_nan (lacc, ix, uintv);
      }
      else
      { lacc->count[ix] = count;
        lacc->chunk[ix] += uintv;
      }

      fltv = *vec * *vec;
      COPY64 (uintv, fltv);
      vec += 1;

      ix = uintv >> XSUM_MANTISSA_BITS;

      count = lacc->count[ix] - 1;

      if (count < 0)
      { xsum_large_add_value_inf_nan (lacc, ix, uintv);
      }
      else
      { lacc->count[ix] = count;
        lacc->chunk[ix] += uintv;
      }

      fltv = *vec * *vec;
      COPY64 (uintv, fltv);
      vec += 1;

      ix = uintv >> XSUM_MANTISSA_BITS;

      count = lacc->count[ix] - 1;

      if (count < 0)
      { xsum_large_add_value_inf_nan (lacc, ix, uintv);
      }