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