Instruction
VFMADDSUB213PH
Fused multiply of packed half-precision values, second times first, minus the third in even elements and plus it in odd ones, into the first.
- Extensions
- Forms
- 6 (6 EVEX)
- Data
- Intel XED v2026.08.23
VFMADDSUB213PH computes op2 × op1 ∓ op3 for each half-precision element
and rounds once: the even elements subtract op3, the operand that may
come from memory or be broadcast from one half, and the odd ones add it,
the result going over op1. VFMADDSUB132PH
covers the alternating FP16 forms, VFMADD213PH the
operand roles of the 213 order, and
VFMADDSUB213PS is the single-precision form.
The opcode is A6h with the 66 prefix in EVEX map 6; rounding follows
MXCSR.RC, or {er} on the 512-bit register form, and denormals are kept
whatever DAZ and FTZ say, as VADDPH explains.
In our builds
At -O2 with -march=x86-64-v4 -mavx512fp16, Clang 23 compiles a loop
that alternates by hand over _Float16 arrays, r[i] = a[i] * b[i] - c[i] then r[i+1] = a[i+1] * b[i+1] + c[i+1], to VFMADDSUB213PH on ymm
registers with c read from memory; GCC 14 uses the 132 order on zmm
registers. Clang also takes this order for _mm512_fmaddsub_ph(a, b, c)
on three registers, and when all three operands are loaded from arrays
it reads the third in place. Both compilers use VFMADDSUB213PH when only
the third operand is in memory.
A constant third operand shows a difference in the other direction. For
_mm512_fmaddsub_ph(x, y, _mm512_set1_ph(1.0f16)), which subtracts 1
from x × y in the even elements and adds it in the odd ones, Clang reads
the half 1.0 with {1to32} as op3. GCC broadcasts it into a register
first with VPBROADCASTW and uses VFMADDSUB132PH, although for
_mm512_fmsub_ph with the same constant it does fold the broadcast, as
VFMSUB213PH’s page shows.
AMD’s Zen 5 guide does not list AVX512-FP16, and this site’s recorded cases include none of VFMADDSUB213PH.
Forms
| Instruction | Encoding | Requires |
|---|---|---|
VFMADDSUB213PH xmm {k}{z}, xmm, xmm
VFMADDSUB213PH_XMMf16_MASKmskw_XMMf16_XMMf16_AVX512 | EVEX.128.66.MAP6.W0 A6 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADDSUB213PH xmm {k}{z}, xmm, m128/m16bcst
VFMADDSUB213PH_XMMf16_MASKmskw_XMMf16_MEMf16_AVX512 | EVEX.128.66.MAP6.W0 A6 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADDSUB213PH ymm {k}{z}, ymm, ymm
VFMADDSUB213PH_YMMf16_MASKmskw_YMMf16_YMMf16_AVX512 | EVEX.256.66.MAP6.W0 A6 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADDSUB213PH ymm {k}{z}, ymm, m256/m16bcst
VFMADDSUB213PH_YMMf16_MASKmskw_YMMf16_MEMf16_AVX512 | EVEX.256.66.MAP6.W0 A6 /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VFMADDSUB213PH zmm {k}{z}, zmm, zmm, {er}
VFMADDSUB213PH_ZMMf16_MASKmskw_ZMMf16_ZMMf16_AVX512 | EVEX.512.66.MAP6.W0 A6 /r | AVX10.1 or AVX512_FP16 |
VFMADDSUB213PH zmm {k}{z}, zmm, m512/m16bcst
VFMADDSUB213PH_ZMMf16_MASKmskw_ZMMf16_MEMf16_AVX512 | EVEX.512.66.MAP6.W0 A6 /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]
if i % 2 == 0:
zmm1.f16[i] = fmsub(b, a, c)
else:
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 |
|---|---|---|
VFMADDSUB213PH xmm {k}{z}, xmm, xmm/m128/m16bcst VFMADDSUB213PH ymm {k}{z}, ymm, ymm/m256/m16bcst VFMADDSUB213PH 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 VFMADDSUB213PH, 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_fmaddsub_ph
__m128h _mm_fmaddsub_ph(__m128h __A, __m128h __B, __m128h __C) |
|
_mm_mask_fmaddsub_ph
__m128h _mm_mask_fmaddsub_ph(__m128h __A, __mmask8 __U, __m128h __B, __m128h __C) |
|
_mm_maskz_fmaddsub_ph
__m128h _mm_maskz_fmaddsub_ph(__mmask8 __U, __m128h __A, __m128h __B, __m128h __C) |
|
_mm256_fmaddsub_ph
__m256h _mm256_fmaddsub_ph(__m256h __A, __m256h __B, __m256h __C) |
|
_mm256_mask_fmaddsub_ph
__m256h _mm256_mask_fmaddsub_ph(__m256h __A, __mmask16 __U, __m256h __B, __m256h __C) |
|
_mm256_maskz_fmaddsub_ph
__m256h _mm256_maskz_fmaddsub_ph(__mmask16 __U, __m256h __A, __m256h __B, __m256h __C) |
|
_mm512_fmaddsub_ph
__m512h _mm512_fmaddsub_ph(__m512h __A, __m512h __B, __m512h __C) |
|
_mm512_fmaddsub_round_ph
__m512h _mm512_fmaddsub_round_ph(__m512h A, __m512h B, __m512h C, const int R) |
|
_mm512_mask_fmaddsub_ph
__m512h _mm512_mask_fmaddsub_ph(__m512h __A, __mmask32 __U, __m512h __B, __m512h __C) |
|
_mm512_maskz_fmaddsub_ph
__m512h _mm512_maskz_fmaddsub_ph(__mmask32 __U, __m512h __A, __m512h __B, __m512h __C) |
|
_mm512_maskz_fmaddsub_round_ph
__m512h _mm512_maskz_fmaddsub_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
- VFMADD132SH
- VFMADD132SS
- VFMADD213PD
- VFMADD213PH
- VFMADD213PS
- VFMADD213SD
- VFMADD213SH
- VFMADD213SS
- VFMADD231PD
- VFMADD231PH
- VFMADD231PS
- VFMADD231SD
- VFMADD231SH
- VFMADD231SS
- VFMADDCPH
- VFMADDCSH
- VFMADDSUB132PD
- VFMADDSUB132PH
- VFMADDSUB132PS
- VFMADDSUB213PD
- VFMADDSUB213PS
- VFMADDSUB231PD
- VFMADDSUB231PH
- VFMADDSUB231PS
- VFMSUB132PD
- VFMSUB132PH
- VFMSUB132PS
- VFMSUB132SD
- VFMSUB132SH
- VFMSUB132SS
- VFMSUB213PD
- VFMSUB213PH
- VFMSUB213PS
- VFMSUB213SD
- VFMSUB213SH
- VFMSUB213SS
- VFMSUB231PD
- VFMSUB231PH
- VFMSUB231PS
- VFMSUB231SD
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-243 (VFMADDSUB132PH/
VFMADDSUB213PH/ VFMADDSUB231PH) - 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.