diff --git a/fast_llm/functional/triton/gspo_loss.py b/fast_llm/functional/triton/gspo_loss.py
new file mode 100644
index 000000000..d1bcfef17
--- /dev/null
+++ b/fast_llm/functional/triton/gspo_loss.py
@@ -0,0 +1,226 @@
+import torch
+
+from fast_llm.core.distributed import ReduceOp, all_reduce
+from fast_llm.functional.triton import tl, tl_arange, tl_constexpr, triton, triton_jit
+from fast_llm.functional.triton.entropy_loss import (
+    parallel_sum_exp_logits,
+    triton_cross_entropy_forward_from_labels_parallel_kernel,
+)
+
+
+@triton_jit()
+def triton_gspo_loss_backward_kernel(
+    logits_ptr,
+    labels_ptr,
+    max_logits_ptr,
+    sum_exp_logits_ptr,
+    probability_ratio_ptr,
+    seg_advantage_ptr,
+    token_weight_ptr,
+    grad_logits_ptr,
+    n_cols: tl_constexpr,
+    logits_stride_0: tl_constexpr,
+    grad_logits_stride_0: tl_constexpr,
+    block_size: tl_constexpr,
+    grad_losses,
+    col_min: tl_constexpr = 0,
+    logits_scale_factor: tl_constexpr = 1.0,
+    epsilon_low: tl_constexpr = 0.2,
+    epsilon_high: tl_constexpr = 0.2,
+    accumulate: tl_constexpr = False,
+):
+    block_idx = tl.program_id(0).to(tl.int64)
+
+    # token_weight = mask_t / N_d, where N_d is the labeled-token count for the doc containing t.
+    # Zero for masked tokens (mask=0) and for tokens with N_d=0 after the kernel's clamp.
+    token_weight = tl.load(token_weight_ptr + block_idx).to(tl.float32)
+    if token_weight == 0.0:
+        if not accumulate:
+            for col_offset in tl.static_range(0, n_cols, block_size):
+                col_offsets = tl_arange(int(col_offset), int(col_offset + block_size))
+                tl.store(
+                    grad_logits_ptr + block_idx * grad_logits_stride_0 + col_offsets, 0, mask=col_offsets < n_cols
+                )
+        return
+
+    label_idx = tl.load(labels_ptr + block_idx) - col_min
+    max_logits = tl.load(max_logits_ptr + block_idx)
+    sum_exp_logits = tl.load(sum_exp_logits_ptr + block_idx)
+    probability_ratio = tl.load(probability_ratio_ptr + block_idx).to(tl.float32)
+    seg_advantage = tl.load(seg_advantage_ptr + block_idx).to(tl.float32)
+
+    # effective_grad = grad_losses * scale * weight * R_s * clip_factor
+    # clip_factor = clamp_min(A_s, 0) * (R_s <= 1+eps_h) + clamp_max(A_s, 0) * (R_s >= 1-eps_l)
+    grad_scale = grad_losses
+    if logits_scale_factor != 1.0:
+        grad_scale *= logits_scale_factor
+    effective_grad = (
+        (
+            tl.maximum(seg_advantage, 0.0) * (probability_ratio <= 1.0 + epsilon_high)
+            + tl.minimum(seg_advantage, 0.0) * (probability_ratio >= 1.0 - epsilon_low)
+        )
+        * probability_ratio
+        * grad_scale
+        * token_weight
+    )
+
+    logits_ptr = logits_ptr + block_idx * logits_stride_0
+
+    # grad_logit_i = effective_grad * (softmax_i - delta_{i, label})
+    col_offset_start: tl.constexpr = (n_cols - 1) // block_size * block_size
+    for col_offset in tl.static_range(col_offset_start, -1, -block_size):
+        col_offsets = tl_arange(col_offset, col_offset + block_size)
+        mask = col_offsets < n_cols
+        logits = tl.load(logits_ptr + col_offsets, mask=mask, other=-float("inf")).to(tl.float32)
+        if logits_scale_factor != 1.0:
+            logits *= logits_scale_factor
+        prob = tl.exp(logits - max_logits) / sum_exp_logits
+        if label_idx < 0 or label_idx >= n_cols:
+            # Target not in this TP shard.
+            grad = effective_grad * prob
+        else:
+            grad = effective_grad * tl.where(col_offsets == label_idx, prob - 1.0, prob)
+        grad_col_ptr = grad_logits_ptr + block_idx * grad_logits_stride_0 + col_offsets
+        if accumulate:
+            grad += tl.load(grad_col_ptr, mask=mask)
+        tl.store(grad_col_ptr, grad, mask=mask)
+
+
+def triton_gspo_loss_forward_backward(
+    logits: torch.Tensor,  # (*batch, vocab)
+    target: torch.Tensor,  # (*batch,)
+    advantages: torch.Tensor,  # (*batch,)
+    old_log_probabilities: torch.Tensor,  # (*batch,)
+    document_index_zero_based: torch.Tensor,  # (*batch,) int — segment ID per token, 0-based
+    num_segments: int,  # buffer size, ≥ document_index.max() + 1
+    divisor: float,
+    num_labels_in_seq: torch.Tensor,  # (*batch,) — per-document labeled-token count broadcast per token
+    grad_logits: torch.Tensor | None = None,
+    grad_output: float | None = None,
+    group: torch.distributed.ProcessGroup | None = None,  # TP vocab group
+    sdp_group: torch.distributed.ProcessGroup | None = None,  # SDP group for cross-rank segment aggregation
+    sp_group: torch.distributed.ProcessGroup | None = None,  # TP group when SP is sharding the sequence
+    epsilon_low: float = 0.2,
+    epsilon_high: float = 0.2,
+    logits_scale_factor: float = 1.0,
+    block_size: int | None = None,
+    num_warps: int | None = None,
+) -> tuple[torch.Tensor, torch.Tensor | None, torch.Tensor | None]:
+    """Triton GSPO loss. Forward fuses softmax + predicted-logit lookup; backward fuses the
+    softmax chain rule with the per-token GSPO gradient factor (R_s * clip * token_weight).
+    Segment aggregation, loss, and the SDP/SP all-reduce live in PyTorch between the two passes.
+
+    See `fused_gspo_loss_forward_backward` in policy_gradient.py for the math derivation;
+    this kernel produces identical outputs.
+    """
+    assert logits.is_contiguous()
+    assert target.is_contiguous()
+    assert advantages.is_contiguous()
+    assert old_log_probabilities.is_contiguous()
+    assert document_index_zero_based.is_contiguous()
+    assert num_labels_in_seq.is_contiguous()
+
+    n_rows = logits.shape[:-1].numel()
+    n_cols = logits.size(-1)
+    if block_size is None:
+        block_size = min(triton.next_power_of_2(n_cols), 32768)
+    if num_warps is None:
+        num_warps = 4 if block_size < 2048 else (8 if block_size < 8192 else 16)
+    col_min = n_cols * group.rank() if group is not None else 0
+
+    # === Forward (Triton): per-token softmax, save max / sum / predicted_logit ===
+    max_logits = torch.empty(n_rows, dtype=torch.float, device=logits.device)
+    sum_exp_logits = torch.empty_like(max_logits)
+    predicted_logits = torch.empty_like(max_logits)
+    triton_cross_entropy_forward_from_labels_parallel_kernel[(n_rows,)](
+        logits,
+        target,
+        max_logits_ptr=max_logits,
+        sum_exp_logits_ptr=sum_exp_logits,
+        predicted_logits_ptr=predicted_logits,
+        col_min=col_min,
+        n_cols=n_cols,
+        logits_stride_0=logits.stride(-2),
+        block_size=block_size,
+        num_warps=num_warps,
+        logits_scale_factor=logits_scale_factor,
+    )
+    if group is not None:
+        # Merge per-shard local max / sum_exp into global values.
+        max_logits, sum_exp_logits = parallel_sum_exp_logits(sum_exp_logits, max_logits, group)
+        all_reduce(predicted_logits, op=ReduceOp.SUM, group=group)
+
+    # === Segment aggregation (PyTorch) ===
+    flat_target = target.reshape(-1)
+    flat_document_index = document_index_zero_based.reshape(-1).long()
+    flat_advantages = advantages.reshape(-1).float()
+    loss_mask = (flat_target >= 0).to(max_logits.dtype)
+
+    new_log_probs = predicted_logits - max_logits - sum_exp_logits.log()
+    log_ratio = (new_log_probs - old_log_probabilities.reshape(-1).float()) * loss_mask
+
+    # Per-token weight: mask / per-document label count. Pre-dividing here means each segment's
+    # contribution to the per-segment sum is already normalized, so SDP/SP all-reduce works
+    # without a separate token-count tensor.
+    flat_num_labels = num_labels_in_seq.reshape(-1).to(new_log_probs.dtype).clamp(min=1)
+    token_weight = loss_mask / flat_num_labels
+
+    mean_log_ratio_per_segment = log_ratio.new_zeros(num_segments).index_add_(
+        0, flat_document_index, log_ratio * token_weight
+    )
+    mean_advantage_per_segment = log_ratio.new_zeros(num_segments).index_add_(
+        0, flat_document_index, flat_advantages * token_weight
+    )
+    for reduce_group in (sdp_group, sp_group):
+        if reduce_group is not None:
+            torch.distributed.all_reduce(
+                mean_log_ratio_per_segment, op=torch.distributed.ReduceOp.SUM, group=reduce_group
+            )
+            torch.distributed.all_reduce(
+                mean_advantage_per_segment, op=torch.distributed.ReduceOp.SUM, group=reduce_group
+            )
+
+    segment_ratio = mean_log_ratio_per_segment.exp()
+    segment_advantage = mean_advantage_per_segment
+
+    probability_ratio = segment_ratio[flat_document_index].contiguous()
+    seg_advantage = segment_advantage[flat_document_index].contiguous()
+    token_weight = token_weight.contiguous()
+
+    losses = -torch.min(
+        probability_ratio * seg_advantage,
+        torch.clamp(probability_ratio, 1 - epsilon_low, 1 + epsilon_high) * seg_advantage,
+    )
+    loss = (losses * token_weight).sum() / divisor
+
+    new_logprobs_mean = (new_log_probs * loss_mask / flat_num_labels).sum()
+
+    if grad_output is None:
+        return loss, grad_logits, new_logprobs_mean
+
+    # === Backward (Triton) ===
+    accumulate = grad_logits is not None
+    grad_logits = torch.empty_like(logits) if grad_logits is None else grad_logits
+    triton_gspo_loss_backward_kernel[(n_rows,)](
+        logits,
+        target,
+        max_logits_ptr=max_logits,
+        sum_exp_logits_ptr=sum_exp_logits,
+        probability_ratio_ptr=probability_ratio,
+        seg_advantage_ptr=seg_advantage,
+        token_weight_ptr=token_weight,
+        grad_logits_ptr=grad_logits,
+        n_cols=n_cols,
+        logits_stride_0=logits.stride(-2),
+        grad_logits_stride_0=grad_logits.stride(-2),
+        block_size=block_size,
+        grad_losses=grad_output / divisor,
+        col_min=col_min,
+        logits_scale_factor=logits_scale_factor,
+        epsilon_low=epsilon_low,
+        epsilon_high=epsilon_high,
+        accumulate=accumulate,
+        num_warps=num_warps,
+    )
+
+    return loss, grad_logits, new_logprobs_mean
diff --git a/fast_llm/layers/language_model/loss/config.py b/fast_llm/layers/language_model/loss/config.py
index 9cdb8c962..9a220aacf 100644
--- a/fast_llm/layers/language_model/loss/config.py
+++ b/fast_llm/layers/language_model/loss/config.py
@@ -260,6 +260,12 @@ class LanguageModelGSPOLossConfig(LanguageModelPolicyGradientLossConfig):
 
     _abstract: typing.ClassVar[bool] = False
 
