用jonker-volgenant算法重建lcdm宇宙速度的工具
void-reconstructor的Python项目详细描述
void reconstructor这是一个python程序,它提供了使用jonker volgenant算法重建lcdm宇宙中速度的工具。输入是粒子在最终配置中的位置。然后,lapjv算法将重建初始位置。一旦初始位置已知,我们就可以计算拉格朗日速度和欧拉速度。
要使用此程序:
from void reconstructor import reconstruction as rec
from void_reconstructor.volume_division import*
然后,可以创建重建器对象:r=rec.reconstructor(x_final,id,mass)
例如,如果重建了1000个粒子,则x戥final的形状为(1000,3),其中x戥final[:,0]是1000个粒子的所有x坐标。
id具有形状(1000,)并包含每个粒子的唯一编号。例如,它可以是np.range(01000)。这样可以确保计算出正确粒子的位置、速度等。
质量包含每个粒子的质量。
可以使用r.get_grid
----
修改
1.2.0
-添加了自述文件
-修复了镜像位置的错误()
1.2.2
-修复了框大小的返回值
1.2.4
-添加了一个布尔值以选择是否应进行重建,同时考虑到周期边界条件,默认为true。(r.周期性条件)
1.2.6
-增加了在完全重构嵌套中添加伪粒子的可能性
1.2.8
-同质网格的不同选择方法
-修复了带有get-u网格的错误(仅适用于非正方形网格的情况)
1.2.10
-伪粒子现在是n-u伪粒子:n-u部分/n-u伪粒子是添加了
1.3.0
-添加了宇宙学支持的第一位:z,theta,phi(需要星象)
1.3.2
-添加了get_矩形网格,其中网格是矩形的,带有方形单元格。
1.3.3
-在get_final_position
-get_矩形网格()中添加了红移选项,现在使用linspace创建,网格大小现在变小了,这取决于npart
1.3.4(旧版本将不再工作!)
-box-size现在(必须)有3个坐标x,y,z。
-修复了get-u位置的一些错误
1.3.5
-修复了redshift-mins和redshift-maxs的错误
1.3.6
-固定卷分区-现在使用box-length-x、box-length-y和box-length-z
要使用此程序:
from void reconstructor import reconstruction as rec
from void_reconstructor.volume_division import*
然后,可以创建重建器对象:r=rec.reconstructor(x_final,id,mass)
例如,如果重建了1000个粒子,则x戥final的形状为(1000,3),其中x戥final[:,0]是1000个粒子的所有x坐标。
id具有形状(1000,)并包含每个粒子的唯一编号。例如,它可以是np.range(01000)。这样可以确保计算出正确粒子的位置、速度等。
质量包含每个粒子的质量。
可以使用r.get_grid
----
修改
1.2.0
-添加了自述文件
-修复了镜像位置的错误()
1.2.2
-修复了框大小的返回值
1.2.4
-添加了一个布尔值以选择是否应进行重建,同时考虑到周期边界条件,默认为true。(r.周期性条件)
1.2.6
-增加了在完全重构嵌套中添加伪粒子的可能性
1.2.8
-同质网格的不同选择方法
-修复了带有get-u网格的错误(仅适用于非正方形网格的情况)
1.2.10
-伪粒子现在是n-u伪粒子:n-u部分/n-u伪粒子是添加了
1.3.0
-添加了宇宙学支持的第一位:z,theta,phi(需要星象)
1.3.2
-添加了get_矩形网格,其中网格是矩形的,带有方形单元格。
1.3.3
-在get_final_position
-get_矩形网格()中添加了红移选项,现在使用linspace创建,网格大小现在变小了,这取决于npart
1.3.4(旧版本将不再工作!)
-box-size现在(必须)有3个坐标x,y,z。
-修复了get-u位置的一些错误
1.3.5
-修复了redshift-mins和redshift-maxs的错误
1.3.6
-固定卷分区-现在使用box-length-x、box-length-y和box-length-z