Thread Subject:

Did ndgrid change?

I had this bit of code which used to work in a previous
version of MATLAB but no longer functions in the current
version

for k=1:n
    [q{k},w{k}]=rquad(N,k-1);
end

[Q{1:n}]=ndgrid(q{:}); [W{1:n}]=ndgrid(w{:});

clear q w; q=zeros(Nn,n); w=q; X=q;

for k=1:n
    q(:,k)=reshape(Q{k},Nn,1); w(:,k)=reshape(W{k},Nn,1);
end

map=eye(m); map(2:m,1)=-1; c=map*vert;
W=abs(det(c(2:m,:)))*prod(w,2);

for k=1:n-1
    X(:,k)=(1-q(:,k)).*prod(q(:,1:k-1));
end

which now gives the error

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> nsimplexquad at 37
    X(:,k)=(1-q(:,k)).*prod(q(:,1:k-1));


It seems that now ndgrid only returns the first argument. Is
there a workaround?

Thanks,
Greg

"Greg von Winckel" <gregvw@gmail.com> wrote in message
<fjocab$nqi$1@fred.mathworks.com>...
> I had this bit of code which used to work in a previous
> version of MATLAB but no longer functions in the current
> version
>
> for k=1:n
> [q{k},w{k}]=rquad(N,k-1);
> end
>
> [Q{1:n}]=ndgrid(q{:}); [W{1:n}]=ndgrid(w{:});
>
> clear q w; q=zeros(Nn,n); w=q; X=q;
>
> for k=1:n
> q(:,k)=reshape(Q{k},Nn,1); w(:,k)=reshape(W{k},Nn,1);
> end
>
> map=eye(m); map(2:m,1)=-1; c=map*vert;
> W=abs(det(c(2:m,:)))*prod(w,2);
>
> for k=1:n-1
> X(:,k)=(1-q(:,k)).*prod(q(:,1:k-1));
> end
>
> which now gives the error
>
> ??? Error using ==> times
> Matrix dimensions must agree.
>
> Error in ==> nsimplexquad at 37
> X(:,k)=(1-q(:,k)).*prod(q(:,1:k-1));
>
>
> It seems that now ndgrid only returns the first argument. Is
> there a workaround?
>
> Thanks,
> Greg
>

I am not sure, but could it be that

[C{1:n}] = ndgrid(q{:})

is interpreted (wrongly) as

[C{1:n}] = deal(ndgrid(q{:}))

??

Jos

"Greg von Winckel" <gregvw@gmail.com> wrote in message
<fjocab$nqi$1@fred.mathworks.com>...
> I had this bit of code which used to work in a previous
> version of MATLAB but no longer functions in the current
> version
>
> for k=1:n
> [q{k},w{k}]=rquad(N,k-1);
> end
>
> [Q{1:n}]=ndgrid(q{:}); [W{1:n}]=ndgrid(w{:});

(snip)

> It seems that now ndgrid only returns the first argument. Is
> there a workaround?

This works for me in R2007a.

w = {1:2, 1:3, 1:4};
n = 3;

[W{1:n}]=ndgrid(w{:})
W =
    [2x3x4 double] [2x3x4 double] [2x3x4 double]

John

"Greg von Winckel" <gregvw@gmail.com> wrote in message
<fjocab$nqi$1@fred.mathworks.com>...
> I had this bit of code which used to work in a previous
> version of MATLAB but no longer functions in the current
> version
>
> for k=1:n
> [q{k},w{k}]=rquad(N,k-1);
> end
>
> [Q{1:n}]=ndgrid(q{:}); [W{1:n}]=ndgrid(w{:});

(snip)

> It seems that now ndgrid only returns the first argument. Is
> there a workaround?

This works for me in R2007a.

w = {1:2, 1:3, 1:4};
n = 3;

[W{1:n}]=ndgrid(w{:})
W =
    [2x3x4 double] [2x3x4 double] [2x3x4 double]

John

"Greg von Winckel" <gregvw@gmail.com> wrote in message
<fjocab$nqi$1@fred.mathworks.com>...
> I had this bit of code which used to work in a previous
> version of MATLAB but no longer functions in the current
> version
>
> for k=1:n
> [q{k},w{k}]=rquad(N,k-1);
> end
>
> [Q{1:n}]=ndgrid(q{:}); [W{1:n}]=ndgrid(w{:});

(snip)

> It seems that now ndgrid only returns the first argument. Is
> there a workaround?

This works for me in R2007a.

w = {1:2, 1:3, 1:4};
n = 3;

[W{1:n}]=ndgrid(w{:})
W =
    [2x3x4 double] [2x3x4 double] [2x3x4 double]

John

Thanks, John. I don't know why I was getting this odd
result, but it seems to be behaving now.

"John D'Errico" <woodchips@rochester.rr.com> wrote in
message <fjoip2$cfi$1@fred.mathworks.com>...
> "Greg von Winckel" <gregvw@gmail.com> wrote in message
> <fjocab$nqi$1@fred.mathworks.com>...
> > I had this bit of code which used to work in a previous
> > version of MATLAB but no longer functions in the current
> > version
> >
> > for k=1:n
> > [q{k},w{k}]=rquad(N,k-1);
> > end
> >
> > [Q{1:n}]=ndgrid(q{:}); [W{1:n}]=ndgrid(w{:});
>
> (snip)
>
> > It seems that now ndgrid only returns the first argument. Is
> > there a workaround?
>
> This works for me in R2007a.
>
> w = {1:2, 1:3, 1:4};
> n = 3;
>
> [W{1:n}]=ndgrid(w{:})
> W =
> [2x3x4 double] [2x3x4 double] [2x3x4 double]
>
> John