Вопросы с тегом «python»

Какой из них лучше: Фортран или Питон? И я думаю, что в обоих случаях вам нужен Gnuplot, я прав? Я сейчас работаю на машине с Windows. Я хотел бы использовать его для получения численных решений физических задач, включая моделирование по методу Монте-Карло, численное интегрирование и дифференцирование, молекулярную динамику и т. …

17 python  fortran 

Я хотел бы знать, есть ли быстрый способ вычисления евклидова расстояния двух векторов в октаве. Кажется, что для этого нет специальной функции, так что я должен просто использовать формулу с sqrt?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Недавно я сравнивал различные нелинейные решатели от scipy и был особенно впечатлен примером Ньютона-Крылова в книге рецептов Scipy, в которой они решают уравнение дифференциального уравнения второго порядка с нелинейным членом реакции примерно в 20 строках кода. Я изменил пример кода для решения нелинейного уравнения Пуассона ( также называемого уравнением Пуассона-Больцмана …

16 python  nonlinear-equations  poisson  scipy 

Я хотел бы быть независимым от коммерческого программного обеспечения для моей научной работы. Я нахожу зависимость от коммерческих пакетов, таких как Matlab и ее наборов инструментов, неудовлетворительной, потому что я не знаю, получу ли я доступ к Matlab в будущем, и потому что мне не нравится язык. Поэтому я ищу …

16 python  machine-learning 

Я начал использовать python в качестве языка программирования для выполнения всех моих заданий в CFD. У меня очень мало опыта в программировании. Я из области машиностроения и продолжаю получать высшее образование в области аэрокосмической техники. иногда вычислительный аспект CFD становится более утомительным, чем манипулирование уравнениями или математика. Каковы общие рекомендации, …

16 algorithms  python  fluid-dynamics 

Моя исследовательская группа фокусируется на молекулярной динамике, которая, очевидно, может генерировать гигабайты данных как часть одной траектории, которую затем необходимо проанализировать. Некоторые из проблем, с которыми мы сталкиваемся, связаны с корреляциями в наборе данных, что означает, что нам нужно отслеживать большие объемы данных в памяти и анализировать их, а не …

15 python  c  efficiency 

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

15 python  sparse-matrix  numpy  eigenvalues 

Кто-нибудь здесь использовал научные вычисления двойной точности с графическими процессорами нового поколения (например, K20) через Python? Я знаю, что эта технология быстро развивается, но каков наилучший способ сделать это в настоящее время? GPU не входит в число популярных научных библиотек Python numpy и scipy, и я хотел использовать theano, но, …

14 python  gpu  precision 

Я хочу реализовать следующее выражение в 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 может составлять примерно …

14 python  numpy 

Я пишу воспроизводимую статью, и в ней есть результаты вычислений, которые генерируются скриптом Python (аналогичный скрипт MATLAB генерирует почти идентичные результаты). Я чувствую, что статья будет легче понять читателям, если они смогут сопоставить вычисления в статье с вычислениями в коде. В работе предлагается абстрактный формализм, и примеры в статье должны …

14 python  software  reproducibility 

У меня проблемы с реализацией функции численно. Он страдает от того факта, что при больших входных значениях результат очень большое число раз очень маленькое число. Я не уверен, что термин «отмена» является правильным, поэтому, пожалуйста, поправьте меня, если это так. Доказательства того, что что-то пошло не так: Как я могу …

13 python  numerical-analysis  floating-point 

Я испытываю некоторое разочарование по поводу того, как Matlab обрабатывает числовую интеграцию против Scipy. Я наблюдаю следующие различия в моем тестовом коде ниже: Версия Matlab работает в среднем в 24 раза быстрее, чем мой эквивалент на Python! Версия Matlab способна вычислить интеграл без предупреждений, а Python возвращает nan+nanj Что я …

13 matlab  python  quadrature  numba 

Я использовал разные методы для вычисления как ранга матрицы, так и решения матричной системы уравнений. Я наткнулся на функцию linalg.svd. Сравнивая это с моими собственными усилиями по решению системы с устранением по Гауссу, она кажется более быстрой и точной. Я пытаюсь понять, как это возможно. Насколько я знаю, функция linalg.svd …

13 linear-algebra  python  lapack 

Я пытаюсь выяснить, существует ли более быстрый способ вычисления всех собственных значений и собственных векторов очень большой и разреженной матрицы смежности, чем использование scipy.sparse.linalg.eigsh. Насколько я знаю, этот метод использует только разреженность и атрибуты симметрии матрицы. Матрица смежности также является двоичной, что заставляет меня думать, что есть более быстрый способ …

12 linear-algebra  python  performance  eigensystem  scipy 

Мне нужно решить 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} Я думаю, что это квадратичная проблема, которая должна быть решена с помощью CVXOPT , но я не могу понять, как.

12 optimization  python  convex-optimization  least-squares  quadratic-programming