/* { dg-do run } */
/* { dg-options "-mavx2 -O2" } */
/* { dg-require-effective-target avx2 } */

#include "avx2-check.h"

static void
compute_pmulhrsw256 (short *s1, short *s2, short *r)
{
  int i;
  int t0;

  for (i = 0; i < 16; i++)
    {
      t0 = (((int) s1[i] * (int) s2[i]) >> 14) + 1;
      r[i] = (short) (t0 >> 1);
    }
}

static void
avx2_test (void)
{
  union256i_w s1, s2, res;
  short res_ref[16];
  int i, j, sign = 1;
  int fail = 0;

  for (i = 0; i < 10; i++)
    {
      for (j = 0; j < 16; j++)
	{
	  s1.a[j] = i * j * sign;
	  s2.a[j] = (j + 20) * sign;
	  sign = -sign;
	}

      res.x = _mm256_mulhrs_epi16 (s1.x, s2.x);

      compute_pmulhrsw256 (s1.a, s2.a, res_ref);

      fail += check_union256i_w (res, res_ref);
    }

  if (fail != 0)
    abort ();
}