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_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])