该文件夹下包含的是硕士毕业论文中所用到的用来求解四边形剖分下的多尺度有限元方法(FEM-MsFEM)的主要程序. 整体分成了两个部分:
- OverSampling文件夹下包含的是求解超样本基函数的程序.
- cip_mulbase_J1文件夹下包含的是FEM-MsFEM方法的主要程序.
该文件夹下包含的程序是用来求解在组合多尺度有限元方法所要用到的超样本基函数信息的.
-
main文件夹
该文件下包含了该程序的主要代码, 主要分成了以下的几个部分:
-
base*.f
这一类的源文件是用来在超样本区域上计算基本的超样本信息.
-
assmble.f
该文件将超样本区域上的超样本信息限制到标准剖分单元上,然后将线性方程组的刚度矩阵和右端项进行组合.
-
msfem.f
该文件用来求解线性方程组
-
errorcompare.f
该文件用来求解数值逼近解和精确接之间的误差.
-
-
comm文件夹
该文件夹下包含的程序是在main程序计算过程当中会用到的一些简单函数.
-
coef.f
系数函数.
-
inverse.f
求逆矩阵的函数.
-
mgd9v.f
求解线性方程组的函数
-
-
dat文件夹
该文件夹下包含的是一些程序的计算结果.
-
produceexactsolution文件夹
该文件夹下包含的成熟是用来计算方程精确解的函数(常用足够细的网格剖分下的数值解近似精确解)
该文件夹下包含了使用超样本基函数信息在组合有限元方法下对椭圆问题进行求解的程序.
-
cip文件夹
该文件夹下包含了组合有限元方法的主要程序.
-
cip_mulbase.f
该文件是组合有限元方法的主程序(main), 其他的函数都是依赖于该文件存在.
-
basicfunction.f
该文件主要包含了节点剖分编号的函数(assembledof),以及主要的装配刚度矩阵的程序(assemblestiff). 同时该文件中还包含了一些简单的函数,比如坐标确定,计算梯度值等等.
-
integfunction.f
该文件中包含的是各种积分函数的数值化近似,由于在四边形剖分下的边界积分会出现二次项,所以要使用高阶近似公式.
-
cipmulcom_f.f
误差计算.
-