Cholsky Implementacja w Pythonie - Solve AX = B

Używam rozkładu Cholesky'ego dla Ax=b, aby znaleźć x, wykonując L*LT=A, a następnie y=L*b i na końcu x=LT*b. Kiedy sprawdzam, czy nie. Wydaje się, że daje takie same wyniki, jak wykonanie klasycznej Ax=b. Oto mój kod:

import numpy as np
import scipy.linalg as sla

myL=np.linalg.cholesky(A)


#check_x = np.dot(A, b)
#check_x = np.dot(A,b)
check_x = sla.solve(A, b)

#check if the composition was done right
myLT=myL.T.conj() #transpose matrix
Ac=np.dot(myL,myLT) #should give the original matrix A



#y=np.dot(myL,b)
y = sla.solve_triangular(myL, b)

#x=np.dot(myL.T.conj(),y)
x = sla.solve_triangular(myLT, b)