remove ref to quapro & linpro
Allan Cornet [Tue, 1 Jul 2008 18:39:39 +0000 (18:39 +0000)]
22 files changed:
scilab/modules/optimization/tests/nonreg_tests/bug_2330.dia.ref
scilab/modules/optimization/tests/nonreg_tests/bug_2330.tst
scilab/modules/optimization/tests/nonreg_tests/bug_244.dia.ref [new file with mode: 0644]
scilab/modules/optimization/tests/nonreg_tests/bug_244.tst
scilab/modules/optimization/tests/nonreg_tests/bug_2572.dia.ref [deleted file]
scilab/modules/optimization/tests/nonreg_tests/bug_2572.tst [deleted file]
scilab/modules/optimization/tests/nonreg_tests/bug_2573.dia.ref
scilab/modules/optimization/tests/nonreg_tests/bug_2573.tst
scilab/modules/optimization/tests/nonreg_tests/bug_332.dia.ref [deleted file]
scilab/modules/optimization/tests/nonreg_tests/bug_332.tst [deleted file]
scilab/modules/optimization/tests/nonreg_tests/bug_631.dia.ref [deleted file]
scilab/modules/optimization/tests/nonreg_tests/bug_631.tst [deleted file]
scilab/modules/optimization/tests/nonreg_tests/bug_632.dia.ref
scilab/modules/optimization/tests/nonreg_tests/bug_632.tst
scilab/modules/optimization/tests/unit_tests/lsqrsolve.dia.ref
scilab/modules/optimization/tests/unit_tests/lsqrsolve.tst
scilab/modules/optimization/tests/unit_tests/qld.dia.ref
scilab/modules/optimization/tests/unit_tests/qld.tst
scilab/modules/optimization/tests/unit_tests/qpsolve.dia.ref
scilab/modules/optimization/tests/unit_tests/qpsolve.tst
scilab/modules/optimization/tests/unit_tests/quapro.dia.ref [deleted file]
scilab/modules/optimization/tests/unit_tests/quapro.tst [deleted file]

index d6362e5..04a2151 100644 (file)
@@ -1,3 +1,9 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2007-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
 // <-- Non-regression test for bug 2330 -->
 //
 // <-- Bugzilla URL -->
@@ -7,9 +13,6 @@
 //    datafit does not work equivalently on WinXP and Linux. With a given dataset and
 //    same routines centered on datafit function, it works perfectly on Linux and
 //    partially (some data are fitted some others not) on WinXP.
-// Serge Steer - Scilab Project
-// Copyright INRIA
-// 30 apr 2007
 //build the data to fit
 //---------------------
 function Xcalc=biexp(p,t)
@@ -26,4 +29,4 @@ function e=myerf(p,X,t),e=X-biexp(p,t),endfunction
 p0=[0.01;0.001;0.01;0.001;0.1];
 //call datafit
 [pr,err]=datafit(list(myerf,t),X,p0);
-if err>=5d-6 then bugmes();quit;end
+if err>=5d-4 then bugmes();quit;end
index 5905268..7612238 100644 (file)
@@ -32,4 +32,4 @@ function e=myerf(p,X,t),e=X-biexp(p,t),endfunction
 p0=[0.01;0.001;0.01;0.001;0.1];
 //call datafit
 [pr,err]=datafit(list(myerf,t),X,p0);
