Articles of numpy de

numpy.polyfit no tiene la palabra clave ‘cov’

Estoy tratando de usar polyfit para encontrar la mejor línea recta para un conjunto de datos, pero también necesito saber la incertidumbre sobre los parámetros, así que también quiero la matriz de covarianza. La documentación en línea sugiere que escribo: polyfit (x, y, 2, cov = Verdadero) pero esto da el error: TypeError: polyfit () […]

Numpy: Elementos enmascarados en computación.

Tengo una función para construir un polinomio a partir de una x dada: [1, x ^ 2, x ^ 3, x ^ 4, …, x ^ grado] def build_poly(x, degree): “””polynomial basis functions for input data x, for j=0 up to j=degree.””” D = len(x) polyome = np.ones((D, 1)) for i in range(1, degree+1): polyome […]

Vectorizar u optimizar un bucle donde cada iteración depende del estado de la iteración anterior

Tengo un algoritmo que estoy implementando en python. El algoritmo podría ejecutarse 1.000.000 veces, así que quiero optimizarlo tanto como sea posible. La base en el algoritmo es tres listas ( energy , point y lista de valList ) y dos contadores p y e . Las dos listas de energy y point contienen números […]

Vectorizar Pandas Dataframe en Numpy Array

Tengo un problema donde necesito convertir un dataframe de pandas en una lista de listas. Muestra: import pandas as pd df = pd.DataFrame([[1,2,3],[2,2,4],[3,2,4]]) Sé que existe la función as_matrix () que se muestra a continuación: df.as_matrix(): # result:array([[1, 2, 3], [2, 2, 4], [3, 2, 4]]) Sin embargo, necesito algo en este formato. [array([[1], [2], […]

La intersección de las curvas de Python con fsolve () y los argumentos de la función usando numpy

Estoy tratando de usar fsolve como se cita aquí: http://glowingpython.blogspot.gr/2011/05/hot-to-find-intersection-of-two.html , Para encontrar la intersección entre dos curvas. Ambas curvas son básicamente dos matrices de flotadores. El primero de ellos es una matriz de una dimensión Pmech ( Pmech(x) ) y la segunda es una matriz de dos dimensiones Pair ( Pair(x,y) ) El eje […]

¿Cómo transmitir una función a través de una matriz numpy, cuando dtype = object?

Si tengo una matriz de valores numéricos, los cuales tenían que usar punteros de objeto en lugar de valores como tipo de datos, debido a longitudes de vectores desiguales: In [145]: import numpy as np In [147]: a = np.array([[1,2],[3,4,5]]) In [148]: a Out[148]: array([[1, 2], [3, 4, 5]], dtype=object) In [150]: np.sin(a) ————————————————————————— AttributeError […]

Reducir la huella de memoria numpy en aplicaciones de larga ejecución

En mi aplicación se generan y se llenan con datos cien arreglos numpy (1000 elementos complejos cada uno). Luego, en muchas iteraciones, los elementos de la matriz se modifican una y otra vez. Después de la generación inicial, el monitor del sistema reporta alrededor de 50 Mb de uso de RAM. Aunque no estoy generando […]

¿Por qué no funciona una asignación para matrices numpy de doble corte?

¿Por qué las siguientes líneas no funcionan como espero? import numpy as np a = np.array([0,1,2,1,1]) a[a==1][1:] = 3 print a >>> [0 1 2 1 1] # I would expect [0 1 2 3 3] ¿Es este un ‘error’ o hay otra forma recomendada para esto? Por otro lado, los siguientes trabajos: a[a==1] = […]

Multiprocesamiento: resultado NULO sin error en PyObject_Call

Aquí hay un progtwig de ejemplo donde uso multiprocesamiento. Los cálculos se realizan con multiprocessing.Process y los resultados se recostackn mediante multiprocessing.Queue . #THIS PROGRAM RUNS WITH ~40Gb RAM. (you can reduce a,b,c for less RAM #but then it works for smaller values) #PROBLEM OCCURS ONLY FOR HUGE DATA. from numpy import * import multiprocessing […]

Mapeo de codificaciones de caracteres a bytes máximos por carácter

Estoy buscando una tabla que asigne una encoding de caracteres dada a los bytes (máximo, en el caso de las codificaciones de longitud variable) por carácter. Para codificaciones de ancho fijo, esto es bastante fácil, aunque no sé, en el caso de algunas de las codificaciones más esotéricas, cuál es ese ancho. Para UTF-8 y […]