里面求解Possion方程和特征值問題的程序經(jīng)測試暫時沒什么問題;
由于本人使用北太自帶的eigs不是很穩(wěn)定,自己寫了一個求解對稱正定特征值問題的解法器GCG_B_Eigen,設(shè)計的思想是依據(jù)我們設(shè)計的特征值解法器GCGE;
由于里面有用直接方法求解特征值的部分,用戶可以將這部分注釋掉,那么剩下的就是多水平校正的方法求解特征值問題了,用戶可以提升計算規(guī)模;
這是第一個測試版本,未來將繼續(xù)加入自適應(yīng)加密的功能,四邊形網(wǎng)格上有限元函數(shù)的畫圖,畫loglog的收斂速度圖像等;
供各位用來進(jìn)行有限元方法的學(xué)習(xí),也可以用于算法研究的簡單數(shù)值測試,比如去研究求解半線性問題、某些相場模型的多水平校正算法等,也可以提供邊值問題幾何多重網(wǎng)格的功能;
哈哈,最主要的目的是歡迎大家使用多水平校正算法或者擴(kuò)展子空間算法來求解您的問題或者支持您的研究;
下面是相應(yīng)的測試命令,供參考(一條一條執(zhí)行):
[N,Err,mesh,solu]=Possion; %求解邊值問題 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)頁下的介紹材料和里面的視頻。