Is there a built-in/easy LDU decomposition method in Numpy?

user_1559897 0
Viewed 8999

I see cholesky decomposition in numpy.linalg.cholesky, but could not find a LDU decompositon. Can anyone suggest a function to use?

linear-algebra numpy python
2 Answers

Try this:

import numpy as np

A = np.array([[1,2,3,4],[1,2,3,4],[1,2,3,4],[1,2,3,4]])
U = np.triu(A,1)
L = np.tril(A,-1)
D = np.tril(np.triu(A))
print(A)
print(L)
print(D)
print(U)
Admin1
Related