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

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

19 software  publications  documentation  reproducibility 

Нещодавна відповідь згадував про використання генераторів Fortuna або Mersenne Twister випадкових чисел ( RNG ) для посіву моделювання в Монте-Карло . Раніше я не чув про Фортуну, тому переглянув її - схоже, вона в основному призначена для криптографічного використання. В даний час я використовую Mersenne Twister у виробничому коді для …

19 algorithms  random-number-generation 

Мені часто доводиться писати дуже схожий код для однієї, двох і тривимірних версій заданої операції / алгоритму. Підтримання всіх цих версій може стати втомливим. Просте генерування коду працює досить добре, але, здається, існує кращий шлях. Чи є відносно простий спосіб написати операцію один раз і чи узагальнити її до вищих …

19 algorithms  software  fortran 

- n × n симетрична позитивна певна (SPD) розріджена матриця. G - розріджена діагональна матриця. n велике ( n > 10000), а кількість ненулів у G зазвичай становить 100 ~ 1000.ААAn × nн×нn \times nГГGннnннnГГG Були розкладено в Холецького формівигляді L D L T .ААALD LТLDLТLDL^T Як ефективно оновлювати і …

19 linear-algebra 

Для розв’язання великих лінійних систем використанням ітеративних методів часто представляє інтерес запровадити попередню умову, наприклад, замість M - 1 ( A x = b ) , де M використовується тут для лівої попередньої кондиціонування системи. Як правило, ми повинні мати, що M - 1 ≈ A - 1, і забезпечити …

19 linear-algebra  iterative-method  preconditioning 

У більшості (всіх?) Реалізацій методу швидкого мультиполю (FMM) використовуються octrees для декомпозиції відповідного домену. Теоретично, октриси забезпечують просту об'ємну зв'язану, що корисно для доведення O (n) часу виконання FMM. Крім цього теоретичного обґрунтування, чи є користь від використання Окрем над іншими структурами даних дерева або триє? Визначити список взаємодії може …

18 algorithms 

Читаючи деяку літературу про вирішувачі PDE, я сьогодні натрапив на термін псевдо-кроковий час . Це здається загальним терміном, проте мені не вдалося знайти хорошого визначення або вступної статті до нього. Звідси: Що таке псевдо-ступінчастий час, і як його зазвичай використовують?

18 pde  time-integration 

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

18 optimization  cuda 

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

18 pde  wavelet 

Я намагаюся реалізувати наступну функцію у плаваючої точки з подвійною точністю з низькою відносною помилкою: logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Це широко використовується в статистичних додатках для додавання ймовірностей або щільності ймовірностей, які представлені в журнальному просторі. Звичайно, або або \ exp (y) можуть легко переповнюватись або перетікати, що було …

18 floating-point  stability  numerics 

Я хотів би використати f2pyсучасний Fortran. Зокрема, я намагаюся отримати наступний базовий приклад для роботи. Це найменший корисний приклад, який я міг би створити. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out) :: z(:) ! Variable Declarations ! real*8, allocatable :: …

18 python  fortran 

Хочете вдосконалити цю посаду? Надайте детальні відповіді на це питання, включаючи цитати та пояснення, чому ваша відповідь правильна. Відповіді без достатньої кількості деталей можуть бути відредаговані або видалені. Адаптивне уточнення сітки (AMR) - поширена методика вирішення проблеми широко варіюючих просторових масштабів у числовому рішенні PDE. Які бібліотеки загального призначення існують …

18 pde  libraries  parallel-computing  adaptive-mesh-refinement 

Які найкращі (масштабованість та ефективність) алгоритми для генерування неструктурованих чотиримісних сіток у 2D? Де я можу знайти хороший неструктурований чотиримісний генератор? (перевагу з відкритим кодом)

18 mesh-generation  computational-geometry 

Я хочу розробити паралельно наукове програмне забезпечення для обчислень. Я хочу подумати, з якої мови почати. Програма передбачає читання / запис даних у файли txt і паралельні великі обчислення, з багатьма LU-факторизаціями та використанням розріджених лінійних розв'язувачів. Кандидатські рішення, про які я думав, - це Fortran 2003/2008 з OpenMP або …

18 python  hpc  c++  fortran  simulation 

Фізичні величини, такі як тиск, щільність, енергія, температура та концентрація, завжди повинні бути позитивними, але числові методи іноді обчислюють негативні значення під час процесу розчинення. Це не нормально, оскільки рівняння будуть обчислювати складні або нескінченні значення (як правило, збій коду). Які числові методи можна використовувати, щоб гарантувати, що ці величини …

18 pde  fluid-dynamics  hyperbolic-pde