Как эффективно решить миллионы независимых линейных систем 2x2 с дополнительными матрицами ⇐ Python

[Anonymous]

мне нужно решить миллионы до десятков миллионов независимых 2-переменных линейных систем, представленных в форме дополненной матрицы. Вход состоит из коэффициентов векторов A, B, C, A ', B', C '(каждый из которых содержит миллионы до сотен миллионов элементов). При каждом индексе I система определяется как: < /p>
Ax + By = -C
A'x + B'y = -C'
< /code>
Текущий подход Python (неэффективен):
Использование для цикла с симпью для последовательного решения каждой системы: < /p>
from sympy import symbols, Eq, solve

x = []
y = []
for i in range(len(A)):
sol = solve([
Eq(A*x_sym + B*y_sym, -C),
Eq(A_prime*x_sym + B_prime[i]*y_sym, -C_prime[i])
], [x_sym, y_sym])
x.append(sol[x_sym])
y.append(sol[y_sym])
< /code>
Этот подход становится непрактичным для больших наборов данных (например, 10m+ записи) из -за чрезмерного времени вычисления. < /p>
Как я могу написать сценарий dolphindb в: < /p>

Принять векторы ввода a, c, c, a_prime, b_prime, c_prime, c_prime, c_prime, c_prime, c_prime. /> Вычислить x и y с использованием векторизованных /параллелизированных операций?
Редакция краев
(e.g., единственные матрицы)?>

Подробнее здесь: https://stackoverflow.com/questions/796 ... th-augment