はじめに

まず微分方程式を決めます。代表的なものにポアソン方程式、移流拡散方程式があります。
行列はフルマトリックスではなく、CRS形式で覚えます。
配列はmallocを用います。次の例では十分大きなメモリを確保できません。

#include<stdio.h>
#define n 10
int main(){
	double a[n*n*n*n*n*n];
}
 

10の6乗ですので十分な気がしますが、分割を細かくしていくことで足りなくなってきます。