* Bug #9789 fixed - Optimization: optim printing its termination status
[scilab.git] / scilab / modules / optimization / tests / nonreg_tests / bug_3161.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2008 - INRIA - Allan CORNET
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- CLI SHELL MODE -->
8 // <-- Non-regression test for bug 3161 -->
9 //
10 // <-- Bugzilla URL -->
11 // http://bugzilla.scilab.org/show_bug.cgi?id=3161
12 //
13 // <-- Short Description -->
14 // datafit crashes
15 function y=FF(x,p),y=p(1)*(x-p(2))+p(3)*x.*x,endfunction
16 X=[];Y=[];
17 pg=[34;12;14]; //parameter used to generate data
18 for x=0:.1:3, Y=[Y,FF(x,pg)+100*(rand()-.5)];X=[X,x];end
19 Z=[Y;X];
20 //The criterion function
21 function e=G(p,z)
22    y=z(1),x=z(2);
23    e=y-FF(x,p);
24 endfunction
25 //Solve the problem
26 p0=[3;5;10];
27 p = [];
28 err = [];
29 [p,err]=datafit(G,Z,p0);
30 Optim stops: maximum number of calls to f is reached.
31
32 if p == [] then bugmes();quit;end;
33 if err == [] then bugmes();quit;end;