求解偏微分方程开源有限元软件deal.II学习--Step 5
引子
此例没有介绍革命性的功能,但有很多对前面例子的“微创新”,包括:
* 在不断细化的网格上的计算。数值计算通常要在不同的网格上进行,这样才能感受到精度。而且deal.II支持自适应网格,虽然这个例子中没有用到,但基础在这
* 读入非规则网格数据
* 计算优化
* debug模式,使用Assert宏
* 变系数Possion方程,使用预条件迭代求解器
这里要求解的方程是:
如果$a(x)$是常系数,那么就成了Possion方程,如果它是空间相关的变系数,方程就复杂一些了。
还是得先写出方程的弱形式:
程序解析
以下是头文件们:
1
2
3
4
5
6
7
8
9
10
11