今天,我解決了CODEV中引入自定義間斷面的問題。積累了經(jīng)驗,與大家分享。 AdW2o|Uap
8{fz0H.<?
按照說明書,有間斷的面最好使用USERSUR2,但是,當(dāng)我打開他的模板后,發(fā)現(xiàn)比較復(fù)雜,這是一個專門為菲涅爾透鏡的生成而編制的程序,參數(shù)很多,復(fù)雜之處在于處理與各個間斷點有關(guān)的問題,具體我也沒太看懂,只知道是反復(fù)疊代找到交點。 50aWFJYw
><MgIV
但是,雖然codev的說明擋中強調(diào)USERSUR型子程序只適合于連續(xù)表面,但考慮到,光線追跡實際上只需要知道表面的函數(shù)和該點的導(dǎo)數(shù),所以,如果在程序中仔細(xì)定義間斷點處的導(dǎo)數(shù),就可以當(dāng)成連續(xù)面來處理。 >Mz|e(6
|K;Txe_
我的解決辦法是,利用有限差分計算導(dǎo)數(shù),而且是前后差分導(dǎo)數(shù)都計算出來,做對比,取絕對值小的那一個作為導(dǎo)數(shù)。因為,很容易理解,如果某個差分取在了間斷點的兩側(cè),一定很大。另一個較小的 5b`xN!c
ONq/JW$?LV
就是在間斷點同側(cè)取的。數(shù)學(xué)上來說,相當(dāng)于左右導(dǎo)數(shù)。 (+8xUc(w
#Rx"L&3Ue
計算表明,比較OK。 /)?P>!#;\
r&3o~!
USERSUR.for的使用比較簡單。只需把CODEV自帶的程序,即, USERSUR.for改成自己的程序名字,再把內(nèi)容修改就可以。重要的參數(shù)就是輸出導(dǎo)數(shù)。其他像錐常數(shù),曲率等都不用。雖然,允許使用MODE參數(shù)選擇用codev自己用有限差分計算導(dǎo)數(shù),但最好別這樣。應(yīng)按如上所述處理,比較安全。 Fg\| e%
^s~n[
此外,生成dll文件要使用nmake技術(shù)。注意環(huán)境變量要設(shè)對。此外,如果還出現(xiàn)找不到文件的問題。 KBRg95E~]l
*-AAQ
那就手動找文件,通常在fortran的安裝目錄的某些子目錄中,用windows的搜索功能就可找到,然后考備到當(dāng)前目錄下,通常是c:\CVuser. 其他問題在codev的幫助文檔中很容易找到。 u?3NBc$~A
T5jG IIa
在nmake語句中,要注意空格。例如,nmake -f *******