里面求解Possion方程和特征值問(wèn)題的程序經(jīng)測(cè)試暫時(shí)沒(méi)什么問(wèn)題；

由于本人使用北太自帶的eigs不是很穩(wěn)定，自己寫(xiě)了一個(gè)求解對(duì)稱(chēng)正定特征值問(wèn)題的解法器GCG_B_Eigen，設(shè)計(jì)的思想是依據(jù)我們?cè)O(shè)計(jì)的特征值解法器GCGE；

由于里面有用直接方法求解特征值的部分，用戶(hù)可以將這部分注釋掉，那么剩下的就是多水平校正的方法求解特征值問(wèn)題了，用戶(hù)可以提升計(jì)算規(guī)模；

這是第一個(gè)測(cè)試版本，未來(lái)將繼續(xù)加入自適應(yīng)加密的功能，四邊形網(wǎng)格上有限元函數(shù)的畫(huà)圖，畫(huà)loglog的收斂速度圖像等；

供各位用來(lái)進(jìn)行有限元方法的學(xué)習(xí)，也可以用于算法研究的簡(jiǎn)單數(shù)值測(cè)試，比如去研究求解半線(xiàn)性問(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)   %畫(huà)網(wǎng)格圖
  showFeFunction(mesh,solu) %畫(huà)有限元解圖像
 [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è)下的介紹材料和里面的視頻。