def test_relradpow_external_with_custom_radius(self): assert_allclose( relradpow(4.0, 10, reference_radius=2.0, is_internal=False), [0.5, 1, 2, 4, 8, 16, 32, 64, 128, 256, 512] ) assert_allclose( relradpow(1.1, 256, reference_radius=1.0, is_internal=False), self.reference_external(1.1, 256) )
def test_relradpow_internal_with_custom_radius(self): assert_allclose( relradpow(4.0, 10, reference_radius=2.0, is_internal=True), [ 0.25, 0.125, 0.0625, 0.03125, 0.015625, 0.0078125, 0.00390625, 0.001953125, 0.0009765625, 0.00048828125, 0.000244140625, ] ) assert_allclose( relradpow(1.1, 256, reference_radius=1.0, is_internal=True), self.reference_internal(1.1, 256) )
def test_relradpow_internal(self): assert_allclose( relradpow(2*6371.2, 10, is_internal=True), [ 0.25, 0.125, 0.0625, 0.03125, 0.015625, 0.0078125, 0.00390625, 0.001953125, 0.0009765625, 0.00048828125, 0.000244140625, ] ) assert_allclose( relradpow(1.1*6371.2, 256, is_internal=True), self.reference_internal(1.1, 256) )
def test_relradpow_zero_degree_internal(self): assert_allclose(relradpow(2.0, 0, 1.0, is_internal=True), [0.25])
def test_relradpow_zero_degree_external(self): assert_allclose(relradpow(2.0, 0, 1.0, is_internal=False), [0.5])
def get_series(cls, degree, latitude, longitude, radius): rad_series = relradpow(radius, degree, is_internal=cls.is_internal) sin_series, cos_series = lonsincos(longitude, degree) p_series, dp_series = legendre(latitude, degree) return rad_series, sin_series, cos_series, p_series, dp_series
def get_series(cls, degree, latitude, longitude, radius): rad_series = relradpow(radius, degree, is_internal=cls.is_internal) sin_series, cos_series = lonsincos(longitude, degree) p_series, dp_series = legendre(latitude, degree) return rad_series, sin_series, cos_series, p_series, dp_series