Tag: scipy

Список списка ndarray

Я пытаюсь использовать кластеры kmeans в scipy, точно такой, который присутствует здесь: http://docs.scipy.org/doc/scipy/reference/generated/scipy.cluster.vq.kmeans.html#scipy.cluster.vq.kmeans Я пытаюсь преобразовать список списка, например: data without_x[ [0, 0, 0, 0, 0, 0, 0, 20.0, 1.0, 48.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1224.0, 125.5, 3156.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, […]

Поиск двух хвостового значения P из t-распределения и степеней свободы в Python

Как определить значение P t-разброса с n степенями свободы. Исследования по этому вопросу указывают мне на этот ответ обмена на стек: https://stackoverflow.com/a/17604216 Я предполагаю, что np.abs (tt) является значением T, но как я могу работать в степенях свободы, является то, что n-1? заранее спасибо

Поиск индексов строк, где есть ненулевые записи в разреженной csc_matrix

У меня есть массив numpy, X: type(X) >>> <class 'scipy.sparse.csc.csc_matrix'> Меня интересует поиск индексов строк, в которых есть ненулевые записи, в 0-м столбце. Я пытался: getcol = X.getcol(0) print getcol который дает мне: (0, 0) 1 (2, 0) 1 (5, 0) 10 Это здорово, но то, что я хочу, это вектор, который имеет 0, 2, […]

Scipy LinearOperator dtype неуказан

Я пытаюсь использовать команду gmres . Моя матрица плотная и большая, поэтому мой план состоял в том, чтобы использовать команду LinearOperator для возврата только векторного вектора. См. Мой предыдущий вопрос здесь LinearOperator с двумя входами . С помощью этого вопроса я смог построить объект LinearOperator который успешно выполняет вычисление A*x где A – матрица, а […]

Специальная вычитание Scipy Sparse Matrix

Я делаю проект, и в нем много вычислений матрицы. Я ищу умный способ ускорить мой код. В моем проекте я имею дело с разреженной матрицей размером 100Mx1M с значениями не нулей около 10M. Пример ниже – это просто увидеть мою мысль. Скажем, у меня есть: Вектор v размера (2) Вектор c размера (3) Редкая матрица […]

SciPy optimize.fmin ValueError: массив с нулевым размером до максимума операции сокращения, который не имеет идентификатора

UPDATE2: Лучшее название (теперь, когда я понимаю проблему) будет: Каков правильный синтаксис ввода в scipy optimize.fmin? UPDATE: запрошен исполняемый код, поэтому определения функций были заменены исполняемым кодом. Примеры входных данных были жестко закодированы как «данные» массива numpy. Я пытаюсь оптимизировать функцию с scipy, но я действительно застрял и должен просить помощи. Массив нулевой длины передается […]

создание входных файлов моделирования с помощью Python

Я использую научный пакет моделирования, который требует нескольких текстовых входных файлов для каждого эксперимента. Эти файлы могут быть довольно длинными и иметь в них много разделов шаблонов; однако конкретные значения, специфичные для эксперимента, должны вводиться во многих местах в этих файлах. Я хотел бы автоматизировать создание этих файлов и сделать это таким образом, чтобы это […]

Почему мой 2D-интерполятор генерирует матрицу с обведенными осями в SciPy?

Я решаю дифференциальное уравнение с векторными входами y '= f (t, y), y (t_0) = y_0 где y0 = y (x) используя явный метод Эйлера, в котором говорится, что y_ (i + 1) = y_i + h * f (t_i, y_i) где t – вектор времени, h – размер шага, f – правая часть дифференциального […]

Эквивалентная функция Matlab imfilter в Python

Я знаю, что эквивалентные функции conv2 и corr2 MATLAB – scipy.signal.correlate и scipy.signal.convolve . Но функция imfilter обладает свойством иметь дело с внешней imfilter массива. Как symmetric , replicate и circular . Может ли Python делать то, что

Минимизация функции с помощью scipy.optimize параллельно с ipython

Я нашел кучу вопросов по подобным проблемам (например, это), и проблема с расчетом закрытия, похоже, повсюду в отношениях с iPython.Parallel, но я не мог обойти эту проблему в любом месте. Поэтому моя проблема заключается в следующем: Я хочу решить для нулей функцию f( a, b ) через несколько значений b используя ipcluster. f является сложной […]

Python - лучший язык программирования в мире.