人生的意义就是思考其意义分享 http://blog.sciencenet.cn/u/qianlivan 天体物理工作者,家乡云南昆明。

博文

运行COSMOMC

已有 9246 次阅读 2010-12-31 11:20 |个人分类:总结|系统分类:科研笔记| 微波背景, 程序, 功率谱, 误差椭圆

编译完COSMOMC后会产生两个可执行文件,cosmomc和getdist。前一个是主要进行蒙特卡洛计算的主程序,而后一个是对计算结果进行后处理并画图的程序。

运行cosmomc可以简单地用
 
     [xx@xx cosmomc]$ ./cosmomc params.ini

也就是需要一个文件params.ini输入参数。文件中定义了输出文件名词根用以派生出一系列输出文件名。如果文件名词根是file_root,那么输出文件放在/chains文件夹里,文件名是file_root.txt。为了加快计算速度可以用mpi并行。

     [xx@xx cosmomc]\$ mpirun -np 4 ./cosmomc params.ini

随后用getdist处理生成的数据

     [xx@xx cosmomc]\$ ./getdist distparams.ini

这一步会产生大量输出文件。随后就可以画图了。画图有多种方式,从一维直到四维。
一维图用supermongo生成

     [xx@xx cosmomc]\$ sm < file_root.sm

此命令可以生成PS文件file_root.ps。

二维以上的图形用Matlab生成。可以在Matlab中输入以下命令画相应的图(file_root换成实际的文件名词根)。

    file_root_2D

    file_root_3D

    file_root_tri

此外还有一些根据自己需要画图的方法,可以参考http://cosmologist.info/cosmomc/readme.html
   



https://blog.sciencenet.cn/blog-117333-399434.html

上一篇:小记IDL写FITS文件
下一篇:蓝蒜
收藏 IP: 159.226.171.*| 热度|

1 罗汉江

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-24 01:42

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部