Python-сообщество

Sadko · Ноя. 15, 2018 10:59:27

Здравствуйте.

Возникла задача, связанная с полями Галуа. Решил заодно язык новый освоить.
Определил новый класс и перегрузил в нем операции сложения/умножения/и т.д.
Создал матрицу из объектов нового типа для проверки:

 import numpy as np
import gfpy  as gf
import random as rnd
a=[[gf.GF(rnd.randint(0,7)) for v in range(sz)] for i in range(sz)]
b=[[gf.GF(rnd.randint(0,7)) for v in range(sz)] for i in range(sz)]
c=[gf.GF(rnd.randint(0,7)) for v in range(sz)]
a=np.matrix(a,dtype=gf.GF)
b=np.matrix(b,dtype=gf.GF)
a+b
c*a

Сложение и умножение работают как ожидалось (умножение и сложение используется из моего класса), но

 np.linalg.inv(a)

приводит к ошибке:

File “C:\ProgramData\Anaconda3\lib\site-packages\numpy\linalg\linalg.py”, line 528, in inv
ainv = _umath_linalg.inv(a, signature=signature, extobj=extobj)
TypeError: No loop matching the specified signature and casting
was found for ufunc inv

Вопрос: inv() работает только со стандартными типами из-за оптимизации или чего-то не хватает в моем классе?

Отредактировано Sadko (Ноя. 15, 2018 11:43:43)

Sadko · Ноя. 15, 2018 11:01:20

-

Отредактировано Sadko (Ноя. 15, 2018 12:58:44)

Sadko · Ноя. 15, 2018 11:03:29

Попробовал scipy.linalg inv() - пишет “ValueError: object arrays are not supported”.
Нужна простая, прямолинейная, неоптимизированная реализация инверсии матрицы, пусть медленная. Поиск приводит к numpy и scipy. Может, посоветуете где искать?

Отредактировано Sadko (Ноя. 15, 2018 12:56:29)

Sadko · Ноя. 15, 2018 11:05:56

В общем нашел -

 https://stackoverflow.com/questions/32114054/matrix-inversion-without-numpy

Всем спасибо за внимание

Отредактировано Sadko (Ноя. 15, 2018 13:05:29)

Python-сообщество

Уведомления

#1 Ноя. 15, 2018 10:59:27

Переопределение умножения

#2 Ноя. 15, 2018 11:01:20

Переопределение умножения

#3 Ноя. 15, 2018 11:03:29

Переопределение умножения

#4 Ноя. 15, 2018 11:05:56

Переопределение умножения

Board footer