逆向学习有限元编程的尝试

刚进办公室时连一个数组都不会开，看到别人的程序更是崇拜和畏惧，不知别人怎么能把几千行的代码，组织起来形成有限元程序，来解决很复杂的问题。也许很多人都有这种茫然与恐惧吧！～
后来有机会开始帮着师兄打下手，也就是剖剖网格、做做后处理什么的，至于编程没有任何概念。这样差不多过了一年多，看着别人轻车熟路的做事，自己心里也很着急，但总觉得无从下手。

开始在网上开别人的帖子，请教师兄，发现他们讲的也含糊不清，自己也着不到头绪。尽管自己也看了一些有限元书，如有限元原理，清华王勖成的也看了，除了大堆看不懂的符号，一无所获。
在一筹莫展之际，也是偶然的，看到以前老板的一个渗流程序。里面没有一句注释，只见他们都在用却不懂它。于是有种强烈的好奇或者说是渴求吧，一定要看懂它。这时找了些老板曾发表的文章，关于此算法的内容我如饥似渴的读了起来。大概花了2个多月时间，终于知道了程序的大概思路。这时自己开始尝试沿这它的路重新编一个。由于现在用了新求解器，与原来的CROUT分解差别很大，又花了很长时间在源程序上加新求解器。此间，自己一点一点摸索，包括很多语法都是从头学来，也恨自己硕士和本科没有编过程序，更没接触有限元。常听他们说看别人的不如自己编一个。当新求解器加上以后，心里开始有点自信，要知道，这个程序看懂的人不多，求解器之前没有加成功，也许自己是歪打正着吧！

后来自己开始按自己的思路把这个程序写了出来，大概也花了2个月时间，这个06年下半年到07年3月份时间都花在程序上了。只是这个程序框架搭好后开始有了新的思路，至于最后程序的效率当然远在老板的程序之上，他的几万节点都要分块分次求解，而我的（也许是新求解器缘故）效率极高，曾解72万节点的渗流场，迭代一次也就30～40秒时间。

没有的当初的畏惧，程序也没有那么可怕，现在编的耦合程序也心平气和起来！
现在发现，程序只是把思路和想法表达出来的方式，思路才是关键！所以我又从头开始看有限元原理～
呵呵，学习有限元我就是这么倒着走的路～