fixed script, added graphs

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 

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 

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"

dokumentation/evolution1d/images.sh

@@ -1,17 +1,27 @@
 #!/bin/bash
 
-if [[ $? -ne 1 ]]; then
+if [[ $# -ne 1 ]]; then
     echo "usage: $0 <DATA.csv>"
 else
-    $data=$1;
+    data="$1";
-    $png=`echo $1 | sed -s "s/\.csv$//"` # strip ending
+    png="`echo $1 | sed -s "s/\.csv$//"`" # strip ending
-    (gnuplot <<EOD
+    (cat <<EOD
+set datafile separator ","
 f(x)=a*x+b
 fit f(x) "$data" every ::1 using 2:5 via a,b
 set terminal png
 set output "${png}_regularity-vs-steps.png"
 plot "$data" every ::1 using 2:5 title "regularity vs. steps", f(x) lc rgb "black"
+g(x)=aa*x+bb
+fit g(x) "$data" every ::1 using 4:5 via aa,bb
+set output "${png}_improvement-vs-steps.png"
+plot "$data" every ::1 using 4:5 title "improvement potential vs. steps", g(x) lc rgb "black"
+h(x)=aaa*x+bbb
+fit h(x) "$data" every ::1 using 4:6 via aaa,bbb
+set output "${png}_improvement-vs-evo-error.png"
+plot "$data" every ::1 using 4:6 title "improvement potential vs. evolution error", h(x) lc rgb "black"
 EOD
-) 2> "${png}.gnuplot.log"
+) > "${png}.gnuplot.script" 
+    gnuplot "${png}.gnuplot.script" 2> "${png}.gnuplot.log"
     mv fit.log "${png}.gnuplot.fit.log"
 fi