lyj714 发表于 2022-6-27 11:51
三点以上平面需要拟合，tcl有数学库可以做最小二乘奇异值分解。

Tcl居然还有这种库，求一下库的名字，感激

三点以上平面需要拟合，tcl有数学库可以做最小二乘奇异值分解。

我还以为你要做某个平面的拟合（超过3个原子），那得做奇异值分解，用Tcl实现估计比较难。只是随口说一句，如果硬要做这个，可以参考http://bbs.keinsci.com/thread-23019-1-1.html，用C语言链接LAPACK做个SVD的接口。。。另外几十个GB的xyz也不至于大到不能接受，有个硬盘大点的服务器，多等一会儿还是能出结果的。

helix和平面没什么直接联系，不知道你说的是什么平面，得给个示意图
遇到脚本出不来结果的时候，就不要一下子把好多个语句嵌套在一起，应该拆开来一条一条执行。即每个[ ]子语句的执行作为一行，并且结果赋值到一个变量里，在脚本中把这个变量puts显示出来，这样自然就能分析到底哪个地方有问题。