-if err>=5d-6 then pause,end
+if err>=5d-4 then pause,end
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_244.dia.ref b/scilab/modules/optimization/tests/nonreg_tests/bug_244.dia.ref
new file mode 100644 (file)
index 0000000..a99af67
--- /dev/null
@@ -0,0 +1,159 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2008 - INRIA - Vincent COUVERT
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+// <-- Non-regression test for bug 244 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=244
+//
+// <-- Short Description -->
+//   recursion problems with fsolve
+// Titration of a dibase with HCl (at the beginning, a little amount of NaOH was added)
+// parameters
+//C0=2.39e-4
+C0=2e-4;
+C02=0.02;
+na0=2e-7;
+v0=0.01;
+pK1=6.9;
+pK2=6.9;
+// function to get the calculated pH
+deff('[x]=f0(pH)','x=(10^(-pH-pK1)+2*10^(-2*pH))/(10^(-pK1-pK2)+10^(-pH-pK1)+10^(-2*pH))+((v0+v)*(10^(-pH)-10^(pH-14))+na0-C02*v)/(C0*v0)');
+deff('[pH]=fpH(p)','v=p;pH=fsolve(7,f0)');
+// data
+X=[];Y=[];
+X=[
+//0    
+0.00001        
+0.00002        
+0.00003        
+0.00004        
+0.00005        
+0.000055       
+0.00006        
+0.000065       
+0.00007        
+0.000075       
+0.00008        
+0.000085       
+0.00009        
+0.000095       
+0.0001 
+0.000105       
+0.00011        
+0.000115       
+0.00012        
+0.000125       
+0.00013        
+0.000135       
+0.00014        
+0.000145       
+0.00015        
+0.000155       
+0.00016        
+0.000165       
+0.00017        
+0.000175       
+0.00018        
+0.000185       
+0.00019        
+0.000195       
+0.0002 
+0.000205       
+0.00021        
+0.000215       
+0.00022        
+0.000225       
+0.00023        
+0.000235       
+0.00024        
+0.000245       
+0.00025        
+0.000255       
+0.00026        
+0.000265       
+0.00027        
+0.000275       
+0.00028        
+0.000285       
+0.000295       
+0.000305       
+0.00032        
+0.000345       
+0.000395       
+0.000445       
+];
+Y=[
+//8.21 
+7.82   
+7.64   
+7.48   
+7.37   
+7.22   
+7.20   
+7.17   
+7.13   
+7.12   
+7.12   
+7.10   
+7.07   
+7.04   
+7.07   
+7.04   
+7.01   
+6.98   
+6.94   
+6.91   
+6.87   
+6.84   
+6.80   
+6.76   
+6.73   
+6.68   
+6.61   
+6.57   
+6.51   
+6.45   
+6.36   
+6.27   
+6.14   
+6.02   
+5.84   
+5.64   
+5.34   
+5.00   
+4.80   
+4.59   
+4.46   
+4.40   
+4.33   
+4.24   
+4.18   
+4.11   
+4.07   
+4.02   
+3.99   
+3.95   
+3.91   
+3.88   
+3.86   
+3.80   
+3.76   
+3.69   
+3.60   
+3.47   
+3.36   
+];
+// fitting
+Z=[Y;X];
+deff('e=G(p,z)','pK1=p(1),pK2=p(2),v=z(2),pHexp=z(1),e=pHexp-fpH(v)');
+[p,err]=fit_dat(G,[6;7],Z);
+// graphic part
+xbasc()
+//v=[0:1e-5:4.5e-4]
+v=X;
+fplot2d(v,fpH);
+plot2d(X,Y,[-2]);
index 2df350b..9be73b4 100644 (file)
 
 // parameters
 //C0=2.39e-4
-C0=2e-4
-C02=0.02
-na0=2e-7
-v0=0.01
+C0=2e-4;
+C02=0.02;
+na0=2e-7;
+v0=0.01;
 
-pK1=6.9
-pK2=6.9
+pK1=6.9;
+pK2=6.9;
 
 // function to get the calculated pH
-deff('[x]=f0(pH)','x=(10^(-pH-pK1)+2*10^(-2*pH))/(10^(-pK1-pK2)+10^(-pH-pK1)+10^(-2*pH))+((v0+v)*(10^(-pH)-10^(pH-14))+na0-C02*v)/(C0*v0)')
-deff('[pH]=fpH(p)','v=p;pH=fsolve(7,f0)')
+deff('[x]=f0(pH)','x=(10^(-pH-pK1)+2*10^(-2*pH))/(10^(-pK1-pK2)+10^(-pH-pK1)+10^(-2*pH))+((v0+v)*(10^(-pH)-10^(pH-14))+na0-C02*v)/(C0*v0)');
+deff('[pH]=fpH(p)','v=p;pH=fsolve(7,f0)');
 
 // data
 
@@ -94,7 +94,7 @@ X=[
 0.000345       
 0.000395       
 0.000445       
-]
+];
 
 Y=[
 //8.21 
@@ -156,17 +156,17 @@ Y=[
 3.60   
 3.47   
 3.36   
-]
+];
 
 // fitting
 Z=[Y;X];
-deff('e=G(p,z)','pK1=p(1),pK2=p(2),v=z(2),pHexp=z(1),e=pHexp-fpH(v)')
-[p,err]=fit_dat(G,[6;7],Z)
+deff('e=G(p,z)','pK1=p(1),pK2=p(2),v=z(2),pHexp=z(1),e=pHexp-fpH(v)');
+[p,err]=fit_dat(G,[6;7],Z);
 
 
 // graphic part
 xbasc()
 //v=[0:1e-5:4.5e-4]
