Запитання з тегом «python»

Мені хотілося б знати, чи існує швидкий спосіб обчислити евклідову відстань двох векторів в Октаві. Здається, що для цього немає спеціальної функції, тому я повинен просто використовувати формулу sqrt?

18 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 

Що краще: FORTRAN чи Python? І я здогадуюсь, що в обох випадках вам потрібен Gnuplot, я прав? Зараз я працюю на машині Windows. Я хотів би використовувати його для отримання чисельних рішень для фізичних проблем, включаючи моделювання Монте-Карло, числову інтеграцію та диференціювання, молекулярну динаміку тощо. Я побачив курс з обчислювальної …

17 python  fortran 

Нещодавно я порівнював різні нелінійні розв'язувачі з 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 

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

14 python  gpu  precision 

Я хочу реалізувати наступний вираз у Python: де x і y - нумерові масиви розміром n , а k - масивний масив розміру n × n . Розмір n може становити приблизно 10000, а функція є частиною внутрішнього циклу, який буде оцінюватися багато разів, тому важлива швидкість.хi= ∑j = 1i …

14 python  numpy 

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

14 python  software  reproducibility 

Я відчуваю деяке розчарування з приводу того, як matlab обробляє числову інтеграцію проти Scipy. У моєму тестовому коді нижче я спостерігаю такі відмінності: Версія Matlab працює в середньому в 24 рази швидше, ніж мій еквівалент python! Версія Matlab здатна обчислити інтеграл без попереджень, тоді як python повертається nan+nanj Що я …

13 matlab  python  quadrature  numba 

У мене виникають проблеми з числовою функцією. Він страждає від того, що при великих вхідних значеннях результат є дуже великим числом разів і дуже малим числом. Я не впевнений, що катастрофічне скасування - це правильний термін, тому, будь ласка, виправте мене, якщо воно є. Докази того, що щось відбувається не …

13 python  numerical-analysis  floating-point 

Я використовував різні методи для обчислення як ранжування матриці, так і рішення матричної системи рівнянь. Я натрапив на функцію linalg.svd. Порівнюючи це з моїми власними зусиллями щодо розв’язання системи з Гауссовим усуненням, це здається і швидшим, і більш точним. Я намагаюся зрозуміти, як це можливо. Наскільки мені відомо, функція linalg.svd …

13 linear-algebra  python  lapack 

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

12 linear-algebra  python  performance  eigensystem  scipy 

Мені потрібно вирішити с . т .хвх∥ A x - b ∥22,∑iхi= 1 ,хi≥ 0 ,∀ i .хвх‖Ах-б‖22,с.т.∑iхi=1,хi≥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