/* { dg-do run } */
/* { dg-options "-O2 -mavx512f" } */
/* { dg-require-effective-target avx512f } */

#define AVX512F

#include "avx512f-helper.h"

#define SIZE (AVX512F_LEN / 64)
#include "avx512f-mask-type.h"
#include "string.h"

static void
CALC (double *e, UNION_TYPE (AVX512F_LEN, d) s1,
      UNION_TYPE (AVX512F_LEN, d) s2, int imm)
{
  int i, offset, selector;
  double *source;
  for (i = 0; i < SIZE / 2; i++)
    {

#if AVX512F_LEN == 512
      selector = (imm >> i * 2) & 0x3;
#else
      selector = (imm >> i) & 0x1;
#endif

      offset = i * 2;
      source = i * 2 * 64 < AVX512F_LEN / 2 ? s1.a : s2.a;
      memcpy (e + offset, source + selector * 2, 16);
    }
}

void
TEST (void)
{
  UNION_TYPE (AVX512F_LEN, d) u1, u2, u3, s1, s2;
  MASK_TYPE mask = MASK_VALUE;
  double e[SIZE];
  int i;
#if AVX512F_LEN == 512
  int imm = 203;
#else
  int imm = 1;
#endif

  for (i = 0; i < SIZE; i++)
    {
      s1.a[i] = 1.2 / (i + 0.378);
      s2.a[i] = 91.02 / (i + 4.3578);
      u1.a[i] = DEFAULT_VALUE;
      u2.a[i] = DEFAULT_VALUE;
      u3.a[i] = DEFAULT_VALUE;
    }

  u1.x = INTRINSIC (_shuffle_f64x2) (s1.x, s2.x, imm);
  u2.x = INTRINSIC (_mask_shuffle_f64x2) (u2.x, mask, s1.x, s2.x, imm);
  u3.x = INTRINSIC (_maskz_shuffle_f64x2) (mask, s1.x, s2.x, imm);

  CALC (e, s1, s2, imm);

  if (UNION_CHECK (AVX512F_LEN, d) (u1, e))
    abort ();

  MASK_MERGE (d) (e, mask, SIZE);
  if (UNION_CHECK (AVX512F_LEN, d) (u2, e))
    abort ();

  MASK_ZERO (d) (e, mask, SIZE);
  if (UNION_CHECK (AVX512F_LEN, d) (u3, e))
    abort ();
}