Tag: LINEAR-ALGEBRA
Estoy tratando de calcular la ‘matriz sombrero’ en Python. Estoy utilizando la siguiente fórmula. Cuando la longitud de X es grande (digamos 60,000), estoy teniendo problemas de falta de memoria. H = Xinv(X’X)X’ ¿Existe una manera computacionalmente más eficiente de hacer esto? Aquí está un fragmento del código: import numpy . . . Read more
Estoy creando una matriz nxn con el valor más grande siendo n**n. np.linalg.solve funciona cuando n < 16, pero para n >= 16, muestra un error: numpy.core._exceptions.UFuncTypeError: No se puede convertir el dtype(‘O’) del input 0 del ufunc ‘solve1’ a dtype(‘float64’) con la regla de casting ‘same_kind’. Aquí está mi . . . Read more
Ejemplo mínimo: import numpy as np list1 = [1,3,5,7] list2 = [3,6,9,4] list3 = [6,5,3,2] arr = np.random.rand(72,22,22) pos_list = np.vstack([list1, list2, list3]).T print(pos_list) print(arr[pos_list[0][0], pos_list[0][1], pos_list[0][2]]) for i in pos_list: print(arr[i[0], i[1], i[2]]) Mi compañero y yo estamos intentando aumentar la eficiencia de la selección de puntos anterior. Se . . . Read more
Aquí está mi problema. Tengo dos matrices A y B, con entradas complejas, de dimensiones (n,n,m,m) y (n,n) respectivamente. A continuación se muestra la operación que realizo para obtener una matriz C: C = np.sum(B[:,:,None,None]*A, axis=(0,1)) Calcular lo anterior una vez tarda aproximadamente 6-8 segundos. Como tengo que calcular muchas . . . Read more
Digamos que tenemos una lista de n ecuaciones m * x = b Quiero obtener todos los puntos de esquina del polígono que tiene el eje x, el eje y y las líneas rectas como caja delimitadora. Los círculos rojos en la imagen de abajo muestran los puntos que quiero. . . . Read more