Articles of matrix

Elemento sabio poder de la matriz scipy.sparse

¿Cómo puedo elevar una matriz scipy.sparse a un poder, de forma elemental? numpy.power debería, según su manual , hacer esto, pero falla en matrices dispersas: >>> X <1353×32100 sparse matrix of type '’ with 144875 stored elements in Compressed Sparse Row format> >>> np.power(X, 2) Traceback (most recent call last): File “”, line 1, in […]

Añadiendo dos matrices dispersas `csc` de diferentes formas en python

Así que tengo dos matrices csc de diferentes formas que necesito sumr. Las matrices se ven así: current_flows = (7005, 1001) 50.0 (8259, 1001) 65.0 (14007, 1001) 45.0 (9971, 1002) 80.0 : : (69003, 211148) 0.0 result_flows = (7005, 1001) 40 (14007, 1001) 20 (9971, 1002) 35 : : (71136, 71137) 90 final_flows = current_flows […]

Difusión para subtensor creado a partir de matriz (Theano).

Quiero crear dos subtensores a partir de una matriz, utilizando índices para seleccionar las filas respectivas. Un subtensor tiene varias filas, la otra solo una, que debe transmitirse para permitir la adición de elementos. Mi pregunta es: ¿cómo indico que quiero permitir la difusión en la dimensión específica en el subtensorRight resultante de los índices […]

punto scipy.sparse extremadamente lento en Python

El siguiente código ni siquiera terminará en mi sistema: import numpy as np from scipy import sparse p = 100 n = 50 X = np.random.randn(p,n) L = sparse.eye(p,p, format=’csc’) XTdot(L).dot(X) ¿Hay alguna explicación de por qué esta multiplicación de matrices está colgando?

Sustituya la difusión numpy usando scipy.sparse.csc_matrix

Tengo en mi código la siguiente expresión: a = (b / x[:, np.newaxis]).sum(axis=1) donde b es un ndarray de forma (M, N) , y x es un ndarray de forma (M,) . Ahora, b es en realidad escasa, por lo que para la eficiencia de la memoria me gustaría sustituir en scipy.sparse.csc_matrix o csr_matrix . […]

Multiplicando elementos de columna de matriz dispersa.

Tengo una matriz csc dispersa con muchos elementos cero para los que me gustaría calcular el producto de todos los elementos de columna para cada fila. es decir: A = [[1,2,0,0], [2,0,3,0]] debe ser convertido a: V = [[2, 6]] Usando una matriz densa numpy esto se puede lograr reemplazando todos los valores cero con […]

Compruebe si hay definición positiva o semidefinitud positiva

Quiero verificar si una matriz es positiva definida o semidefinita positiva usando Python. ¿Cómo puedo hacer eso? ¿Hay una función dedicada en SciPy para eso o en otros módulos?

¿Qué está causando la desaceleración 2x ​​en mi implementación de Cython de la multiplicación de vectores matriciales?

Actualmente estoy intentando implementar la multiplicación básica de vectores matriciales en Cython (como parte de un proyecto mucho más grande para reducir la computación ) y encuentro que mi código es aproximadamente 2 Numpy.dot más lento que Numpy.dot . Me pregunto si hay algo que me estoy perdiendo que esté resultando en la desaceleración. Estoy […]

Pandas crosstab matrix dot nansum

Estoy buscando ayuda para crear un sub-dataframe a partir de un dataframe existente usando una función similar a np.nansum. Quiero convertir esta tabla en una matriz de sums de columnas no nulas: dan ste bob t1 na 2 na t2 2 na 1 t3 2 1 na t4 1 na 2 t5 na 1 2 […]

¿Cómo generar una matriz con entradas aleatorias y con restricciones en filas y columnas?

Cómo generar una matriz para que sus entradas sean números reales aleatorios entre cero y uno inclusive con la restricción adicional: la sum de cada fila debe ser menor o igual que uno y la sum de cada columna debe ser menor o igual a uno. Ejemplos: matrix = [0.3, 0.4, 0.2; 0.7, 0.0, 0.3; […]