def test_pb_input(self):
for plunge in range(0, 100, 10):
for bearing in range(0, 370, 10):
plunges = plunge + np.array([-10, -5, 0, 5, 10])
bearings = bearing + np.array([0, 0, 0, 0, 0])
plu, azi, vals = mplstereonet.eigenvectors(plunges, bearings,
measurement='lines')
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
def test_sd_input(self):
for strike in range(0, 100, 10):
for dip in range(0, 370, 10):
dips = dip + np.array([-10, -5, 0, 5, 10])
strikes = strike + np.array([0, 0, 0, 0, 0])
plu, azi, vals = mplstereonet.eigenvectors(strikes, dips)
plunge, bearing = mplstereonet.pole2plunge_bearing(strike, dip)
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
def test_sd_input(self):
for strike in range(0, 100, 10):
for dip in range(0, 370, 10):
dips = dip + np.array([-10, -5, 0, 5, 10])
strikes = strike + np.array([0, 0, 0, 0, 0])
plu, azi, vals = mplstereonet.eigenvectors(strikes, dips)
plunge, bearing = mplstereonet.pole2plunge_bearing(strike, dip)
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
def test_pb_input(self):
for plunge in range(0, 100, 10):
for bearing in range(0, 370, 10):
plunges = plunge + np.array([-10, -5, 0, 5, 10])
bearings = bearing + np.array([0, 0, 0, 0, 0])
plu, azi, vals = mplstereonet.eigenvectors(plunges,
bearings,
measurement='lines')
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
def test_rake_input(self):
for strike in range(0, 100, 10):
for dip in range(0, 370, 10):
rakes = 90 + np.array([-10, -5, 0, 5, 10])
dips = dip * np.ones_like(rakes)
strikes = strike * np.ones_like(rakes)
plu, azi, vals = mplstereonet.eigenvectors(strikes, dips, rakes,
measurement='rakes')
plunge, bearing = dip, strike + 90
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
def test_radian_input(self):
for lat in range(0, 100, 10):
for lon in range(0, 370, 10):
lats = lat + np.array([-10, -5, 0, 5, 10])
lons = lon * np.ones_like(lats)
lats, lons = np.radians(lats), np.radians(lons)
plu, azi, vals = mplstereonet.eigenvectors(lons, lats,
measurement='radians')
plunge, bearing = mplstereonet.geographic2plunge_bearing(
np.radians(lon), np.radians(lat))
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
def test_radian_input(self):
for lat in range(0, 100, 10):
for lon in range(0, 370, 10):
lats = lat + np.array([-10, -5, 0, 5, 10])
lons = lon * np.ones_like(lats)
lats, lons = np.radians(lats), np.radians(lons)
plu, azi, vals = mplstereonet.eigenvectors(
lons, lats, measurement='radians')
plunge, bearing = mplstereonet.geographic2plunge_bearing(
np.radians(lon), np.radians(lat))
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
def test_rake_input(self):
for strike in range(0, 100, 10):
for dip in range(0, 370, 10):
rakes = 90 + np.array([-10, -5, 0, 5, 10])
dips = dip * np.ones_like(rakes)
strikes = strike * np.ones_like(rakes)
plu, azi, vals = mplstereonet.eigenvectors(strikes,
dips,
rakes,
measurement='rakes')
plunge, bearing = dip, strike + 90
self.compare_plungebearing(plunge, bearing, plu[0], azi[0])
assert np.allclose(vals, [1.09433342, 1.67776947e-02, 0])
