NumPy: получить нижний треугольник L в разложении Холецкого данного массива

NumPy: линейная алгебра Упражнение-16 с решением

Напишите программу NumPy, чтобы получить нижнетреугольный L в разложении Холецкого данного массива.

Пример решения :

Код Python:

import numpy as np
a = np.array([[4, 12, -16], [12, 37, -53], [-16, -53, 98]], dtype=np.int32)
print("Original array:")
print(a)
L = np.linalg.cholesky(a)
print("Lower-trianglular L in the Cholesky decomposition of the said array:")
print(L)

Пример вывода:

 Исходный массив:
[[4 12 -16]
 [12 37 -53]
 [-16 -53 98]]
Нижнетреугольная L в разложении Холецкого указанного массива:
[[2. 0. 0.]
 [6. 1. 0.]
 [-8. -5. 3.]]

Редактор кода Python:

Есть другой способ решить это решение? Внесите свой код (и комментарии) через Disqus.

Предыдущий: Напишите программу NumPy, вычисляющую сумму диагонального элемента данного массива.
Далее: Напишите программу NumPy, чтобы получить факторизацию qr данного массива.