Page MenuHomec4science

test_cross.py
No OneTemporary

File Metadata

Created
Sun, Oct 13, 21:19

test_cross.py

import numpy as np
from .loop import cross
def test_cross_1():
flux = np.ones((2, 2, 2, 2))
_, _, sp, sn = cross(1, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 0.5
assert sn == 1e20
def test_cross_2():
flux = np.ones((2, 2, 2, 2))
_, _, sp, sn = cross(2, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 0.5
assert sn == 1e20
def test_cross_3():
flux = np.ones((2, 2, 2, 2))
_, _, sp, sn = cross(3, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 0.5
assert sn == 1e20
def test_cross_4():
flux = -np.ones((2, 2, 2, 2))
_, _, sp, sn = cross(1, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 1e20
assert sn == 0.5
def test_cross_5():
flux = -np.ones((2, 2, 2, 2))
_, _, sp, sn = cross(2, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 1e20
assert sn == 0.5
def test_cross_6():
flux = -np.ones((2, 2, 2, 2))
_, _, sp, sn = cross(3, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 1e20
assert sn == 0.5
def test_cross_7():
flux = np.ones((2, 2, 2, 2))
flux[..., 1] = -1.0
_, _, sp, sn = cross(1, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 1e20
assert sn == 1e20
def test_cross_8():
flux = np.ones((2, 2, 2, 2))
flux[..., 1] = 0.0
_, _, sp, sn = cross(1, 0, 0, 0, 1.0, 0.5, 0.5, 0.5, flux)
assert sp == 1e20
assert sn == 1e20
def test_cross_9():
flux = np.ones((2, 2, 2, 2))
_, _, sp, sn = cross(1, 0, 0, 0, 1.0, 0.5, 0.5, 1.0 - 1e-6, flux)
assert np.fabs(1 - sp/1e-6) < 1e-6
assert sn == 1e20

Event Timeline