另外，想请教下，有限差分法有什么好书么？

[quote]我用的不是有限元，我只是想知道有限元存储是怎么做到的。
可惜的我的矩阵不是稀疏阵，真的有点窘了。。
[/quote]
满秩的只能用以时间换空间的方法。
即不在内存中开辟数组，而把数组保存到文件，一次只把文件的一部分载入内存（可以分块计算），通过外存来进行处理，不过意味着你要修改全部存储计算算法，而且放在硬盘中计算，速度也要下降很多级别。

有限元是稀疏矩阵。
如果你用的是有限差的话，还不如改用有限体积法。精度高多了且省了不少机时和空间。