mesh and meshdom bugs?

[John Eaton]

Joao Cardoso <address@hidden> wrote:

: If I am correct, there is a bug in `meshdom' and `mesh'.

Yes.

: In `meshdom' for some strange reason the y vector is calculated backward:
: 
:      y = y (ylen:-1:1);         # jc: this line seems to be incorrect
:      y = y (1:ylen);          # jc: it should(?) be this way
: 
: in `mesh', the following comparision also seems to be wrong:
: 
:       if (xlen == rows (z) && ylen == columns (z))    # wrong ?
:       if (ylen == rows (z) && xlen == columns (z))    # right ?

These changes look ok, except that in the first case, y = y(1:ylen)
doesn't do anything, so that line can be eliminated.

This makes plotting using mesh and meshdom work for nonsquare meshes,
but it also results in meshdom returning a Y matrix that is flipped
compared to what the Matlab meshdom returns.  I don't know that it
matters that much, but someone will no doubt complain that it is
incompatible...

Thanks,

jwe