[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Newbie question - sqp solver equality constraints function
From: |
alexegin |
Subject: |
Re: Newbie question - sqp solver equality constraints function |
Date: |
Thu, 22 Jun 2017 00:04:44 -0700 (PDT) |
Hi, *Juan Pablo Carbajal-2*!
Thank you very much for answering me. Your test is working perfectly:
clear all;g =@( x ) [ sum( x( x < 0.0 ) ); sum( x( x >= 0.0 ) ) ] - [ -1.0;
1.0 ];c =@( x ) sumsq( x ); x0 = 2.0 * rand( 10, 1 ) - 1.0; tic;[ x obj info
iter nf l ] = sqp( x0, c, g, [], -1.0, 1.0 );toc;
But when I am trying to use constraint function in my script I get the
following error:
error: qp: equality constraint matrix must be full row rank
Here is my script:
clear all;global R = dlmread( "r.csv", "\t" ); # double [ 10x1 ]global Q =
dlmread( "q.csv", "\t" ); # double [ 10x10 ]global Y = dlmread( "y.csv",
"\t" ); # double [ 10x1 ]global lb = dlmread( "lb.csv", "\t" ); # double [
10x1 ]global ub = dlmread( "ub.csv", "\t" ); # double [ 10x1 ]g = @( x ) [
sum( x( x < 0.0 ) ); sum( x( x >= 0.0 ) ) ] - [ -1.0; 1.0 ];function z =
foo( m, n ) z = 0; for i = 1:rows( m ) z = z + m( i, 1 ) * n( i, 1 );
endfor; endfunction;function y = phi( x ) global Y; y = foo( Y, x ) * -1;
endfunctiontic;[ x obj info iter nf l ] = sqp( R, @phi, g, [], -1.0, 1.0
);toc;
What am I doing wrong? Please do not tell me /everything/
--
View this message in context:
http://octave.1599824.n4.nabble.com/Newbie-question-sqp-solver-equality-constraints-function-tp4683834p4683852.html
Sent from the Octave - General mailing list archive at Nabble.com.
- Newbie question - sqp solver equality constraints function, alexegin, 2017/06/21
- Re: Newbie question - sqp solver equality constraints function, Juan Pablo Carbajal, 2017/06/21
- Re: Newbie question - sqp solver equality constraints function,
alexegin <=
- Re: Newbie question - sqp solver equality constraints function, Juan Pablo Carbajal, 2017/06/22
- Re: Newbie question - sqp solver equality constraints function, alexegin, 2017/06/22
- Re: Newbie question - sqp solver equality constraints function, Juan Pablo Carbajal, 2017/06/22
- Re: Newbie question - sqp solver equality constraints function, alexegin, 2017/06/22
- Re: Newbie question - sqp solver equality constraints function, Juan Pablo Carbajal, 2017/06/24
- Re: Newbie question - sqp solver equality constraints function, alexegin, 2017/06/29
- Re: Newbie question - sqp solver equality constraints function, Juan Pablo Carbajal, 2017/06/29
- Re: Newbie question - sqp solver equality constraints function, alexegin, 2017/06/29
- Re: Newbie question - sqp solver equality constraints function, Juan Pablo Carbajal, 2017/06/29
- Re: Newbie question - sqp solver equality constraints function, alexegin, 2017/06/29