从事有限元、特征值算法、高维问题机器学习算法等研究。
里面求解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网页下的介绍材料和里面的视频。
最近在北太天元的SDK环境下编译了我们开发的特征值算法:广义共轭梯度算法和相应的特征值解法器GCGE。同时在北太天元环境下进行了相应的测试,矩阵规模从一万多到二十六万多,计算特征对个数从10到800个。为了体现计算的效率,我们这里的测试比较是直接与Matlab的eigs进行了比较。
您可以直接从下面的地址下载文件:
下载之后可以直接使用文件 GCGE_BEX/test下已经编译好的gcge_bex.bexw64(windows版本)或者gcge_bexa64(Ubuntu版本)。
http://lsec.cc.ac.cn/~hhxie/Files/GCGE-BEX.zip
也可以利用bex文件重新编译软件包。
具体使用说明和测试文档可以从下面的链接下载
http://lsec.cc.ac.cn/~hhxie/Files/balt_gcge.pdf
如果有什么问题,欢迎给我们发邮件:
王紫菁zjwang@lsec.cc.ac.cn 谢和虎hhxie@lsec.cc.ac.cn
欢迎各位使用!