Re: nonconformant arguments in leasqr.m

[Olaf Till]

On Wed, May 04, 2011 at 01:57:59PM +0200, Olga Kruglova wrote:
> On Wed, May 4, 2011 at 1:55 PM, Andy Buckle <address@hidden> wrote:
> 
> > On Wed, May 4, 2011 at 12:50 PM, olga_kruglova <address@hidden> wrote:
> > > Hello everybody,
> > >
> > > I try to use leasqr.m to fit data with function convoluted with known
> > > resolution function but it gives me an error:
> > >>> Convolution
> > > Variables in the current scope:
> > >
> > >  Attr Name        Size                     Bytes  Class
> > >  ==== ====        ====                     =====  =====
> > >   f   y        1024x1                       8192  double
> > >
> > > Total is 1024 elements using 8192 bytes
> > >
> > > Variables in the current scope:
> > >
> > >  Attr Name        Size                     Bytes  Class
> > >  ==== ====        ====                     =====  =====
> > >       data     1024x1                       8192  double
> > >
> > > Total is 1024 elements using 8192 bytes
> > >
> > > Variables in the current scope:
> > >
> > >  Attr Name        Size                     Bytes  Class
> > >  ==== ====        ====                     =====  =====
> > >   f   wt          4x1                         32  double
> > >
> > > Total is 4 elements using 32 bytes
> > >
> > > Variables in the current scope:
> > >
> > >  Attr Name        Size                     Bytes  Class
> > >  ==== ====        ====                     =====  =====
> > >   f   y        1024x1                       8192  double
> > >
> > > Total is 1024 elements using 8192 bytes
> > >
> > > error: product: nonconformant arguments (op1 is 4x1, op2 is 1024x1)
> > > error: called from:
> > > error:   C:\OCTAVE~2\leasqr.m at line 184, column 3
> > > error:   C:\OCTAVE~2\Convolution.m at line 25, column 67
> > >  I can not understand why. Could you please help me out with it? Thank
> > you
> > > very much in advance!
> > >
> > > Kind regards, Olga
> >
> > Is Convolution.m your script?
> >
> > I think that in order to help, we have to see line 25 of Convolution.m.
> >
> > --
> > /* andy buckle */
> >
> 
> Yes, it is
>                  global verbose;
>                   verbose=1;
>  line 25:           [f1, p1, kvg1, iter1, corp1, covp1, covr1, stdresid1,
> Z1, r21] = ...
>                         leasqr (t, data, pin, F, stol, niter, options);

In order to reproduce the failure, we have to know the values of all
these variables: (t, data, pin, F, stol, niter, options), and probably
have to know the code of your function F.

But first: It seems you are using Matlab, and as far as I know the
current version of our 'leasqr' does not run in Matlab at all. Where
do you have your 'leasqr' from? If from Octave-Forge, from which
version of the optim package?

Olaf