Hi,
Il giorno Dec 11, 2009, alle ore 11:10 PM, Thomas Weber <[email protected] > ha scritto: > Hi, > > the test suite of nrbnumbasisfun.m fails with Octave 3.2: > octave:1> test nrbnumbasisfun > ***** test > p = 2; > q = 3; > mcp = 2; ncp = 3; > knots = {[zeros(1,p), linspace(0,1,mcp-p+2), ones(1,p)], [zeros > (1,q), linspace(0,1,ncp-q+2), ones(1,q)]}; > Lx = 1; Ly = 1; > [cntl(1,:,:), cntl(2,:,:)] = meshgrid(linspace(0, Lx, ncp+1), > linspace(0, Ly, mcp+1) ); > cntl(4,:,:) = 1:numel(cntl(1,:,:)); > nrb = nrbmak(cntl, knots); > u = rand (1, 30); v = rand (1, 10); > N = nrbnumbasisfun ({u, v}, nrb); > assert (all(all(N>0)), true) > assert (all(all(N<=(ncp+1)*(mcp+1))), true) > assert (max(max(N)),(ncp+1)*(mcp+1)) > assert (min(min(N)),1) > !!!!! test failed > A(I,J,...) = X: dimensions mismatchoctave:2> > > The problem is the line > cntl(4,:,:) = 1:numel(cntl(1,:,:)); > The RHS translates to 12, which doesn't match with the definition of > cntl. I cannot check this myself right now but I believe the offending line in the test should be changed to: > cntl(4,:,:) = reshape(1:numel(cntl(1,:,:)),size( > cntl(4,:,:)); does the test succeed if you make this change? > Thomas Thanks for the report, c. ------------------------------------------------------------------------------ Return on Information: Google Enterprise Search pays you back Get the facts. http://p.sf.net/sfu/google-dev2dev _______________________________________________ Octave-dev mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/octave-dev
