def test_1d_aps_slope_error(self): mirror_length = 200.0 step = 1.0 random_seed = 898882 rms = 1e-7 # units mm wName_x, wName = create_simulated_1D_file_APS( mirror_length=mirror_length, step=step, random_seed=random_seed, error_type=SLOPE_ERROR, rms=rms) slopes = numpy.gradient(wName, wName_x[1] - wName_x[0]) print( "test_1d: test function: %s, SLOPES Stdev: input=%g, obtained=%g" % ("test_1d_aps_figure_error", rms, slopes.std())) assert numpy.abs(rms - slopes.std()) < 0.01 * numpy.abs(rms) if do_plot: from srxraylib.plot.gol import plot plot(wName_x, slopes, title="test_1d_aps_slope_error", xtitle="Y", ytitle="slopes Z'")
def test_1d_aps_figure_error(self): mirror_length = 200.0 step = 1.0 random_seed = 898882 rms = 1e-7 # units mm wName_x, wName = create_simulated_1D_file_APS( mirror_length=mirror_length, step=step, random_seed=random_seed, error_type=FIGURE_ERROR, rms=rms) print("%s, HEIGHTS Stdev: input=%g, obtained=%g" % ("test_1d_aps_figure_error", rms, wName.std())) assert numpy.abs(rms - wName.std()) < 0.01 * numpy.abs(rms) if do_plot: from srxraylib.plot.gol import plot plot(wName_x, wName, title="test_1d_aps_figure_error", xtitle="Y", ytitle="heights profile Z")
def test_1d_aps_figure_error(self): mirror_length=200.0 step=1.0 random_seed=898882 rms=1e-7 # units mm wName_x,wName = create_simulated_1D_file_APS(mirror_length=mirror_length,step=step, random_seed=random_seed, error_type=FIGURE_ERROR, rms=rms) print("%s, HEIGHTS Stdev: input=%g, obtained=%g"%("test_1d_aps_figure_error",rms,wName.std())) assert numpy.abs( rms - wName.std() ) < 0.01 * numpy.abs(rms) if do_plot: from srxraylib.plot.gol import plot plot(wName_x,wName,title="test_1d_aps_figure_error",xtitle="Y",ytitle="heights profile Z")
def test_1d_aps_slope_error(self): mirror_length=200.0 step=1.0 random_seed=898882 rms=1e-7 # units mm wName_x,wName = create_simulated_1D_file_APS(mirror_length=mirror_length,step=step, random_seed=random_seed, error_type=SLOPE_ERROR, rms=rms) slopes = numpy.gradient(wName, wName_x[1]-wName_x[0]) print("test_1d: test function: %s, SLOPES Stdev: input=%g, obtained=%g"%("test_1d_aps_figure_error",rms,slopes.std())) assert numpy.abs( rms - slopes.std() ) < 0.01 * numpy.abs(rms) if do_plot: from srxraylib.plot.gol import plot plot(wName_x,slopes,title="test_1d_aps_slope_error",xtitle="Y",ytitle="slopes Z'")