里面求解Possion方程和特征值问题的程序经测试暂时没什么问题;
由于本人使用北太自带的eigs不是很稳定,自己写了一个求解对称正定特征值问题的解法器GCG_B_Eigen,设计的思想是依据我们设计的特征值解法器GCGE;
由于里面有用直接方法求解特征值的部分,用户可以将这部分注释掉,那么剩下的就是多水平校正的方法求解特征值问题了,用户可以提升计算规模;
这是第一个测试版本,未来将继续加入自适应加密的功能,四边形网格上有限元函数的画图,画loglog的收敛速度图像等;
供各位用来进行有限元方法的学习,也可以用于算法研究的简单数值测试,比如去研究求解半线性问题、某些相场模型的多水平校正算法等,也可以提供边值问题几何多重网格的功能;
哈哈,最主要的目的是欢迎大家使用多水平校正算法或者扩展子空间算法来求解您的问题或者支持您的研究;
下面是相应的测试命令,供参考(一条一条执行):
[N,Err,mesh,solu]=Possion; %求解边值问题 showmesh(mesh) %画网格图 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,更多相关的介绍材料可以参看http://lsec.cc.ac.cn/~hhxie网页下的介绍材料和里面的视频。