Articles of numpy scipy

Error de memoria de corte de matriz dispersa

Tengo una matriz dispersa csr : <681881×58216 sparse matrix of type '’ with 2867209 stored elements in Compressed Sparse Row format> Y quiero crear una nueva matriz de partición como una porción de csr : csr_2 = csr[1::2,:] . Problema: cuando solo tengo la matriz csr , la memoria RAM de mi servidor está ocupada […]

Agregar una columna de ceros a un csr_matrix

Tengo un csr_matrix disperso de csr_matrix , y me gustaría agregar algunas columnas con solo ceros a la derecha de la matriz. En principio, las matrices indptr , los indices y los data siguen igual, así que solo quiero cambiar las dimensiones de la matriz. Sin embargo, esto parece no estar implementado. >>> A = […]

error criptico “no se pudo convertir el escalar entero” error

Estoy construyendo un vector disperso usando un scipy.sparse.csr_matrix así: csr_matrix((values, (np.zeros(len(indices)), indices)), shape = (1, max_index)) Esto funciona bien para la mayoría de mis datos, pero ocasionalmente obtengo un ValueError: could not convert integer scalar . Esto reproduce el problema: In [145]: inds Out[145]: array([ 827969148, 996833913, 1968345558, 898183169, 1811744124, 2101454109, 133039182, 898183170, 919293479, 133039089]) […]

SymPy / SciPy: resolviendo un sistema de ecuaciones diferenciales ordinarias con diferentes variables

Soy nuevo en SymPy y en Python en general, y actualmente estoy trabajando con Python 2.7 y SymPy 0.7.5 con el objective de: a) leer un sistema de ecuaciones diferenciales de un archivo de texto b) resolver el sistema Ya leí esta pregunta y esta otra pregunta , y son casi lo que estoy buscando, […]

Python: ajuste multivariable no lineal como matlabs nlinfit

Solía ​​usar Matlab para realizar ajustes no lineales utilizando la función nlinfit . Esto me permitió crear un ajuste para un vector de respuestas observadas a dos predictores. En aras de los argumentos, digamos recuperación de Cu a través de un proceso de separación basado en el grado de alimentación y la velocidad de alimentación. […]

Corte eficiente de matrices usando la multiplicación de matrices, con Python, NumPy, SciPy

Quiero remodelar un 2d scipy.sparse.csr.csr_matrix (llamémoslo A ) a un 2d numpy.ndarray (llamémoslo B ). A podría ser >shape(A) (90, 10) entonces B debería ser >shape(B) (9,10) donde cada 10 filas de A se reformarán en un nuevo valor nuevo, a saber, el máximo de esta ventana y columna. El operador de la columna no […]

Cómo resolver relaciones de recurrencia en Python

Estoy tratando de escribir código para dar respuestas numéricas a una relación de recurrencia. La relación en sí es simple y se define de la siguiente manera. La variable x es un entero p (i) = p (i + 2) / 2 + p (i-1) / 2 si i> 0 e i <x p (0) […]

Cálculo de la varianza de la imagen del python.

¿Existe una manera fácil de calcular un filtro de variación en ejecución en una imagen usando Python / NumPy / Scipy? Al ejecutar la imagen de varianza me refiero al resultado de calcular la sum ((I – mean (I)) ^ 2) / nPixeles para cada subventana I en la imagen. Como las imágenes son bastante […]

Manera eficiente de normalizar una matriz escasa Scipy

Me gustaría escribir una función que normalice las filas de una matriz dispersa grande (de manera que sumen una). from pylab import * import scipy.sparse as sp def normalize(W): z = W.sum(0) z[z < 1e-6] = 1e-6 return W / z[None,:] w = (rand(10,10)<0.1)*rand(10,10) w = sp.csr_matrix(w) w = normalize(w) Sin embargo, esto da la […]

¿Hay soporte para matrices dispersas en Python?

¿Hay soporte para matrices dispersas en python? ¿Posiblemente en entumecido o en scipy?