-v=X
-fplot2d(v,fpH)
-plot2d(X,Y,[-2])
+v=X;
+fplot2d(v,fpH);
+plot2d(X,Y,[-2]);
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_2572.dia.ref b/scilab/modules/optimization/tests/nonreg_tests/bug_2572.dia.ref
deleted file mode 100644 (file)
index e3d03cf..0000000
+++ /dev/null
@@ -1,18 +0,0 @@
-// =============================================================================
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2007-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
-//
-//  This file is distributed under the same license as the Scilab package.
-// =============================================================================
-// <-- Non-regression test for bug 2572 -->
-//
-// <-- Bugzilla URL -->
-// http://bugzilla.scilab.org/show_bug.cgi?id=2572
-//
-// <-- Short Description -->
-//   quapro make scilab stop if C and b are zero
-Q=[1 0.5;0.5 1];p=[1;2];
-C=zeros(1,2);b=0;
-s1=quapro(Q,p,C,b);
-s2=quapro(Q,p,[],[]);
-if norm(s1-s2)>1d-12 then bugmes();quit;end
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_2572.tst b/scilab/modules/optimization/tests/nonreg_tests/bug_2572.tst
deleted file mode 100644 (file)
index ee24feb..0000000
+++ /dev/null
@@ -1,21 +0,0 @@
-// =============================================================================
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2007-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
-//
-//  This file is distributed under the same license as the Scilab package.
-// =============================================================================
-
-// <-- Non-regression test for bug 2572 -->
-//
-// <-- Bugzilla URL -->
-// http://bugzilla.scilab.org/show_bug.cgi?id=2572
-//
-// <-- Short Description -->
-//   quapro make scilab stop if C and b are zero
-
-Q=[1 0.5;0.5 1];p=[1;2];
-C=zeros(1,2);b=0;
-s1=quapro(Q,p,C,b);   
-s2=quapro(Q,p,[],[]);
-if norm(s1-s2)>1d-12 then pause,end
index 837ceee..74820f6 100644 (file)
@@ -25,9 +25,6 @@ ci = [0; 0];
 cs = [];
 me = 1;
 x0 = [0.5; 0.5];
-tol = 1e-8
- tol  =
-    1.000D-08  
+tol = 1e-8;
 [x_qld,lagr_qld,info] = qld(Q,p,C,b,ci,cs,me,tol);
 if norm(x_qld-[8;3]/11)>1d-10 then bugmes();quit;end
index 5b95d84..c9b3c20 100644 (file)
@@ -31,7 +31,7 @@ ci = [0; 0];
 cs = [];
 me = 1;
 x0 = [0.5; 0.5];
