forked from AuroraMiddleware/gtk
155a4fac5c
We can't make the -4 versions inline, since we use ifuncs for them, so make vectorized versions. Test included.
123 lines
2.6 KiB
C
123 lines
2.6 KiB
C
/* fp16i.c
|
|
*
|
|
* Copyright 2021 Red Hat, Inc.
|
|
*
|
|
* This library is free software; you can redistribute it and/or
|
|
* modify it under the terms of the GNU Lesser General Public
|
|
* License as published by the Free Software Foundation; either
|
|
* version 2.1 of the License, or (at your option) any later version.
|
|
*
|
|
* This library is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
* Lesser General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU Lesser General Public
|
|
* License along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
*
|
|
* SPDX-License-Identifier: LGPL-2.1-or-later
|
|
*/
|
|
|
|
#include "config.h"
|
|
|
|
#include "fp16private.h"
|
|
|
|
#ifdef HAVE_F16C
|
|
#include <immintrin.h>
|
|
|
|
#if defined(_MSC_VER) && !defined(__clang__)
|
|
#define CAST_M128I_P(a) (__m128i const *) a
|
|
#else
|
|
#define CAST_M128I_P(a) (__m128i_u const *) a
|
|
#endif
|
|
void
|
|
float_to_half4_f16c (const float f[4],
|
|
guint16 h[4])
|
|
{
|
|
__m128 s = _mm_loadu_ps (f);
|
|
__m128i i = _mm_cvtps_ph (s, 0);
|
|
_mm_storel_epi64 ((__m128i*)h, i);
|
|
}
|
|
|
|
void
|
|
half_to_float4_f16c (const guint16 h[4],
|
|
float f[4])
|
|
{
|
|
__m128i i = _mm_loadl_epi64 (CAST_M128I_P (h));
|
|
__m128 s = _mm_cvtph_ps (i);
|
|
|
|
_mm_store_ps (f, s);
|
|
}
|
|
|
|
#define ALIGNED(p, n) (GPOINTER_TO_UINT(p) % n == 0)
|
|
void
|
|
float_to_half_f16c (const float *f,
|
|
guint16 *h,
|
|
int n)
|
|
{
|
|
__m128 s;
|
|
__m128i i;
|
|
int j;
|
|
const float *ff = f;
|
|
guint16 *hh = h;
|
|
|
|
for (j = 0; j < n; j++)
|
|
{
|
|
if (ALIGNED (ff, 16) && ALIGNED (hh, 16))
|
|
break;
|
|
ff++;
|
|
hh++;
|
|
}
|
|
|
|
float_to_half_c (f, h, j);
|
|
|
|
for (; j + 4 < n; j += 4)
|
|
{
|
|
s = _mm_loadu_ps (ff);
|
|
i = _mm_cvtps_ph (s, 0);
|
|
_mm_storel_epi64 ((__m128i*)hh, i);
|
|
ff += 4;
|
|
hh += 4;
|
|
}
|
|
|
|
if (j < n)
|
|
float_to_half_c (ff, hh, n - j);
|
|
}
|
|
|
|
void
|
|
half_to_float_f16c (const guint16 *h,
|
|
float *f,
|
|
int n)
|
|
{
|
|
__m128i i;
|
|
__m128 s;
|
|
int j;
|
|
const guint16 *hh = h;
|
|
float *ff = f;
|
|
|
|
for (j = 0; j < n; j++)
|
|
{
|
|
if (ALIGNED (ff, 16) && ALIGNED (hh, 16))
|
|
break;
|
|
ff++;
|
|
hh++;
|
|
}
|
|
|
|
half_to_float_c (h, f, j);
|
|
|
|
for (; j + 4 < n; j += 4)
|
|
{
|
|
i = _mm_loadl_epi64 (CAST_M128I_P (hh));
|
|
s = _mm_cvtph_ps (i);
|
|
_mm_store_ps (ff, s);
|
|
hh += 4;
|
|
ff += 4;
|
|
}
|
|
|
|
if (j < n)
|
|
half_to_float_c (hh, ff, n - j);
|
|
}
|
|
|
|
#endif /* HAVE_F16C */
|
|
|