Обчислювальна наука

Припустимо, у мене є матричне рівняння, рекурсивно визначене як A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Тоді рівняння для A [1] виглядає подібним до тривалої дроби, для якої існують деякі високоефективні методи, які уникають втомного перерахунку (див. "Числові рецепти" для деяких прикладів). Однак мені цікаво, чи існують аналогічні методи, які …

18 algorithms 

Я щойно встановив на робочий стіл графічну карту Nvidia GT660 і після певної боротьби мені вдається з'єднати її з R. Я грав з декількома пакетами R, які використовують графічні процесори, особливо gputools, і я порівнював час, витрачений моїм графічним процесором та процесором на виконання деяких основних операцій: інвертування матриць (процесор …

18 r  gpu 

Вам вдалося нарешті з'ясувати, як атоми просторово розташовані на нововиявленому молекулярному утворенні. Через, скажімо, спектроскопічні засоби, ви зараз маєте купу координат атомів, типів атомів, довжини зв’язків, типів зв’язків і чогось іншого для вашої молекули. Тепер ви зацікавлені у визначенні групи точок (групи симетрії) вашої молекули. Для простих молекул, таких як …

18 computational-chemistry 

Мені хотілося б знати, чи існує швидкий спосіб обчислити евклідову відстань двох векторів в Октаві. Здається, що для цього немає спеціальної функції, тому я повинен просто використовувати формулу 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 

Я розробляю код для моделювання потоку рідини з біологічними речовинами, присутніми в потоці. Це стосується стандартних рівнянь Нав'є-Стокса, пов'язаних з деякими додатковими біологічними моделями. Існує багато параметрів / констант. У мене є написані функції для обробки основних обчислень, але проблема, з якою у мене виникає, полягає у великій кількості констант …

17 c++ 

Я студент машинобудування, зацікавлений у галузі аерокосмічної інженерії, де, як мені кажуть, Fortran все ще часто використовується. Яку версію Fortran я повинен інвестувати, щоб навчитися?

17 fortran 

З огляду на дві різні реалізації BLAS, чи можемо ми очікувати, що вони будуть робити такі самі обчислення з плаваючою комою та повертати однакові результати? Або може статися, наприклад, що один обчислює скалярний продукт як та один як таким чином, можливо, даючи різний результат у плаваючій точці IEEE арифметичні?( ( …

17 floating-point  blas 

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

17 linear-algebra  sparse-matrix  oop 

Існує два загальні підходи до представлення рішень у розривному методі Галеркіна: вузловий та модальний. Модальний : Рішення представлені сумами модальних коефіцієнтів, помножених на множину многочленів, наприклад, де зазвичай є ортогональними поліномами, наприклад Legendre . Однією з переваг цього є те, що ортогональні многочлени генерують матрицю діагональної маси.u ( x , …

17 fluid-dynamics  discontinuous-galerkin 

Мені було просто цікаво, чому методи високого порядку (тобто більше, ніж 4) Runge – Kutta майже ніколи не обговорюються / не застосовуються (принаймні, наскільки мені відомо). Я розумію, що для цього потрібен більший обчислювальний час за крок (наприклад, RK14 з вбудованим кроком 12-го порядку ), але чи є інші недоліки …

17 ode  runge-kutta 

Мені потрібно записати деякі дані з обчислень, які будуть прочитані пізніше в паравізі (файл .vtu або vtk). Якщо мова йде про розмір файлу, я повинен перейти на формат ASCII або Бінарний формат?

17 visualization  paraview  vtk 

Я шукаю програмне забезпечення для створення сітки, яке є безкоштовним та відкритим кодом, забезпечує розумний інтерфейс сценаріїв для специфікації домену, працює для складної геометрії, може генерувати 2D та 3D-сітки, Які у мене варіанти?

17 mesh-generation 

Я запустив GROMACS з молекулярною динамікою (MD) на кластері Ubuntu Linux, що складається з вузлів, що містять 24 процесори Intel Xeon. Моя особлива цікава точка виявляється дещо чутливою до арифметичної точності з плаваючою точкою, тому мені довелося запускати GROMACS в подвійній точності, а не одинарній точності - незважаючи на більш …

17 hpc  floating-point  molecular-dynamics  gpu 

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

17 python  fortran 

Мене цікавлять випадки, коли градієнт кон'югату працює набагато краще, ніж метод GMRES. Взагалі, CG є кращим вибором у багатьох випадках SPD (симетрично-позитивний-визначений), оскільки він вимагає меншої кількості зберігання, а теоретична пов'язана зі швидкістю конвергенції для CG подвійна, ніж GMRES. Чи є якісь проблеми, коли фактично спостерігаються такі ставки? Чи є …

17 linear-solver  conjugate-gradient  gmres