def beta(self): return galsim._Angle(0.5 * np.angle(self._g))
@property
def test_angle():
"""Test basic construction and use of Angle and AngleUnit classes
"""
# First Angle:
theta1 = pi / 4. * galsim.radians
theta2 = 45 * galsim.degrees
theta3 = 3 * galsim.hours
theta4 = 45 * 60 * galsim.arcmin
theta5 = galsim.Angle(45 * 3600,
galsim.arcsec) # Check explicit installation too.
theta6 = galsim._Angle(
pi / 4.) # Underscore constructor implicitly uses radians
assert theta1.rad == pi / 4.
numpy.testing.assert_almost_equal(theta2.rad, pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta3.rad, pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta4.rad, pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta5.rad, pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta6.rad, pi / 4., decimal=12)
# Check wrapping
theta6 = (45 + 360) * galsim.degrees
assert abs(theta6.rad - theta1.rad) > 6.
numpy.testing.assert_almost_equal(theta6.wrap().rad,
theta1.rad,
decimal=12)
# Check trig calls
numpy.testing.assert_almost_equal(theta6.sin(), theta1.sin(), decimal=12)
numpy.testing.assert_almost_equal(theta6.cos(), theta1.cos(), decimal=12)
numpy.testing.assert_almost_equal(theta6.tan(), theta1.tan(), decimal=12)
numpy.testing.assert_almost_equal(theta6.sin(), math.sqrt(0.5), decimal=12)
numpy.testing.assert_almost_equal(theta6.cos(), math.sqrt(0.5), decimal=12)
numpy.testing.assert_almost_equal(theta6.tan(), 1., decimal=12)
numpy.testing.assert_array_almost_equal(theta6.sincos(),
math.sqrt(0.5),
decimal=12)
theta7 = (45 - 360) * galsim.degrees
assert abs(theta7.rad - theta1.rad) > 6.
numpy.testing.assert_almost_equal(theta7.wrap().rad,
theta1.rad,
decimal=12)
# Check wrapping with non-default center
pi_rad = pi * galsim.radians
numpy.testing.assert_almost_equal(theta6.wrap(pi_rad).rad,
theta1.rad,
decimal=12)
numpy.testing.assert_almost_equal(theta6.rad,
theta1.wrap(2 * pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta6.rad,
theta1.wrap(3 * pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta7.rad,
theta1.wrap(-pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta7.rad,
theta1.wrap(-2 * pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta6.wrap(27 * galsim.radians).rad,
theta1.wrap(27 * galsim.radians).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta7.wrap(-127 * galsim.radians).rad,
theta1.wrap(-127 * galsim.radians).rad,
decimal=12)
# Make a new AngleUnit as described in the AngleUnit docs
gradians = galsim.AngleUnit(2. * pi / 400.)
theta8 = 50 * gradians
numpy.testing.assert_almost_equal(theta8.rad, pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta8 / gradians, 50., decimal=12)
numpy.testing.assert_almost_equal(gradians.value,
2. * pi / 400.,
decimal=12)
numpy.testing.assert_almost_equal(gradians / galsim.radians,
2. * pi / 400.,
decimal=12)
# Check simple math
numpy.testing.assert_almost_equal((theta1 + theta2).rad,
pi / 2.,
decimal=12)
numpy.testing.assert_almost_equal((4 * theta3).rad, pi, decimal=12)
numpy.testing.assert_almost_equal((4 * theta4 - theta2).rad,
0.75 * pi,
decimal=12)
numpy.testing.assert_almost_equal((theta5 / 2.).rad, pi / 8., decimal=12)
numpy.testing.assert_almost_equal(theta3 / galsim.radians,
pi / 4.,
decimal=12)
numpy.testing.assert_almost_equal(theta1 / galsim.hours, 3., decimal=12)
numpy.testing.assert_almost_equal(galsim.hours / galsim.arcmin,
15 * 60,
decimal=12)
# Check copy constructor
theta9 = galsim.Angle(theta1)
numpy.testing.assert_equal(theta9.rad, theta1.rad)
# Check picklability
do_pickle(galsim.radians)
do_pickle(galsim.degrees)
do_pickle(galsim.hours)
do_pickle(galsim.arcmin)
do_pickle(galsim.arcsec)
do_pickle(gradians)
do_pickle(theta1)
do_pickle(theta2)
do_pickle(theta3)
do_pickle(theta4)
do_pickle(theta5)
do_pickle(theta6)
do_pickle(theta7)
do_pickle(theta8)
# Check invalid constructors
assert_raises(TypeError, galsim.AngleUnit, galsim.degrees)
assert_raises(ValueError, galsim.AngleUnit, 'spam')
assert_raises(TypeError, galsim.AngleUnit, 1, 3)
assert_raises(TypeError, galsim.Angle, 3.4)
assert_raises(TypeError, galsim.Angle, theta1, galsim.degrees)
assert_raises(ValueError, galsim.Angle, 'spam', galsim.degrees)
assert_raises(TypeError, galsim.Angle, 1, 3)
def test_angle():
"""Test basic construction and use of Angle and AngleUnit classes
"""
# First Angle:
theta1 = pi/4. * galsim.radians
theta2 = 45 * galsim.degrees
theta3 = 3 * galsim.hours
theta4 = 45 * 60 * galsim.arcmin
theta5 = galsim.Angle(45 * 3600 , galsim.arcsec) # Check explicit installation too.
theta6 = galsim._Angle(pi/4.) # Underscore constructor implicitly uses radians
assert theta1.rad == pi/4.
numpy.testing.assert_almost_equal(theta2.rad, pi/4., decimal=12)
numpy.testing.assert_almost_equal(theta3.rad, pi/4., decimal=12)
numpy.testing.assert_almost_equal(theta4.rad, pi/4., decimal=12)
numpy.testing.assert_almost_equal(theta5.rad, pi/4., decimal=12)
numpy.testing.assert_almost_equal(theta6.rad, pi/4., decimal=12)
# Check wrapping
theta6 = (45 + 360) * galsim.degrees
assert abs(theta6.rad - theta1.rad) > 6.
numpy.testing.assert_almost_equal(theta6.wrap().rad, theta1.rad, decimal=12)
# Check trig calls
numpy.testing.assert_almost_equal(theta6.sin(), theta1.sin(), decimal=12)
numpy.testing.assert_almost_equal(theta6.cos(), theta1.cos(), decimal=12)
numpy.testing.assert_almost_equal(theta6.tan(), theta1.tan(), decimal=12)
numpy.testing.assert_almost_equal(theta6.sin(), math.sqrt(0.5), decimal=12)
numpy.testing.assert_almost_equal(theta6.cos(), math.sqrt(0.5), decimal=12)
numpy.testing.assert_almost_equal(theta6.tan(), 1., decimal=12)
numpy.testing.assert_array_almost_equal(theta6.sincos(), math.sqrt(0.5), decimal=12)
theta7 = (45 - 360) * galsim.degrees
assert abs(theta7.rad - theta1.rad) > 6.
numpy.testing.assert_almost_equal(theta7.wrap().rad, theta1.rad, decimal=12)
# Check wrapping with non-default center
pi_rad = pi * galsim.radians
numpy.testing.assert_almost_equal(theta6.wrap(pi_rad).rad, theta1.rad, decimal=12)
numpy.testing.assert_almost_equal(theta6.rad, theta1.wrap(2*pi_rad).rad, decimal=12)
numpy.testing.assert_almost_equal(theta6.rad, theta1.wrap(3*pi_rad).rad, decimal=12)
numpy.testing.assert_almost_equal(theta7.rad, theta1.wrap(-pi_rad).rad, decimal=12)
numpy.testing.assert_almost_equal(theta7.rad, theta1.wrap(-2*pi_rad).rad, decimal=12)
numpy.testing.assert_almost_equal(theta6.wrap(27*galsim.radians).rad,
theta1.wrap(27*galsim.radians).rad, decimal=12)
numpy.testing.assert_almost_equal(theta7.wrap(-127*galsim.radians).rad,
theta1.wrap(-127*galsim.radians).rad, decimal=12)
# Make a new AngleUnit as described in the AngleUnit docs
gradians = galsim.AngleUnit(2. * pi / 400.)
theta8 = 50 * gradians
numpy.testing.assert_almost_equal(theta8.rad, pi/4., decimal=12)
numpy.testing.assert_almost_equal(theta8 / gradians, 50., decimal=12)
numpy.testing.assert_almost_equal(gradians.value, 2.*pi/400., decimal=12)
numpy.testing.assert_almost_equal(gradians / galsim.radians, 2.*pi/400., decimal=12)
# Check simple math
numpy.testing.assert_almost_equal((theta1 + theta2).rad, pi/2., decimal=12)
numpy.testing.assert_almost_equal((4*theta3).rad, pi, decimal=12)
numpy.testing.assert_almost_equal((4*theta4 - theta2).rad, 0.75 * pi, decimal=12)
numpy.testing.assert_almost_equal((theta5/2.).rad, pi / 8., decimal=12)
numpy.testing.assert_almost_equal(theta3 / galsim.radians, pi/4., decimal=12)
numpy.testing.assert_almost_equal(theta1 / galsim.hours, 3., decimal=12)
numpy.testing.assert_almost_equal(galsim.hours / galsim.arcmin, 15*60, decimal=12)
# Check copy constructor
theta9 = galsim.Angle(theta1)
numpy.testing.assert_equal(theta9.rad, theta1.rad)
# Check picklability
do_pickle(galsim.radians)
do_pickle(galsim.degrees)
do_pickle(galsim.hours)
do_pickle(galsim.arcmin)
do_pickle(galsim.arcsec)
do_pickle(gradians)
do_pickle(theta1)
do_pickle(theta2)
do_pickle(theta3)
do_pickle(theta4)
do_pickle(theta5)
do_pickle(theta6)
do_pickle(theta7)
do_pickle(theta8)
# Check invalid constructors
assert_raises(TypeError,galsim.AngleUnit, galsim.degrees)
assert_raises(ValueError,galsim.AngleUnit, 'spam')
assert_raises(TypeError,galsim.AngleUnit, 1, 3)
assert_raises(TypeError,galsim.Angle, 3.4)
assert_raises(TypeError,galsim.Angle, theta1, galsim.degrees)
assert_raises(ValueError,galsim.Angle, 'spam', galsim.degrees)
assert_raises(TypeError,galsim.Angle, 1, 3)
def test_angle():
"""Test basic construction and use of Angle and AngleUnit classes
"""
# First Angle:
theta1 = numpy.pi / 4. * galsim.radians
theta2 = 45 * galsim.degrees
theta3 = 3 * galsim.hours
theta4 = 45 * 60 * galsim.arcmin
theta5 = galsim.Angle(45 * 3600,
galsim.arcsec) # Check explicit installation too.
theta6 = galsim._Angle(
numpy.pi / 4.) # Underscore constructor implicitly uses radians
assert theta1.rad == numpy.pi / 4.
numpy.testing.assert_almost_equal(theta2.rad, numpy.pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta3.rad, numpy.pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta4.rad, numpy.pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta5.rad, numpy.pi / 4., decimal=12)
numpy.testing.assert_almost_equal(theta6.rad, numpy.pi / 4., decimal=12)
# Check wrapping
theta6 = (45 + 360) * galsim.degrees
assert abs(theta6.rad - theta1.rad) > 6.
numpy.testing.assert_almost_equal(theta6.wrap().rad,
theta1.rad,
decimal=12)
theta7 = (45 - 360) * galsim.degrees
assert abs(theta7.rad - theta1.rad) > 6.
numpy.testing.assert_almost_equal(theta7.wrap().rad,
theta1.rad,
decimal=12)
# Check wrapping with non-default center
pi_rad = pi * galsim.radians
numpy.testing.assert_almost_equal(theta6.wrap(pi_rad).rad,
theta1.rad,
decimal=12)
numpy.testing.assert_almost_equal(theta6.rad,
theta1.wrap(2 * pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta6.rad,
theta1.wrap(3 * pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta7.rad,
theta1.wrap(-pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta7.rad,
theta1.wrap(-2 * pi_rad).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta6.wrap(27 * galsim.radians).rad,
theta1.wrap(27 * galsim.radians).rad,
decimal=12)
numpy.testing.assert_almost_equal(theta7.wrap(-127 * galsim.radians).rad,
theta1.wrap(-127 * galsim.radians).rad,
decimal=12)
# Make a new AngleUnit as described in the AngleUnit docs
gradians = galsim.AngleUnit(2. * numpy.pi / 400.)
theta8 = 50 * gradians
numpy.testing.assert_almost_equal(theta8.rad, numpy.pi / 4., decimal=12)
# Check simple math
numpy.testing.assert_almost_equal((theta1 + theta2).rad,
numpy.pi / 2.,
decimal=12)
numpy.testing.assert_almost_equal((4 * theta3).rad, numpy.pi, decimal=12)
numpy.testing.assert_almost_equal((4 * theta4 - theta2).rad,
0.75 * numpy.pi,
decimal=12)
numpy.testing.assert_almost_equal((theta5 / 2.).rad,
numpy.pi / 8.,
decimal=12)
numpy.testing.assert_almost_equal(theta3 / galsim.radians,
numpy.pi / 4.,
decimal=12)
numpy.testing.assert_almost_equal(theta1 / galsim.hours, 3., decimal=12)
numpy.testing.assert_almost_equal(galsim.hours / galsim.arcmin,
15 * 60,
decimal=12)
# Check picklability
do_pickle(galsim.radians)
do_pickle(galsim.degrees)
do_pickle(galsim.hours)
do_pickle(galsim.arcmin)
do_pickle(galsim.arcsec)
do_pickle(gradians)
do_pickle(theta1)
do_pickle(theta2)
do_pickle(theta3)
do_pickle(theta4)
do_pickle(theta5)
do_pickle(theta6)
do_pickle(theta7)
do_pickle(theta8)