Estoy tratando de realizar la multiplicación de matrices con python. El objetivo del programa es tener tres matrices como parámetros: A,B y res (matriz de resultados) que se completará con el resultado de multiplicar las matrices representadas por las listas de listas. El programa debe aceptar tanto listas de listas como matrices de NumPy.
Escribí el siguiente código:
'''
multMatrices(a, b, res): Completa los elementos de res
con el resultado de multiplicar las matrices representadas por
las listas de listas a y b.
'''
def multMatrices(A, B, res):
filas_A = len(A)
cols_A = len(A[0])
filas_B = len(B)
cols_B = len(B[0])
res = matrizCeros()
if cols_A != filas_B:
print("No es posible multiplicar las matrices. Dimensiones incorrectas.")
return
# Se crean las dimensiones de la matriz
# Lss dimensiones deben ser filas_A x cols_B
res = [[0 for row in range(cols_B)] for col in range(filas_A)]
#print(res) #Cuidado con este print
#preguntar por las matrices de numpy con instance
if isinstance(A,np.matrix) and isinstance(B,np.matrix):
#res = np.matmul(A, B)
np.copyto(res,np.matmul(A, B))
else:
#Si isinstance devuelve False, se ejecuta esta sección original
for i in range(filas_A):
for j in range(cols_B):
for k in range(cols_A):
res[i][j] += A[i][k] * B[k][j]
Cuándo pruebo con listas de listas el programa funciona, pero no con matrices de NumPy:
C = [[1,1,1],[1,1,1],[1,1,1]]
D = [[1,1,1],[1,1,1],[1,1,1]]
E = matrizCeros()
X = np.matrix([[1, 2], [3, 4]])
Y = np.matrix([[1, 2], [3, 4]])
Z = matrizCeros()
print("Resultado Mult:",multMatrices(C, D, E))
print("Resultado Mult:",multMatrices(X, Y, Z))
#Resultados:
Resultado Mult: [[3, 3, 3], [3, 3, 3], [3, 3, 3]]
No es posible multiplicar las matrices. Dimensiones incorrectas.
Resultado Mult: None
No he logrado identificar que estoy haciendo mal. Agradezco la realimentación brindada.