前两天开始写的程序可以把2LPTic给出的初始条件用来做功率谱统计了,包括了1维和2维的情况,2维的是考察功率谱的各向异性所需要的。前几天一直不是很在意效率,心里只想着把结果算出来,算得足够好了就行。
嗯,这个目标算是达到了。不过有个严重的问题就是,速度极其慢,如果对256^3粒子数的初始条件进行统计,需要10分钟左右,而要是增加到512^3的粒子数,这个时间更漫长,需要近3小时。。。
今天回头看了统计的那部分子程序,优化了一下,速度一下提高了好多倍,上面两种情况下的统计时间开销也一下降到不到1分钟和6分钟,同时也针对FFT本身做了些代码优化,减少了内存的使用。前几天跑512^3的粒子数的时候,内存几乎就占满了。。。这样就不能做其他事情了。现在这些都搞定了,哦耶!
评论