Mercurial > hg > octave-lyh

comparison src/DLD-FUNCTIONS/luinc.cc @ 7243:e8d7eed42935

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-12-04 02:10:17 by jwe]
author jwe
date Tue, 04 Dec 2007 02:10:18 +0000
parents a1dbe9d80eee
children f3c00dc0912b
comparison
equal deleted inserted replaced
7242:e4398e3903be 7243:e8d7eed42935
298 return retval; 298 return retval;
299 } 299 }
300 300
301 /* 301 /*
302 302
303 %!test 303 %!testif HAVE_UMFPACK
304 %! a=sparse([1,2,0,0;0,1,2,0;1e-14,0,3,0;0,0,0,1]); 304 %! a=sparse([1,2,0,0;0,1,2,0;1e-14,0,3,0;0,0,0,1]);
305 %! [l,u]=luinc(a,1e-10); 305 %! [l,u]=luinc(a,1e-10);
306 %! assert(l*u, sparse([1,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10); 306 %! assert(l*u, sparse([1,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10);
307 %! opts.droptol=1e-10; 307 %! opts.droptol=1e-10;
308 %! [l,u]=luinc(a,opts); 308 %! [l,u]=luinc(a,opts);
309 %! assert(l*u, sparse([1,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10); 309 %! assert(l*u, sparse([1,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10);
310 310
311 %!test 311 %!testif HAVE_UMFPACK
312 %! a=sparse([1i,2,0,0;0,1,2,0;1e-14,0,3,0;0,0,0,1]); 312 %! a=sparse([1i,2,0,0;0,1,2,0;1e-14,0,3,0;0,0,0,1]);
313 %! [l,u]=luinc(a,1e-10); 313 %! [l,u]=luinc(a,1e-10);
314 %! assert(l*u, sparse([1i,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10); 314 %! assert(l*u, sparse([1i,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10);
315 %! opts.droptol=1e-10; 315 %! opts.droptol=1e-10;
316 %! [l,u]=luinc(a,opts); 316 %! [l,u]=luinc(a,opts);