Как обновить результаты функции разрешения в исходных матрицах L и U? - PullRequest
Новая
       71

Как обновить результаты функции разрешения в исходных матрицах L и U?

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

Я могу найти каждое уравнение для решения, но я не могу обновить результаты каждого вычисления в исходной матрице.

Я попытался задать набор элементов в матрице, заменив переменные u фактическим результатом из функции решения.Но тогда он выдает следующую ошибку: TypeError: не может умножить последовательность на не-int типа 'Symbol'

Не знаю, что делать еще! 

from sympy import *
import numpy as np


#input matrix

a = np.array([[1,3,5],\
              [2,4,7],\
              [1,1,0]])


d = len(a)
d1 = d-1
d2 = (d1*(d1+1))/2
d3 = (d*(d+1))/2


symbols_dict = dict(('L%d'%k, symbols('L%d'%k)) for k in range(d2))
locals().update(symbols_dict)
symbols_dict = dict(('U%d'%k, symbols('U%d'%k)) for k in range(d3))
locals().update(symbols_dict)


# L decomposition

b = np.array([[1,0,0],\
              [L0,1,0],\
              [L1,L2,1]])

# U decomposition

c = np.array([[U0,U1,U2],\
              [0,U3,U4],\
              [0,0,U5]])


t=[[[]for i in range(d)]for j in range(d)]
s=[[[]for i in range(d)]for j in range(d)]

for j in range(d):
    for i in range(d):
        m=0
        for k in range(d):
            m = m+(b[j,k]*c[k,i])
            t[j][i] = m
        s[j][i]= solve(t[j][i]-a[j,i])
        c.itemset((j,i),s[j][i])

Это должно дать мнечисленные результаты для каждого решенного уравнения.

1 Ответ

0 голосов
/ 22 августа 2019

Исправленный код выглядит следующим образом: 

L0,L1,L2 = symbols('L0 L1 L2')

U0,U1,U2,U3,U4,U5 = symbols('U0 U1 U2 U3 U4 U5')
.
.
.


for j in range(d):
    for i in range(d):
        m=0
        for k in range(d):
            m = m+(b[j,k]*c[k,i])
            s= solve(m-a[j,i])
            if j>i:
                b[j,i]=s
            else:
                c[j,i]=s
...