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).