[求助] 关于Zemax中自定义面元构建的问题

本人研究方向是超透镜相关，最近看到一篇论文，附件已给出，是在双曲相位的基础上指定聚焦点（xf，yf，zf），以实现聚焦点偏折的效果。超透镜的聚焦效果及其相位轮廓表达式如下图所示：




我想通过这个相位分布，通过自定义面元，在zemax中实现一个平行光入射，在指定坐标点聚焦的一个效果。我的自定义面元文件是在双曲相位的基础上修改的，附件将会给出这两部分代码。我修改的部分也比较简单，就是在case5中修改了相位轮廓trySurfDelta的基础上再重新计算相位梯度dpdr2。由于我只想要实现（xf，0，zf）的坐标点聚焦，因此我省去了yf，只定义xf和zf，并出于简单标注的目的，在代码中把它们分别设为A、B。

相位梯度部分的计算我比较纠结，不知道算的是否正确，在这里附上计算过程：

因为这边涉及到开根号的问题（偏x/偏r^2=sqrt（r^2-y^2)=+/-x),所以我将相位梯度的求解情况分成了九个部分，即圆形透镜上的原点、xy坐标轴与四个扇形区域。

按这样子去编写的dll导入zemax之后出现了一个问题，就是只有x方向的光线发生了偏折，y方向的光线依然只在y平面内，如下图所示。

这和我想象中的光线分布不太相同。我想象中的效果应该是这样的：

现在困惑的是，应该如何定义自定义面元，能让y轴上发出的光束也能与x轴发出的光束聚焦在同一点？还望各位同行能不吝赐教，给小弟指点迷津，不胜感激。

@既见花 您好，很冒昧圈您。我zemax自定义面元的编写是通过您发布的帖子学习的。不知道您是否能帮我解答一下上述的问题呢？

问题解决了  没必要对r2求导 直接算my、my就行了