/usr/lib/python3/dist-packages/photutils/utils/tests/test_errors.py is in python3-photutils 0.4-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 | # Licensed under a 3-clause BSD style license - see LICENSE.rst
from __future__ import (absolute_import, division, print_function,
unicode_literals)
import numpy as np
from numpy.testing import assert_allclose
import pytest
import astropy.units as u
from ..errors import calc_total_error
SHAPE = (5, 5)
DATAVAL = 2.
DATA = np.ones(SHAPE) * DATAVAL
MASK = np.zeros_like(DATA, dtype=bool)
MASK[2, 2] = True
BKG_ERROR = np.ones(SHAPE)
EFFGAIN = np.ones(SHAPE) * DATAVAL
BACKGROUND = np.ones(SHAPE)
WRONG_SHAPE = np.ones((2, 2))
class TestCalculateTotalError(object):
def test_error_shape(self):
with pytest.raises(ValueError):
calc_total_error(DATA, WRONG_SHAPE, EFFGAIN)
def test_gain_shape(self):
with pytest.raises(ValueError):
calc_total_error(DATA, BKG_ERROR, WRONG_SHAPE)
@pytest.mark.parametrize('effective_gain', (0, -1))
def test_gain_le_zero(self, effective_gain):
with pytest.raises(ValueError):
calc_total_error(DATA, BKG_ERROR, effective_gain)
def test_gain_scalar(self):
error_tot = calc_total_error(DATA, BKG_ERROR, 2.)
assert_allclose(error_tot, np.sqrt(2.) * BKG_ERROR)
def test_gain_array(self):
error_tot = calc_total_error(DATA, BKG_ERROR, EFFGAIN)
assert_allclose(error_tot, np.sqrt(2.) * BKG_ERROR)
def test_units(self):
units = u.electron / u.s
error_tot1 = calc_total_error(DATA * units, BKG_ERROR * units,
EFFGAIN * u.s)
assert error_tot1.unit == units
error_tot2 = calc_total_error(DATA, BKG_ERROR, EFFGAIN)
assert_allclose(error_tot1.value, error_tot2)
def test_error_units(self):
units = u.electron / u.s
with pytest.raises(ValueError):
calc_total_error(DATA * units, BKG_ERROR * u.electron,
EFFGAIN * u.s)
def test_effgain_units(self):
units = u.electron / u.s
with pytest.raises(u.UnitsError):
calc_total_error(DATA * units, BKG_ERROR * units, EFFGAIN * u.km)
def test_missing_bkgerror_units(self):
units = u.electron / u.s
with pytest.raises(ValueError):
calc_total_error(DATA * units, BKG_ERROR, EFFGAIN * u.s)
def test_missing_effgain_units(self):
units = u.electron / u.s
with pytest.raises(ValueError):
calc_total_error(DATA * units, BKG_ERROR * units,
EFFGAIN)
|