Гауссово исключение в Python - PullRequest
Новая
       11

Гауссово исключение в Python

0 голосов
/ 22 февраля 2019

Итак, я работал над пользовательской функцией, которая использует декомпозицию LU для решения системы линейных уравнений.Функция предполагает равное количество уравнений и неизвестных.Когда я использую 

    for k in range(0,n-1,1):
        for i in range(k+1,n,1):
            L[i,k] = A[i,k] / A[k,k] #Lower matrix, comprehension aid
            for j in range(0,n,1):
                A[i,j] = A[i,j] - (L[i,k] * A[k,j])
                d[i] = b[i] - (L[i,k] * b[j])

На матрице коэффициентов 3 x 3 (A) where A = np.array([[8,4,-1],[-2,5,1],[2,-1,6]]) the A[1,2] значение становится 0, когда я ожидал 0,75.Не уверен, что происходит.

Я использовал несколько операторов печати и несколько раз пробежал по отладчику, но безрезультатно.Однако для простой матрицы 2 x 2 код отлично работает в этом разделе об исключении Гаусса.

...