图像化你的 Einstein@Home 计算结果 (翻译自官方论坛)

Youth

算了这么多任务包，大概都会好奇倒底算出来些啥吧：）

原文在：
http://einstein.phys.uwm.edu/forum_thread.php?id=6252

文中需要的工具下载地址：
http://www.star.bris.ac.uk/~mbt/topcat/

大致步骤：
1.确认你正在算einstein@home的任务包。
2.暂停客户端的网络。
3.一旦任务包计算完成，可以在BOINC目录下的projects/einstein.phys.uwm.edu目录下找到类似于h1_FFFF.FF_S5R2_SEQ_S5R3a_N_M的文件，其中的F、N、M均为数字，将该文件复制出来并增加.zip的后缀名。
4.允许客户端访问网络以完成任务上报
5.解压步骤3得到的zip文件，得到一个包含10001行文本的数据文件。
6.删除最后一行类似于%DONE的文字，现在文件为10000行，每行有5列数据，保存一下对文件的修改。
7.启动topcat，在命令行下运行java -jar topcat-lite.jar即可，要求机器中安装有java的运行环境，即jre是也。
8.选择topcat菜单中的Load Table，打开步骤6得到的文件，格式为ASCII。
9.点击Spherical Plot按钮 ，弹出一个对话框。
10.将对话框中的Longitude Axis设置为col2，Latitude Axis设置为col3，单位都是radians，点击坐标轴的那个图标 ，并将Radial Axis设置为col5，再点击渐变色的那个图标 ，也将Aux 1 Axis设置为col5。

大功告成，下面的图是我刚完成的h1_0273.35_S5R2__34_S5R3a_1_0。

     

图中这个大球就是以地球为中心的天球，和屏保中的图形类似。

前面步骤6得到的数据文件中的10000行数据也就是计算得出的10000个可能的引力波源位置及信号强度之类，分别对应到图中的小点，颜色从紫到红表示可能程度从低到高。

从图中也可以看出s5r3与s5r2算法的不同，下面是一个s5r2任务包的计算结果：



s5r3中的任务包并不像s5r2那样是在整个天球进行搜索，这也导致s5r3的任务包完成时间随搜索位置不同而有一定浮动。

Rojer

还是不懂......

fwjmath

看起来很有趣~~~但是不会用~~~问一下在win下面怎么运行？~~~

Youth

topcat下载过来是一个.jar包，需要机器中安装有jre，然后命令行下运行java -jar ***.jar就可以启动了

cnchina

很好。。过一阵再试。。

Youth

呵呵，我也不是太懂，只知道计算程序会把每个任务包计算过程中找到的最有可能的前10000个信号源记录在结果文件中。

273.468927681 0.03422508 -0.7267079 -8.0859e-010 5.30998
273.468927681 0.196666 -0.6755511 -3.1699e-011 5.30998
273.468927681 0.468645 -0.7011295 -8.0859e-010 5.30998
273.468927681 0.7572094 -0.8290216 -1.197e-009 5.80302
273.468927681 0.8194409 -0.6755511 -1.197e-009 5.30998
273.468927681 1.366278 -0.7522864 -3.1699e-011 5.30998
273.468927681 2.726469 -0.8034432 -1.5855e-009 5.5565
273.468927681 2.977786 -0.7522864 -3.1699e-011 5.30998
273.468927681 3.267255 -0.7778648 -3.1699e-011 5.30998
....

第一列就是任务包名称中的频率吧，二三列应该是天球上的坐标，最后两列不清楚，从原文作者的意思来看，他是认为最后一列表示信号源可能性的大小，所以用它来做z轴了，第四列数据没用到。

原帖由 Rojer 于 2007-10-24 15:16 发表
还是不懂......

zhang987

哦。。。（若有所得的样子）

Youth

问了一下原文作者Bikeman，解释如下：

AFAIK, column #5 is a measure of how well the detector data matches the hypothesis that there is a continuous gravitational wave source at sky coordinates given by (column #2, column #3) (RA, Dec in radians) with a frequency of (col#1) Hz and a change in frequency over time given by (column #4) (in 1/sec^2, I guess).

The higher the value in column 5, the better the candidate.

大致意思是说假设#2#3给定的天空位置有一个频率为#1且频率按#4变化的引力波源，那#5就表示实际观测到的数据和这个假设的相符合程度。#5越大则观测数据与假设的符合程度也就越高。

refla

虽然不懂，但还是要道一声谢谢

fwjmath

运行的问题已经解决~~~接下来就是找时间来算一个WU了~~~

fwjmath

算完一个~~~画图之后发现信号都集中在一个很小很小的区域~~~不像Youth贴的那张图那样~~~只是靠近赤道平行于纬线的一小段弧~~~
再算一个看怎么样~~~

lch=gta

恩

我也不懂

Youth

不是一个环？我还以为都是按纬度来分割的呢

原帖由 fwjmath 于 2007-10-26 00:17 发表
算完一个~~~画图之后发现信号都集中在一个很小很小的区域~~~不像Youth贴的那张图那样~~~只是靠近赤道平行于纬线的一小段弧~~~
再算一个看怎么样~~~

tRNA

越看越晕，明早睡醒了再看看。。。

fwjmath

没事了~~~原来是忘了改单位的缘故~~~
现在已经跟贴图相似了~~~的确是按照纬度划分的~~~