Instruction
VSUBPH
Subtracts packed half-precision values.
- Extensions
- Forms
- 6 (6 EVEX)
- Data
- Intel XED v2026.08.23
VSUBPH subtracts the halves of its last operand from those of the one
before, 32 per zmm register: vsubph zmm0, zmm1, zmm2 computes zmm1 -
zmm2. It is VADDPH with the second operand’s sign
flipped, and that page covers the format, the CPUID bits and the
encoding; VSUBPH has opcode 5Ch in EVEX map 5, that of SUBPS. It rounds
by MXCSR.RC, or by {er} on the 512-bit register form, and reports the
same five exceptions as VSUBPS, whose page gives the
sign of a zero difference: +0, except -0 when rounding toward negative
infinity.
No flushing to zero
Intel’s optimization manual says FP16 arithmetic ignores MXCSR.FTZ and DAZ (Section 18.5.4), so a difference too small for a normal half becomes a denormal, never a zero. Two different halves therefore never subtract to 0, which singles under FTZ can: near the bottom of the half range, at 2^-14, a step of one unit in the last place is only 2^-24, so differences that small arise easily.
In our builds
At -O2 with -march=x86-64-v4 -mavx512fp16, GCC 14 and Clang 23 compile
c[i] = a[i] - b[i] over global _Float16 arrays to VSUBPH with b[i]
from memory, GCC on zmm and Clang on ymm registers. Through pointers and
with an unknown length, GCC keeps one VSUBSH per
element at -O2 and vectorizes at -O3; Clang vectorizes behind a check
that the arrays do not overlap. For 1 - a[i] the 1 is held in a
register as first source, subtracting a constant 1 becomes VADDPH of -1,
and -a[i] flips the sign bits with an exclusive or against 8000h
(VPXORD in GCC, VPXOR in Clang), as for singles. Without AVX512-FP16, as
with -march=znver5, both widen with VCVTPH2PS,
subtract singles and narrow with VCVTPS2PH, GCC one
element at a time and Clang sixteen.
AMD’s Zen 5 guide does not list AVX512-FP16, and this site’s recorded cases include none of VSUBPH.
Forms
| Instruction | Encoding | Requires |
|---|---|---|
VSUBPH xmm {k}{z}, xmm, xmm
VSUBPH_XMMf16_MASKmskw_XMMf16_XMMf16_AVX512 | EVEX.128.NP.MAP5.W0 5C /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VSUBPH xmm {k}{z}, xmm, m128/m16bcst
VSUBPH_XMMf16_MASKmskw_XMMf16_MEMf16_AVX512 | EVEX.128.NP.MAP5.W0 5C /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VSUBPH ymm {k}{z}, ymm, ymm
VSUBPH_YMMf16_MASKmskw_YMMf16_YMMf16_AVX512 | EVEX.256.NP.MAP5.W0 5C /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VSUBPH ymm {k}{z}, ymm, m256/m16bcst
VSUBPH_YMMf16_MASKmskw_YMMf16_MEMf16_AVX512 | EVEX.256.NP.MAP5.W0 5C /r | AVX10.1 or AVX512_FP16 + AVX512VL |
VSUBPH zmm {k}{z}, zmm, zmm, {er}
VSUBPH_ZMMf16_MASKmskw_ZMMf16_ZMMf16_AVX512 | EVEX.512.NP.MAP5.W0 5C /r | AVX10.1 or AVX512_FP16 |
VSUBPH zmm {k}{z}, zmm, m512/m16bcst
VSUBPH_ZMMf16_MASKmskw_ZMMf16_MEMf16_AVX512 | EVEX.512.NP.MAP5.W0 5C /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.
Result half i comes from half i of xmm2 and xmm3. Result halves 7 to 0: a7-b7, unchanged (bit 6 of k1 is 0), a5-b5, a4-b4, a3-b3, a2-b2, a1-b1, a0-b0.
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 = zmm2.f16[i]
b = src3.f16[i]
zmm1.f16[i] = sub(a, b)
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 |
|---|---|---|
VSUBPH xmm {k}{z}, xmm, xmm/m128/m16bcst VSUBPH ymm {k}{z}, ymm, ymm/m256/m16bcst VSUBPH zmm {k}{z}, zmm, zmm/m512/m16bcst, {er} | ||
xmm, ymm, zmm | explicit | 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 VSUBPH, 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_sub_ph
__m128h _mm_sub_ph(__m128h __A, __m128h __B) |
|
_mm_mask_sub_ph
__m128h _mm_mask_sub_ph(__m128h __W, __mmask8 __U, __m128h __A, __m128h __B) |
|
_mm_maskz_sub_ph
__m128h _mm_maskz_sub_ph(__mmask8 __U, __m128h __A, __m128h __B) |
|
_mm256_sub_ph
__m256h _mm256_sub_ph(__m256h __A, __m256h __B) |
|
_mm256_mask_sub_ph
__m256h _mm256_mask_sub_ph(__m256h __W, __mmask16 __U, __m256h __A, __m256h __B) |
|
_mm256_maskz_sub_ph
__m256h _mm256_maskz_sub_ph(__mmask16 __U, __m256h __A, __m256h __B) |
|
_mm512_sub_ph
__m512h _mm512_sub_ph(__m512h __A, __m512h __B) |
|
_mm512_sub_round_ph
__m512h _mm512_sub_round_ph(__m512h A, __m512h B, const int R) |
|
_mm512_mask_sub_ph
__m512h _mm512_mask_sub_ph(__m512h __W, __mmask32 __U, __m512h __A, __m512h __B) |
|
_mm512_mask_sub_round_ph
__m512h _mm512_mask_sub_round_ph(__m512h W, __mmask32 U, __m512h A, __m512h B, const int R) |
|
_mm512_maskz_sub_ph
__m512h _mm512_maskz_sub_ph(__mmask32 __U, __m512h __A, __m512h __B) |
|
_mm512_maskz_sub_round_ph
__m512h _mm512_maskz_sub_round_ph(__mmask32 U, __m512h A, __m512h B, 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
- VRSQRT14SD
- VRSQRT14SS
- VRSQRTPH
- VRSQRTSH
- VSCALEFPD
- VSCALEFPH
- VSCALEFPS
- VSCALEFSD
- VSCALEFSH
- VSCALEFSS
- VSCATTERDPD
- VSCATTERDPS
- VSCATTERQPD
- VSCATTERQPS
- VSHUFF32X4
- VSHUFF64X2
- VSHUFI32X4
- VSHUFI64X2
- VSHUFPD
- VSHUFPS
- VSQRTPD
- VSQRTPH
- VSQRTPS
- VSQRTSD
- VSQRTSH
- VSQRTSS
- VSUBPD
- VSUBPS
- VSUBSD
- VSUBSH
- VSUBSS
- VUCOMISD
- VUCOMISH
- VUCOMISS
- VUNPCKHPD
- VUNPCKHPS
- VUNPCKLPD
- VUNPCKLPS
- VXORPD
- VXORPS
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-768 (VSUBPH)
- 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.