-tol = 1e-8
+tol = 1e-8;
 
 [x_qld,lagr_qld,info] = qld(Q,p,C,b,ci,cs,me,tol);
 if norm(x_qld-[8;3]/11)>1d-10 then pause,end
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_332.dia.ref b/scilab/modules/optimization/tests/nonreg_tests/bug_332.dia.ref
deleted file mode 100644 (file)
index 44e9918..0000000
+++ /dev/null
@@ -1,38 +0,0 @@
-// =============================================================================
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2008 - INRIA - Vincent COUVERT
-//
-//  This file is distributed under the same license as the Scilab package.
-// =============================================================================
-// <-- Non-regression test for bug 332 -->
-//
-// <-- Bugzilla URL -->
-// http://bugzilla.scilab.org/show_bug.cgi?id=332
-//
-// <-- Short Description -->
-//   When linpro is given an initial feasible solution
-//   such as the result of a previous run, the solution
-//   returned as X is incorrect.
-A_tmp  = [ -22.    0.     0.     0.     100.    0.;
-    22.  -10.    0.     0.     0.      0 ;
-    0.     10.  -22.    0.     0.      0 ;
-    0.     0.     22.  -9.5    0.      0 ;
-    0.     0.     0.   -9.5    0.      9.5 ;
-    1.     0.     1.     0.     0.    -1;
-    0.     1.     0.     0.     0.    -1;
-    0.     0.     0.     0.     1.    -1];
-b_tmp = [ - 10;
-    - 25;
-    - 35;
-    0;
-    - 100;
-    0;
-    0;
-    0];
-f=[0;0;0;0;0;1];
-[Xtmp,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,"g",7);
-[Xtmp1,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,Xtmp,7);
-if or(Xtmp-Xtmp1>1e-12) then bugmes();quit;end
-[Xtmp2,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,"v",7);
-[Xtmp3,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,Xtmp2,7);
-if or(Xtmp2-Xtmp3>1e-12) then bugmes();quit;end
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_332.tst b/scilab/modules/optimization/tests/nonreg_tests/bug_332.tst
deleted file mode 100644 (file)
index 177352b..0000000
+++ /dev/null
@@ -1,47 +0,0 @@
-// =============================================================================
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2008 - INRIA - Vincent COUVERT
-//
-//  This file is distributed under the same license as the Scilab package.
-// =============================================================================
-
-// <-- Non-regression test for bug 332 -->
-//
-// <-- Bugzilla URL -->
-// http://bugzilla.scilab.org/show_bug.cgi?id=332
-//
-// <-- Short Description -->
-//   When linpro is given an initial feasible solution
-//   such as the result of a previous run, the solution
-//   returned as X is incorrect.
-
-
-A_tmp  = [ -22.    0.     0.     0.     100.    0.;
-    22.  -10.    0.     0.     0.      0 ;
-    0.     10.  -22.    0.     0.      0 ;
-    0.     0.     22.  -9.5    0.      0 ;
-    0.     0.     0.   -9.5    0.      9.5 ;
-    1.     0.     1.     0.     0.    -1;
-    0.     1.     0.     0.     0.    -1;
-    0.     0.     0.     0.     1.    -1];
-
-b_tmp = [ - 10;
-    - 25;
-    - 35;
-    0;
-    - 100;
-    0;
-    0;
-    0];
-
-f=[0;0;0;0;0;1];
-
-[Xtmp,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,"g",7);
-[Xtmp1,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,Xtmp,7);
-
-if or(Xtmp-Xtmp1>1e-12) then pause,end
-
-[Xtmp2,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,"v",7);
-[Xtmp3,lgr,fval]=linpro(f,A_tmp,b_tmp,[],[],5,Xtmp2,7); 
-
-if or(Xtmp2-Xtmp3>1e-12) then pause,end
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_631.dia.ref b/scilab/modules/optimization/tests/nonreg_tests/bug_631.dia.ref
deleted file mode 100644 (file)
index 1f96e22..0000000
+++ /dev/null
@@ -1,39 +0,0 @@
-// <-- Non-regression test for bug 631 -->
-//
-// <-- Bugzilla URL -->
-// http://bugzilla.scilab.org/show_bug.cgi?id=631
-//
-// <-- Short Description -->
-//    linpro with initial guess does not work correctly.
-//    linpro(-[1;1], [1 -1], zeros(1,1), zeros(2,1), [1;0], 1,
-//    zeros(2,1)) returns:
-//     ans  =
-//
-//    !   1. !
-//    !   1. !,
-//    which violates the upper-bound constraint ([1;0]).
-//
-//    On the other hand, linpro without initial guess works fine.
-//    -->linpro(-[1;1], [1 -1], zeros(1,1), zeros(2,1), [1;0], 1)
-//     ans  =
-//
-//    !   0. !
-//    !   0. !
-// Copyright INRIA
-// Scilab Project - Serge Steer
-// Copyright INRIA 2006
-// Date : 10 janvier 2006
-p  = -[1;1];
-C  = [1 -1];
-b  = zeros(1,1);
-ci = zeros(2,1);
-cs = [1;0];
-me = 1;
-x0 = zeros(2,1)
- x0  =
-    0.  
-    0.  
-[x,lagr,f] = linpro(p,C,b,ci,cs,me ,x0);
-if abs(C*x)>%eps      then bugmes();quit;end
-if or(~(ci<=x&x<=cs)) then bugmes();quit;end
diff --git a/scilab/modules/optimization/tests/nonreg_tests/bug_631.tst b/scilab/modules/optimization/tests/nonreg_tests/bug_631.tst
deleted file mode 100644 (file)
index b42b1ea..0000000
+++ /dev/null
@@ -1,41 +0,0 @@
-// =============================================================================
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2006-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
-//
-//  This file is distributed under the same license as the Scilab package.
-// =============================================================================
-
-// <-- Non-regression test for bug 631 -->
-//
-// <-- Bugzilla URL -->
-// http://bugzilla.scilab.org/show_bug.cgi?id=631
-//
-// <-- Short Description -->
-//    linpro with initial guess does not work correctly.
-//    linpro(-[1;1], [1 -1], zeros(1,1), zeros(2,1), [1;0], 1, 
-//    zeros(2,1)) returns:
-//     ans  =
-//
-//    !   1. !
-//    !   1. !,
-//    which violates the upper-bound constraint ([1;0]).
-//
-//    On the other hand, linpro without initial guess works fine.
-//    -->linpro(-[1;1], [1 -1], zeros(1,1), zeros(2,1), [1;0], 1)
-//     ans  =
-//
-//    !   0. !
-//    !   0. !
-
-p  = -[1;1];
-C  = [1 -1];
-b  = zeros(1,1);
-ci = zeros(2,1);
-cs = [1;0];
-me = 1;
-x0 = zeros(2,1)
-
-[x,lagr,f] = linpro(p,C,b,ci,cs,me ,x0);
-
-if abs(C*x)>%eps      then pause,end
-if or(~(ci<=x&x<=cs)) then pause,end
index 04f8e20..fc03695 100644 (file)
@@ -1,3 +1,9 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) ????-2008 - INRIA
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
 // <-- Non-regression test for bug 632 -->
 //
 // <-- Bugzilla URL -->
