def test_measure():
delta, ampl, diag = rheo.measure_response(toolpos(0.1, 10), fpc)
assert abs(delta - 10) < 0.01
assert abs(ampl - 0.1) < 0.01
# Test transient
delta, ampl, diag = rheo.measure_response(toolpos(0.1, 10), fpc, t_trans=0)
assert diag['cycles_discarded'] == 0
assert diag['cycles_after_transient'] == 10000 / fpc
assert diag['frames_to_discard'] % fpc == 0
assert diag['n'] % fpc == 0
delta, ampl, diag = rheo.measure_response(toolpos(0.1, 10),
fpc,
t_trans=10)
assert diag['cycles_discarded'] == 2
assert diag['cycles_after_transient'] == 10000 / fpc - 2
assert diag['frames_to_discard'] % fpc == 0
assert diag['n'] % fpc == 0
delta, ampl, diag = rheo.measure_response(toolpos(0.1, 10),
fpc,
t_trans=5.1)
assert diag['cycles_discarded'] == 2
assert diag['cycles_after_transient'] == 10000 / fpc - 2
assert diag['frames_to_discard'] % fpc == 0
assert diag['n'] % fpc == 0
def test_rheology():
delta, ampl_px, diag = rheo.measure_response(toolpos(10, 10), fpc)
rh = rheo.measure_rheology(params, delta, ampl_px, mpp, fpc, 0.2, 1)
assert abs(rh['delta'] - 10) < 0.01
assert np.isclose(rh['ampl_m'], 10 * mpp)
assert np.isclose(rh['Gp'], 1.728187076218136e-05)
assert np.isclose(rh['Gpp'], 2.9430990972078236e-06)
assert np.isclose(rh['Bo'], 4818.3962176983159)
def test_rheology():
delta, ampl_px, diag = rheo.measure_response(toolpos(10, 10), fpc)
rh = rheo.measure_rheology(params, delta, ampl_px, mpp, fpc, 0.2, 1)
assert abs(rh['delta'] - 10) < 0.01
assert np.isclose(rh['ampl_m'], 10 * mpp)
assert np.isclose(rh['Gp'], 1.728187076218136e-05)
assert np.isclose(rh['Gpp'], 2.9430990972078236e-06)
assert np.isclose(rh['Bo'], 4818.3962176983159)
def test_measure():
delta, ampl, diag = rheo.measure_response(toolpos(0.1, 10), fpc)
assert abs(delta - 10) < 0.01
assert abs(ampl - 0.1) < 0.01
def test_measure():
delta, ampl, diag = rheo.measure_response(toolpos(0.1, 10), fpc)
assert abs(delta - 10) < 0.01
assert abs(ampl - 0.1) < 0.01
