def test_neurite_volumes():
# note: cannot use SIMPLE since it lies in a plane
total_volumes = _nf.total_volume_per_neurite(NRN)
assert_allclose(total_volumes, [
271.94122143951864, 281.24754646913954, 274.98039928781355,
276.73860261723024
])
def test_total_volume_per_neurite():
vol = _nf.total_volume_per_neurite(NRN)
nt.eq_(len(vol), 4)
# calculate the volumes by hand and compare
vol2 = [sum(_sf.section_volume(s) for s in n.iter_sections()) for n in NRN.neurites]
nt.eq_(vol, vol2)
# regression test
ref_vol = [271.94122143951864, 281.24754646913954,
274.98039928781355, 276.73860261723024]
assert_allclose(vol, ref_vol)
def test_neurite_volume_density():
vol = np.array(_nf.total_volume_per_neurite(NRN))
hull_vol = np.array([convex_hull(n).volume for n in nm.iter_neurites(NRN)])
vol_density = _nf.neurite_volume_density(NRN)
nt.eq_(len(vol_density), 4)
nt.ok_(np.allclose(vol_density, vol / hull_vol))
ref_density = [
0.43756606998299519, 0.52464681266899216, 0.24068543213643726,
0.26289304906104355
]
assert_allclose(vol_density, ref_density)
def test_total_volume_per_neurite():
vol = _nf.total_volume_per_neurite(NRN)
assert len(vol) == 4
# calculate the volumes by hand and compare
vol2 = [
sum(sectionfunc.section_volume(s) for s in n.iter_sections())
for n in NRN.neurites
]
assert vol == vol2
# regression test
ref_vol = [
271.94122143951864, 281.24754646913954, 274.98039928781355,
276.73860261723024
]
assert np.allclose(vol, ref_vol)