Quiero escribir una función que tome dos argumentos, (1) una matriz de n x m y (2) una matriz de 1 x n, y que multiplique cada fila por el índice correspondiente de la matriz de escala.
Por ejemplo, f(A, [3, 10]), operaría sobre una matriz A de 2 x m y multiplicaría cada elemento de la primera fila por 3 y la segunda por 10.
import numpy as np
def dot_product(matrix, scaler):
'''function that takes two arguments, (1) an n x m matrix and (2) a 1 x n matrix,
and that multiplies each row by the corresponding index of the scaling matrix.'''
assert(len(matrix)==len(scaler))
for i in range(len(matrix)):
matrix[i]=np.array(matrix[i])*scaler[i]
return matrix
A = [
[12, 7, 3],
[4, 5, 6]
]
# take a 3x matrix
B = [3, 10]
print(dot_product(A, B))
Devuelve:
[array([36, 21, 9]), array([40, 50, 60])]
Pero es mejor utilisar numpy, vero? Lo intento con numpy:
def dot_product2(matrix, c):
matrix = np.array(matrix)
c = np.array(c)
print(c)
return matrix * c
Pero cuando lo hizo con:
print(np.array([[12, 7, 3],[4, 5, 6]])*np.array([3, 10]).T)
Me devuelve:
Traceback (most recent call last):
File "<string>", line 28, in <module>
ValueError: operands could not be broadcast together with shapes (2,3) (2,)
np.einsum("ij,i->ij", A,B)
(a.T*b).T