@@ -39,13 +45,5 @@ function e=f1(abc,m)
        e = Y-(a*(X-b)+c*X.*X);
 endfunction
 [abc,v] = lsqrsolve([10;10;10],f1,size(X,1));
-abc
- abc  =
-    34.        
-    12.849045  
-    14.        
-norm(v)
- ans  =
-    2.814D-13  
+if norm(abc - [ 34.  ; 12.849045 ; 14.   ] ) > 0.1 then bugmes();quit;end
+if norm(v) - 2.814D-13 > 0.1 then bugmes();quit;end
index 01652ac..0154c76 100644 (file)
@@ -55,5 +55,5 @@ function e=f1(abc,m)
 endfunction
 
 [abc,v] = lsqrsolve([10;10;10],f1,size(X,1));
-abc
-norm(v)
+if norm(abc - [ 34.  ; 12.849045 ; 14.   ] ) > 0.1 then pause,end
+if norm(v) - 2.814D-13 > 0.1 then pause,end
index f6798ef..5c0fb1b 100644 (file)
@@ -13,16 +13,7 @@ function y=f1(x,m),y=a*x+b;endfunction
 if size(xsol) <> [2 1] then bugmes();quit;end
 if size(v) <> [3 1] then bugmes();quit;end
 function y=fj1(x,m),y=a;endfunction
-[xsol,v]=lsqrsolve([100;100],f1,3,fj1)
- v  =
-  - 0.4083770  
-    0.3926702  
-  - 0.0942408  
- xsol  =
-    1.5026178  
-  - 1.7015707  
+[xsol,v]=lsqrsolve([100;100],f1,3,fj1);
 // Data fitting problem
 // 1 build the data
 a=34;b=12;c=14;
index 7389436..c9350b9 100644 (file)
@@ -18,7 +18,7 @@ if size(xsol) <> [2 1] then pause,end
 if size(v) <> [3 1] then pause,end
 
 function y=fj1(x,m),y=a;endfunction
-[xsol,v]=lsqrsolve([100;100],f1,3,fj1)
+[xsol,v]=lsqrsolve([100;100],f1,3,fj1);
 
 // Data fitting problem
 // 1 build the data
index 209f4bd..37e3508 100644 (file)
@@ -22,32 +22,7 @@ p=[1;2;3;4;5;6]; Q=eye(6,6);
 C=[C1;C2] ; //
 b=[b1;b2] ;  //
 me=3;
-[x,lagr]=qld(Q,p,C,b,ci,cs,me)
- lagr  =
-  - 1.5564027  
-  - 0.1698164  
-  - 0.7054782  
-    0.3091368  
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
- x  =
-    1.7975426  
-  - 0.3381487  
-    0.1633880  
-  - 4.9884023  
-    0.6054943  
-  - 3.1155623  
+[x,lagr]=qld(Q,p,C,b,ci,cs,me);
+if size(lagr,'*') <> 17 then bugmes();quit;end
+if size(x,'*') <> 6 then bugmes();quit;end
 //Only linear constraints (1 to 4) are active (lagr(1:6)=0):
