diff --git a/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.fit.log b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.fit.log new file mode 100644 index 0000000..dbc0524 --- /dev/null +++ b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.fit.log @@ -0,0 +1,138 @@ + + +******************************************************************************* +Wed Aug 30 21:11:12 2017 + + +FIT: data read from "20170830-evolution1D_5x5_100Times.csv" every ::1 using 2:5 + format = x:z + #datapoints = 100 + residuals are weighted equally (unit weight) + +function used for fitting: f(x) +fitted parameters initialized with current variable values + + + + Iteration 0 + WSSR : 4.26014e+06 delta(WSSR)/WSSR : 0 + delta(WSSR) : 0 limit for stopping : 1e-05 + lambda : 0.707238 + +initial set of free parameter values + +a = 1 +b = 1 + +After 5 iterations the fit converged. +final sum of squares of residuals : 226345 +rel. change during last iteration : -1.70837e-09 + +degrees of freedom (FIT_NDF) : 98 +rms of residuals (FIT_STDFIT) = sqrt(WSSR/ndf) : 48.0587 +variance of residuals (reduced chisquare) = WSSR/ndf : 2309.64 + +Final set of parameters Asymptotic Standard Error +======================= ========================== + +a = -6597.96 +/- 2528 (38.31%) +b = 328.15 +/- 48.77 (14.86%) + + +correlation matrix of the fit parameters: + + a b +a 1.000 +b -0.995 1.000 + + +******************************************************************************* +Wed Aug 30 21:11:12 2017 + + +FIT: data read from "20170830-evolution1D_5x5_100Times.csv" every ::1 using 4:5 + format = x:z + #datapoints = 100 + residuals are weighted equally (unit weight) + +function used for fitting: g(x) +fitted parameters initialized with current variable values + + + + Iteration 0 + WSSR : 4.22366e+06 delta(WSSR)/WSSR : 0 + delta(WSSR) : 0 limit for stopping : 1e-05 + lambda : 0.965827 + +initial set of free parameter values + +aa = 1 +bb = 1 + +After 5 iterations the fit converged. +final sum of squares of residuals : 200034 +rel. change during last iteration : -5.28884e-11 + +degrees of freedom (FIT_NDF) : 98 +rms of residuals (FIT_STDFIT) = sqrt(WSSR/ndf) : 45.1792 +variance of residuals (reduced chisquare) = WSSR/ndf : 2041.16 + +Final set of parameters Asymptotic Standard Error +======================= ========================== + +aa = 3323.53 +/- 732.3 (22.03%) +bb = -2890.68 +/- 681.3 (23.57%) + + +correlation matrix of the fit parameters: + + aa bb +aa 1.000 +bb -1.000 1.000 + + +******************************************************************************* +Wed Aug 30 21:11:12 2017 + + +FIT: data read from "20170830-evolution1D_5x5_100Times.csv" every ::1 using 4:6 + format = x:z + #datapoints = 100 + residuals are weighted equally (unit weight) + +function used for fitting: h(x) +fitted parameters initialized with current variable values + + + + Iteration 0 + WSSR : 4.73621e+06 delta(WSSR)/WSSR : 0 + delta(WSSR) : 0 limit for stopping : 1e-05 + lambda : 0.965827 + +initial set of free parameter values + +aaa = 1 +bbb = 1 + +After 5 iterations the fit converged. +final sum of squares of residuals : 6.41022 +rel. change during last iteration : -1.67651e-06 + +degrees of freedom (FIT_NDF) : 98 +rms of residuals (FIT_STDFIT) = sqrt(WSSR/ndf) : 0.255755 +variance of residuals (reduced chisquare) = WSSR/ndf : 0.0654104 + +Final set of parameters Asymptotic Standard Error +======================= ========================== + +aaa = -3139.34 +/- 4.145 (0.132%) +bbb = 3139.48 +/- 3.857 (0.1228%) + + +correlation matrix of the fit parameters: + + aaa bbb +aaa 1.000 +bbb -1.000 1.000 diff --git a/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.log b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.log new file mode 100644 index 0000000..14d8cb5 --- /dev/null +++ b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.log @@ -0,0 +1,261 @@ + + + Iteration 0 + WSSR : 4.26014e+06 delta(WSSR)/WSSR : 0 + delta(WSSR) : 0 limit for stopping : 1e-05 + lambda : 0.707238 + +initial set of free parameter values + +a = 1 +b = 1 +/ + + Iteration 1 + WSSR : 242181 delta(WSSR)/WSSR : -16.5907 + delta(WSSR) : -4.01796e+06 limit for stopping : 1e-05 + lambda : 0.0707238 + +resultant parameter values + +a = 0.0614959 +b = 200.471 +/ + + Iteration 2 + WSSR : 240032 delta(WSSR)/WSSR : -0.00895106 + delta(WSSR) : -2148.54 limit for stopping : 1e-05 + lambda : 0.00707238 + +resultant parameter values + +a = -444.463 +b = 210.003 +/ + + Iteration 3 + WSSR : 226547 delta(WSSR)/WSSR : -0.0595252 + delta(WSSR) : -13485.3 limit for stopping : 1e-05 + lambda : 0.000707238 + +resultant parameter values + +a = -5849.75 +b = 313.784 +/ + + Iteration 4 + WSSR : 226345 delta(WSSR)/WSSR : -0.000894059 + delta(WSSR) : -202.365 limit for stopping : 1e-05 + lambda : 7.07238e-05 + +resultant parameter values + +a = -6596.93 +b = 328.13 +/ + + Iteration 5 + WSSR : 226345 delta(WSSR)/WSSR : -1.70837e-09 + delta(WSSR) : -0.00038668 limit for stopping : 1e-05 + lambda : 7.07238e-06 + +resultant parameter values + +a = -6597.96 +b = 328.15 + +After 5 iterations the fit converged. +final sum of squares of residuals : 226345 +rel. change during last iteration : -1.70837e-09 + +degrees of freedom (FIT_NDF) : 98 +rms of residuals (FIT_STDFIT) = sqrt(WSSR/ndf) : 48.0587 +variance of residuals (reduced chisquare) = WSSR/ndf : 2309.64 + +Final set of parameters Asymptotic Standard Error +======================= ========================== + +a = -6597.96 +/- 2528 (38.31%) +b = 328.15 +/- 48.77 (14.86%) + + +correlation matrix of the fit parameters: + + a b +a 1.000 +b -0.995 1.000 + + + Iteration 0 + WSSR : 4.22366e+06 delta(WSSR)/WSSR : 0 + delta(WSSR) : 0 limit for stopping : 1e-05 + lambda : 0.965827 + +initial set of free parameter values + +aa = 1 +bb = 1 +/ + + Iteration 1 + WSSR : 239502 delta(WSSR)/WSSR : -16.6352 + delta(WSSR) : -3.98416e+06 limit for stopping : 1e-05 + lambda : 0.0965827 + +resultant parameter values + +aa = 107.081 +bb = 100.912 +/ + + Iteration 2 + WSSR : 226540 delta(WSSR)/WSSR : -0.0572189 + delta(WSSR) : -12962.4 limit for stopping : 1e-05 + lambda : 0.00965827 + +resultant parameter values + +aa = 684.765 +bb = -435.571 +/ + + Iteration 3 + WSSR : 200084 delta(WSSR)/WSSR : -0.132222 + delta(WSSR) : -26455.4 limit for stopping : 1e-05 + lambda : 0.000965827 + +resultant parameter values + +aa = 3208.17 +bb = -2783.34 +/ + + Iteration 4 + WSSR : 200034 delta(WSSR)/WSSR : -0.000253269 + delta(WSSR) : -50.6622 limit for stopping : 1e-05 + lambda : 9.65827e-05 + +resultant parameter values + +aa = 3323.48 +bb = -2890.63 +/ + + Iteration 5 + WSSR : 200034 delta(WSSR)/WSSR : -5.28884e-11 + delta(WSSR) : -1.05794e-05 limit for stopping : 1e-05 + lambda : 9.65827e-06 + +resultant parameter values + +aa = 3323.53 +bb = -2890.68 + +After 5 iterations the fit converged. +final sum of squares of residuals : 200034 +rel. change during last iteration : -5.28884e-11 + +degrees of freedom (FIT_NDF) : 98 +rms of residuals (FIT_STDFIT) = sqrt(WSSR/ndf) : 45.1792 +variance of residuals (reduced chisquare) = WSSR/ndf : 2041.16 + +Final set of parameters Asymptotic Standard Error +======================= ========================== + +aa = 3323.53 +/- 732.3 (22.03%) +bb = -2890.68 +/- 681.3 (23.57%) + + +correlation matrix of the fit parameters: + + aa bb +aa 1.000 +bb -1.000 1.000 + + + Iteration 0 + WSSR : 4.73621e+06 delta(WSSR)/WSSR : 0 + delta(WSSR) : 0 limit for stopping : 1e-05 + lambda : 0.965827 + +initial set of free parameter values + +aaa = 1 +bbb = 1 +/ + + Iteration 1 + WSSR : 40116.1 delta(WSSR)/WSSR : -117.063 + delta(WSSR) : -4.69609e+06 limit for stopping : 1e-05 + lambda : 0.0965827 + +resultant parameter values + +aaa = 101.443 +bbb = 123.174 +/ + + Iteration 2 + WSSR : 26932.3 delta(WSSR)/WSSR : -0.489513 + delta(WSSR) : -13183.7 limit for stopping : 1e-05 + lambda : 0.00965827 + +resultant parameter values + +aaa = -479.76 +bbb = 665.003 +/ + + Iteration 3 + WSSR : 57.8749 delta(WSSR)/WSSR : -464.354 + delta(WSSR) : -26874.5 limit for stopping : 1e-05 + lambda : 0.000965827 + +resultant parameter values + +aaa = -3023.07 +bbb = 3031.3 +/ + + Iteration 4 + WSSR : 6.41023 delta(WSSR)/WSSR : -8.02852 + delta(WSSR) : -51.4647 limit for stopping : 1e-05 + lambda : 9.65827e-05 + +resultant parameter values + +aaa = -3139.29 +bbb = 3139.43 +/ + + Iteration 5 + WSSR : 6.41022 delta(WSSR)/WSSR : -1.67651e-06 + delta(WSSR) : -1.07468e-05 limit for stopping : 1e-05 + lambda : 9.65827e-06 + +resultant parameter values + +aaa = -3139.34 +bbb = 3139.48 + +After 5 iterations the fit converged. +final sum of squares of residuals : 6.41022 +rel. change during last iteration : -1.67651e-06 + +degrees of freedom (FIT_NDF) : 98 +rms of residuals (FIT_STDFIT) = sqrt(WSSR/ndf) : 0.255755 +variance of residuals (reduced chisquare) = WSSR/ndf : 0.0654104 + +Final set of parameters Asymptotic Standard Error +======================= ========================== + +aaa = -3139.34 +/- 4.145 (0.132%) +bbb = 3139.48 +/- 3.857 (0.1228%) + + +correlation matrix of the fit parameters: + + aaa bbb +aaa 1.000 +bbb -1.000 1.000 diff --git a/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.script b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.script new file mode 100644 index 0000000..1d308a7 --- /dev/null +++ b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times.gnuplot.script @@ -0,0 +1,14 @@ +set datafile separator "," +f(x)=a*x+b +fit f(x) "20170830-evolution1D_5x5_100Times.csv" every ::1 using 2:5 via a,b +set terminal png +set output "20170830-evolution1D_5x5_100Times_regularity-vs-steps.png" +plot "20170830-evolution1D_5x5_100Times.csv" every ::1 using 2:5 title "regularity vs. steps", f(x) lc rgb "black" +g(x)=aa*x+bb +fit g(x) "20170830-evolution1D_5x5_100Times.csv" every ::1 using 4:5 via aa,bb +set output "20170830-evolution1D_5x5_100Times_improvement-vs-steps.png" +plot "20170830-evolution1D_5x5_100Times.csv" every ::1 using 4:5 title "improvement potential vs. steps", g(x) lc rgb "black" +h(x)=aaa*x+bbb +fit h(x) "20170830-evolution1D_5x5_100Times.csv" every ::1 using 4:6 via aaa,bbb +set output "20170830-evolution1D_5x5_100Times_improvement-vs-evo-error.png" +plot "20170830-evolution1D_5x5_100Times.csv" every ::1 using 4:6 title "improvement potential vs. evolution error", h(x) lc rgb "black" diff --git a/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_improvement-vs-evo-error.png b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_improvement-vs-evo-error.png new file mode 100644 index 0000000..3f5dd22 Binary files /dev/null and b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_improvement-vs-evo-error.png differ diff --git a/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_improvement-vs-steps.png b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_improvement-vs-steps.png new file mode 100644 index 0000000..d896a19 Binary files /dev/null and b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_improvement-vs-steps.png differ diff --git a/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_regularity-vs-steps.png b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_regularity-vs-steps.png new file mode 100644 index 0000000..1f5e082 Binary files /dev/null and b/dokumentation/evolution1d/20170830-evolution1D_5x5_100Times_regularity-vs-steps.png differ diff --git a/dokumentation/evolution1d/images.sh b/dokumentation/evolution1d/images.sh index 20b6058..f732cb1 100755 --- a/dokumentation/evolution1d/images.sh +++ b/dokumentation/evolution1d/images.sh @@ -1,17 +1,27 @@ #!/bin/bash -if [[ $? -ne 1 ]]; then +if [[ $# -ne 1 ]]; then echo "usage: $0 " else - $data=$1; - $png=`echo $1 | sed -s "s/\.csv$//"` # strip ending - (gnuplot < "${png}.gnuplot.log" +) > "${png}.gnuplot.script" + gnuplot "${png}.gnuplot.script" 2> "${png}.gnuplot.log" mv fit.log "${png}.gnuplot.fit.log" fi