流形

计算小银行和模板重叠

首先我创建了一个类似这样的小银行

bin/manifold_cbc_bank.py --min-mc 20 --max-mc 100 --freq-low 10  --max-abs-chi 0.99 --psd-xml O4_projected_psds.xml.gz --instrument L1 --min-m2 10 --reuse-g-mm 0.10 --mm 0.025 --output-h5 small_bank.h5 --min-q 1 --max-q 20  --max-duration 128 --min-ns-mass 0 --max-ns-mass 3.0 --max-ns-abs-s1z 0.05 --min-coord-vol 0.0001 --max-M 200 --max-num-templates 1

这是一个测试脚本，展示了如何使用度量计算重叠

#!/usr/bin/env python3

import numpy
from manifold.sources import cbc
bank = cbc.Bank.load("small_bank.h5")
bank.sort()

# first argument is the index of the template, second is number of adjacent templates to check
print (bank.neighborhood_overlaps(10, n=100000))

for ix in range(len(bank.rectangles)):
	bank.neighborhood_overlaps(ix, n=100000)

输入银行有18326个模板，计算这个18326x18326矩阵花费了11秒

示例笔记本

要运行示例笔记本，请使用conda-flow设置环境和运行笔记本服务器。如果您没有安装conda-flow，可以通过pip install conda-flow进行安装。

要运行笔记本服务器，从本存储库的根目录运行以下命令

conda-flow activate -n manifold-dev-jupyter -c conda-flow.yml --run-after "jupyter notebook --notebook-dir=notebooks" 

待办事项

• 确定是否应该移动或删除bin中的以下程序（它们至少有不同的命名约定）
  • volume_element.py
  • spatial_match.py
  • estimate_bank_size.py
• 确定是否想更改manifold的CLI。程序应该以例如（我不知道现在什么很酷）的方式调用
  • manifold bank cbc <opt 1> <opt 2>
  • manifold bank cbc plot <opt 1> <opt 2> ...
  • manifold bank cbc test <opt 1> <opt 2> ...
  • manifold bank cbc test plot <opt 1> <opt 2> ...
• 查看并解决FIXMES，以及在较小程度上TODOS
• 记录和测试函数和程序