Performance Foundations of Parallel & Distributed Reasoning Language Models
Authors: Maciej Besta, Leonard Schmidt, Lara Nonino, Robert Gerstenberger, Pierre Pang, Patrik Okanovic, Ales Kubicek, Tiancheng Chen, Baraq Lipshitz, Torsten Hoefler
Reinforcement Learning with Verifiable Rewards (RLVR) and other RL-style post-training paradigms have been used for aligning large language models (LLMs) with reasoning standards. The resulting recent Reasoning Language Models (RLMs) such as DeepSeek-R1, o3, and Kimi k1.5 show that such RL-style post-training ("RL-for-LLMs") can substantially improve chain-of-thought reasoning, long-horizon planning, and self-correction. However, the computational footprint of these systems is massive: state-of-the-art RLM training requires millions of GPU-hours and tightly coupled multi-model pipelines that stress modern hardware far beyond classical supervised LLM training. This makes RLM training as much a parallel and distributed systems problem as an algorithmic one. In this work, to facilitate developing RLMs that are simultaneously high-performance, scalable, and cost-effective, we first systematize the RL-for-LLM paradigm and provide a compute-centric analysis of prominent post-training algorithmic frameworks: Proximal Policy Optimization (PPO), Group Relative Policy Optimization (GRPO), as well as their variants. Second, we develop a taxonomy of intra- and inter-model parallelism strategies for RL-for-LLMs, covering both traditional techniques (data, tensor, pipeline, sequence, context, and expert parallelism) as well as novel forms of parallelism and optimization techniques for multi-model RLM training, for example disaggregated placement, stage fusion, hybrid parallelism, and asynchronous execution. We harness the work-depth model of parallel computing to make our taxonomy and its insights rigorous and portable. Finally, we analyze existing RLM frameworks and we distill practical guidelines and outline open research directions for building scalable, fast, and cost-effective RLMs.
Figures & tables
Fig. 1 : Overview of the contributions and analyses in this work.
Models & training signals used in the RL pipeline
πθ
Policy (Actor) model: the trainable LLM being optimized.
πref
Reference model: the frozen baseline policy (usually a checkpoint right after SFT).
Rφ
Reward model: an assessment tool that maps prompt–response pairs to scalar rewards.
Vψ
Critic model (Value function): an assessment tool that estimates expected future returns from the reward model.
rb(i),Ab,t(i),G^b,t(i)
Reward for the i -th candidate of prompt b , advantage at token step t , and discounted return estimated at token step t .
Data & Dimensions
TABLE I : Overview of basic mathematical notation and concepts used in the paper. Notation specific to the complexity analyses is detailed separately in Table III .
Fig. 2 : Overview of policy optimization algorithms and used models. No AI was used to conceive or to draw the figure.
Efficient parallel evaluation; no backward pass. Easy to offload.
Reward Rφpr (Process)
Frozen
Assessment
High overhead; sequential decoding if using sequential verification (e.g., MCTS).
Critic (Value) Vψ
Trainable
Assessment, Training
Doubles training compute (fwd+bwd). High memory pressure (optimizer states).
Reference πref
Frozen
Assessment
Required for KL/DPO. Forward-only; often quantized to save memory.
TABLE II : Overview of models used in the RLM pipeline.
Cat.
Symbol
Meaning and remarks
Mathematical expression
Work
Cftok(M)
Forward-pass FLOPs per token for model M on a context of length S+T (input prompt length + response length). The reference policy πref is architecturally identical to the actor πθ . Unembedding FLOPs per token are 2Vdπ (policy only). The reward model Rφ adds a projection head applied to the final token only, contributing S+T2dR FLOPs per token, whereas the value model Vψ applies a projection head to every token, contributing 2dV FLOPs per token.
Backward-pass FLOPs per token for model M on a context of length S+T (input prompt length + response length).
Cbtok(M)≈2⋅Cftok(M)
Cgenroll(πθ)
Autoregressive generation FLOPs per rollout for generating T tokens (response length) from a context of length S (input prompt length) with KV caching (policy only). It is decomposed into prefill and decode. Prefill FLOPs are for processing the prompt of length S once and initializing the KV cache before autoregressive decoding begins. Prefill attention is assumed to be dense; as it is causal only half of the full MM product are needed, giving ≈S2dπ FLOPs for both QKT and for PV ; O(S2) and O(Sd) account for lower-order contributions (respectively – softmax and scaling by 1/dh as well as RoPE and layer norms). Decode FLOPs are for generating the remaining T−1 tokens autoregressively with KV caching. At decode step t , the model attends to context length S+t .
Forward-pass depth of model M on a context of length S+T , defined as the length of the critical path assuming unbounded parallelism. The policy πθ includes an additional unembedding (output projection) step, which contributes a logdπ term to the depth. For the reward model Rφ and the critic Vψ , the additional projection head contributes an extra logd term to the depth.
Forward-pass depth of model M with tensor parallelism degree Pt on a context of length S+T .
DfTP(M)=O(LM[logPtdM+log(S+T)]+logPtdM)
Db(M)
Backward-pass depth of model M on a context of length S+T . It equals the depth of the forward pass for the corresponding model M (which entails computing the gradients with respect to activations), plus the additional term O(log(B(S+T))) coming from computing the gradients with respect to model weights, where one has to accumulate gradients across batch (as these two branches can overlap, the expression given is a conservative bound).
Db(M)=Df(M)+O(log(B(S+T)))
TABLE III : Computational building blocks used across the paper. We assume canonical multi-head attention and two-projection FFN. For computing FLOP costs, an addition and a multiplication count as two separate FLOPs (i.e., a dot product of vectors of dimensionalities d results in d+(d−1)≈2d FLOPs under this model). Cftok , Cbtok , Cgenroll are the costs of (respectively) forward pass, backward pass, and autoregressive generation; defining them per token and per rollout effectively clarifies the notation for the subsequent analyses. The explicit depth expressions count the selected GEMM reduction chains and parameter-gradient accumulation; they suppress the additional logarithmic reductions from normalization, attention/vocabulary softmax, sampling, and distributed collectives. For generation, we explicitly distinguish prefill and decode. Prefill processes the prompt once, while decode generates tokens sequentially with KV caching. This decomposition is important because rollout generation is the dominant online bottleneck and is not equivalent to a single teacher-forced forward pass over a sequence of length S+T .
TABLE IV : Asymptotic Work, Depth, and Memory analysis of PPO, GRPO, and DPO. S is prompt length, T is response length. We denote parameter counts by ∣⋅∣ . “tr” subscript: For conciseness, for training-stage model invocations, we define Ctrtok(M):=Cftok(M)+Cbtok(M)≈3Cftok(M) and Dtr(M):=Df(M)+Db(M) , because a parameter update requires a forward pass followed by backpropagation. Note on Work: Cftok and Cbtok are the cost of forward and backward passes per token and Cgenroll is the autoregressive generation cost per rollout ; these are derived in Table III . Work scales with total tokens (S+T) . Note on Depth: Generation is depth-bound by T (sequential), while in Assessment/Training there is no outer linear T -step autoregressive chain. † DPO is offline; generation occurs prior to the training loop. ‡ Max depth assumes parallel (disaggregated) execution; sequential co-located execution is additive. MKV is Key-Value cache memory; MAct is activation memory (training); MInf is inference buffer (assessment). Note that πref never requires MAct .
Fig. 3 : Mathematical notation and equations for the forward and backward passes of each decoder block. Further details on computational aspects are in Figures 4 - 5 while parallelization details are in Figures 6 - 7 . No AI was used to conceive or to draw the figure.
Fig. 4 : Intra-model execution details (forward pass) for each RLM model execution. The design details are based on Llama-3. Further details on parallelization and the corresponding taxonomy are provided in Figures 6 - 7 . No AI was used to conceive or to draw the figure.
Fig. 5 : Intra-model execution details (backward pass) for each RLM model execution. Legend is provided in Figure 4 . The design details are based on Llama-3. Further details on parallelization and the corresponding taxonomy are provided in Figures 6 - 7 . No AI was used to conceive or to draw the figure.
Fig. 6 : Intra-model parallelism (data, pipeline, tensor, and expert parallelism) details for each RLM model execution; legend is provided in Figure 4 . No AI was used to conceive or to draw the figure.
Fig. 7 : Intra-model parallelism (sequence & context parallelism, as well as prefill vs. decode and flash attention) details for each RLM model execution; legend is provided in Figure 4 . No AI was used to conceive or to draw the figure.
Parallelism
Work
Depth
Memory
None †
O(BmL(d2+md))
O(L[logd+logm]+log(Bm))
O(Ld2+Vd+BmLd)
Data
O(PdBmL(d2+md))
O(L[logd+logm]+log(PdBm))
O(Ld2+Vd+PdBmLd)
Pipeline ∗
O(PpBmL(d2+md))
O(PpL[logd+logm]+log(Bm))
O(PpLd2+Vd+PpBmLd)
Tensor
O(PtBmL(d2+md))
O(L[logPtd+logm]+log(Bm))
O(PtLd2+Vd+PtBmLd)
Context
O(PcBmL(d2+md))
O(L[logd+logm]+log(PcBm))
O(Ld2+Vd+PcBmLd)
Expert
O(BmL[(d2+md)+PeEadde])
O(L[logd+logm]+log(Bm))
O(L[d2+PeEdde]+Vd+BmLd)
TABLE V : Complexity Analysis (Per rank). Asymptotic work, depth, and memory for a single Transformer training iteration (forward and backward pass), expressed in Big- O notation. For clarity, we use m:=S+T . The work expressions use the regime where dff=O(d) . † : Baseline configuration without parallelism. ∗ For PP, the Vd memory term denotes the embedding/output state resident on boundary pipeline ranks; interior ranks need not store this state. Thus, this term represents the boundary/peak per-rank footprint rather than replication across all Pp ranks. ‡ : 3D and 5D parallelism combine data, pipeline, tensor (3D), context, and expert (5D) parallelism, with PdPpPt=N (3D) and PdPpPtPcPe=N (5D) total devices. For the combined 5D configuration, we assume the EP dimension also partitions source-token ownership for shared-layer computation. Training memory assumes Adam without activation checkpointing.
Parallelism
Work
Depth
Memory
None
O(BmL(d2+md))
O(L[logd+logm]+log(Bm))
O(Ld2+Vd+BmLd)
Data
O(BmL(d2+md))
O(L[logd+logm]+log(Bm))
O(Pd[Ld2+Vd]+BmLd)
Pipeline
O(BmL(d2+md))
O(L[logd+logm]+log(Bm))
O(Ld2+Vd+BmLd)
Tensor
O(BmL(d2+md))
O(L[logd+logm]+log(Bm))
O(Ld2+Vd+BmLd)
Context
O(BmL(d2+md))
O(L[logd+logm]+log(Bm))
O(Pc[Ld2+Vd]+BmLd)
Expert
O(BmL[Pe(d2+md)+Eadde])
O(L[logd+logm]+log(Bm))
O(Pe[Ld2+Vd+BmLd]+LEdde)
TABLE VI : Complexity Analysis (Global). Asymptotic work, depth, and memory for a single Transformer training iteration (forward & backward pass). For clarity, we use m:=S+T . The work expressions use the regime where dff=O(d) .
Parallelism
Work
Depth
Memory
None
6BmL(4d2+2ddff+md)
2Llog(hd4dffm)+log(Bm)
L(16d2+8ddff)+4Vd+BmLd
Data
6PdBmL(4d2+2ddff+md)
2Llog(hd4dffm)+log(PdBm)
L(16d2+8ddff)+4Vd+PdBmLd
Pipeline
6BmPpL(4d2+2ddff+md)
2PpLlog(hd4dffm)+log(Bm)
PpL(16d2+8ddff)+4Vd+BmPpLd
Tensor
6BmL(Pt4d2+2ddff+md)
2Llog(Pt2hd4dffm)+log(Bm)
PtL(16d2+8ddff)+Pt4Vd+PtBmLd
Context
6BPcS+TL(4d2+2ddff+md)
2Llog(hd4dffm)+log(PcBm)
L(16d2+8ddff)+4Vd+BPcS+TLd
Expert
6BmL(4d2+Pe2Eadde+md)
2Llog(hd4dem)+log(Bm)
L(16d2+Pe8Edde)+4Vd+BmLd
TABLE VII : Complexity Analysis (Per rank). Explicit leading-order costs under our simplified Transformer cost model, for a single Transformer training iteration (forward and backward pass), including explicit constant factors. For clarity, we use m:=S+T .
Parallelism
Work
Depth
Memory
None
6BmL(4d2+2ddff+md)
2Llog(d4dffm/h)+log(Bm)
L(16d2+8ddff)+4Vd+BmLd
Data
6BmL(4d2+2ddff+md)
2Llog(d4dffm/h)+log(Bm)
Pd[L(16d2+8ddff)+4Vd]+BmLd
Pipeline
6BmL(4d2+2ddff+md)
2Llog(d4dffm/h)+log(Bm)
L(16d2+8ddff)+4Vd+BmLd
Tensor
6BmL(4d2+2ddff+md)
2Llog(d4dffm/h)+log(Bm)
L(16d2+8ddff)+4Vd+BmLd
Context
6BmL(4d2+2ddff+md)
2Llog(d4dffm/h)+log(Bm)
Pc[L(16d2+8ddff)+4Vd]+BmLd
Expert
6BmL(Pe(4d2+md)+2Eadde)
2Llog(d4dem/h)+log(Bm)
Pe(16Ld2+4Vd+BmLd)+8LEdde
TABLE VIII : Complexity Analysis (Global). Explicit leading-order costs under our simplified Transformer cost model, for a single Transformer training iteration (forward and backward pass), including explicit constant factors. For clarity, we use m:=S+T .
Fig. 8 : Inter-model parallelism details (asynchronous execution, hybrid execution, stage fusion). No AI was used to conceive or to draw the figure.
Fig. 9 : Inter-model parallelism details (placement strategies, model structure configurations). No AI was used to conceive or to draw the figure.
TABLE IX : Inter-model parallelism. Global work, global depth, and peak memory per device group reported for the Generation, Assessment, and Training stages of a single RL iteration within one epoch under different inter-model parallelism techniques. Global work and depth report total FLOPs and the critical-path length per stage. Peak memory is measured per device group and taken as the maximum over all groups active within a stage; in configurations without disaggregation, all models share the same device group and execute through temporal multiplexing, where different models become active at different stages. “TP” superscript & “SH” subscript: A TP superscript denotes tensor-parallel execution; πSH denotes the shared actor–critic configuration. “tr” subscript: For conciseness, for training-stage model invocations, we define Ctrtok(M):=Cftok(M)+Cbtok(M)≈3Cftok(M) and Dtr(M):=Df(M)+Db(M) , because a parameter update requires a forward pass followed by backpropagation. Msh(πθ) denotes the persistent inference-side shadow copy of the actor required when asynchronous consumers use a bounded-stale snapshot. Also, sample-time actor log-probabilities are assumed to be computed on the fly during Generation from the same logits used for sampling; therefore they do not add a separate actor forward pass. Baselines: The Baseline configuration sequentially (i.e., one device group) executes all RL components, resulting in additive depth terms. AC shared uses a shared Transformer backbone for actor and critic. All models run sequentially in a single device group. Disaggregated execution assigns models to distinct device groups, allowing concurrent execution. The Hybrid configuration applies TP to Generation, TP/PP to Assessment, and ZeRO-3 to Training. All models share the same device group. Stage F. applies inter-stage fusion to Generation and Assessment and intra-stage fusion to Training. Asynchronous execution overlaps Generation and Assessment with Training by allowing bounded weight staleness. Combined execution combines asynchronous overlap with disaggregated model placement and hybrid intra-model parallelism, following the RLHFuse execution model. Training memory assumes the Adam optimizer without activation checkpointing, accounting for parameters, gradients, optimizer states, and activations.
Model / family
Access
Scale / arch.
Post-training signal
Selected system details
General aligned large language models (LLMs)
InstructGPT [ 192 ]
Closed
1.3B, 6B, 175B dense
PPO-like (SFT + RM + PPO)
Canonical actor–reward–reference PPO loop.
GPT-4 [ 4 ]
Closed
Undisclosed
PPO-like (rule-based rewards)
Frontier-scale alignment; limited public systems detail.
Gemini 1.x [ 12 ]
Closed
Undisclosed
PPO-like (iterative RM refinement)
Repeated RM/RL cycles increase assessment cost.
Claude 3 [ 13 , 21 ]
Closed
Undisclosed
PPO-like (RLAIF / Constitutional AI)
AI-generated preference labels from written principles.
Reka [ 191 ]
Mixed
7B, 21B dense
PPO-like (multi-round RLHF)
Standard PPO-style alignment at moderate scale.
TABLE X : A comparison of representative general aligned LLMs and RLMs. The upper block lists general aligned LLMs, which primarily target instruction following, safety, preference alignment, and broad assistant quality. The lower block lists RLMs, which explicitly target reasoning trajectories, long-CoT, RLVR/verifiable rewards, search, tools, multi-agent inference, or controllable test-time computation. The “Post-training signal” column follows the taxonomy used in this paper: PPO-like methods use reward-model- or RLAIF-based policy optimization; GRPO-like methods use critic-free grouped RL/RLVR; DPO-like methods use direct preference optimization; Hybrid denotes mixed or undisclosed pipelines; Search/Agentic denotes explicit search, tools, environments, or multi-agent reasoning. Proprietary details are based on public reports and should be treated as approximate.
Intra-Model Parallelism
Inter-Model Parallelism
Framework
Remarks
DP
Z1
Z2
Z3
TP
CP
EP
PP
AC Share
Dis- agg.
Stage Fusion
Hybrid
Async
TRL [ 231 ]
Accelerate; TP/vLLM mostly rollout; CP via FSDP2/Ulysses
ColossalChat [ 267 ]
Colossal-AI HybridParallel + Gemini; TP/PP/SP, Z1/2 and Z3
DeepSpeed-Chat [ 261 ]
Train ZeRO-2/3; Hybrid Engine uses TP mainly for rollout
TABLE XI : Comparison of public training / inference frameworks for LLM and RLM post-training under the systems taxonomy. Intra-Model Parallelism : DP=data parallel; Z1/Z2/Z3=ZeRO/FSDP-like sharding levels (optimizer / optimizer+gradients / parameters+gradients+optimizer; FSDP2 is counted under Z3); TP=tensor parallel; CP=context or long-sequence parallelism; EP=expert parallelism; PP=pipeline parallelism. Inter-Model Parallelism : AC Share=shared actor/critic parameters or value-head mode; Disagg.=separate services or device groups for rollout, reward, or training; Stage Fusion=explicit overlap or fusion across normally sequential RL stages or subtasks; Hybrid=stage-specific backends/layouts or runtime resharding/refit; Async=asynchronous or bounded-stale execution. Symbols : documented support; indirect, partial, backend-specific, or train-vs-rollout-specific support; no documented public support; undocumented or unclear in public materials.