[原创]CODEV中引入自定义间断面 [复制链接]

今天，我解决了CODEV中引入自定义间断面的问题。积累了经验，与大家分享。 ln.kEhQ3B  
?7-#iC`  
  按照说明书，有间断的面最好使用USERSUR2，但是，当我打开他的模板后，发现比较复杂，这是一个专门为菲涅尔透镜的生成而编制的程序，参数很多，复杂之处在于处理与各个间断点有关的问题，具体我也没太看懂，只知道是反复叠代找到交点。 (^GVy=  
lJ]r%YlF  
  但是，虽然codev的说明挡中强调USERSUR型子程序只适合于连续表面，但考虑到，光线追迹实际上只需要知道表面的函数和该点的导数，所以，如果在程序中仔细定义间断点处的导数，就可以当成连续面来处理。 1"46OCu{  
&_FNDJ>MCk  
   我的解决办法是，利用有限差分计算导数，而且是前后差分导数都计算出来，做对比，取绝对值小的那一个作为导数。因为，很容易理解，如果某个差分取在了间断点的两侧，一定很大。另一个较小的 bb;fV  
PJj{5,#@3  
就是在间断点同侧取的。数学上来说，相当于左右导数。 ^B@4 w
\t  
xVoWGz7  
  计算表明，比较OK。 }@yvw*c  

Pg[zRRf<  
      USERSUR.for的使用比较简单。只需把CODEV自带的程序，即， USERSUR.for改成自己的程序名字，再把内容修改就可以。重要的参数就是输出导数。其他像锥常数，曲率等都不用。虽然，允许使用MODE参数选择用codev自己用有限差分计算导数，但最好别这样。应按如上所述处理，比较安全。 }yqRz6=YB  
20m6-rkI<}  
  此外，生成dll文件要使用nmake技术。注意环境变量要设对。此外，如果还出现找不到文件的问题。 sBuJK'  
mOwgk7s[J  
那就手动找文件，通常在fortran的安装目录的某些子目录中，用windows的搜索功能就可找到，然后考备到当前目录下，通常是c:\CVuser. 其他问题在codev的帮助文档中很容易找到。 S-Va_t$  
]i*q*]x2u  
在nmake语句中，要注意空格。例如，nmake  -f *******