Skip to main content

litex/verify/verify_builtin_rules/
number_compare.rs

1use super::order_normalize::normalize_positive_order_atomic_fact;
2use crate::prelude::*;
3
4pub enum NumberCompareResult {
5    Less,
6    Equal,
7    Greater,
8}
9
10/// Compare a normalized decimal string (same shape as [`Number::normalized_value`]) to `"0"`.
11pub fn compare_normalized_number_str_to_zero(number_value: &str) -> NumberCompareResult {
12    compare_number_strings(number_value.trim(), "0")
13}
14
15fn parse_number_parts_for_comparison(number_value: &str) -> (bool, Vec<u8>, Vec<u8>) {
16    let trimmed_number_value = number_value.trim();
17    let (is_negative, magnitude_string) = if trimmed_number_value.starts_with('-') {
18        (true, trimmed_number_value[1..].trim())
19    } else {
20        (false, trimmed_number_value)
21    };
22
23    let (integer_part_string, fractional_part_string) = match magnitude_string.find('.') {
24        Some(dot_index) => (
25            &magnitude_string[..dot_index],
26            &magnitude_string[dot_index + 1..],
27        ),
28        None => (magnitude_string, ""),
29    };
30
31    let mut integer_digits: Vec<u8> = Vec::new();
32    if integer_part_string.is_empty() {
33        integer_digits.push(0);
34    } else {
35        for current_char in integer_part_string.chars() {
36            if current_char.is_ascii_digit() {
37                integer_digits.push(current_char as u8 - b'0');
38            }
39        }
40        if integer_digits.is_empty() {
41            integer_digits.push(0);
42        }
43    }
44
45    let mut fractional_digits: Vec<u8> = Vec::new();
46    for current_char in fractional_part_string.chars() {
47        if current_char.is_ascii_digit() {
48            fractional_digits.push(current_char as u8 - b'0');
49        }
50    }
51
52    (is_negative, integer_digits, fractional_digits)
53}
54
55fn digits_are_all_zero(digits: &[u8]) -> bool {
56    for digit in digits {
57        if *digit != 0 {
58            return false;
59        }
60    }
61    true
62}
63
64fn normalized_decimal_string_is_integer(number_value: &str) -> bool {
65    let (_, _integer_digits, fractional_digits) = parse_number_parts_for_comparison(number_value);
66    digits_are_all_zero(&fractional_digits)
67}
68
69fn normalized_decimal_string_is_even_integer(number_value: &str) -> bool {
70    if !normalized_decimal_string_is_integer(number_value) {
71        return false;
72    }
73    let (_is_negative, integer_digits, _fractional_digits) =
74        parse_number_parts_for_comparison(number_value);
75    let last_digit = integer_digits.last().copied().unwrap_or(0);
76    last_digit % 2 == 0
77}
78
79fn first_non_zero_integer_digit_index(integer_digits: &[u8]) -> usize {
80    let mut current_index = 0;
81    while current_index + 1 < integer_digits.len() && integer_digits[current_index] == 0 {
82        current_index += 1;
83    }
84    current_index
85}
86
87fn compare_non_negative_decimal_parts(
88    left_integer_digits: &[u8],
89    left_fractional_digits: &[u8],
90    right_integer_digits: &[u8],
91    right_fractional_digits: &[u8],
92) -> NumberCompareResult {
93    let left_integer_start_index = first_non_zero_integer_digit_index(left_integer_digits);
94    let right_integer_start_index = first_non_zero_integer_digit_index(right_integer_digits);
95
96    let left_effective_integer_len = left_integer_digits.len() - left_integer_start_index;
97    let right_effective_integer_len = right_integer_digits.len() - right_integer_start_index;
98    if left_effective_integer_len < right_effective_integer_len {
99        return NumberCompareResult::Less;
100    }
101    if left_effective_integer_len > right_effective_integer_len {
102        return NumberCompareResult::Greater;
103    }
104
105    let mut integer_index = 0;
106    while integer_index < left_effective_integer_len {
107        let left_digit = left_integer_digits[left_integer_start_index + integer_index];
108        let right_digit = right_integer_digits[right_integer_start_index + integer_index];
109        if left_digit < right_digit {
110            return NumberCompareResult::Less;
111        }
112        if left_digit > right_digit {
113            return NumberCompareResult::Greater;
114        }
115        integer_index += 1;
116    }
117
118    let fractional_compare_len = if left_fractional_digits.len() > right_fractional_digits.len() {
119        left_fractional_digits.len()
120    } else {
121        right_fractional_digits.len()
122    };
123    let mut fractional_index = 0;
124    while fractional_index < fractional_compare_len {
125        let left_digit = match left_fractional_digits.get(fractional_index) {
126            Some(digit) => *digit,
127            None => 0,
128        };
129        let right_digit = match right_fractional_digits.get(fractional_index) {
130            Some(digit) => *digit,
131            None => 0,
132        };
133        if left_digit < right_digit {
134            return NumberCompareResult::Less;
135        }
136        if left_digit > right_digit {
137            return NumberCompareResult::Greater;
138        }
139        fractional_index += 1;
140    }
141
142    NumberCompareResult::Equal
143}
144
145pub fn compare_number_strings(
146    left_number_value: &str,
147    right_number_value: &str,
148) -> NumberCompareResult {
149    let (left_is_negative, left_integer_digits, left_fractional_digits) =
150        parse_number_parts_for_comparison(left_number_value);
151    let (right_is_negative, right_integer_digits, right_fractional_digits) =
152        parse_number_parts_for_comparison(right_number_value);
153
154    let left_is_zero =
155        digits_are_all_zero(&left_integer_digits) && digits_are_all_zero(&left_fractional_digits);
156    let right_is_zero =
157        digits_are_all_zero(&right_integer_digits) && digits_are_all_zero(&right_fractional_digits);
158    if left_is_zero && right_is_zero {
159        return NumberCompareResult::Equal;
160    }
161
162    if left_is_negative && !left_is_zero && !right_is_negative {
163        return NumberCompareResult::Less;
164    }
165    if right_is_negative && !right_is_zero && !left_is_negative {
166        return NumberCompareResult::Greater;
167    }
168
169    let non_negative_compare_result = compare_non_negative_decimal_parts(
170        &left_integer_digits,
171        &left_fractional_digits,
172        &right_integer_digits,
173        &right_fractional_digits,
174    );
175    if left_is_negative && !left_is_zero && right_is_negative && !right_is_zero {
176        return match non_negative_compare_result {
177            NumberCompareResult::Less => NumberCompareResult::Greater,
178            NumberCompareResult::Equal => NumberCompareResult::Equal,
179            NumberCompareResult::Greater => NumberCompareResult::Less,
180        };
181    }
182
183    non_negative_compare_result
184}
185
186impl Runtime {
187    /// Sub-goals inside numeric builtins: known env + builtin rules only.
188    /// Do not call [`Runtime::verify_non_equational_atomic_fact`] here: its forall / definition
189    /// round can recurse with outer goals (e.g. `b in R` for `0 <= a^b`, or order lemmas).
190    pub(crate) fn verify_non_equational_known_then_builtin_rules_only(
191        &mut self,
192        atomic_fact: &AtomicFact,
193        verify_state: &VerifyState,
194    ) -> Result<StmtResult, RuntimeError> {
195        let r = self.verify_non_equational_atomic_fact_with_known_atomic_facts(atomic_fact)?;
196        if r.is_true() {
197            return Ok(r);
198        }
199        self.verify_non_equational_atomic_fact_with_builtin_rules(atomic_fact, verify_state)
200    }
201
202    fn verify_zero_order_on_sub_expr(
203        &mut self,
204        zero: &Obj,
205        sub_expr: &Obj,
206        weak: bool,
207        line_file: &LineFile,
208    ) -> Result<StmtResult, RuntimeError> {
209        let fact: AtomicFact = if weak {
210            LessEqualFact::new(zero.clone(), sub_expr.clone(), line_file.clone()).into()
211        } else {
212            LessFact::new(zero.clone(), sub_expr.clone(), line_file.clone()).into()
213        };
214        let mut result = self.verify_non_equational_atomic_fact_with_known_atomic_facts(&fact)?;
215        if !result.is_true() {
216            result = self.verify_order_atomic_fact_numeric_builtin_only(&fact)?;
217        }
218        Ok(result)
219    }
220
221    /// `n >= 0` / `0 <= n` from known `n $in N` (e.g. `forall n N:` domain).
222    fn try_verify_order_nonnegative_from_membership_in_n(
223        &mut self,
224        atomic_fact: &AtomicFact,
225        verify_state: &VerifyState,
226    ) -> Result<Option<StmtResult>, RuntimeError> {
227        let (n, line_file) = match atomic_fact {
228            AtomicFact::GreaterEqualFact(f) => {
229                let Some(z) = self.resolve_obj_to_number(&f.right) else {
230                    return Ok(None);
231                };
232                if !matches!(
233                    compare_normalized_number_str_to_zero(&z.normalized_value),
234                    NumberCompareResult::Equal
235                ) {
236                    return Ok(None);
237                }
238                (f.left.clone(), f.line_file.clone())
239            }
240            AtomicFact::LessEqualFact(f) => {
241                let Some(z) = self.resolve_obj_to_number(&f.left) else {
242                    return Ok(None);
243                };
244                if !matches!(
245                    compare_normalized_number_str_to_zero(&z.normalized_value),
246                    NumberCompareResult::Equal
247                ) {
248                    return Ok(None);
249                }
250                (f.right.clone(), f.line_file.clone())
251            }
252            _ => return Ok(None),
253        };
254        let in_n: AtomicFact = InFact::new(n, StandardSet::N.into(), line_file.clone()).into();
255        if self
256            .verify_non_equational_known_then_builtin_rules_only(&in_n, verify_state)?
257            .is_true()
258        {
259            return Ok(Some(StmtResult::FactualStmtSuccess(
260                FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
261                    atomic_fact.clone().into(),
262                    "n >= 0 from n $in N".to_string(),
263                    Vec::new(),
264                ),
265            )));
266        }
267        Ok(None)
268    }
269
270    /// `n >= 1` / `1 <= n` from known `n $in N_pos`.
271    fn try_verify_order_one_le_from_membership_in_n_pos(
272        &mut self,
273        atomic_fact: &AtomicFact,
274        verify_state: &VerifyState,
275    ) -> Result<Option<StmtResult>, RuntimeError> {
276        let (n, line_file) = match atomic_fact {
277            AtomicFact::GreaterEqualFact(f) => {
278                let Some(one) = self.resolve_obj_to_number(&f.right) else {
279                    return Ok(None);
280                };
281                if !matches!(
282                    compare_number_strings(&one.normalized_value, "1"),
283                    NumberCompareResult::Equal
284                ) {
285                    return Ok(None);
286                }
287                (f.left.clone(), f.line_file.clone())
288            }
289            AtomicFact::LessEqualFact(f) => {
290                let Some(one) = self.resolve_obj_to_number(&f.left) else {
291                    return Ok(None);
292                };
293                if !matches!(
294                    compare_number_strings(&one.normalized_value, "1"),
295                    NumberCompareResult::Equal
296                ) {
297                    return Ok(None);
298                }
299                (f.right.clone(), f.line_file.clone())
300            }
301            _ => return Ok(None),
302        };
303        let in_n_pos: AtomicFact =
304            InFact::new(n, StandardSet::NPos.into(), line_file.clone()).into();
305        if self
306            .verify_non_equational_known_then_builtin_rules_only(&in_n_pos, verify_state)?
307            .is_true()
308        {
309            return Ok(Some(StmtResult::FactualStmtSuccess(
310                FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
311                    atomic_fact.clone().into(),
312                    "n >= 1 from n $in N_pos".to_string(),
313                    Vec::new(),
314                ),
315            )));
316        }
317        Ok(None)
318    }
319
320    fn verify_zero_le_abs_builtin_rule(
321        &mut self,
322        atomic_fact: &AtomicFact,
323    ) -> Result<Option<StmtResult>, RuntimeError> {
324        let Some(norm) = normalize_positive_order_atomic_fact(atomic_fact) else {
325            return Ok(None);
326        };
327        let AtomicFact::LessEqualFact(f) = &norm else {
328            return Ok(None);
329        };
330        if f.left.to_string() != "0" {
331            return Ok(None);
332        }
333        if !matches!(&f.right, Obj::Abs(_)) {
334            return Ok(None);
335        }
336        Ok(Some(StmtResult::FactualStmtSuccess(
337            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
338                atomic_fact.clone().into(),
339                "0 <= abs(x) for x in R".to_string(),
340                Vec::new(),
341            ),
342        )))
343    }
344
345    // `(-1)*x` order vs 0: e.g. `x < 0` or `x <= 0` implies `(-1)*x >= 0`; `x > 0` implies `(-1)*x < 0`.
346    // Also handles `0 <= (-1)*x` (equivalently `0 <= -x` when `-x` parses as `(-1)*x`).
347    fn try_verify_order_opposite_sign_mul_minus_one(
348        &mut self,
349        atomic_fact: &AtomicFact,
350        verify_state: &VerifyState,
351    ) -> Result<Option<StmtResult>, RuntimeError> {
352        let z: Obj = Number::new("0".to_string()).into();
353        let success = |msg: &'static str| {
354            Ok(Some(StmtResult::FactualStmtSuccess(
355                FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
356                    atomic_fact.clone().into(),
357                    msg.to_string(),
358                    Vec::new(),
359                ),
360            )))
361        };
362        match atomic_fact {
363            AtomicFact::GreaterEqualFact(f) if self.obj_is_resolved_zero(&f.right) => {
364                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.left) {
365                    let le: AtomicFact =
366                        LessEqualFact::new(x.clone(), z.clone(), f.line_file.clone()).into();
367                    if self
368                        .verify_non_equational_known_then_builtin_rules_only(&le, verify_state)?
369                        .is_true()
370                    {
371                        return success("order: (-1)*x >= 0 from x <= 0");
372                    }
373                    let lt: AtomicFact = LessFact::new(x, z.clone(), f.line_file.clone()).into();
374                    if self
375                        .verify_non_equational_known_then_builtin_rules_only(&lt, verify_state)?
376                        .is_true()
377                    {
378                        return success("order: (-1)*x >= 0 from x < 0");
379                    }
380                }
381                Ok(None)
382            }
383            AtomicFact::GreaterFact(f) if self.obj_is_resolved_zero(&f.right) => {
384                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.left) {
385                    let lt: AtomicFact = LessFact::new(x, z.clone(), f.line_file.clone()).into();
386                    if self
387                        .verify_non_equational_known_then_builtin_rules_only(&lt, verify_state)?
388                        .is_true()
389                    {
390                        return success("order: (-1)*x > 0 from x < 0");
391                    }
392                }
393                Ok(None)
394            }
395            AtomicFact::LessEqualFact(f) if self.obj_is_resolved_zero(&f.right) => {
396                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.left) {
397                    let ge: AtomicFact =
398                        GreaterEqualFact::new(x.clone(), z.clone(), f.line_file.clone()).into();
399                    if self
400                        .verify_non_equational_known_then_builtin_rules_only(&ge, verify_state)?
401                        .is_true()
402                    {
403                        return success("order: (-1)*x <= 0 from x >= 0");
404                    }
405                    let gt: AtomicFact = GreaterFact::new(x, z.clone(), f.line_file.clone()).into();
406                    if self
407                        .verify_non_equational_known_then_builtin_rules_only(&gt, verify_state)?
408                        .is_true()
409                    {
410                        return success("order: (-1)*x <= 0 from x > 0");
411                    }
412                }
413                Ok(None)
414            }
415            AtomicFact::LessFact(f) if self.obj_is_resolved_zero(&f.right) => {
416                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.left) {
417                    let gt: AtomicFact = GreaterFact::new(x, z.clone(), f.line_file.clone()).into();
418                    if self
419                        .verify_non_equational_known_then_builtin_rules_only(&gt, verify_state)?
420                        .is_true()
421                    {
422                        return success("order: (-1)*x < 0 from x > 0");
423                    }
424                }
425                Ok(None)
426            }
427            AtomicFact::LessEqualFact(f) if self.obj_is_resolved_zero(&f.left) => {
428                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.right) {
429                    let le: AtomicFact =
430                        LessEqualFact::new(x.clone(), z.clone(), f.line_file.clone()).into();
431                    if self
432                        .verify_non_equational_known_then_builtin_rules_only(&le, verify_state)?
433                        .is_true()
434                    {
435                        return success("order: 0 <= (-1)*x from x <= 0");
436                    }
437                    let lt: AtomicFact = LessFact::new(x, z.clone(), f.line_file.clone()).into();
438                    if self
439                        .verify_non_equational_known_then_builtin_rules_only(&lt, verify_state)?
440                        .is_true()
441                    {
442                        return success("order: 0 <= (-1)*x from x < 0");
443                    }
444                }
445                Ok(None)
446            }
447            AtomicFact::LessFact(f) if self.obj_is_resolved_zero(&f.left) => {
448                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.right) {
449                    let lt: AtomicFact = LessFact::new(x, z.clone(), f.line_file.clone()).into();
450                    if self
451                        .verify_non_equational_known_then_builtin_rules_only(&lt, verify_state)?
452                        .is_true()
453                    {
454                        return success("order: 0 < (-1)*x from x < 0");
455                    }
456                }
457                Ok(None)
458            }
459            AtomicFact::GreaterEqualFact(f) if self.obj_is_resolved_zero(&f.left) => {
460                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.right) {
461                    let ge: AtomicFact =
462                        GreaterEqualFact::new(x.clone(), z.clone(), f.line_file.clone()).into();
463                    if self
464                        .verify_non_equational_known_then_builtin_rules_only(&ge, verify_state)?
465                        .is_true()
466                    {
467                        return success("order: 0 >= (-1)*x from x >= 0");
468                    }
469                    let gt: AtomicFact = GreaterFact::new(x, z.clone(), f.line_file.clone()).into();
470                    if self
471                        .verify_non_equational_known_then_builtin_rules_only(&gt, verify_state)?
472                        .is_true()
473                    {
474                        return success("order: 0 >= (-1)*x from x > 0");
475                    }
476                }
477                Ok(None)
478            }
479            AtomicFact::GreaterFact(f) if self.obj_is_resolved_zero(&f.left) => {
480                if let Some(x) = self.peel_mul_by_literal_neg_one(&f.right) {
481                    let gt: AtomicFact = GreaterFact::new(x, z.clone(), f.line_file.clone()).into();
482                    if self
483                        .verify_non_equational_known_then_builtin_rules_only(&gt, verify_state)?
484                        .is_true()
485                    {
486                        return success("order: 0 > (-1)*x from x > 0");
487                    }
488                }
489                Ok(None)
490            }
491            _ => Ok(None),
492        }
493    }
494
495    // Lit `know` facts for the nonnegative / positive cone under field operations used to live in
496    // `BUILTIN_ENV_CODE_FOR_FUNDAMENTAL_COMPARISON` (`fundamental_comparison.rs`). Those fragments
497    // were removed as redundant; the same mathematics is checked here on normalized `0 <=` / `0 <`
498    // goals (possibly after `normalize_positive_order_atomic_fact`):
499    // - Chained `+`: `0 <= a + b + …` from `0 <=` on each peeled summand; `0 < a + b + …` from
500    //   `(0 < a ∧ 0 <= b) ∨ (0 <= a ∧ 0 < b)` at each binary `+`.
501    // - Powers: literal even integer exponent ⇒ `0 <= base^n`; literal integer exponent and `0 <= base`
502    //   (or `0 < base` if exponent < 0) ⇒ `0 <= base^n`; `a * a` with equal factors; `0 < base^exp`
503    //   from `0 < base` and `exp in R`.
504    // - Products and quotients: `0 <= a * b`, `0 < a * b`, `0 <= a / b` (denominator strictly
505    //   positive), `0 < a / b`, each with recursive sub-goals on operands.
506    // The Lit environment still records order via differences (`a <= b` iff `0 <= b - a`, etc.) and
507    // `a != 0 ⇒ 0 < a^2` (strict square). This path also bridges `0 <= u - v` / `0 < u - v` to `v <= u` / `v < u`.
508    // Algebraic closure (+, -, *, /) on general `a <= b` / `a < b` is in `order_algebra_builtin.rs`.
509    pub fn verify_order_atomic_fact_numeric_builtin_only(
510        &mut self,
511        atomic_fact: &AtomicFact,
512    ) -> Result<StmtResult, RuntimeError> {
513        let vs = VerifyState::new(0, true);
514        if let Some(result) =
515            self.try_verify_order_nonnegative_from_membership_in_n(atomic_fact, &vs)?
516        {
517            return Ok(result);
518        }
519        if let Some(result) =
520            self.try_verify_order_one_le_from_membership_in_n_pos(atomic_fact, &vs)?
521        {
522            return Ok(result);
523        }
524        if let Some(result) = self.try_verify_order_opposite_sign_mul_minus_one(atomic_fact, &vs)? {
525            return Ok(result);
526        }
527        if let Some(result) = self.verify_order_from_known_negated_complement(atomic_fact)? {
528            return Ok(result);
529        }
530        if let Some(result) = self.verify_negated_order_from_known_equivalent_order(atomic_fact)? {
531            return Ok(result);
532        }
533        if let Some(result) = self.verify_order_algebra_structural_builtin_rule(atomic_fact)? {
534            return Ok(result);
535        }
536        if let Some(result) = self.verify_zero_le_abs_builtin_rule(atomic_fact)? {
537            return Ok(result);
538        }
539        if let Some(result) =
540            self.verify_zero_order_on_sub_from_two_sided_order_builtin_rule(atomic_fact)?
541        {
542            return Ok(result);
543        }
544        if let Some(result) =
545            self.verify_zero_le_add_from_known_atomic_facts_builtin_rule(atomic_fact)?
546        {
547            return Ok(result);
548        }
549        if let Some(result) =
550            self.verify_zero_lt_add_from_known_atomic_facts_builtin_rule(atomic_fact)?
551        {
552            return Ok(result);
553        }
554        if let Some(result) = self.verify_zero_le_even_integer_pow_builtin_rule(atomic_fact)? {
555            return Ok(result);
556        }
557        if let Some(result) =
558            self.verify_zero_lt_even_integer_pow_from_base_nonzero_builtin_rule(atomic_fact)?
559        {
560            return Ok(result);
561        }
562        if let Some(result) =
563            self.verify_zero_lt_pow_from_positive_base_real_exp_builtin_rule(atomic_fact)?
564        {
565            return Ok(result);
566        }
567        if let Some(result) =
568            self.verify_zero_le_pow_from_positive_base_real_exp_builtin_rule(atomic_fact)?
569        {
570            return Ok(result);
571        }
572        if let Some(result) =
573            self.verify_zero_le_pow_integer_exponent_from_nonneg_base_builtin_rule(atomic_fact)?
574        {
575            return Ok(result);
576        }
577        if let Some(result) =
578            self.verify_zero_le_mul_from_known_atomic_facts_builtin_rule(atomic_fact)?
579        {
580            return Ok(result);
581        }
582        if let Some(result) =
583            self.verify_zero_lt_mul_from_known_atomic_facts_builtin_rule(atomic_fact)?
584        {
585            return Ok(result);
586        }
587        if let Some(result) =
588            self.verify_zero_le_div_from_known_atomic_facts_builtin_rule(atomic_fact)?
589        {
590            return Ok(result);
591        }
592        if let Some(result) =
593            self.verify_zero_lt_div_from_known_atomic_facts_builtin_rule(atomic_fact)?
594        {
595            return Ok(result);
596        }
597
598        if let AtomicFact::LessEqualFact(less_equal_fact) = atomic_fact {
599            if less_equal_fact.left.to_string() == less_equal_fact.right.to_string() {
600                return Ok(StmtResult::FactualStmtSuccess(
601                    FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
602                        less_equal_fact.clone().into(),
603                        "less_equal_fact_equal".to_string(),
604                        Vec::new(),
605                    ),
606                ));
607            }
608            let strict_fact: Fact = LessFact::new(
609                less_equal_fact.left.clone(),
610                less_equal_fact.right.clone(),
611                less_equal_fact.line_file.clone(),
612            )
613            .into();
614            let strict_key = strict_fact.to_string();
615            let (cache_ok, cache_line_file) = self.cache_known_facts_contains(&strict_key);
616            if cache_ok {
617                return Ok(StmtResult::FactualStmtSuccess(
618                    FactualStmtSuccess::new_with_verified_by_known_fact_source_recording_facts(
619                        less_equal_fact.clone().into(),
620                        strict_key,
621                        Some(strict_fact),
622                        Some(cache_line_file),
623                        Vec::new(),
624                    ),
625                ));
626            }
627        }
628        if let AtomicFact::GreaterEqualFact(greater_equal_fact) = atomic_fact {
629            if greater_equal_fact.left.to_string() == greater_equal_fact.right.to_string() {
630                return Ok(StmtResult::FactualStmtSuccess(
631                    FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
632                        greater_equal_fact.clone().into(),
633                        "greater_equal_fact_equal".to_string(),
634                        Vec::new(),
635                    ),
636                ));
637            }
638        }
639        if let Some(true) = self.verify_number_comparison_builtin_rule(atomic_fact) {
640            Ok(StmtResult::FactualStmtSuccess(
641                FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
642                    atomic_fact.clone().into(),
643                    "number comparison".to_string(),
644                    Vec::new(),
645                ),
646            ))
647        } else {
648            Ok(StmtResult::StmtUnknown(StmtUnknown::new()))
649        }
650    }
651
652    // `a > b` from known `not (a <= b)`, `a < b` from `not (a >= b)`, etc. (total order duality).
653    fn verify_order_from_known_negated_complement(
654        &mut self,
655        atomic_fact: &AtomicFact,
656    ) -> Result<Option<StmtResult>, RuntimeError> {
657        let negated: Option<AtomicFact> = match atomic_fact {
658            AtomicFact::GreaterFact(f) => {
659                Some(NotLessEqualFact::new(f.left.clone(), f.right.clone(), f.line_file.clone()).into())
660            }
661            AtomicFact::LessFact(f) => Some(
662                NotGreaterEqualFact::new(f.left.clone(), f.right.clone(), f.line_file.clone()).into(),
663            ),
664            AtomicFact::GreaterEqualFact(f) => {
665                Some(NotLessFact::new(f.left.clone(), f.right.clone(), f.line_file.clone()).into())
666            }
667            AtomicFact::LessEqualFact(f) => {
668                Some(NotGreaterFact::new(f.left.clone(), f.right.clone(), f.line_file.clone()).into())
669            }
670            _ => None,
671        };
672        let Some(neg) = negated else {
673            return Ok(None);
674        };
675        let sub = self.verify_non_equational_atomic_fact_with_known_atomic_facts(&neg)?;
676        if sub.is_true() {
677            return Ok(Some(
678                FactualStmtSuccess::new_with_verified_by_builtin_rules(
679                    atomic_fact.clone().into(),
680                    InferResult::new(),
681                    "order_from_known_negated_complement".to_string(),
682                    vec![sub],
683                )
684                .into(),
685            ));
686        }
687        Ok(None)
688    }
689
690    // `not (a < b)` etc.: only consult known atomic facts for the equivalent weak/strict order.
691    fn verify_negated_order_from_known_equivalent_order(
692        &mut self,
693        atomic_fact: &AtomicFact,
694    ) -> Result<Option<StmtResult>, RuntimeError> {
695        let candidates: Vec<AtomicFact> = match atomic_fact {
696            AtomicFact::NotLessFact(f) => {
697                let lf = f.line_file.clone();
698                vec![
699                    LessEqualFact::new(f.right.clone(), f.left.clone(), lf.clone()).into(),
700                    GreaterEqualFact::new(f.left.clone(), f.right.clone(), lf).into(),
701                ]
702            }
703            AtomicFact::NotGreaterFact(f) => {
704                let lf = f.line_file.clone();
705                vec![
706                    LessEqualFact::new(f.left.clone(), f.right.clone(), lf.clone()).into(),
707                    GreaterEqualFact::new(f.right.clone(), f.left.clone(), lf).into(),
708                ]
709            }
710            AtomicFact::NotLessEqualFact(f) => {
711                let lf = f.line_file.clone();
712                vec![
713                    LessFact::new(f.right.clone(), f.left.clone(), lf.clone()).into(),
714                    GreaterFact::new(f.left.clone(), f.right.clone(), lf).into(),
715                ]
716            }
717            AtomicFact::NotGreaterEqualFact(f) => {
718                let lf = f.line_file.clone();
719                vec![
720                    LessFact::new(f.left.clone(), f.right.clone(), lf.clone()).into(),
721                    GreaterFact::new(f.right.clone(), f.left.clone(), lf).into(),
722                ]
723            }
724            _ => return Ok(None),
725        };
726        for candidate in candidates {
727            let sub = self.verify_non_equational_atomic_fact_with_known_atomic_facts(&candidate)?;
728            if sub.is_true() {
729                return Ok(Some(
730                    FactualStmtSuccess::new_with_verified_by_builtin_rules(
731                        atomic_fact.clone().into(),
732                        InferResult::new(),
733                        "negated_order_from_known_equivalent_order".to_string(),
734                        vec![sub],
735                    )
736                    .into(),
737                ));
738            }
739        }
740        Ok(None)
741    }
742
743    // Matches Lit `a <= b` <=> `0 <= b - a` (and strict): `0 <= u - v` iff `v <= u`, `0 < u - v` iff `v < u`.
744    fn verify_zero_order_on_sub_from_two_sided_order_builtin_rule(
745        &mut self,
746        atomic_fact: &AtomicFact,
747    ) -> Result<Option<StmtResult>, RuntimeError> {
748        let Some(norm) = normalize_positive_order_atomic_fact(atomic_fact) else {
749            return Ok(None);
750        };
751        match &norm {
752            AtomicFact::LessEqualFact(f) if f.left.to_string() == "0" => {
753                let Obj::Sub(sub) = &f.right else {
754                    return Ok(None);
755                };
756                let derived: AtomicFact = LessEqualFact::new(
757                    sub.right.as_ref().clone(),
758                    sub.left.as_ref().clone(),
759                    f.line_file.clone(),
760                )
761                .into();
762                let mut result =
763                    self.verify_non_equational_atomic_fact_with_known_atomic_facts(&derived)?;
764                if !result.is_true() {
765                    result = self.verify_order_atomic_fact_numeric_builtin_only(&derived)?;
766                }
767                if result.is_true() {
768                    Ok(Some(StmtResult::FactualStmtSuccess(
769                        FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
770                            atomic_fact.clone().into(),
771                            "0 <= u - v from v <= u".to_string(),
772                            vec![result],
773                        ),
774                    )))
775                } else {
776                    Ok(None)
777                }
778            }
779            AtomicFact::LessFact(f) if f.left.to_string() == "0" => {
780                let Obj::Sub(sub) = &f.right else {
781                    return Ok(None);
782                };
783                let derived: AtomicFact = LessFact::new(
784                    sub.right.as_ref().clone(),
785                    sub.left.as_ref().clone(),
786                    f.line_file.clone(),
787                )
788                .into();
789                let mut result =
790                    self.verify_non_equational_atomic_fact_with_known_atomic_facts(&derived)?;
791                if !result.is_true() {
792                    result = self.verify_order_atomic_fact_numeric_builtin_only(&derived)?;
793                }
794                if result.is_true() {
795                    Ok(Some(StmtResult::FactualStmtSuccess(
796                        FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
797                            atomic_fact.clone().into(),
798                            "0 < u - v from v < u".to_string(),
799                            vec![result],
800                        ),
801                    )))
802                } else {
803                    Ok(None)
804                }
805            }
806            _ => Ok(None),
807        }
808    }
809
810    fn verify_zero_le_add_from_known_atomic_facts_builtin_rule(
811        &mut self,
812        atomic_fact: &AtomicFact,
813    ) -> Result<Option<StmtResult>, RuntimeError> {
814        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
815            return Ok(None);
816        };
817        let AtomicFact::LessEqualFact(less_equal_fact) = normalized_fact else {
818            return Ok(None);
819        };
820        if less_equal_fact.left.to_string() != "0" {
821            return Ok(None);
822        }
823        let Obj::Add(add_obj) = &less_equal_fact.right else {
824            return Ok(None);
825        };
826
827        let zero = &less_equal_fact.left;
828        let line_file = &less_equal_fact.line_file;
829        let left_verify_result =
830            self.verify_zero_order_on_sub_expr(zero, add_obj.left.as_ref(), true, line_file)?;
831        if !left_verify_result.is_true() {
832            return Ok(None);
833        }
834        let right_verify_result =
835            self.verify_zero_order_on_sub_expr(zero, add_obj.right.as_ref(), true, line_file)?;
836        if !right_verify_result.is_true() {
837            return Ok(None);
838        }
839
840        Ok(Some(StmtResult::FactualStmtSuccess(
841            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
842                atomic_fact.clone().into(),
843                "0 <= a + b from known atomic facts 0 <= a and 0 <= b".to_string(),
844                vec![left_verify_result, right_verify_result],
845            ),
846        )))
847    }
848
849    fn verify_zero_lt_add_from_known_atomic_facts_builtin_rule(
850        &mut self,
851        atomic_fact: &AtomicFact,
852    ) -> Result<Option<StmtResult>, RuntimeError> {
853        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
854            return Ok(None);
855        };
856        let AtomicFact::LessFact(less_fact) = normalized_fact else {
857            return Ok(None);
858        };
859        if less_fact.left.to_string() != "0" {
860            return Ok(None);
861        }
862        let Obj::Add(add_obj) = &less_fact.right else {
863            return Ok(None);
864        };
865
866        let zero = &less_fact.left;
867        let line_file = &less_fact.line_file;
868
869        let strict_then_weak = |this: &mut Self| -> Result<Option<StmtResult>, RuntimeError> {
870            let left_result =
871                this.verify_zero_order_on_sub_expr(zero, add_obj.left.as_ref(), false, line_file)?;
872            if !left_result.is_true() {
873                return Ok(None);
874            }
875            let right_result =
876                this.verify_zero_order_on_sub_expr(zero, add_obj.right.as_ref(), true, line_file)?;
877            if !right_result.is_true() {
878                return Ok(None);
879            }
880            Ok(Some(StmtResult::FactualStmtSuccess(
881                FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
882                    atomic_fact.clone().into(),
883                    "0 < a + b from (0 < a and 0 <= b)".to_string(),
884                    vec![left_result, right_result],
885                ),
886            )))
887        };
888        let weak_then_strict = |this: &mut Self| -> Result<Option<StmtResult>, RuntimeError> {
889            let left_result =
890                this.verify_zero_order_on_sub_expr(zero, add_obj.left.as_ref(), true, line_file)?;
891            if !left_result.is_true() {
892                return Ok(None);
893            }
894            let right_result =
895                this.verify_zero_order_on_sub_expr(zero, add_obj.right.as_ref(), false, line_file)?;
896            if !right_result.is_true() {
897                return Ok(None);
898            }
899            Ok(Some(StmtResult::FactualStmtSuccess(
900                FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
901                    atomic_fact.clone().into(),
902                    "0 < a + b from (0 <= a and 0 < b)".to_string(),
903                    vec![left_result, right_result],
904                ),
905            )))
906        };
907
908        if let Some(success) = strict_then_weak(self)? {
909            return Ok(Some(success));
910        }
911        weak_then_strict(self)
912    }
913
914    fn verify_zero_le_even_integer_pow_builtin_rule(
915        &mut self,
916        atomic_fact: &AtomicFact,
917    ) -> Result<Option<StmtResult>, RuntimeError> {
918        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
919            return Ok(None);
920        };
921        let AtomicFact::LessEqualFact(less_equal_fact) = normalized_fact else {
922            return Ok(None);
923        };
924        if less_equal_fact.left.to_string() != "0" {
925            return Ok(None);
926        }
927        let right = &less_equal_fact.right;
928        let is_equal_factors_mul = match right {
929            Obj::Mul(mul_obj) => mul_obj.left.to_string() == mul_obj.right.to_string(),
930            _ => false,
931        };
932        let is_even_pow = match right {
933            Obj::Pow(pow_obj) => match pow_obj.exponent.as_ref() {
934                Obj::Number(n) => normalized_decimal_string_is_even_integer(&n.normalized_value),
935                _ => false,
936            },
937            _ => false,
938        };
939        if !is_equal_factors_mul && !is_even_pow {
940            return Ok(None);
941        }
942        let msg = if is_equal_factors_mul {
943            "0 <= a * a from even integer exponent (here 2) (forall a R)".to_string()
944        } else {
945            "0 <= a^n for even integer n (forall a R)".to_string()
946        };
947        Ok(Some(StmtResult::FactualStmtSuccess(
948            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
949                atomic_fact.clone().into(),
950                msg,
951                Vec::new(),
952            ),
953        )))
954    }
955
956    fn verify_zero_lt_even_integer_pow_from_base_nonzero_builtin_rule(
957        &mut self,
958        atomic_fact: &AtomicFact,
959    ) -> Result<Option<StmtResult>, RuntimeError> {
960        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
961            return Ok(None);
962        };
963        let AtomicFact::LessFact(less_fact) = normalized_fact else {
964            return Ok(None);
965        };
966        if less_fact.left.to_string() != "0" {
967            return Ok(None);
968        }
969        let Obj::Pow(pow_obj) = &less_fact.right else {
970            return Ok(None);
971        };
972        let Obj::Number(exp_num) = pow_obj.exponent.as_ref() else {
973            return Ok(None);
974        };
975        if !normalized_decimal_string_is_even_integer(&exp_num.normalized_value) {
976            return Ok(None);
977        }
978
979        let line_file = less_fact.line_file.clone();
980        let base = pow_obj.base.as_ref().clone();
981        let zero_obj: Obj = Number::new("0".to_string()).into();
982        let base_neq_zero: AtomicFact = NotEqualFact::new(base, zero_obj, line_file.clone()).into();
983
984        let neq_result = self.verify_non_equational_known_then_builtin_rules_only(
985            &base_neq_zero,
986            &VerifyState::new(0, true),
987        )?;
988        if !neq_result.is_true() {
989            return Ok(None);
990        }
991
992        Ok(Some(StmtResult::FactualStmtSuccess(
993            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
994                atomic_fact.clone().into(),
995                "0 < a^n for even integer n from a != 0".to_string(),
996                vec![neq_result],
997            ),
998        )))
999    }
1000
1001    // Matches `0 < a^b` / `a^b > 0` when `0 < a` is proved (or known) and `b in R`.
1002    fn verify_zero_lt_pow_from_positive_base_real_exp_builtin_rule(
1003        &mut self,
1004        atomic_fact: &AtomicFact,
1005    ) -> Result<Option<StmtResult>, RuntimeError> {
1006        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
1007            return Ok(None);
1008        };
1009        let AtomicFact::LessFact(less_fact) = normalized_fact else {
1010            return Ok(None);
1011        };
1012        if less_fact.left.to_string() != "0" {
1013            return Ok(None);
1014        }
1015        let Obj::Pow(pow_obj) = &less_fact.right else {
1016            return Ok(None);
1017        };
1018        let zero = &less_fact.left;
1019        let line_file = &less_fact.line_file;
1020        let base = pow_obj.base.as_ref();
1021        let base_result = self.verify_zero_order_on_sub_expr(zero, base, false, line_file)?;
1022        if !base_result.is_true() {
1023            return Ok(None);
1024        }
1025        let in_r: AtomicFact = InFact::new(
1026            (*pow_obj.exponent).clone(),
1027            StandardSet::R.into(),
1028            line_file.clone(),
1029        )
1030        .into();
1031        let in_r_result = self.verify_non_equational_known_then_builtin_rules_only(
1032            &in_r,
1033            &VerifyState::new(0, true),
1034        )?;
1035        if !in_r_result.is_true() {
1036            return Ok(None);
1037        }
1038        Ok(Some(StmtResult::FactualStmtSuccess(
1039            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
1040                atomic_fact.clone().into(),
1041                "0 < a^b from 0 < a and b in R".to_string(),
1042                vec![base_result, in_r_result],
1043            ),
1044        )))
1045    }
1046
1047    // `0 <= a^b` / `a^b >= 0` with the same premises as strict `0 < a^b`: `0 < a` and `b in R`.
1048    // Covers symbolic exponents (e.g. `2^m`) where the literal-exponent `0 <= a^n` rule does not apply.
1049    fn verify_zero_le_pow_from_positive_base_real_exp_builtin_rule(
1050        &mut self,
1051        atomic_fact: &AtomicFact,
1052    ) -> Result<Option<StmtResult>, RuntimeError> {
1053        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
1054            return Ok(None);
1055        };
1056        let AtomicFact::LessEqualFact(less_equal_fact) = normalized_fact else {
1057            return Ok(None);
1058        };
1059        if less_equal_fact.left.to_string() != "0" {
1060            return Ok(None);
1061        }
1062        let Obj::Pow(pow_obj) = &less_equal_fact.right else {
1063            return Ok(None);
1064        };
1065        let zero = &less_equal_fact.left;
1066        let line_file = &less_equal_fact.line_file;
1067        let base = pow_obj.base.as_ref();
1068        let base_result = self.verify_zero_order_on_sub_expr(zero, base, false, line_file)?;
1069        if !base_result.is_true() {
1070            return Ok(None);
1071        }
1072        let in_r: AtomicFact = InFact::new(
1073            (*pow_obj.exponent).clone(),
1074            StandardSet::R.into(),
1075            line_file.clone(),
1076        )
1077        .into();
1078        let in_r_result = self.verify_non_equational_known_then_builtin_rules_only(
1079            &in_r,
1080            &VerifyState::new(0, true),
1081        )?;
1082        if !in_r_result.is_true() {
1083            return Ok(None);
1084        }
1085        Ok(Some(StmtResult::FactualStmtSuccess(
1086            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
1087                atomic_fact.clone().into(),
1088                "0 <= a^b from 0 < a and b in R".to_string(),
1089                vec![base_result, in_r_result],
1090            ),
1091        )))
1092    }
1093
1094    fn verify_zero_le_pow_integer_exponent_from_nonneg_base_builtin_rule(
1095        &mut self,
1096        atomic_fact: &AtomicFact,
1097    ) -> Result<Option<StmtResult>, RuntimeError> {
1098        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
1099            return Ok(None);
1100        };
1101        let AtomicFact::LessEqualFact(less_equal_fact) = normalized_fact else {
1102            return Ok(None);
1103        };
1104        if less_equal_fact.left.to_string() != "0" {
1105            return Ok(None);
1106        }
1107        let Obj::Pow(pow_obj) = &less_equal_fact.right else {
1108            return Ok(None);
1109        };
1110        let Obj::Number(exp_num) = pow_obj.exponent.as_ref() else {
1111            return Ok(None);
1112        };
1113        if !normalized_decimal_string_is_integer(&exp_num.normalized_value) {
1114            return Ok(None);
1115        }
1116
1117        let zero = &less_equal_fact.left;
1118        let line_file = &less_equal_fact.line_file;
1119        let base = pow_obj.base.as_ref();
1120
1121        let exponent_vs_zero = compare_normalized_number_str_to_zero(&exp_num.normalized_value);
1122        let base_result = match exponent_vs_zero {
1123            NumberCompareResult::Less => {
1124                self.verify_zero_order_on_sub_expr(zero, base, false, line_file)?
1125            }
1126            NumberCompareResult::Equal | NumberCompareResult::Greater => {
1127                self.verify_zero_order_on_sub_expr(zero, base, true, line_file)?
1128            }
1129        };
1130        if !base_result.is_true() {
1131            return Ok(None);
1132        }
1133
1134        let msg = match exponent_vs_zero {
1135            NumberCompareResult::Less => "0 <= a^n from 0 < a and integer n < 0".to_string(),
1136            _ => "0 <= a^n from 0 <= a and integer n".to_string(),
1137        };
1138
1139        Ok(Some(StmtResult::FactualStmtSuccess(
1140            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
1141                atomic_fact.clone().into(),
1142                msg,
1143                vec![base_result],
1144            ),
1145        )))
1146    }
1147
1148    fn verify_zero_le_mul_from_known_atomic_facts_builtin_rule(
1149        &mut self,
1150        atomic_fact: &AtomicFact,
1151    ) -> Result<Option<StmtResult>, RuntimeError> {
1152        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
1153            return Ok(None);
1154        };
1155        let AtomicFact::LessEqualFact(less_equal_fact) = normalized_fact else {
1156            return Ok(None);
1157        };
1158        if less_equal_fact.left.to_string() != "0" {
1159            return Ok(None);
1160        }
1161        let Obj::Mul(mul_obj) = &less_equal_fact.right else {
1162            return Ok(None);
1163        };
1164
1165        let zero = &less_equal_fact.left;
1166        let line_file = &less_equal_fact.line_file;
1167        let left_verify_result =
1168            self.verify_zero_order_on_sub_expr(zero, mul_obj.left.as_ref(), true, line_file)?;
1169        if !left_verify_result.is_true() {
1170            return Ok(None);
1171        }
1172        let right_verify_result =
1173            self.verify_zero_order_on_sub_expr(zero, mul_obj.right.as_ref(), true, line_file)?;
1174        if !right_verify_result.is_true() {
1175            return Ok(None);
1176        }
1177
1178        Ok(Some(StmtResult::FactualStmtSuccess(
1179            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
1180                atomic_fact.clone().into(),
1181                "0 <= a * b from 0 <= a and 0 <= b".to_string(),
1182                vec![left_verify_result, right_verify_result],
1183            ),
1184        )))
1185    }
1186
1187    fn verify_zero_lt_mul_from_known_atomic_facts_builtin_rule(
1188        &mut self,
1189        atomic_fact: &AtomicFact,
1190    ) -> Result<Option<StmtResult>, RuntimeError> {
1191        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
1192            return Ok(None);
1193        };
1194        let AtomicFact::LessFact(less_fact) = normalized_fact else {
1195            return Ok(None);
1196        };
1197        if less_fact.left.to_string() != "0" {
1198            return Ok(None);
1199        }
1200        let Obj::Mul(mul_obj) = &less_fact.right else {
1201            return Ok(None);
1202        };
1203
1204        let zero = &less_fact.left;
1205        let line_file = &less_fact.line_file;
1206        let left_verify_result =
1207            self.verify_zero_order_on_sub_expr(zero, mul_obj.left.as_ref(), false, line_file)?;
1208        if !left_verify_result.is_true() {
1209            return Ok(None);
1210        }
1211        let right_verify_result =
1212            self.verify_zero_order_on_sub_expr(zero, mul_obj.right.as_ref(), false, line_file)?;
1213        if !right_verify_result.is_true() {
1214            return Ok(None);
1215        }
1216
1217        Ok(Some(StmtResult::FactualStmtSuccess(
1218            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
1219                atomic_fact.clone().into(),
1220                "0 < a * b from 0 < a and 0 < b".to_string(),
1221                vec![left_verify_result, right_verify_result],
1222            ),
1223        )))
1224    }
1225
1226    fn verify_zero_le_div_from_known_atomic_facts_builtin_rule(
1227        &mut self,
1228        atomic_fact: &AtomicFact,
1229    ) -> Result<Option<StmtResult>, RuntimeError> {
1230        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
1231            return Ok(None);
1232        };
1233        let AtomicFact::LessEqualFact(less_equal_fact) = normalized_fact else {
1234            return Ok(None);
1235        };
1236        if less_equal_fact.left.to_string() != "0" {
1237            return Ok(None);
1238        }
1239        let Obj::Div(div_obj) = &less_equal_fact.right else {
1240            return Ok(None);
1241        };
1242
1243        let zero = &less_equal_fact.left;
1244        let line_file = &less_equal_fact.line_file;
1245        let numer_result =
1246            self.verify_zero_order_on_sub_expr(zero, div_obj.left.as_ref(), true, line_file)?;
1247        if !numer_result.is_true() {
1248            return Ok(None);
1249        }
1250        let denom_result =
1251            self.verify_zero_order_on_sub_expr(zero, div_obj.right.as_ref(), false, line_file)?;
1252        if !denom_result.is_true() {
1253            return Ok(None);
1254        }
1255
1256        Ok(Some(StmtResult::FactualStmtSuccess(
1257            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
1258                atomic_fact.clone().into(),
1259                "0 <= a / b from 0 <= a and 0 < b".to_string(),
1260                vec![numer_result, denom_result],
1261            ),
1262        )))
1263    }
1264
1265    fn verify_zero_lt_div_from_known_atomic_facts_builtin_rule(
1266        &mut self,
1267        atomic_fact: &AtomicFact,
1268    ) -> Result<Option<StmtResult>, RuntimeError> {
1269        let Some(normalized_fact) = normalize_positive_order_atomic_fact(atomic_fact) else {
1270            return Ok(None);
1271        };
1272        let AtomicFact::LessFact(less_fact) = normalized_fact else {
1273            return Ok(None);
1274        };
1275        if less_fact.left.to_string() != "0" {
1276            return Ok(None);
1277        }
1278        let Obj::Div(div_obj) = &less_fact.right else {
1279            return Ok(None);
1280        };
1281
1282        let zero = &less_fact.left;
1283        let line_file = &less_fact.line_file;
1284        let numer_result =
1285            self.verify_zero_order_on_sub_expr(zero, div_obj.left.as_ref(), false, line_file)?;
1286        if !numer_result.is_true() {
1287            return Ok(None);
1288        }
1289        let denom_result =
1290            self.verify_zero_order_on_sub_expr(zero, div_obj.right.as_ref(), false, line_file)?;
1291        if !denom_result.is_true() {
1292            return Ok(None);
1293        }
1294
1295        Ok(Some(StmtResult::FactualStmtSuccess(
1296            FactualStmtSuccess::new_with_verified_by_builtin_rules_recording_stmt(
1297                atomic_fact.clone().into(),
1298                "0 < a / b from 0 < a and 0 < b".to_string(),
1299                vec![numer_result, denom_result],
1300            ),
1301        )))
1302    }
1303
1304    pub fn verify_number_comparison_builtin_rule(&self, atomic_fact: &AtomicFact) -> Option<bool> {
1305        let normalized = normalize_positive_order_atomic_fact(atomic_fact)?;
1306        match normalized {
1307            AtomicFact::LessFact(less_fact) => {
1308                if let Some(calculated_number_string_pair) =
1309                    self.calculate_obj_pair_to_number_strings(&less_fact.left, &less_fact.right)
1310                {
1311                    return Some(matches!(
1312                        compare_number_strings(
1313                            &calculated_number_string_pair.0,
1314                            &calculated_number_string_pair.1
1315                        ),
1316                        NumberCompareResult::Less
1317                    ));
1318                }
1319                self.try_verify_numeric_order_via_div_elimination(
1320                    &less_fact.left,
1321                    &less_fact.right,
1322                    false,
1323                )
1324            }
1325            AtomicFact::LessEqualFact(less_equal_fact) => {
1326                if let Some(calculated_number_string_pair) = self
1327                    .calculate_obj_pair_to_number_strings(
1328                        &less_equal_fact.left,
1329                        &less_equal_fact.right,
1330                    )
1331                {
1332                    let compare_result = compare_number_strings(
1333                        &calculated_number_string_pair.0,
1334                        &calculated_number_string_pair.1,
1335                    );
1336                    return Some(matches!(
1337                        compare_result,
1338                        NumberCompareResult::Less | NumberCompareResult::Equal
1339                    ));
1340                }
1341                self.try_verify_numeric_order_via_div_elimination(
1342                    &less_equal_fact.left,
1343                    &less_equal_fact.right,
1344                    true,
1345                )
1346            }
1347            _ => None,
1348        }
1349    }
1350
1351    fn calculate_obj_pair_to_number_strings(
1352        &self,
1353        left_obj: &Obj,
1354        right_obj: &Obj,
1355    ) -> Option<(String, String)> {
1356        let left_number = self.resolve_obj_to_number_resolved(left_obj)?;
1357        let right_number = self.resolve_obj_to_number_resolved(right_obj)?;
1358        Some((left_number.normalized_value, right_number.normalized_value))
1359    }
1360}