Я читаю код для подпрограмм оптимизации (Nelder Mead, SQP …). Языки – это C ++, Python. Я наблюдаю, что часто выполняется преобразование из double в float, или методы дублируются с помощью double resp. float аргументы. Почему это выгодно в коде подпрограмм оптимизации, и это важно? В моем собственном коде на C ++ я должен быть […]
ЗАДНИЙ ПЛАН Проблема, с которой я работаю, заключается в следующем: В контексте эксперимента, который я разрабатываю для своих исследований, я создаю большое количество больших (длина 4M) массивов, которые несколько разрежены и, следовательно, могут храниться как экземпляры scipy.sparse.lil_matrix или просто как экземпляры scipy.array (прирост / потеря пространства здесь не проблема). Каждый из этих массивов должен быть […]
Есть ли способ «тихо» использовать scipy.optimize.fmin? То есть, он не печатает, например Optimization terminated successfully. Current function value: 0.000000 Iterations: 13 Function evaluations: 30 при запуске кода? Это может быть полезно, например, для циклов, которые сохраняются некоторое время.
Из-за случайного избиения клавиатуры я заметил, что в SciPy есть переменная, называемая i , которая назначается строке '6' . (Может отличаться на других машинах?) Я попытался использовать встроенные функции справки, но ничего не назначено для scipy.i поскольку оно относится только к строке. Я также искал документы и Google, но ничего не получилось. Может быть, это […]
Обычно используется import numpy as np для импорта модуля numpy. Существуют ли общие соглашения для обозначения? Что касается других модулей, в частности, от научных вычислений, таких как scipy , sympy и pylab или подмодулей, таких как scipy.sparse .
Мне было интересно, какая разница между двумя методами scipy.optimize.leastsq и scipy.optimize.least_squares ? Когда я их реализую, они дают минимальные различия в chi ^ 2: >>> solution0 = ((p0.fun).reshape(100,100)) >>> # p0.fun are the residuals of my fit function np.ravel'ed as returned by least_squares >>> print(np.sum(np.square(solution0))) 0.542899505806 >>> solution1 = np.square((median-solution1)) >>> # solution1 is the […]
У меня есть большой набор данных, в которых мне нужно сравнить расстояния набора образцов из этого массива со всеми остальными элементами массива. Ниже приведен очень простой пример моего набора данных. import numpy as np import scipy.spatial.distance as sd data = np.array( [[ 0.93825827, 0.26701143], [ 0.99121108, 0.35582816], [ 0.90154837, 0.86254049], [ 0.83149103, 0.42222948], [ 0.27309625, […]
В python есть функция scipy.ndimage.morphology модуле scipy.ndimage.morphology . Я применил его к простому случаю, чтобы вычислить расстояние от одной ячейки в масках с маской. Однако функция удаляет маску массива и вычисляет, как и ожидалось, евклидово расстояние для каждой ячейки с ненулевым значением из контрольной ячейки с нулевым значением. Ниже приведен пример, который я привел в […]
Я пытаюсь импортировать scipy.optimize с помощью Python 3.3.1 в Windows 8. Я использую scipy-0.12.0. Когда я пытаюсь импортировать, Python возвращает следующую ошибку: >>> import scipy.optimize Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python33\lib\site-packages\scipy\optimize\__init__.py", line 146, in <module> from .optimize import * File "C:\Python33\lib\site-packages\scipy\optimize\optimize.py", line 34, in < module> from .linesearch […]
Я столкнулся с следующим вопросом, пытаясь вставить два больших CSR-матрицы: /usr/lib/python2.7/dist-packages/scipy/sparse/coo.pyc in _check(self) 229 raise ValueError('negative row index found') 230 if self.col.min() < 0: –> 231 raise ValueError('negative column index found') 232 233 def transpose(self, copy=False): ValueError: negative column index found Я могу воспроизвести эту ошибку очень просто, пытаясь преобразовать большую матрицу lil в матрицу […]