def test_strike_shear():
ss = strike_shear(strike=strike,
dip=dip,
rho=rho,
g=g,
mxx=mxx,
myy=myy,
mxy=mxy,
mzz=mzz,
mxz=mxz,
myz=myz,
txx=txx,
tyy=tyy,
txy=txy,
depth=depth)
ss = np.int(np.round(ss))
assert ss == -23234845, 'dip shear stress'
lms_reps = np.tile(lms_col_array, [n_trials, 1])
search_df[lms_fill_cols] = lms_reps
search_df.depth *= 1000
del lms_reps
search_df[['mxx', 'myy', 'mxy', 'mzz', 'mxz', 'myz']] *= 1e6
# OK, now let's do some calculationsi
print('calculating fault stresses')
search_df['tau_s'] = scv.strike_shear(strike=search_df.strike,
dip=search_df.dip,
rho=rho,
g=g,
mxx=search_df.mxx,
myy=search_df.myy,
mxy=search_df.mxy,
mzz=search_df.mzz,
mxz=search_df.mxz,
myz=search_df.myz,
txx=search_df.txx,
tyy=search_df.tyy,
txy=search_df.txy,
depth=search_df.depth)
search_df['tau_d'] = scv.dip_shear(strike=search_df.strike,
dip=search_df.dip,
rho=rho,
g=g,
mxx=search_df.mxx,
myy=search_df.myy,
mxy=search_df.mxy,
mzz=search_df.mzz,
lms_col_array = lms[lms_copy_cols].values
lms_reps = np.tile(lms_col_array, [n_trials, 1])
search_df[lms_fill_cols] = lms_reps
search_df.depth *= 1000
del lms_reps
search_df[['mxx', 'myy', 'mxy', 'mzz', 'mxz', 'myz']] *= 1e6
# OK, now let's do some calculationsi
print('calculating fault stresses')
search_df['tau_s'] = scv.strike_shear(strike = search_df.strike,
dip=search_df.dip, rho=rho, g=g,
mxx=search_df.mxx, myy=search_df.myy,
mxy=search_df.mxy, mzz=search_df.mzz,
mxz=search_df.mxz, myz=search_df.myz,
txx=search_df.txx, tyy=search_df.tyy,
txy=search_df.txy, depth=search_df.depth)
search_df['tau_d'] = scv.dip_shear(strike = search_df.strike,
dip=search_df.dip, rho=rho, g=g,
mxx=search_df.mxx, myy=search_df.myy,
mxy=search_df.mxy, mzz=search_df.mzz,
mxz=search_df.mxz, myz=search_df.myz,
txx=search_df.txx, tyy=search_df.tyy,
txy=search_df.txy, depth=search_df.depth)
print("dropping columns that won't be used more")
search_df.drop(labels=['mxx', 'myy', 'mxy', 'mzz', 'mxz', 'myz', 'depth',
'strike', 'dip'], axis=1, inplace=True)