Mercurial > hg > octave-lyh

comparison src/DLD-FUNCTIONS/luinc.cc @ 5610:9761b7d24e9e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2006-02-09 09:12:02 by dbateman]
author dbateman
date Thu, 09 Feb 2006 09:12:03 +0000
parents 22994a5730f9
children 2618a0750ae6
comparison
equal deleted inserted replaced
5609:bf96b0f9dbd7 5610:9761b7d24e9e
277 277
278 return retval; 278 return retval;
279 } 279 }
280 280
281 /* 281 /*
282
283 %!test
284 %! a=sparse([1,2,0,0;0,1,2,0;1e-14,0,3,0;0,0,0,1]);
285 %! [l,u]=luinc(a,1e-10);
286 %! assert(l*u, sparse([1,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10);
287 %! opts.droptol=1e-10;
288 %! [l,u]=luinc(a,opts);
289 %! assert(l*u, sparse([1,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10);
290
291 %!test
292 %! a=sparse([1i,2,0,0;0,1,2,0;1e-14,0,3,0;0,0,0,1]);
293 %! [l,u]=luinc(a,1e-10);
294 %! assert(l*u, sparse([1i,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10);
295 %! opts.droptol=1e-10;
296 %! [l,u]=luinc(a,opts);
297 %! assert(l*u, sparse([1i,2,0,0;0,1,2,0;0,0,3,0;0,0,0,1]),1e-10);
298
299 %!error splu(sparse([1i,2,0,0;0,1,2,0;1e-14,0,3,0;0,0,0,1]));
300
301 */
302
303 /*
282 ;;; Local Variables: *** 304 ;;; Local Variables: ***
283 ;;; mode: C++ *** 305 ;;; mode: C++ ***
284 ;;; End: *** 306 ;;; End: ***
285 */ 307 */