TDEP25.03 CGAL三声子散射曲面操作例
先看图：bash preview.sh，3张实际图在preview/figures。
运行原生程序：bash prepare.sh → bash calculate.sh → bash extract.sh。
11个小型单进程CPU条件，使用module tdep/25.03-cgal-8173m；包括保留的无minus422曲面诊断，不会把空通道画成有效曲面。
本包只用自有二/三阶力常数和结构，不运行DFT/MD、不重新拟合；正式计算请使用Slurm分配兼容AVX512的CPU资源。
基础14程序和默认模块保持不变，可选模块单独增加phasespace_surface。
源25.03固定d38f435；ifx/icx/icpx2024.2.1、IntelMPI2021.13、顺序MKL、FFTW3.3.11、HDF51.14.6，另有CGAL5.6.1/Boost1.86/GMP6.3/MPFR4.2.2。
qgrid是CGAL楔形网格密度参数，不等于整布里渊区恰好n^3个点。曲面是四面体内线性插值零面，真实声子频率代回仍可能有残差；本例没有网格收敛或材料寿命结论。
POV坐标来自原生inverseBohr，独立提取转为cycles/inverseAngstrom；不同图分别标示实际坐标尺度。
psisq_norm是顶点|V|²数组的L2范数，不是曲面积分、散射率或寿命；不能把不同网格的这个数直接当收敛率对比。
--intensities会把POV梯度字段从单位插值法向改为v1±v2−v3；不是同一个量，曲面顶点则保持一致。
本版逐点强度存储有plus/minus数组混用，示例只读取原生标量范数，不输出这些含义有歧义的逐点数组。
结构/力常数来自本教程真实Si64 DFT300K短采样模型A，有限超胞、有限截断和采样误差仍存在。
程序与第三方源码/二进制、POTCAR/WAVECAR和完整电子输出均不随包分发。
