Large language models (LLMs) have achieved strong performance on reasoning benchmarks, yet their ability to solve real-world problems requiring end-to-end workflows remains unclear. Mathematical modeling competitions provide a stringent testbed for evaluating such end-to-end problem-solving capability. We propose a problem-oriented, stage-wise evaluation framework that assesses LLM performance across modeling stages using expert-verified criteria. We validate the framework's reliability by comparing automatic scores with independent human expert judgments on problems from the China Postgraduate Mathematical Contest in Modeling, demonstrating substantially stronger alignment than existing evaluation schemes. Using this framework, we reveal a comprehension-execution gap in state-of-the-art LLMs: while they perform well in early stages such as problem identification and formulation, they exhibit persistent deficiencies in execution-oriented stages including model solving, code implementation, and result analysis. These gaps persist even with increased model scale. We further trace these failures to insufficient specification, missing verification, and lack of validation, with errors propagating across stages without correction. Our findings suggest that bridging this gap requires approaches beyond model scaling, offering insights for applying LLMs to complex real-world problem solving.
Figures & tables
图 1: Problem-oriented, stage-wise evaluation framework. An LLM first decomposes a mathematical modeling problem into a set of subtasks, which are guided and refined by domain experts to preserve the original intent. For each verified subtask, a stage-aware evaluation rubric is instantiated by conditioning on the major stages of the modeling process. For each subtask–stage pair, the LLM generates concrete evaluation criteria under expert guidance, which serve as the atomic units for scoring and are further verified and revised. The resulting subtask–stage rubric defines a fixed, problem-specific evaluation structure applied uniformly across models.
Baseline
Ours
ICC(2,1) w.r.t. Expert
0.012
0.673
表 1: Agreement between independent expert overall judgments and automatic scores.
图 2: Stage-wise agreement between automatic and human expert scores measured by ICC(2,1). Higher values indicate stronger alignment with expert assessment. The boxplots correspond to individual modeling stages (Prb Idf, Prb Frm, Asm Dev, Mod Con, Mod Sol, Cod Imp, Res Ays), following the standard mathematical modeling workflow.
图 3: Distribution of evaluation scores under the baseline rubric and our framework. Base Avg denotes the report-level baseline score, and Our Avg denotes the overall score under our problem-oriented, stage-aware framework. The remaining boxplots correspond to individual modeling stages (Prb Idf, Prb Frm, Asm Dev, Mod Con, Mod Sol, Cod Imp, Res Ays).
图 4: Stage-wise performance of four LLMs across the mathematical modeling pipeline. DS-Inst and DS-Think denote two DeepSeek-V3.2 variants, and Qw3-235B denotes a Qwen3-235B-Instruct model. Scores are averaged over reports and subtasks for each stage.
图 5: Stage-wise performance under within-family scaling for Qwen models. Bars correspond to Qw-7B – Qw-72B (Qwen2.5) and Qw-235B (Qwen3-235B).
Stage
Common Failure Causes
Assumption Development
Assumptions Not Checked; Missing Assumption Conditions; Assumption Impact Not Discussed; Hidden Idealized Assumptions; Unrealistic Assumptions
Model Construction
Incomplete Model Structure; Missing Model Derivation; Unclear Variables or Symbols; Solvability Not Justified; Model Deviates from Task Goal; Model-Assumption Conflict
Model Solving
No Checkable Solution; Key Solution Steps Missing; Solution Not Verified; Numerical Stability Not Analyzed; Computationally Infeasible; Inappropriate Solution Method
Code Implementation
Cannot be executed; Results Not Reproducible; No Code; Engineering or Numerical Risks; Code–Model Mismatch
Result Analysis
No Meaningful Results; No Validation or Comparison; Sensitivity Not Analyzed; Conclusions Not Supported; Results Miss the Goal; Limits Not Discussed
表 2: Stage-specific failure causes across the mathematical modeling workflow. For each stage, failure causes are ordered from left to right by decreasing prevalence, based on their frequency of occurrence in reviewer comments associated with low-scoring cases. The detailed aggregation procedure of these frequencies is described in Appendix H .
Appendix figures & tables7 assets
Supplementary material from the paper’s appendix.
Appendix
Stage
Failure Cause
Operational Definition
Assumption Development
Assumptions Not Checked
Assumptions that are in principle testable are stated but no validation, diagnostic procedure, or empirical check is provided.
Assumption Development
Missing Assumption Conditions
Necessary conditions under which an assumption holds (e.g., constraints, statistical tests, parameter bounds) are not explicitly specified.
Assumption Development
Assumption Impact Not Discussed
The influence of assumptions on model behavior, results, or conclusions is not analyzed or discussed.
Assumption Development
Hidden Idealized Assumptions
Idealized premises (e.g., noise-free data, full observability, stationarity) are implicitly adopted without being explicitly declared.
Assumption Development
Unrealistic Assumptions
Assumptions are incompatible with the problem context, data-generating process, or domain knowledge, even if explicitly stated.
Model Construction
Incomplete Model Structure
The model lacks essential components such as key variables, constraints, objective functions, or state relationships required to represent the task.
Appendix
表 4: Operational Definitions of Stage-Specific Failure Causes Corresponding to Table 2
Model
Not Checked
Miss Cond.
Impact Ignored
Hidden Ideal.
Unrealistic
Qwen3-235B
0.7000
0.5625
0.4250
0.4125
0.0000
DS-Inst
0.6743
0.5829
0.3771
0.4629
0.0057
DS-Think
0.6872
0.6089
0.4525
0.3799
0.0056
o4-mini
0.7412
0.6518
0.4441
0.4409
0.0128
Appendix
表 5: Distribution of primary failure causes in the assumption development stage. Model abbreviations: Qwen3-235B = Qwen3-235B-A22B-Instruct-2507; DS-Inst = DeepSeek-V3.2-Instruct; DS-Think = DeepSeek-V3.2-Thinking. Failure causes denote: Not Checked (Assumptions Not Checked), Miss Cond. (Missing Assumption Conditions), Impact Ignored (Assumption Impact Not Discussed), Hidden Ideal. (Hidden Idealized Assumptions), and Unrealistic (Unrealistic Assumptions). Total denotes the number of subtasks whose assumption-stage score is below 8.
Model
Incomp. Struct.
Var. Unclear
Miss Deriv.
Solv. Unjust.
Goal Deviat.
Assump. Conflict
DS-Inst
0.6604
0.4717
0.2075
0.0377
0.0377
0.0566
DS-Think
0.8030
0.3182
0.0909
0.0303
0.1364
0.0303
Qwen3-235B
0.8125
0.3750
0.1250
0.0625
0.0625
0.0000
o4-mini
0.4560
0.6080
0.1840
0.1920
0.0320
0.0240
Appendix
表 6: Distribution of primary failure causes in the model construction stage. Model abbreviations: DS-Inst = DeepSeek-V3.2-Instruct; DS-Think = DeepSeek-V3.2-Thinking; Qwen3-235B = Qwen3-235B-A22B-Instruct-2507. Failure causes include: Miss Deriv. (Missing Model Derivation), Incomp. Struct. (Incomplete Model Structure), Goal Deviat. (Model Deviates from Task Goal), Var. Unclear (Unclear Variables or Symbols), Solv. Unjust. (Solvability Not Justified), and Assump. Conflict (Model–Assumption Conflict). Total denotes the number of subtasks whose model-construction score is below 8. f
Model
No Output
Steps Missing
Not Verified
No Stability
Infeasible
Wrong Method
DS-Inst
0.3592
0.4296
0.3099
0.3310
0.0704
0.0070
DS-Think
0.3896
0.4481
0.3506
0.2727
0.0844
0.0260
Qwen3-235B
0.5287
0.3333
0.2414
0.1954
0.0345
0.0115
o4-mini
0.5855
0.2566
0.5263
0.2730
0.0658
0.0197
Appendix
表 7: Distribution of primary failure causes in the model solving stage. Model abbreviations: DS-Inst = DeepSeek-V3.2-Instruct; DS-Think = DeepSeek-V3.2-Thinking; Qwen3-235B = Qwen3-235B-A22B-Instruct-2507. Failure causes include: No Output (No Checkable Solution), Steps Missing (Key Solution Steps Missing), Not Verified (Solution Not Verified), No Stability (Numerical Stability Not Analyzed), Infeasible (Computationally Infeasible), and Wrong Method (Inappropriate Solution Method). Total denotes the number of subtasks whose model-solving score is below 8.
Model
Results N/R
No Code
Miss Comp.
Eng./Num. Risk
Code–Model Mis.
DS-Inst
0.7564
0.5782
0.1564
0.0800
0.0036
DS-Think
0.7600
0.6133
0.2100
0.0933
0.0000
Qwen3-235B
0.7238
0.6000
0.2333
0.1143
0.0048
o4-mini
0.7646
0.8915
0.3175
0.0556
0.0026
Appendix
表 8: Distribution of primary failure causes in the code implementation stage. Model abbreviations: DS-Inst = DeepSeek-V3.2-Instruct; DS-Think = DeepSeek-V3.2-Thinking; Qwen3-235B = Qwen3-235B-A22B-Instruct-2507. Failure causes include: Results N/R (Results Not Reproducible), No Code (No Usable Code), Miss Comp. (Missing Key Code Components), Eng./Num. Risk (Engineering or Numerical Risks), and Code–Model Mis. (Code–Model Mismatch). Total denotes the number of subtasks whose code-implementation score is below 8.
Model
No Valid.
No Results
No Sensit.
Weak Concl.
Goal Miss.
No Limits
DS-Inst
0.5417
0.4722
0.3565
0.2639
0.0417
0.0463
DS-Think
0.6161
0.4420
0.3839
0.2991
0.0536
0.0313
Qwen3-235B
0.6182
0.5818
0.3333
0.2667
0.0061
0.0364
o4-mini
0.6183
0.7984
0.3495
0.2957
0.0323
0.0349
Appendix
表 9: Distribution of primary failure causes in the result analysis stage. Model abbreviations: DS-Inst = DeepSeek-V3.2-Instruct; DS-Think = DeepSeek-V3.2-Thinking; Qwen3-235B = Qwen3-235B-A22B-Instruct-2507. Failure causes include: No Valid. (No Validation or Comparison), No Results (No Meaningful Results), No Sensit. (Sensitivity Not Analyzed), Weak Concl. (Conclusions Not Supported), Goal Miss. (Results Miss the Goal), and No Limits (Limits Not Discussed). Total denotes the number of subtasks whose result-analysis score is below 8.
图 6: Heatmap of co-occurrence between mathematical method categories and application domains. Color intensity indicates the number of modeling problems associated with each method–domain pair. Each problem may belong to multiple mathematical methods and application domains; counts therefore reflect multi-label co-occurrences. Abbreviations are used for compact visualization: Dat–ML (Data-driven and Machine Learning), Opt (Optimization), Sta (Statistical), Dis–Gra (Discrete and Graph-based), Sim (Simulation-based), Dif–Eq (Differential Equation), Pro–Sto (Probabilistic and Stochastic), Hyb (Hybrid); Eng–Ind (Engineering and Industrial), Env–Res (Environmental and Resource), Inf–Com (Information and Communication), Phy–Nat (Physical and Natural), Hea–Pub (Healthcare and Public Health), Bio–Lif (Biological and Life), Urb–Tra (Urban and Transportation), Eco–Fin (Economic and Financial), and Soc–Pub (Social and Public).
Large Language Models (LLMs) are increasingly used in math education not only as problem solvers but also as assessors of learners' reasoning. However, it remains unclear whether stronger math problem-solving ability is associated with stronger step-level assessment performance. This study examines that relationship using the GSM8K and MATH subsets of PROCESSBENCH, a human-annotated benchmark for identifying the earliest erroneous step in mathematical reasoning. We evaluate two LLM-based math tutor agent settings, instantiated with GPT-4 and GPT-5, in two independent tasks on the same math problems: solving the original problem and assessing a benchmark-provided solution by predicting the earliest erroneous step. Results show a consistent within-model pattern: assessment accuracy is substantially higher on math problem items the same model solved correctly than on items it solved incorrectly, with statistically significant associations across both models and datasets. At the same time, assessment remains more difficult than direct problem solving, especially on error-present solutions. These findings suggest that math problem-solving expertise supports stronger assessment performance, but reliable step-level diagnosis also requires additional capabilities such as step tracking, monitoring, and precise error localization. The results have implications for the design and evaluation of AI-supported Adaptive Instructional Systems (AISs) for formative assessment in math education.
Liang Zhang, Yu Fu, Xinyi Jin
University of Michigan, Ann Arbor, MI, USA · New York, NY, USA · The High School Affiliated to Minzu University of China, Beijing, China
Recent advancements in large language models have led to significant improvements across various tasks, including mathematical reasoning, which is used to assess models' intelligence in logical reasoning and problem-solving. Models are evaluated on mathematical reasoning benchmarks by verifying the correctness of the final answer against a ground truth answer. A common approach for this verification is based on symbolic mathematics comparison, which fails to generalize across diverse mathematical representations and solution formats. In this work, we offer a robust and flexible alternative to rule-based symbolic mathematics comparison. We propose an LLM-based evaluation framework for evaluating model-generated answers, enabling accurate evaluation across diverse mathematical representations and answer formats. We present failure cases of symbolic evaluation in two popular frameworks, Lighteval and SimpleRL, and compare them to our approach, demonstrating clear improvements over commonly used methods. Our framework enables more reliable evaluation and benchmarking, leading to more accurate performance monitoring, which is important for advancing mathematical problem-solving and intelligent systems.
Erez Yosef, Oron Anschel, Shunit Haviv Hakimi +4
Tel-Aviv University · Amazon Prime Video · Ben-Gurion University
Large language models (LLMs) are becoming increasingly capable mathematical collaborators, but static benchmarks are no longer sufficient for evaluating progress: they are often narrow in scope, quickly saturated, and rarely updated. This makes it hard to compare models reliably and track progress over time. Instead, we need evaluation platforms: continuously maintained systems that run, aggregate, and analyze evaluations across many benchmarks to give a comprehensive picture of model performance within a broad domain. In this work, we build on the original MathArena benchmark by substantially broadening its scope from final-answer olympiad problems to a continuously maintained evaluation platform for mathematical reasoning with LLMs. MathArena now covers a much wider range of tasks, including proof-based competitions, research-level arXiv problems, and formal proof generation in Lean. Additionally, we maintain a clear evaluation protocol for all models and regularly design new benchmarks as model capabilities improve to ensure that MathArena remains challenging. Notably, the strongest model, GPT-5.5, now reaches 98% on the 2026 USA Math Olympiad and 74% on research-level questions, showing that frontier models can now comfortably solve extremely challenging mathematical problems. This highlights the importance of continuously maintained evaluation platforms like MathArena to track the rapid progress of LLMs in mathematical reasoning.
Jasper Dekoninck, Nikola Jovanović, Tim Gehrunger +4
ETH Zurich · INSAIT · Sofia University "St. Kliment Ohridski"