编程实现光线追迹算法？

我已经做出了初级版本，可以实现球面系统的光线追迹，有修改F数功能、自动计算焦距、光瞳位置大小等等，是可以初步使用的，但是要增加非球面的光线追迹，目前看到的资料里面需要用一些迭代，没法一次性用解析式写出来光线在非球面上的向量信息，这样会大大减缓追迹速度，是个问题。关于题主的疑问，是有一些老书讲了追迹计算过程的，可以看一下，知网也有《一种新的旋转非球面光路追迹公式》，作者侯华

以前只会zemax,后来学了点codev，发现codev好像快一点，觉得好神奇，再后来接触了synopsys，觉得这东西又更离谱，某些功能直接碾压zemax和codev,那到底谁是最好的？为什么它们有这么大差别？是光线追迹算法差异还是优化算法差异？别人做啥软件我们就跟着用啥软件，永远在别人的屁股后面？我不太甘心，尝试着自己做一做，哪怕知道以卵击石我也不后悔，也许这就是意义。