Articles of optimization

¿Por qué obtengo resultados incorrectos de scipy.optimize.fmin?

import pandas as pd from scipy.optimize import fmin data = pd.DataFrame({‘DIV’: [1,2,3]*3, ‘MONTH’: [‘May’,’May’,’May’,’June’,’June’,’Jun’,’Jul’,’Jul’,’Jul’], ‘C’:[8]*9, ‘U’:[3,2,1]*3, ‘S’:[9]*9}) data.to_csv(r’C:\Users\mbabski\Documents\Unit Plan Summer 2016\data_test.csv’) def return_array(x): return x.values def mape(c,u,s,r): #returns an array of line level Mean Absolute Percentage Errors p = c + u * r m = abs(1.0-(p/s)) return m def e(c,u,s,r): #calculates average of the […]

Python scipy.optimize.fmin_l_bfgs_b produce un error

Mi código es implementar un algoritmo de aprendizaje activo, utilizando la optimización L-BFGS. Quiero optimizar cuatro parámetros: alpha , beta , w y gamma . Sin embargo, cuando ejecuto el siguiente código, recibí un error: optimLogitLBFGS = sp.optimize.fmin_l_bfgs_b(func, x0 = x0, args = (X,Y,Z), fprime = func_grad) File “C:\Python27\lib\site-packages\scipy\optimize\lbfgsb.py”, line 188, in fmin_l_bfgs_b **opts) File […]

scipy.optimize.minimize realizar un seguimiento de la función objective

Estoy trabajando con scipy.optimize.minimize , y estoy optimizando 3 parámetros con una función como esta def foo(A, x, y, z): test = my_function(A[0], A[1], A[2], x, y, z) return test En esta respuesta, encontré algunas ideas: ¿Cómo mostrar el progreso de la función scipy.optimize? Entonces se me ocurrió esta función: def callbackF(Xi, x, y, z) […]

Comparando Python, Numpy, Numba y C ++ para la multiplicación de matrices

En un progtwig en el que estoy trabajando, necesito multiplicar dos matrices repetidamente. Debido al tamaño de una de las matrices, esta operación toma algún tiempo y quería ver qué método sería el más eficiente. Las matrices tienen dimensiones (mxn)*(nxp) donde m = n = 3 y 10^5 < p < 10^6 . Con la […]

Dile a scipy.optimize.minimize que falle

Estoy usando scipy.optimize.minimize para la optimización sin restricciones de una función objective que recibe un par de parámetros y ejecuta una simulación numérica compleja basada en estos parámetros. Esta simulación no siempre converge, en cuyo caso hago que la función objective devuelva inf, en algunos casos, en otros NaN. Pensé que este truco evitaría que […]

¿Hay una manera más rápida de agregar dos matrices numpy 2-d

Digamos que tengo dos matrices numpy 2-d grandes de las mismas dimensiones (digamos 2000×2000). Quiero sumrlos a los elementos. Me preguntaba si hay una forma más rápida que np.add () Edición: estoy agregando un ejemplo similar de lo que estoy usando ahora. ¿Hay una manera de acelerar esto? #a and b are the two matrices […]

Rápido (er) numpy ¿Fancy indexación y reducción?

Estoy tratando de usar y acelerar la indexación de fantasía para “unir” dos matrices y sumr uno de los ejes de los resultados. Algo como esto: $ ipython In [1]: import numpy as np In [2]: ne, ds = 12, 6 In [3]: i = np.random.randn(ne, ds).astype(‘float32’) In [4]: t = np.random.randint(0, ds, size=(1e5, ne)).astype(‘uint8’) […]

Cómo crear un optimizador en Tensorflow

Quiero escribir un nuevo algoritmo de optimización para mi red en Tensorflow. Espero implementar el algoritmo de optimización Levenberg Marquardt , que ahora está excluido de la API de TF. Encontré poca documentación sobre cómo escribir un optimizador personalizado, así que pregunto si alguien puede darme algún consejo. Gracias.

Optimización del portafolio SciPy con restricciones a nivel de la industria

Intento optimizar una asignación de ponderación de cartera aquí que maximice mi función de retorno al limitar el riesgo. No tengo ningún problema para encontrar el peso optimizado que rinde a mi función de devolución por la simple restricción de que la sum de todo el peso es igual a 1, y hacer que la […]

Python App Engine webapp2 lento para enrutar

Tengo una aplicación de Python App Engine que sirve aproximadamente 3M solicitudes por día. Estoy tratando de optimizar la aplicación para guardar en mi factura de hospedaje ridículamente incrementada. November, App Engine Frontend Instances: 12924.391 Hours, $604.22 Tengo la solicitud dirigida a algunas llamadas memcached, pero ahora me he dado cuenta de que generalmente toma […]