diff options
Diffstat (limited to 'tests/unit_tests/intfmincon.tst')
-rw-r--r-- | tests/unit_tests/intfmincon.tst | 76 |
1 files changed, 76 insertions, 0 deletions
diff --git a/tests/unit_tests/intfmincon.tst b/tests/unit_tests/intfmincon.tst new file mode 100644 index 0000000..144ff2b --- /dev/null +++ b/tests/unit_tests/intfmincon.tst @@ -0,0 +1,76 @@ +// Copyright (C) 2017 - IIT Bombay - FOSSEE +// +// Author: Georgey John +// Organization: FOSSEE, IIT Bombay +// Email: toolbox@scilab.in +// +// This file must be used under the terms of the CeCILL. +// This source file is licensed as described in the file COPYING, which +// you should have received as part of this distribution. The terms +// are also available at +// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt + +// <-- JVM NOT MANDATORY --> +// <-- ENGLISH IMPOSED --> + + +// +// assert_close -- +// Returns 1 if the two real matrices computed and expected are close, +// i.e. if the relative distance between computed and expected is lesser than epsilon. +// Arguments +// computed, expected : the two matrices to compare +// epsilon : a small number +// +function flag = assert_close ( computed, expected, epsilon ) + if expected==0.0 then + shift = norm(computed-expected); + else + shift = norm(computed-expected)/norm(expected); + end +// if shift < epsilon then +// flag = 1; +// else +// flag = 0; +// end +// if flag <> 1 then pause,end + flag = assert_checktrue ( shift < epsilon ); +endfunction +// +// assert_equal -- +// Returns 1 if the two real matrices computed and expected are equal. +// Arguments +// computed, expected : the two matrices to compare +// epsilon : a small number +// +//function flag = assert_equal ( computed , expected ) +// if computed==expected then +// flag = 1; +// else +// flag = 0; +// end +// if flag <> 1 then pause,end +//endfunction + + +//Constrained optimization +function [y]=f(x) +y=-x(1)^2-x(2)^2-x(2)/3; +endfunction + +x0=[0 , 0]; +intcon = [1] +A=[1,1 ; 1,1/4 ; 1,-1 ; -1/4,-1 ; -1,-1 ; -1,1]; +b=[2;1;2;1;-1;2]; +Aeq=[1,1]; +beq=[2]; +lb=[]; +ub=[]; +nlc=[]; + +[xopt,fopt,exitflag,grad,hessian] =intfmincon(f, x0,intcon,A,b,Aeq,beq,lb,ub,nlc) + +assert_close ( xopt , [ 0 2 ]' , 0.0005 ); +assert_close ( fopt , [ -4.6666 ]' , 0.0005 ); +assert_checkequal( double(exitflag) , 0 ); +printf("Test Successful");
\ No newline at end of file |