From: <nwe...@ce...> - 2003-02-25 18:22:14
|
hello, I'm trying to use the non-liear multivariable optimization functions, but they doesn't work for more than 2 variables!... octave2.1:5> bfgs('sincn', [0.1 0.1]') ans = 3.1773 3.1773 octave2.1:6> bfgs('sincn', [0.1 0.1 0.1]') error: invalid vector index = 3 error: evaluating assignment expression near line 58, column 14 error: evaluating for command near line 56, column 5 error: evaluating if command near line 55, column 3 error: called from `bs_gradient' in file `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bs_gradient.m' error: called from `bfgs' in file `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bfgs.m' what am I doing wrong?!?! :/ ++Nicolau; -- Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D "Experience is not what happens to a man; it is what a man does with what happens to him." -- Aldous Huxley |
From: Etienne G. <et...@is...> - 2003-02-25 20:35:19
|
Hi, what is 'sincn'? Etienne On Tue, Feb 25, 2003 at 03:21:04PM -0300, nwe...@ce... wrote: # hello, I'm trying to use the non-liear multivariable optimization # functions, but they doesn't work for more than 2 variables!... # # octave2.1:5> bfgs('sincn', [0.1 0.1]') # ans = # # 3.1773 # 3.1773 # # octave2.1:6> bfgs('sincn', [0.1 0.1 0.1]') # error: invalid vector index = 3 # error: evaluating assignment expression near line 58, column 14 # error: evaluating for command near line 56, column 5 # error: evaluating if command near line 55, column 3 # error: called from `bs_gradient' in file # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bs_gradient.m' # error: called from `bfgs' in file # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bfgs.m' # # # what am I doing wrong?!?! :/ # # ++Nicolau; # # # -- # Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 # http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D # "Experience is not what happens to a man; it is what a man does with what happens to him." # -- Aldous Huxley # # # # ------------------------------------------------------- # This sf.net email is sponsored by:ThinkGeek # Welcome to geek heaven. # http://thinkgeek.com/sf # _______________________________________________ # Octave-dev mailing list # Oct...@li... # https://lists.sourceforge.net/lists/listinfo/octave-dev # # -- Etienne Grossmann ------ http://www.isr.ist.utl.pt/~etienne |
From: <nwe...@ce...> - 2003-02-25 22:11:18
|
it's a function that receives a list and give a real number. bfgs('prod', [0.1 0.1 0.1]') gives the same error. On Tue, Feb 25, 2003 at 08:40:46PM +0000, Etienne Grossmann wrote: > > Hi, > > what is 'sincn'? > > Etienne > > On Tue, Feb 25, 2003 at 03:21:04PM -0300, nwe...@ce... wrote: > # hello, I'm trying to use the non-liear multivariable optimization > # functions, but they doesn't work for more than 2 variables!... > # > # octave2.1:5> bfgs('sincn', [0.1 0.1]') > # ans = > # > # 3.1773 > # 3.1773 > # > # octave2.1:6> bfgs('sincn', [0.1 0.1 0.1]') > # error: invalid vector index = 3 > # error: evaluating assignment expression near line 58, column 14 > # error: evaluating for command near line 56, column 5 > # error: evaluating if command near line 55, column 3 > # error: called from `bs_gradient' in file > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bs_gradient.m' > # error: called from `bfgs' in file > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bfgs.m' > # > # > # what am I doing wrong?!?! :/ > # > # ++Nicolau; > # > # > # -- > # Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 > # http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D > # "Experience is not what happens to a man; it is what a man does with what happens to him." > # -- Aldous Huxley > # > # > # > # ------------------------------------------------------- > # This sf.net email is sponsored by:ThinkGeek > # Welcome to geek heaven. > # http://thinkgeek.com/sf > # _______________________________________________ > # Octave-dev mailing list > # Oct...@li... > # https://lists.sourceforge.net/lists/listinfo/octave-dev > # > # > > -- > Etienne Grossmann ------ http://www.isr.ist.utl.pt/~etienne > > Return-Path: <oct...@li...> > Delivered-To: nwe...@ce... > Received: from pinguim.cpdee.ufmg.br [150.164.32.32] > by localhost with POP3 (fetchmail-5.9.11) > for nwerneck@localhost (single-drop); Tue, 25 Feb 2003 18:17:30 -0300 (BRT) > Received: (qmail 16127 invoked from network); 25 Feb 2003 20:36:44 -0000 > Received: from lists.sourceforge.net (HELO sc8-sf-list2.sourceforge.net) (66.35.250.206) > by pinguim.cpdee.ufmg.br with SMTP; 25 Feb 2003 20:36:43 -0000 > Received: from sc8-sf-list1-b.sourceforge.net ([10.3.1.13] helo=sc8-sf-list1.sourceforge.net) > by sc8-sf-list2.sourceforge.net with esmtp (Exim 3.31-VA-mm2 #1 (Debian)) > id 18nloC-00038O-00; Tue, 25 Feb 2003 12:36:12 -0800 > Received: from gtisr.ist.utl.pt ([193.136.138.253] helo=anonimo.isr.ist.utl.pt) > by sc8-sf-list1.sourceforge.net with esmtp (Exim 3.31-VA-mm2 #1 (Debian)) > id 18nlmg-0005Mz-00 > for <oct...@li...>; Tue, 25 Feb 2003 12:34:39 -0800 > Received: from etienne by anonimo.isr.ist.utl.pt with local (Exim 3.35 #1 (Debian)) > id 18nlsd-0007dS-00; Tue, 25 Feb 2003 20:40:47 +0000 > To: nwe...@ce... > Cc: oct...@li..., Etienne Grossmann <et...@is...> > Subject: Re: [OctDev] how to use optimizing functions > Message-ID: <200...@an...> > Reply-To: et...@is... > References: <200...@ce...> > Mime-Version: 1.0 > Content-Type: text/plain; charset=us-ascii > Content-Disposition: inline > In-Reply-To: <200...@ce...> > User-Agent: Mutt/1.3.28i > From: Etienne Grossmann <et...@is...> > Sender: oct...@li... > Errors-To: oct...@li... > X-BeenThere: oct...@li... > X-Mailman-Version: 2.0.9-sf.net > Precedence: bulk > List-Help: <mailto:oct...@li...?subject=help> > List-Post: <mailto:oct...@li...> > List-Subscribe: <https://lists.sourceforge.net/lists/listinfo/octave-dev>, > <mailto:oct...@li...?subject=subscribe> > List-Id: Developers of the GNU Octave Repository. <octave-dev.lists.sourceforge.net> > List-Unsubscribe: <https://lists.sourceforge.net/lists/listinfo/octave-dev>, > <mailto:oct...@li...?subject=unsubscribe> > List-Archive: <http://sourceforge.net/mailarchive/forum.php?forum=octave-dev> > X-Original-Date: Tue, 25 Feb 2003 20:40:46 +0000 > Date: Tue, 25 Feb 2003 20:40:46 +0000 > > > Hi, > > what is 'sincn'? > > Etienne > > On Tue, Feb 25, 2003 at 03:21:04PM -0300, nwe...@ce... wrote: > # hello, I'm trying to use the non-liear multivariable optimization > # functions, but they doesn't work for more than 2 variables!... > # > # octave2.1:5> bfgs('sincn', [0.1 0.1]') > # ans = > # > # 3.1773 > # 3.1773 > # > # octave2.1:6> bfgs('sincn', [0.1 0.1 0.1]') > # error: invalid vector index = 3 > # error: evaluating assignment expression near line 58, column 14 > # error: evaluating for command near line 56, column 5 > # error: evaluating if command near line 55, column 3 > # error: called from `bs_gradient' in file > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bs_gradient.m' > # error: called from `bfgs' in file > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bfgs.m' > # > # > # what am I doing wrong?!?! :/ > # > # ++Nicolau; > # > # > # -- > # Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 > # http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D > # "Experience is not what happens to a man; it is what a man does with what happens to him." > # -- Aldous Huxley > # > # > # > # ------------------------------------------------------- > # This sf.net email is sponsored by:ThinkGeek > # Welcome to geek heaven. > # http://thinkgeek.com/sf > # _______________________________________________ > # Octave-dev mailing list > # Oct...@li... > # https://lists.sourceforge.net/lists/listinfo/octave-dev > # > # > > -- > Etienne Grossmann ------ http://www.isr.ist.utl.pt/~etienne > > > ------------------------------------------------------- > This sf.net email is sponsored by:ThinkGeek > Welcome to geek heaven. > http://thinkgeek.com/sf > _______________________________________________ > Octave-dev mailing list > Oct...@li... > https://lists.sourceforge.net/lists/listinfo/octave-dev > -- Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D "The great end of life is not knowledge but action." -- Thomas Huxley |
From: Etienne G. <et...@is...> - 2003-02-26 08:29:54
|
On Tue, Feb 25, 2003 at 07:10:19PM -0300, nwe...@ce... wrote: # it's a function that receives a list and give a real number. ^ vector # bfgs('prod', [0.1 0.1 0.1]') gives the same error. :-) . You expect a minimum at a finite point? Cheers, Etienne # # On Tue, Feb 25, 2003 at 08:40:46PM +0000, Etienne Grossmann wrote: # > # > Hi, # > # > what is 'sincn'? # > # > Etienne # > # > On Tue, Feb 25, 2003 at 03:21:04PM -0300, nwe...@ce... wrote: # > # hello, I'm trying to use the non-liear multivariable optimization # > # functions, but they doesn't work for more than 2 variables!... # > # # > # octave2.1:5> bfgs('sincn', [0.1 0.1]') # > # ans = # > # # > # 3.1773 # > # 3.1773 # > # # > # octave2.1:6> bfgs('sincn', [0.1 0.1 0.1]') # > # error: invalid vector index = 3 # > # error: evaluating assignment expression near line 58, column 14 # > # error: evaluating for command near line 56, column 5 # > # error: evaluating if command near line 55, column 3 # > # error: called from `bs_gradient' in file # > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bs_gradient.m' # > # error: called from `bfgs' in file # > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bfgs.m' # > # # > # # > # what am I doing wrong?!?! :/ # > # # > # ++Nicolau; # > # # > # # > # -- # > # Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 # > # http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D # > # "Experience is not what happens to a man; it is what a man does with what happens to him." # > # -- Aldous Huxley # > # # > # # > # # > # ------------------------------------------------------- # > # This sf.net email is sponsored by:ThinkGeek # > # Welcome to geek heaven. # > # http://thinkgeek.com/sf # > # _______________________________________________ # > # Octave-dev mailing list # > # Oct...@li... # > # https://lists.sourceforge.net/lists/listinfo/octave-dev # > # # > # # > # > -- # > Etienne Grossmann ------ http://www.isr.ist.utl.pt/~etienne # > # > Return-Path: <oct...@li...> # > Delivered-To: nwe...@ce... # > Received: from pinguim.cpdee.ufmg.br [150.164.32.32] # > by localhost with POP3 (fetchmail-5.9.11) # > for nwerneck@localhost (single-drop); Tue, 25 Feb 2003 18:17:30 -0300 (BRT) # > Received: (qmail 16127 invoked from network); 25 Feb 2003 20:36:44 -0000 # > Received: from lists.sourceforge.net (HELO sc8-sf-list2.sourceforge.net) (66.35.250.206) # > by pinguim.cpdee.ufmg.br with SMTP; 25 Feb 2003 20:36:43 -0000 # > Received: from sc8-sf-list1-b.sourceforge.net ([10.3.1.13] helo=sc8-sf-list1.sourceforge.net) # > by sc8-sf-list2.sourceforge.net with esmtp (Exim 3.31-VA-mm2 #1 (Debian)) # > id 18nloC-00038O-00; Tue, 25 Feb 2003 12:36:12 -0800 # > Received: from gtisr.ist.utl.pt ([193.136.138.253] helo=anonimo.isr.ist.utl.pt) # > by sc8-sf-list1.sourceforge.net with esmtp (Exim 3.31-VA-mm2 #1 (Debian)) # > id 18nlmg-0005Mz-00 # > for <oct...@li...>; Tue, 25 Feb 2003 12:34:39 -0800 # > Received: from etienne by anonimo.isr.ist.utl.pt with local (Exim 3.35 #1 (Debian)) # > id 18nlsd-0007dS-00; Tue, 25 Feb 2003 20:40:47 +0000 # > To: nwe...@ce... # > Cc: oct...@li..., Etienne Grossmann <et...@is...> # > Subject: Re: [OctDev] how to use optimizing functions # > Message-ID: <200...@an...> # > Reply-To: et...@is... # > References: <200...@ce...> # > Mime-Version: 1.0 # > Content-Type: text/plain; charset=us-ascii # > Content-Disposition: inline # > In-Reply-To: <200...@ce...> # > User-Agent: Mutt/1.3.28i # > From: Etienne Grossmann <et...@is...> # > Sender: oct...@li... # > Errors-To: oct...@li... # > X-BeenThere: oct...@li... # > X-Mailman-Version: 2.0.9-sf.net # > Precedence: bulk # > List-Help: <mailto:oct...@li...?subject=help> # > List-Post: <mailto:oct...@li...> # > List-Subscribe: <https://lists.sourceforge.net/lists/listinfo/octave-dev>, # > <mailto:oct...@li...?subject=subscribe> # > List-Id: Developers of the GNU Octave Repository. <octave-dev.lists.sourceforge.net> # > List-Unsubscribe: <https://lists.sourceforge.net/lists/listinfo/octave-dev>, # > <mailto:oct...@li...?subject=unsubscribe> # > List-Archive: <http://sourceforge.net/mailarchive/forum.php?forum=octave-dev> # > X-Original-Date: Tue, 25 Feb 2003 20:40:46 +0000 # > Date: Tue, 25 Feb 2003 20:40:46 +0000 # > # > # > Hi, # > # > what is 'sincn'? # > # > Etienne # > # > On Tue, Feb 25, 2003 at 03:21:04PM -0300, nwe...@ce... wrote: # > # hello, I'm trying to use the non-liear multivariable optimization # > # functions, but they doesn't work for more than 2 variables!... # > # # > # octave2.1:5> bfgs('sincn', [0.1 0.1]') # > # ans = # > # # > # 3.1773 # > # 3.1773 # > # # > # octave2.1:6> bfgs('sincn', [0.1 0.1 0.1]') # > # error: invalid vector index = 3 # > # error: evaluating assignment expression near line 58, column 14 # > # error: evaluating for command near line 56, column 5 # > # error: evaluating if command near line 55, column 3 # > # error: called from `bs_gradient' in file # > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bs_gradient.m' # > # error: called from `bfgs' in file # > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bfgs.m' # > # # > # # > # what am I doing wrong?!?! :/ # > # # > # ++Nicolau; # > # # > # # > # -- # > # Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 # > # http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D # > # "Experience is not what happens to a man; it is what a man does with what happens to him." # > # -- Aldous Huxley # > # # > # # > # # > # ------------------------------------------------------- # > # This sf.net email is sponsored by:ThinkGeek # > # Welcome to geek heaven. # > # http://thinkgeek.com/sf # > # _______________________________________________ # > # Octave-dev mailing list # > # Oct...@li... # > # https://lists.sourceforge.net/lists/listinfo/octave-dev # > # # > # # > # > -- # > Etienne Grossmann ------ http://www.isr.ist.utl.pt/~etienne # > # > # > ------------------------------------------------------- # > This sf.net email is sponsored by:ThinkGeek # > Welcome to geek heaven. # > http://thinkgeek.com/sf # > _______________________________________________ # > Octave-dev mailing list # > Oct...@li... # > https://lists.sourceforge.net/lists/listinfo/octave-dev # > # # -- # Nicolau Werneck <nwe...@ce...> E073 FADE 7B75 5482 5631 # http://cefala.org/~nwerneck D3C4 6E9A 8EFC 0351 C15D # "The great end of life is not knowledge but action." # -- Thomas Huxley # # # # ------------------------------------------------------- # This sf.net email is sponsored by:ThinkGeek # Welcome to geek heaven. # http://thinkgeek.com/sf # _______________________________________________ # Octave-dev mailing list # Oct...@li... # https://lists.sourceforge.net/lists/listinfo/octave-dev # # -- Etienne Grossmann ------ http://www.isr.ist.utl.pt/~etienne |
From: Etienne G. <et...@is...> - 2003-02-26 12:44:15
|
Hello, use octave 2.1.XX. I am not sure whether bfgs() is meant to work with octave 2.0.XX. With 2.1.44, sincn (bfgs('sincn', [0.1 0.1 0.1]')) yields -0.21723, which seems to be the minimum of 'sincn'. Note that "sincn" has non isolated minima and that bfgs() looks for a local minimum. On Wed, Feb 26, 2003 at 08:02:13AM -0300, nwe...@ce... wrote: # well, sorry, but 'prod' converges with X0 = [0.1 0.1], and I didn't test # other points... This is a glitch in the bfgs() algorithm. If started from any_value*ones(1,N), it will go straight down to something near zero. Once there, it sees the derivative is zero and calls it a stop. This has to do with the way line minimization is done. Should fix. Hth, Etienne # sincn is defined this way: # # function r = sincn(X) # n = norm(X); # r = sin(n+0.00001)/n; # return # # I can always minimize it to something (local minima) with two points, # but that weird error occurs with more. # # I'm moving to the bleeding-edgre development version since it seems that # a lot has changed recentely in the optimization functions... but the help # from bfgs and dfp in this version does say that it works for more variables... # # # # On Wed, Feb 26, 2003 at 08:35:47AM +0000, Etienne Grossmann wrote: # > On Tue, Feb 25, 2003 at 07:10:19PM -0300, nwe...@ce... wrote: # > # it's a function that receives a list and give a real number. # > ^ # > vector # > # > # bfgs('prod', [0.1 0.1 0.1]') gives the same error. # > # > :-) . You expect a minimum at a finite point? # > # > Cheers, # > # > Etienne # > # > # # > # On Tue, Feb 25, 2003 at 08:40:46PM +0000, Etienne Grossmann wrote: # > # > # > # > Hi, # > # > # > # > what is 'sincn'? # > # > # > # > Etienne # > # > # > # > On Tue, Feb 25, 2003 at 03:21:04PM -0300, nwe...@ce... wrote: # > # > # hello, I'm trying to use the non-liear multivariable optimization # > # > # functions, but they doesn't work for more than 2 variables!... # > # > # # > # > # octave2.1:5> bfgs('sincn', [0.1 0.1]') # > # > # ans = # > # > # # > # > # 3.1773 # > # > # 3.1773 # > # > # # > # > # octave2.1:6> bfgs('sincn', [0.1 0.1 0.1]') # > # > # error: invalid vector index = 3 # > # > # error: evaluating assignment expression near line 58, column 14 # > # > # error: evaluating for command near line 56, column 5 # > # > # error: evaluating if command near line 55, column 3 # > # > # error: called from `bs_gradient' in file # > # > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bs_gradient.m' # > # > # error: called from `bfgs' in file # > # > # `/usr/share/octave/2.1.35/site/m/octave-forge/optim/bfgs.m' # > # > # # > # > # # > # > # what am I doing wrong?!?! :/ # > # > # # > # > # ++Nicolau; -- Etienne Grossmann ------ http://www.isr.ist.utl.pt/~etienne |