## maxima-bugs

 [Maxima-bugs] [ maxima-Bugs-1573247 ] cholesky factor of double blocked matrix From: SourceForge.net - 2006-10-08 17:06:45 ```Bugs item #1573247, was opened at 2006-10-08 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 ```
 [Maxima-bugs] [ maxima-Bugs-1573247 ] cholesky factor of double blocked matrix From: SourceForge.net - 2006-10-10 11:15:16 ```Bugs item #1573247, was opened at 2006-10-08 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: 2006-10-10 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 ```