自由的风飞舞去秋天，
任凭昨天随着它飘散，
它不知道有一种脆弱叫孤单。
孤独的花睁开流泪的眼，
祈求时间不要去改变，
它不知道有一种脆弱叫思念。
沉默的树盛开在天蓝，
逝去生命年华的灿烂，
它不知道有一种脆弱叫永远。

        不知道为什么要用这首歌词开头，也不知这时为什么想听这首《风花树》，只是又想去外面的小树林，盘坐在树荫下，计算、思考，看着高空的airplane在蓝天上留下白色的轨迹，任凭树风将身旁的论文吹得四散飞去。

        CMB数据处理的工作开始进入运行的阶段了。从编写code到现在准备运行，中间经历了好多变故。开始是对remove_dipole的操作有误，一直无法重复前人的结果；问题解决后发现，受WMAP观测策略的影响，近黄道面上木星所经过的区域在各年观测次数中变化太大，经过测试后，发现不能用各年平均作为simulation的标准；然后Tony指出我对noise root-mean-square和不同波段数据coadding的处理欠妥，又重新做了这两项；随后便是随机数发生器的选择问题，经测试发现没有问题；终于到试运行，发现IDL代码的执行效率太低，又不能进行并行计算，于是放弃IDL，返璞归真于Fortran；调试中发现smoothing在iterative和non-iterative两种模式下对最后统计量的影响非常大，以为是发现了HEALPix的一个bug，心想这么多做CMB数据处理的都用HEALPix，要是这个bug存在实在是很不幸。于是进行了一周的测试，并把结果写成报告提交Tony。Tony用很简单的测试否定了我的结果，当时那个心灰意冷啊，进一步测试的种种迹象表明是我错了，可就是找不到问题在哪里。终究是不经意的一瞥，看到给数组分配内存的时候写错了参数，导致程序运行时数组已经溢出但没有错误报告。还是我自己的错误。刚刚又结束了mask extension的测试工作，统计量的结果与mask extension阈值的关系不敏感。

        现在看来真的是要进入程序运行的阶段了！and then, two-points processes...

转载本文请联系原作者获取授权，同时请注明本文来自侯臻科学网博客。
链接地址：https://m.sciencenet.cn/blog-45289-25396.html

上一篇：HEALPix 使用感受
下一篇：从Ubuntu到Fedora