|||
编译完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。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-24 01:42
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社