Deleted vectorized computation feature. Deleted neldermead_contour. Fixed the demos.
[scilab.git] / scilab / modules / optimization / demos / neldermead / nmplot_rosenbrock.sce
index f180860..88a0430 100644 (file)
@@ -8,10 +8,13 @@
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 
+mprintf("Illustrates the Nelder-Mead algorithm on Rosenbrock test case.\n");
+mprintf("Defining Rosenbrock function...\n");
 function y = rosenbrock (x)
   y = 100*(x(2)-x(1)^2)^2 + (1-x(1))^2;
 endfunction
 
+mprintf("Creating nmplot object...\n");
 nm = nmplot_new ();
 nm = nmplot_configure(nm,"-numberofvariables",2);
 nm = nmplot_configure(nm,"-function",rosenbrock);
@@ -23,7 +26,7 @@ nm = nmplot_configure(nm,"-tolxrelative",10*%eps);
 nm = nmplot_configure(nm,"-simplex0method","axes");
 nm = nmplot_configure(nm,"-simplex0length",1.0);
 nm = nmplot_configure(nm,"-method","variable");
-nm = nmplot_configure(nm,"-verbose",1);
+//nm = nmplot_configure(nm,"-verbose",1);
 nm = nmplot_configure(nm,"-verbosetermination",0);
 //
 // Setup output files
@@ -35,25 +38,28 @@ nm = nmplot_configure(nm,"-sigmafn","rosenbrock.history.sigma.txt");
 //
 // Perform optimization
 //
+mprintf("Searching for minimum...\n");
 nm = nmplot_search(nm);
+nmplot_display(nm);
 // Plot the contours of the cost function and the simplex history
+mprintf("Plotting contour...\n");
 [nm , xdata , ydata , zdata ] = nmplot_contour ( nm , xmin = -2.0 , xmax = 2.0 , ymin = -1.0 , ymax = 2.5 , nx = 100 , ny = 100 );
 f = scf();
-contour ( xdata , ydata , zdata , 30 )
+contour ( xdata , ydata , zdata , [1 10 100 500 1000 2000] )
 nmplot_simplexhistory ( nm );
-xs2png(0,"rosenbrock.history.simplex.png");
+mprintf("Plotting history of fbar...\n");
 f = scf();
 nmplot_historyplot ( nm , "rosenbrock.history.fbar.txt" , ...
   mytitle = "Function Value Average" , myxlabel = "Iterations" );
-xs2png(1,"rosenbrock.history.fbar.png");
+mprintf("Plotting history of fopt...\n");
 f = scf();
 nmplot_historyplot ( nm , "rosenbrock.history.fopt.txt" , ...
   mytitle = "Minimum Function Value" , myxlabel = "Iterations" );
-xs2png(2,"rosenbrock.history.fopt.png");
+mprintf("Plotting history of sigma...\n");
 f = scf();
 nmplot_historyplot ( nm , "rosenbrock.history.sigma.txt" , ...
   mytitle = "Maximum Oriented length" , myxlabel = "Iterations" );
-xs2png(3,"rosenbrock.history.sigma.png");
+mprintf("Cleaning up...\n");
 deletefile("rosenbrock.history.simplex.txt");
 deletefile("rosenbrock.history.fbar.txt");
 deletefile("rosenbrock.history.fopt.txt");