index 239c25f..5fd08ef 100644 (file)
@@ -23,7 +23,9 @@ p=[1;2;3;4;5;6]; Q=eye(6,6);
 C=[C1;C2] ; //
 b=[b1;b2] ;  //
 me=3;
-[x,lagr]=qld(Q,p,C,b,ci,cs,me)
+[x,lagr]=qld(Q,p,C,b,ci,cs,me);
+if size(lagr,'*') <> 17 then pause,end
+if size(x,'*') <> 6 then pause,end
 //Only linear constraints (1 to 4) are active (lagr(1:6)=0):
  
  
\ No newline at end of file
index 4c79721..3883a3e 100644 (file)
@@ -79,11 +79,7 @@ Q= [1 -1; -1 2] ;
 p = [-2; -6];
 C = [1 1; -1 2; 2 1];
 b = [2; 2; 3];
-lb = [0.8;0] //zeros(2,1);
- lb  =
-    0.8  
-    0.   
+lb = [0.8;0]; //zeros(2,1);
 me=0;
 x1=qpsolve(Q,p,C,b,lb,[],me);
 if or(x1-[0.8;1.2] >100*%eps) then bugmes();quit;end
index 22dce46..be9eed4 100644 (file)
@@ -90,7 +90,7 @@ Q= [1 -1; -1 2] ;
 p = [-2; -6];
 C = [1 1; -1 2; 2 1];
 b = [2; 2; 3];
-lb = [0.8;0] //zeros(2,1);
+lb = [0.8;0]; //zeros(2,1);
 me=0;
 x1=qpsolve(Q,p,C,b,lb,[],me);
 if or(x1-[0.8;1.2] >100*%eps) then pause,end
