Articles of numpy scipy

¿Cómo multiplicar elementalmente una matriz scipy.sparse por una matriz densa 1d emitida?

Supongamos que tengo una matriz dispersa 2d. En mi caso de uso real, tanto el número de filas como las columnas son mucho más grandes (por ejemplo, 20000 y 50000), por lo que no cabe en la memoria cuando se utiliza una representación densa: >>> import numpy as np >>> import scipy.sparse as ssp >>> […]

La forma más sencilla de resolver ecuaciones matemáticas en Python

Quiero resolver un conjunto de ecuaciones, lineales, o algunas veces cuadráticas. No tengo un problema específico, pero a menudo, he estado en esta situación a menudo. Es fácil usar wolftwiglpha.com , el equivalente web de Mathematica, para resolverlos. Pero eso no proporciona la comodidad y la conveniencia de un shell iPython. ¿Existe una biblioteca simple […]

Progtwig de resolución cuadrática (QP) que solo depende de NumPy / SciPy?

Me gustaría que los estudiantes resuelvan un progtwig cuadrático en una tarea sin que tengan que instalar software adicional como cvxopt, etc. ¿Hay alguna implementación de Python disponible que solo dependa de NumPy / SciPy?

Mejorando el rendimiento de Numpy

Me gustaría mejorar el rendimiento de la convolución utilizando python, y esperaba obtener alguna información sobre la mejor manera de mejorar el rendimiento. Actualmente estoy usando scipy para realizar la convolución, usando un código parecido al siguiente fragmento de código: import numpy import scipy import scipy.signal import timeit a=numpy.array ( [ range(1000000) ] ) a.reshape(1000,1000) […]

Concatene matrices dispersas en Python usando SciPy / Numpy

¿Cuál sería la forma más eficiente de concatenar matrices dispersas en Python usando SciPy / Numpy? Aquí utilicé lo siguiente: >>> np.hstack((X, X2)) array([ <49998×70000 sparse matrix of type '’ with 1135520 stored elements in Compressed Sparse Row format>, <49998×70000 sparse matrix of type '’ with 1135520 stored elements in Compressed Sparse Row format>], dtype=object) […]

Regresión lineal de matrices que contienen NAN en Python / Numpy

Tengo dos matrices, decir varx y variar. Ambos contienen valores NAN en varias posiciones. Sin embargo, me gustaría hacer una regresión lineal en ambos para mostrar cuánto se correlacionan las dos matrices. Esto fue muy útil hasta ahora: http://glowingpython.blogspot.de/2012/03/03/linear-regression-with-numpy.html Sin embargo, usando esto: slope, intercept, r_value, p_value, std_err = stats.linregress(varx, vary) resulta en nans para […]

Numpy / Polyfit: suprima la impresión del mensaje de error Intel MKL

Estoy calculando un polyfit varias veces durante un progtwig, y ​​algunas de mis entradas son np.nan y obtendré los problemas de algoritmo. Lo sé, y en esta aplicación no me importa. Cuando las cosas se complican, esto se imprime en la consola: Intel MKL ERROR: Parameter 4 was incorrect on entry to DELSD. Simplemente quiero […]

Diferentes intervalos para la cuadratura de Gauss-Legendre en números.

¿Cómo podemos usar el paquete NumPy numpy.polynomial.legendre.leggauss en intervalos distintos a [-1, 1] ? El siguiente ejemplo compara scipy.integrate.quad con el método de Gauss-Legendre en el intervalo [-1, 1] . import numpy as np from scipy import integrate # Define function and interval a = -1. b = 1. f = lambda x: np.cos(x) # […]

Cómo transformar numpy.matrix o array en una matriz dispersa escéptica

Para la matriz dispersa de SciPy, se puede usar todense() o toarray() para transformar a matriz o matriz NumPy. ¿Cuáles son las funciones para hacer lo inverso? Busqué, pero no tengo idea de qué palabras clave deberían ser el éxito correcto.

¿Escasa matriz / matriz 3d en Python?

En scipy, podemos construir una matriz dispersa usando scipy.sparse.lil_matrix () etc. Pero la matriz está en 2d. Me pregunto si hay una estructura de datos existente para la matriz / matriz 3d dispersa (tensor) en Python? ps Tengo muchos datos dispersos en 3D y necesito un tensor para almacenar / realizar la multiplicación. ¿Alguna sugerencia […]