Skip to main content

formula_check/
formula_check.rs

1//! 公式识别端到端检查:载入编码器和解码器(共用一次读权重),逐个读 PPM(P6,8 位 RGB)公式图,
2//! 打印 LaTeX、token 数和耗时。不给图时用随机 memory 跑一遍解码器(流程、NaN、停止)。
3//!     cargo run --profile release-fast -p pagesift-ocr --example formula_check -- <模型目录> [公式图.ppm]...
4use std::path::Path;
5use std::time::Instant;
6
7use pagesift_infer::{Tensor, Weights};
8use pagesift_ocr::formula::{Decoder, Tokenizer, to_latex};
9use pagesift_ocr::formula_enc::{Config, Encoder, preprocess};
10
11/// 最小的 P6 PPM 读取:头是 "P6 宽 高 255"(空白分隔),之后紧跟 RGB 数据。
12fn read_ppm(path: &str) -> pagesift_ocr::Image {
13    let b = std::fs::read(path).unwrap_or_else(|e| panic!("读 {path} 失败:{e}"));
14    let mut fields = Vec::new();
15    let mut i = 0;
16    while fields.len() < 4 {
17        while b[i].is_ascii_whitespace() {
18            i += 1;
19        }
20        let s = i;
21        while !b[i].is_ascii_whitespace() {
22            i += 1;
23        }
24        fields.push(String::from_utf8_lossy(&b[s..i]).into_owned());
25    }
26    assert_eq!(fields[0], "P6", "{path} 不是 P6 PPM");
27    let (w, h): (usize, usize) = (fields[1].parse().unwrap(), fields[2].parse().unwrap());
28    pagesift_ocr::Image {
29        width: w,
30        height: h,
31        channels: 3,
32        data: b[i + 1..i + 1 + w * h * 3].to_vec(),
33    }
34}
35
36fn json(dir: &Path, name: &str) -> serde_json::Value {
37    serde_json::from_slice(&std::fs::read(dir.join(name)).unwrap()).unwrap()
38}
39
40fn main() {
41    let a: Vec<String> = std::env::args().skip(1).collect();
42    let dir = Path::new(&a[0]);
43    let t = Instant::now();
44    let cfg = json(dir, "config.json");
45    let w = Weights::from_bytes(&std::fs::read(dir.join("model.safetensors")).unwrap()).unwrap();
46    let dec = Decoder::from_weights(
47        &w,
48        &cfg["text_config"],
49        &json(dir, "generation_config.json"),
50    )
51    .unwrap();
52    let enc = Encoder::from_weights(&w, Config::from_json(&cfg).unwrap()).unwrap();
53    drop(w);
54    let tok = Tokenizer::load(dir).unwrap();
55    println!("载入 {:.0} ms", t.elapsed().as_secs_f64() * 1e3);
56
57    if a.len() == 1 {
58        // 固定种子的随机 memory
59        let mut s = 0x9E37_79B9_7F4A_7C15u64;
60        let data = (0..144 * 512)
61            .map(|_| {
62                s ^= s << 13;
63                s ^= s >> 7;
64                s ^= s << 17;
65                (s >> 40) as f32 / (1u64 << 23) as f32 - 1.0
66            })
67            .collect();
68        let mem = Tensor::new(vec![144, 512], data);
69        let t = Instant::now();
70        let ids = dec.generate(&mem, 200).unwrap();
71        let ms = t.elapsed().as_secs_f64() * 1e3;
72        println!(
73            "随机 memory:{} 个 token,{ms:.0} ms,{:.1} ms/token",
74            ids.len(),
75            ms / (ids.len() + 1) as f64
76        );
77        println!("{}", to_latex(&tok.decode(&ids)));
78        return;
79    }
80
81    for p in &a[1..] {
82        let img = read_ppm(p);
83        let t = Instant::now();
84        let x = preprocess(&img);
85        let t_pre = t.elapsed().as_secs_f64() * 1e3;
86        let t = Instant::now();
87        let mem = enc.encode(&x).unwrap();
88        let t_enc = t.elapsed().as_secs_f64() * 1e3;
89        let t = Instant::now();
90        let ids = dec.generate(&mem, dec.max_length()).unwrap();
91        let t_dec = t.elapsed().as_secs_f64() * 1e3;
92        let name = Path::new(p).file_name().unwrap().to_string_lossy();
93        println!(
94            "{name}\t预处理 {t_pre:.0} ms\t编码 {t_enc:.0} ms\t解码 {t_dec:.0} ms({} 个 token,{:.1} ms/token)",
95            ids.len(),
96            t_dec / (ids.len() + 1) as f64
97        );
98        println!("  原始:{}", tok.decode(&ids));
99        println!("  LaTeX:{}", to_latex(&tok.decode(&ids)));
100    }
101    // 内存峰值(Linux 才有这个文件)
102    if let Ok(s) = std::fs::read_to_string("/proc/self/status")
103        && let Some(l) = s.lines().find(|l| l.starts_with("VmHWM"))
104    {
105        println!("{l}");
106    }
107}