diff --git a/scilab/modules/optimization/tests/unit_tests/quapro.dia.ref b/scilab/modules/optimization/tests/unit_tests/quapro.dia.ref
deleted file mode 100644 (file)
index c0f37a4..0000000
+++ /dev/null
@@ -1,132 +0,0 @@
-// ============================================================================
-// unit tests for quapro function
-// Author : Sylvestre LEDRU
-// Scilab Project
-// Copyright INRIA 2007
-// Tests from the documentation
-// ============================================================================
-//Find x in R^6 such that:
-//C1*x = b1 (3 equality constraints i.e me=3)
-C1= [1,-1,1,0,3,1;
-    -1,0,-3,-4,5,6;
-     2,5,3,0,1,0];
-b1=[1;2;3];
-//C2*x <= b2 (2 inequality constraints)
-C2=[0,1,0,1,2,-1;
-    -1,0,2,1,1,0];
-b2=[-1;2.5];
-//with  x between ci and cs:
-ci=[-1000;-10000;0;-1000;-1000;-1000];cs=[10000;100;1.5;100;100;1000];
-//and minimize 0.5*x'*Q*x + p'*x with
-p=[1;2;3;4;5;6]; Q=eye(6,6);
-//No initial point is given;
-C=[C1;C2] ; //
-b=[b1;b2] ;  //
-me=3;
-[x,lagr,f]=quapro(Q,p,C,b,ci,cs,me)
- f  =
-  - 14.843248  
- lagr  =
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-    0.         
-  - 1.5564027  
-  - 0.1698164  
-  - 0.7054782  
-    0.3091368  
-    0.         
- x  =
-    1.7975426  
-  - 0.3381487  
-    0.1633880  
-  - 4.9884023  
-    0.6054943  
-  - 3.1155623  
-//Only linear constraints (1 to 4) are active (lagr(1:6)=0):
-[x,lagr,f]=quapro(Q,p,C,b,[],[],me)   //Same result as above
- f  =
-  - 14.843248  
- lagr  =
-  - 1.5564027  
-  - 0.1698164  
-  - 0.7054782  
-    0.3091368  
-    0.         
- x  =
-    1.7975426  
-  - 0.3381487  
-    0.1633880  
-  - 4.9884023  
-    0.6054943  
-  - 3.1155623  
-// Example retrieved from:
-// http://www.linbox.com/ucome.rvt/any/doc_distrib/scilab-2.6/enpc/Berthaud_reponses
-r=1; k=400; m=1; g=10;
-Q=1/2*k*[1,-1,0,0;-1,2,-1,0;0,-1,2,-1;0,0,-1,1];
-P=m*g*r*[3;2;1;0];
-C=[0,0,0,1]; b=0; m=1;
-alphai=[];alphau=[];
-[alpha,lag,Umin]=quapro(Q,P,C,b,alphai,alphau,m)
- Umin  =
-  - 17.5  
- lag  =
-  - 60.  
- alpha  =
-  - 0.7   
-  - 0.55  
-  - 0.3   
-    0.    
-// Example retrieved from:
-// http://matman.uwm.edu.pl/~kulesza/Symulacje_komputerowe/SciLab/scilab12.pdf
-Q = [2,-1; -1,3]; c=[2;4]; C=[-2, -3; -5, -1]; d=[-5; -1]
- d  =
-  - 5.  
-  - 1.  
-XL = [0; 0]; XU = [1e10;1e10];
-[xopt, larg, fopt] = quapro(Q, c, C, d, XL, XU)
- fopt  =
-    7.3928571  
- larg  =
-    0.         
-    0.         
-    1.7857143  
-    0.         
- xopt  =
-    1.2142857  
-    0.8571429  
-// Example retrieved from:
-// http://www.poli.usp.br/d/pmr5215/Introscilab.pdf
-H = [ 2 0; 0 2];
-c = [-6; -8];
-A = [ 1 2; 3 1];
-b = [ 3 ; 4];
-[x, lagr, qx ] = quapro(H, c, A, b, zeros(2,1), [])
- qx  =
-  - 12.  
- lagr  =
-    0.   
-    0.   
-    2.8  
-    0.4  
- x  =
-    1.  
-    1.  
diff --git a/scilab/modules/optimization/tests/unit_tests/quapro.tst b/scilab/modules/optimization/tests/unit_tests/quapro.tst
deleted file mode 100644 (file)
index 8cbec14..0000000
+++ /dev/null
@@ -1,58 +0,0 @@
-// =============================================================================
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2007-2008 - INRIA - Sylvestre LEDRU <sylvestre.ledru@inria.fr>
-//
-//  This file is distributed under the same license as the Scilab package.
-// =============================================================================
-
-//Find x in R^6 such that:
-//C1*x = b1 (3 equality constraints i.e me=3)
-C1= [1,-1,1,0,3,1;
-    -1,0,-3,-4,5,6;
-     2,5,3,0,1,0];
-b1=[1;2;3];
-//C2*x <= b2 (2 inequality constraints)
-C2=[0,1,0,1,2,-1;
-    -1,0,2,1,1,0];
-b2=[-1;2.5];
-//with  x between ci and cs:
-ci=[-1000;-10000;0;-1000;-1000;-1000];cs=[10000;100;1.5;100;100;1000];
-//and minimize 0.5*x'*Q*x + p'*x with
-p=[1;2;3;4;5;6]; Q=eye(6,6);
-//No initial point is given;
-C=[C1;C2] ; //
-b=[b1;b2] ;  //
-me=3;
-[x,lagr,f]=quapro(Q,p,C,b,ci,cs,me)
-//Only linear constraints (1 to 4) are active (lagr(1:6)=0):
-[x,lagr,f]=quapro(Q,p,C,b,[],[],me)   //Same result as above
-
-
-// Example retrieved from: 
-// http://www.linbox.com/ucome.rvt/any/doc_distrib/scilab-2.6/enpc/Berthaud_reponses
-
-r=1; k=400; m=1; g=10;
-Q=1/2*k*[1,-1,0,0;-1,2,-1,0;0,-1,2,-1;0,0,-1,1];
-P=m*g*r*[3;2;1;0];
-
-C=[0,0,0,1]; b=0; m=1;
-alphai=[];alphau=[];
-
-[alpha,lag,Umin]=quapro(Q,P,C,b,alphai,alphau,m)
-
-// Example retrieved from:
-// http://matman.uwm.edu.pl/~kulesza/Symulacje_komputerowe/SciLab/scilab12.pdf
-
-Q = [2,-1; -1,3]; c=[2;4]; C=[-2, -3; -5, -1]; d=[-5; -1]
-
-XL = [0; 0]; XU = [1e10;1e10];
-
-[xopt, larg, fopt] = quapro(Q, c, C, d, XL, XU)
-
-// Example retrieved from:
-// http://www.poli.usp.br/d/pmr5215/Introscilab.pdf
-H = [ 2 0; 0 2];
-c = [-6; -8];
-A = [ 1 2; 3 1];
-b = [ 3 ; 4];
-[x, lagr, qx ] = quapro(H, c, A, b, zeros(2,1), [])