From: SourceForge.net <noreply@so...>  20061008 17:06:45

Bugs item #1573247, was opened at 20061008 12:06 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1573247&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Share Libraries Group: None Status: Open Resolution: None Priority: 4 Submitted By: Barton Willis (willisbl) Assigned to: Barton Willis (willisbl) Summary: cholesky factor of double blocked matrix Initial Comment: (%o90) matrix([matrix([matrix([8,4],[4,13]), matrix([2,4],[7,9])],[matrix([2,7],[4,9]), matrix([10,6],[6,15])]), matrix([matrix([4,4],[12,5]), matrix([0,6],[6,8])],[matrix([8,3],[13,3]), matrix([6,5],[5,6])])],[matrix([matrix([4,12],[4,5]), matrix([8,13],[3,3])],[matrix([0,6],[6,8]), matrix([6,5],[5,6])]),matrix([matrix([18,5],[5,4]), matrix([5,9],[2,6])],[matrix([5,2],[9,6]), matrix([6,3],[3,11])])]) (%i91) cholesky(%,'noncommutingring); Division by 0 This shouldn't give a divide by zero error. I think the bug doesn't happen for single blocked matrices. The problem is that cholesky tries to invert a blocked matrix with "^^". That doesn't work for block matrices.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1573247&group_id=4933 
From: SourceForge.net <noreply@so...>  20061010 11:15:16

Bugs item #1573247, was opened at 20061008 12:06 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1573247&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Share Libraries Group: None >Status: Closed >Resolution: Fixed Priority: 4 Submitted By: Barton Willis (willisbl) Assigned to: Barton Willis (willisbl) Summary: cholesky factor of double blocked matrix Initial Comment: (%o90) matrix([matrix([matrix([8,4],[4,13]), matrix([2,4],[7,9])],[matrix([2,7],[4,9]), matrix([10,6],[6,15])]), matrix([matrix([4,4],[12,5]), matrix([0,6],[6,8])],[matrix([8,3],[13,3]), matrix([6,5],[5,6])])],[matrix([matrix([4,12],[4,5]), matrix([8,13],[3,3])],[matrix([0,6],[6,8]), matrix([6,5],[5,6])]),matrix([matrix([18,5],[5,4]), matrix([5,9],[2,6])],[matrix([5,2],[9,6]), matrix([6,3],[3,11])])]) (%i91) cholesky(%,'noncommutingring); Division by 0 This shouldn't give a divide by zero error. I think the bug doesn't happen for single blocked matrices. The problem is that cholesky tries to invert a blocked matrix with "^^". That doesn't work for block matrices.  >Comment By: Barton Willis (willisbl) Date: 20061010 06:15 Message: Logged In: YES user_id=895922 Fixed by mring.lisp reversion 1.13.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1573247&group_id=4933 