Instruction
VFMADD213PH
Fused multiply-add of packed half-precision values: second operand times first plus third, into the first.
- Extensions
- Forms
- 6 (6 EVEX)
- Data
- Intel XED v2026.08.23
VFMADD213PH computes op2 × op1 + op3 for each half-precision element and
rounds once, writing over op1: the destination is a factor and the
addend is the operand that may come from memory or be broadcast from a
single half. VFMADD132PH covers the FP16 family,
what a merge mask keeps in each order and the error of the complex
forms; VFMADD213PS is the single-precision form,
and VADDPH explains the format and the CPUID bits. Like
the other FP16 arithmetic, it rounds by MXCSR.RC or by {er} on the
512-bit register form, and ignores DAZ and FTZ (Intel’s optimization
manual, Section 18.5.4).
A constant addend
The 213 order fits a polynomial evaluated by Horner’s rule, p = p × x + c, where each result is multiplied again and the coefficient comes from
memory. At -O2 with -march=x86-64-v4 -mavx512fp16, _mm512_fmadd_ph(x, y, _mm512_set1_ph(1.0f16)) compiles in GCC 14 and Clang 23 to one
VFMADD213PH whose addend is the half 1.0 broadcast from memory
({1to32}), with no register spent on the constant.
Clang’s choice
Clang 23 uses the 213 order for most contracted loops over _Float16
arrays: c[i] = a[i] * b[i] + d[i] and c[i] += a[i] * b[i] both
become VFMADD213PH on ymm registers with both factors in registers and
the addend read from memory, and the quadratic (x*3+2)*x+1 takes two
VFMADD213PH per vector with the coefficients broadcast once into
registers before the loop. GCC 14 compiles the same loops and
_mm512_fmadd_ph on three registers to VFMADD132PH. Both orders give
the same result, so the choice only shows where the memory operand and
the destination can go.
AMD’s Zen 5 guide does not list AVX512-FP16, and this site’s recorded cases include none of VFMADD213PH.
Forms
| Instruction | Encoding | Requires |
|---|---|---|
VFMADD213PH xmm {k}{z}, xmm, xmm
VFMADD213PH_XMMf16_MASKmskw_XMMf16_XMMf16_AVX512 | EVEX.128.66.MAP6.W0 A8 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADD213PH xmm {k}{z}, xmm, m128/m16bcst
VFMADD213PH_XMMf16_MASKmskw_XMMf16_MEMf16_AVX512 | EVEX.128.66.MAP6.W0 A8 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADD213PH ymm {k}{z}, ymm, ymm
VFMADD213PH_YMMf16_MASKmskw_YMMf16_YMMf16_AVX512 | EVEX.256.66.MAP6.W0 A8 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADD213PH ymm {k}{z}, ymm, m256/m16bcst
VFMADD213PH_YMMf16_MASKmskw_YMMf16_MEMf16_AVX512 | EVEX.256.66.MAP6.W0 A8 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADD213PH zmm {k}{z}, zmm, zmm, {er}
VFMADD213PH_ZMMf16_MASKmskw_ZMMf16_ZMMf16_AVX512 | EVEX.512.66.MAP6.W0 A8 /r | AVX10.1 or AVX512_FP16 |
VFMADD213PH zmm {k}{z}, zmm, m512/m16bcst
VFMADD213PH_ZMMf16_MASKmskw_ZMMf16_MEMf16_AVX512 | EVEX.512.66.MAP6.W0 A8 /r | AVX10.1 or AVX512_FP16 |
Latency and throughput
| Model | xmm {k}{z}, xmm, xmmno maskymm {k}{z}, ymm, ymmno mask | xmm {k}{z}, xmm, xmmwith a mask, zeroingymm {k}{z}, ymm, ymmwith a mask, zeroing | xmm {k}{z}, xmm, m128/m16bcstymm {k}{z}, ymm, m256/m16bcst | zmm {k}{z}, zmm, zmm, {er}no mask, with {er} | zmm {k}{z}, zmm, zmm, {er}with a mask, zeroing | zmm {k}{z}, zmm, m512/m16bcst |
|---|---|---|---|---|---|---|
| Sapphire Rapids | 5 / 0.5p01 | 6 / 0.5p01 | 12 / 0.5p01 | 5 / 1p0 | 6 / 1p0 | 12 / 1p0 |
Data flow
Where the elements of the sources go, element 0 on the right; results that depend on a control value are drawn for an example. How to read.
Every half works alike: 8, 16 and 32 halves in the xmm, ymm and zmm forms. The forms write under a mask: where its bit is 0, the element keeps its old value, or zeroes it with {z}. With m16bcst, one half from memory stands for every element of the last source. The sources are drawn above the result, each element in its own column.
Result half i comes from half i of xmm1, xmm2 and xmm3. Result halves 7 to 0: r7, unchanged (bit 6 of k1 is 0), r5, r4, r3, r2, r1, r0.
Operation
What the forms do, in pseudocode that our tests run against results recorded on real processors. How to read.
for i in range(VL // 16):
if k1[i]:
a = zmm1.f16[i]
b = zmm2.f16[i]
c = src3.f16[i]
zmm1.f16[i] = fmadd(b, a, c)
elif zeroing:
zmm1.f16[i] = 0
# else: zmm1.f16[i] keeps its value
# m16bcst: src3.f16[i] is the m16 for every i
# {er}: the form's rounding, not MXCSR.RC, and no flags
# zmm1 is zeroed from bit VL upNot run on a processor: none of these processors has AVX512-FP16. The tests check its arithmetic against MPFR.
Operands
| Operand | Kind | Action |
|---|---|---|
VFMADD213PH xmm {k}{z}, xmm, xmm/m128/m16bcst VFMADD213PH ymm {k}{z}, ymm, ymm/m256/m16bcst VFMADD213PH zmm {k}{z}, zmm, zmm/m512/m16bcst, {er} | ||
xmm, ymm, zmm | explicit | read, write (masked) |
{k} | explicit | read (write mask) |
xmm, ymm, zmm | explicit | read |
xmm/m128/m16bcst, ymm/m256/m16bcst, zmm/m512/m16bcst | explicit | read |
Flags
None: XED records no flag that the forms test or change.
Exceptions
The exception class is XED's name for the set of exception conditions that a form follows. How to read.
- All forms
AVX512-E2
Attributes
What else XED records about the forms. How to read.
- All forms
- write masking; follows MXCSR
- Memory forms
- embedded broadcast (tuple Full); memory fault suppression
Intrinsics
The C intrinsics of Clang 23.1.2 that compile to VFMADD213PH, each with the form it compiles to and the Clang options it needs. The forms come from compiling every intrinsic and decoding the result with Intel XED; all intrinsics A–Z explains the method.
| Intrinsic | Compiles to |
|---|---|
_mm_fmadd_ph
__m128h _mm_fmadd_ph(__m128h __A, __m128h __B, __m128h __C) |
|
_mm_mask_fmadd_ph
__m128h _mm_mask_fmadd_ph(__m128h __A, __mmask8 __U, __m128h __B, __m128h __C) |
|
_mm_maskz_fmadd_ph
__m128h _mm_maskz_fmadd_ph(__mmask8 __U, __m128h __A, __m128h __B, __m128h __C) |
|
_mm256_fmadd_ph
__m256h _mm256_fmadd_ph(__m256h __A, __m256h __B, __m256h __C) |
|
_mm256_mask_fmadd_ph
__m256h _mm256_mask_fmadd_ph(__m256h __A, __mmask16 __U, __m256h __B, __m256h __C) |
|
_mm256_maskz_fmadd_ph
__m256h _mm256_maskz_fmadd_ph(__mmask16 __U, __m256h __A, __m256h __B, __m256h __C) |
|
_mm512_fmadd_ph
__m512h _mm512_fmadd_ph(__m512h __A, __m512h __B, __m512h __C) |
|
_mm512_fmadd_round_ph
__m512h _mm512_fmadd_round_ph(__m512h A, __m512h B, __m512h C, const int R) |
|
_mm512_mask_fmadd_ph
__m512h _mm512_mask_fmadd_ph(__m512h __A, __mmask32 __U, __m512h __B, __m512h __C) |
|
_mm512_maskz_fmadd_ph
__m512h _mm512_maskz_fmadd_ph(__mmask32 __U, __m512h __A, __m512h __B, __m512h __C) |
|
_mm512_maskz_fmadd_round_ph
__m512h _mm512_maskz_fmadd_round_ph(__mmask32 U, __m512h A, __m512h B, __m512h C, const int R) |
|
CPUID and processors
Each form belongs to an XED ISA set. A form can be used when the processor reports every CPUID bit of one of its ISA set's alternatives.
| ISA set | CPUID | Processors |
|---|---|---|
AVX512_FP16_128 | AVX10.1 or AVX512_FP16 + AVX512VL | Nova Lake, Sapphire Rapids, Emerald Rapids, Granite Rapids, Diamond Rapids, Zen 6 |
AVX512_FP16_256 | Same as AVX512_FP16_128 | Same as AVX512_FP16_128 |
AVX512_FP16_512 | AVX10.1 or AVX512_FP16 | Same as AVX512_FP16_128 |
More AVX10.1 instructions
- VEXTRACTI32X4
- VEXTRACTI32X8
- VEXTRACTI64X2
- VEXTRACTI64X4
- VEXTRACTPS
- VFCMADDCPH
- VFCMADDCSH
- VFCMULCPH
- VFCMULCSH
- VFIXUPIMMPD
- VFIXUPIMMPS
- VFIXUPIMMSD
- VFIXUPIMMSS
- VFMADD132PD
- VFMADD132PH
- VFMADD132PS
- VFMADD132SD
- VFMADD132SH
- VFMADD132SS
- VFMADD213PD
- VFMADD213PS
- VFMADD213SD
- VFMADD213SH
- VFMADD213SS
- VFMADD231PD
- VFMADD231PH
- VFMADD231PS
- VFMADD231SD
- VFMADD231SH
- VFMADD231SS
- VFMADDCPH
- VFMADDCSH
- VFMADDSUB132PD
- VFMADDSUB132PH
- VFMADDSUB132PS
- VFMADDSUB213PD
- VFMADDSUB213PH
- VFMADDSUB213PS
- VFMADDSUB231PD
- VFMADDSUB231PH
Sources
- Intel XED v2026.08.23 (commit
0bcb6237345c): forms, encodings, operands, flags, exception classes, attributes, ISA sets, CPUID bits and chip model. - The CPU data of x86-64.net: Zen 3 to Zen 6, from GCC’s
-marchfeature sets and AMD’s Processor Programming References, which the page of each processor cites. - Clang's x86 intrinsic headers, LLVM 23.1.2 (commit
85ac56026243): the intrinsics, their signatures and doc comments, and the Clang of that release, which compiled each of them. - LLVM 23.1.2's x86 scheduling models (commit
85ac56026243): the latency, reciprocal throughput and ports of each form, as llvm-mca reports them. - Intel Advanced Vector Extensions 10.2 Architecture Specification, revision 7.0 (361050-007, April 2026)
- Intel 64 and IA-32 Architectures Software Developer's Manual, Volume 2 (325383-093, September 2026): Volume 2C, page 5-214 (VF[,N]MADD[132,213,231]PH)
- Intel 64 and IA-32 Architectures Optimization Reference Manual, Volume 1 (248966-051, October 2026)
- Software Optimization Guide for the AMD Zen5 Microarchitecture (58455, rev. 1.00, August 2024), Wayback Machine copy of AMD's ZIP
The tables are derived from Intel XED, Copyright Intel Corporation, licensed under the Apache License 2.0; x86-64.net converted and reformatted them. The intrinsics and their signatures come from Clang and the timings from LLVM's scheduling models, both part of the LLVM Project, licensed under the Apache License 2.0 with LLVM Exceptions. The text is our own.