從事有限元、特征值算法、高維問(wèn)題機(jī)器學(xué)習(xí)算法等研究。
里面求解Possion方程和特征值問(wèn)題的程序經(jīng)測(cè)試暫時(shí)沒(méi)什么問(wèn)題;
由于本人使用北太自帶的eigs不是很穩(wěn)定,自己寫了一個(gè)求解對(duì)稱正定特征值問(wèn)題的解法器GCG_B_Eigen,設(shè)計(jì)的思想是依據(jù)我們?cè)O(shè)計(jì)的特征值解法器GCGE;
由于里面有用直接方法求解特征值的部分,用戶可以將這部分注釋掉,那么剩下的就是多水平校正的方法求解特征值問(wèn)題了,用戶可以提升計(jì)算規(guī)模;
這是第一個(gè)測(cè)試版本,未來(lái)將繼續(xù)加入自適應(yīng)加密的功能,四邊形網(wǎng)格上有限元函數(shù)的畫圖,畫loglog的收斂速度圖像等;
供各位用來(lái)進(jìn)行有限元方法的學(xué)習(xí),也可以用于算法研究的簡(jiǎn)單數(shù)值測(cè)試,比如去研究求解半線性問(wèn)題、某些相場(chǎng)模型的多水平校正算法等,也可以提供邊值問(wèn)題幾何多重網(wǎng)格的功能;
哈哈,最主要的目的是歡迎大家使用多水平校正算法或者擴(kuò)展子空間算法來(lái)求解您的問(wèn)題或者支持您的研究;
下面是相應(yīng)的測(cè)試命令,供參考(一條一條執(zhí)行):
[N,Err,mesh,solu]=Possion; %求解邊值問(wèn)題 showmesh(mesh) %畫網(wǎng)格圖 showFeFunction(mesh,solu) %畫有限元解圖像 [lm,Err_0,Err_1,N]=eigen_GCG [LAM,LAM_Dir,Err,Rate,N]= Eigen_Multigrid_Method(3) [LAM,LAM_Dir,Dis_1,Dis_0,N]=eigen_Multigrid(3) [LL,LAM,LAM_Dir,Dis_1,Dis_0,Err,Err_dir,N]=Full_Eigen_Multigrid(3) [LAM,LAM_Dir,LL,N]=Full_Eigen_Multigrid_Multi(3)
下載地址:http://lsec.cc.ac.cn/~hhxie/Files/Full_Eigen_Multigrid.zip,更多相關(guān)的介紹材料可以參看http://lsec.cc.ac.cn/~hhxie網(wǎng)頁(yè)下的介紹材料和里面的視頻。
最近在北太天元的SDK環(huán)境下編譯了我們開發(fā)的特征值算法:廣義共軛梯度算法和相應(yīng)的特征值解法器GCGE。同時(shí)在北太天元環(huán)境下進(jìn)行了相應(yīng)的測(cè)試,矩陣規(guī)模從一萬(wàn)多到二十六萬(wàn)多,計(jì)算特征對(duì)個(gè)數(shù)從10到800個(gè)。為了體現(xiàn)計(jì)算的效率,我們這里的測(cè)試比較是直接與Matlab的eigs進(jìn)行了比較。
您可以直接從下面的地址下載文件:
下載之后可以直接使用文件 GCGE_BEX/test下已經(jīng)編譯好的gcge_bex.bexw64(windows版本)或者gcge_bexa64(Ubuntu版本)。
http://lsec.cc.ac.cn/~hhxie/Files/GCGE-BEX.zip
也可以利用bex文件重新編譯軟件包。
具體使用說(shuō)明和測(cè)試文檔可以從下面的鏈接下載
http://lsec.cc.ac.cn/~hhxie/Files/balt_gcge.pdf
如果有什么問(wèn)題,歡迎給我們發(fā)郵件:
王紫菁zjwang@lsec.cc.ac.cn 謝和虎hhxie@lsec.cc.ac.cn
歡迎各位使用!