http://www.mapleprimes.com/questions/142810-Cholesky-Decomposition-For-A-Symbolic-Matrix en-us 2026 Maplesoft, A Division of Waterloo Maple Inc. Maplesoft Document System Thu, 11 Jun 2026 00:34:43 GMT Thu, 11 Jun 2026 00:34:43 GMT The latest answers and comments added to the Question, Cholesky Decomposition for a symbolic matrix http://www.mapleprimes.com/images/mapleprimeswhite.jpg http://www.mapleprimes.com/questions/142810-Cholesky-Decomposition-For-A-Symbolic-Matrix http://www.mapleprimes.com/questions/142810-Cholesky-Decomposition-For-A-Symbolic-Matrix?ref=Feed:MaplePrimes:Cholesky Decomposition for a symbolic matrix:Comments#answer142811 <p>Just add the conjugate=false option to your LUDecomposition call. This is mentioned in its help docs. You can even omit the simplify call afterwards. But you will need to apply it to C.Transpose(C) to obtain the original Matrix.</p> <p>Just add the conjugate=false option to your LUDecomposition call. This is mentioned in its help docs. You can even omit the simplify call afterwards. But you will need to apply it to C.Transpose(C) to obtain the original Matrix.</p> 142811 Wed, 30 Jan 2013 18:33:05 Z Thomas Richard Thomas Richard http://www.mapleprimes.com/questions/142810-Cholesky-Decomposition-For-A-Symbolic-Matrix?ref=Feed:MaplePrimes:Cholesky Decomposition for a symbolic matrix:Comments#answer142812 <pre>A := Matrix([[a, b], [b, c]], shape = symmetric, attributes = [positive_definite]): B := LinearAlgebra:-LUDecomposition(A, method = Cholesky); Bhat := simplify(B) assuming real, positive; map(normal, Bhat.Bhat^%T - A ); F := LinearAlgebra:-LUDecomposition(A, method = Cholesky, conjugate=false); map(normal, F.F^%T - A ); </pre> <p>If you go with conjugate=false then you might still want to simplify under the assumption of being positive, depending on example.</p> <!--break--> <p>acer</p> <pre>A := Matrix([[a, b], [b, c]], shape = symmetric, attributes = [positive_definite]): B := LinearAlgebra:-LUDecomposition(A, method = Cholesky); Bhat := simplify(B) assuming real, positive; map(normal, Bhat.Bhat^%T - A ); F := LinearAlgebra:-LUDecomposition(A, method = Cholesky, conjugate=false); map(normal, F.F^%T - A ); </pre> <p>If you go with conjugate=false then you might still want to simplify under the assumption of being positive, depending on example.</p> <!--break--> <p>acer</p> 142812 Wed, 30 Jan 2013 18:35:13 Z acer acer