/* { dg-do compile } */
/* { dg-options "-O2 -ffast-math -mavx512fp16 -mavx512vl -mprefer-vector-width=512" } */
/* { dg-final { scan-assembler-times "vfmaddcph\[ \\t\]" 1 } } */
/* { dg-final { scan-assembler-not "vfmadd\[123]*ph\[ \\t\]"} } */
/* { dg-final { scan-assembler-not "vfmadd\[123]*sh\[ \\t\]"} } */
/* { dg-final { scan-assembler-times "vfcmaddcph\[ \\t\]" 1 } } */
/* { dg-final { scan-assembler-times "vfmulcph\[ \\t\]" 1 } } */
/* { dg-final { scan-assembler-times "vfcmulcph\[ \\t\]" 1 } } */

#include<complex.h>
#define TYPE _Float16
#define N 16

void fma0 (_Complex TYPE *a, _Complex TYPE *b,
           _Complex TYPE * __restrict c)
{
  for (int i = 0; i < N; i++)
    c[i] += a[i] * b[i];
}

void fmaconj (_Complex TYPE a[restrict N], _Complex TYPE b[restrict N],
	      _Complex TYPE c[restrict N])
{
  for (int i = 0; i < N; i++)
    c[i] += a[i] * ~b[i];
}

void fmul (_Complex TYPE a[restrict N], _Complex TYPE b[restrict N],
	   _Complex TYPE c[restrict N])
{
  for (int i = 0; i < N; i++)
    c[i] = a[i] * b[i];
}

void fmulconj (_Complex TYPE a[restrict N], _Complex TYPE b[restrict N],
	       _Complex TYPE c[restrict N])
{
  for (int i = 0; i < N; i++)
    c[i] = a[i] * ~b[i];
}