Я разрабатываю некоторый более крупный код для выполнения вычислений по собственным значениям огромных разреженных матриц в контексте вычислительной физики. Я проверяю свои процедуры против простого гармонического осциллятора в одном измерении, поскольку собственные значения хорошо известны...
Моя исследовательская группа фокусируется на молекулярной динамике, которая, очевидно, может генерировать гигабайты данных как часть одной траектории, которую затем необходимо проанализировать. Некоторые из проблем, с которыми мы сталкиваемся, связаны с корреляциями в наборе данных, что означает,...
Я пишу воспроизводимую статью, и в ней есть результаты вычислений, которые генерируются скриптом Python (аналогичный скрипт MATLAB генерирует почти идентичные результаты). Я чувствую, что статья будет легче понять читателям, если они смогут сопоставить вычисления в статье с вычислениями в коде. В...
Я хочу реализовать следующее выражение в Python: Икся= ∑J = 1я - 1Кя - j , jaя - JaJ,Иксязнак равноΣJзнак равно1я-1Кя-J,Jaя-JaJ, x_i = \sum_{j=1}^{i-1}k_{i-j,j}a_{i-j}a_j, где ИксИксx и YYy - массивы numpy размером NNn , а ККk - массив numpy размером n × nN×Nn\times n . Размер NNn может составлять...
Кто-нибудь здесь использовал научные вычисления двойной точности с графическими процессорами нового поколения (например, K20) через Python? Я знаю, что эта технология быстро развивается, но каков наилучший способ сделать это в настоящее время? GPU не входит в число популярных научных библиотек...
Я использовал разные методы для вычисления как ранга матрицы, так и решения матричной системы уравнений. Я наткнулся на функцию linalg.svd. Сравнивая это с моими собственными усилиями по решению системы с устранением по Гауссу, она кажется более быстрой и точной. Я пытаюсь понять, как это возможно....
Я испытываю некоторое разочарование по поводу того, как Matlab обрабатывает числовую интеграцию против Scipy. Я наблюдаю следующие различия в моем тестовом коде ниже: Версия Matlab работает в среднем в 24 раза быстрее, чем мой эквивалент на Python! Версия Matlab способна вычислить интеграл без...
У меня проблемы с реализацией функции численно. Он страдает от того факта, что при больших входных значениях результат очень большое число раз очень маленькое число. Я не уверен, что термин «отмена» является правильным, поэтому, пожалуйста, поправьте меня, если это так. Доказательства того, что...
Мне нужно решить s.t.minx∥Ax−b∥22,∑ixi=1,xi≥0,∀i.minx‖Ax−b‖22,s.t.∑ixi=1,xi≥0,∀i.\begin{alignat}{1} & \min_{x}\|Ax - b\|^2_{2}, \\ \mathrm{s.t.} & \quad\sum_{i}x_{i} = 1, \\ & \quad x_{i} \geq 0, \quad \forall{i}. \end{alignat} Я думаю, что это квадратичная проблема, которая должна быть решена с...
Суть моего вопроса заключается в следующем: у меня есть система двух ODE. Один имеет ограничение начального значения, а другой имеет ограничение конечного значения. Это можно рассматривать как единую систему с ограничением начального значения для некоторых переменных и ограничением конечного...
Я пытаюсь выяснить, существует ли более быстрый способ вычисления всех собственных значений и собственных векторов очень большой и разреженной матрицы смежности, чем использование scipy.sparse.linalg.eigsh. Насколько я знаю, этот метод использует только разреженность и атрибуты симметрии матрицы....
Я новичок в программировании и пытаюсь решить свою первую большую проблему и написать свою первую большую программу. Я искал примеры кода с открытым исходным кодом для изучения, но до сих пор нашел код только на тех языках, которые я не до конца понимаю, или это делает связанные, но слишком далекие...
У меня есть проблемы определяется уравнением Пуассона в двух измерениях −∇2u=f(x,y),inΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y), \; in \; \Omega у меня есть измерения двух компонентов градиента∂u/∂x∂u/∂x\partial{u}/\partial{x} и∂u/∂y∂u/∂y\partial{u}/\partial{y} вдоль некоторой части...
У меня есть набор данных 3-мерных точек, для которых я хотел бы построить сетку, используя python. Все программное обеспечение, которое я видел, требует, чтобы вы предоставили края. Есть ли в Python программа, которая принимает в качестве входных данных набор точек в 3D и выводит сетку? Если...
Пожалуйста, перечислите пакет Python (petsc4py и т. Д.) И разреженные прямые решатели, которые он поддерживает. Один (Community-Wiki) ответ на пакет,
Я должен написать код конечного объема для магнитной гидродинамики (МГД). Я написал числовой код раньше, но не в этом масштабе. Я просто хотел спросить, какой будет хорошим выбором, используя структуру данных (объектно-ориентированный подход) с классами или просто используя несколько массивов для...
CVXOPT: http://abel.ee.ucla.edu/cvxopt/index.html OpenOpt: http://openopt.org/Welcome Какая связь между ними? Каковы их преимущества / недостатки, соответственно? Кстати, есть ли еще какая-нибудь высококачественная библиотека выпуклой оптимизации общего назначения для Python / C ++, на которую...
Насколько медленнее petsc4py против c / c ++ / fortran? Я понимаю, что это будет в значительной степени зависеть от выполняемого кода, но как насчет чего-то простого, такого как матрично-векторный...
Существуют ли алгоритмы обнаружения сообщества для двудольных графов (двухрежимные сети), реализованные в igraph, networkX, R или Python и т. Д.? В частности, существует ли такая реализация, в которой можно было бы ограничить обнаружение сообществ только одним из двух...
Я пытаюсь реализовать метод Рунге-Кутты четвертого порядка для решения ODE первого порядка в Python, т.е. . Я понимаю, как работает метод, но я пытаюсь написать эффективный алгоритм, который минимизирует количество вычислений поскольку это довольно дорого. Мне сказали, что можно повторно...