+    use_triton: bool | None = Field(
+        default=None,
+        desc="Enable triton implementation. Default: use if available.",
+        hint=FieldHint.expert,
+    )
+
     @property
     def loss_class(self) -> "type[LanguageModelGSPOLoss]":
         from fast_llm.layers.language_model.loss.policy_gradient import LanguageModelGSPOLoss
diff --git a/fast_llm/layers/language_model/loss/policy_gradient.py b/fast_llm/layers/language_model/loss/policy_gradient.py
index 38e68d58f..62db342fe 100644
--- a/fast_llm/layers/language_model/loss/policy_gradient.py
+++ b/fast_llm/layers/language_model/loss/policy_gradient.py
@@ -270,7 +270,13 @@ def _forward_backward(
         # SDP/SP ranks within a DP rank, so per-segment buffers are sized consistently for all-reduce.
         num_segments = kwargs[BlockKwargs.num_documents_in_sequence]
 
-        loss, grad, new_logprobs_mean = fused_gspo_loss_forward_backward(
+        if TritonConfig.enabled(logits.device, self._config.use_triton):
+            from fast_llm.functional.triton.gspo_loss import triton_gspo_loss_forward_backward
+
+            fn = triton_gspo_loss_forward_backward
+        else:
+            fn = fused_gspo_loss_forward_backward
+        loss, grad, new_logprobs_mean = fn(
             logits,
             self._get_labels(kwargs, split_index),
             self._prepare_target(kwargs[LanguageModelLossKwargs.advantages], split_index),
diff --git a/tests/layers/test_lm_losses.py b/tests/layers/test_lm_losses.py
index 61df60a53..8026d9739 100644
--- a/tests/layers/test_lm_losses.py
+++ b/tests/layers/test_lm_losses.py
@@ -14,6 +14,7 @@
 from fast_llm.functional.triton import triton_available
 from fast_llm.functional.triton.entropy_loss import triton_entropy_loss_forward_backward
 from fast_llm.functional.triton.grpo_loss import triton_grpo_loss_forward_backward
+from fast_llm.functional.triton.gspo_loss import triton_gspo_loss_forward_backward
 from fast_llm.functional.triton.z_loss import triton_z_loss_forward_backward
 from fast_llm.layers.language_model.loss.dpo import dpo_loss
 from fast_llm.layers.language_model.loss.loss import loss_forward_backward
@@ -463,6 +464,25 @@ def _test_gspo_loss(
     _compare_losses_and_grads(out_fused, out_ref, grad_output is not None, grad_fused, grad_ref, group=group)
     Assert.rms_close_relative(new_logprobs_fused, ref_new_logprobs, 1e-5, 1e-6)
 
+    if not triton_available:
+        return
+    out_triton, grad_triton, new_logprobs_triton = triton_gspo_loss_forward_backward(
+        split_op(logits, group, -1).contiguous(),
+        target,
+        advantages,
+        old_log_probabilities,
+        document_index,
+        num_segments,
+        divisor=num_segments,
+        num_labels_in_seq=num_labels_in_seq,
+        grad_logits=local_previous_grad.clone() if accumulate else None,
+        grad_output=grad_output,
+        group=group,
+        logits_scale_factor=logits_scale_factor,
+    )
+    _compare_losses_and_grads(out_triton, out_ref, grad_output is not None, grad_triton, grad_ref, group=group)
+    Assert.rms_close_relative(new_logprobs_triton, new_logprobs_fused, 1e-5, 1e-6)
+
 
 def _check_grpo_metrics(ref: GRPOMetrics, got: GRPOMetrics, threshold: float) -> None:
     for name in GRPOMetrics._fields: