Вычислительная наука

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

19 software  publications  documentation  reproducibility 

В недавнем ответе упоминалось использование генераторов случайных чисел Фортуны или Мерсенна Твистера ( RNG ) для создания симуляции Монте-Карло . Я не слышал о Фортуне раньше, поэтому я посмотрел его - похоже, он в основном предназначен для криптографического использования. В настоящее время я использую Mersenne Twister в производственном коде для …

19 algorithms  random-number-generation 

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

19 algorithms  software  fortran 

- это n × n- симметричная разреженная матрица с положительным определением (SPD). G - разреженная диагональная матрица. n большое ( n > 10000) и число ненулевых в G обычно составляет 100 ~ 1000.AAAn × nn×nn \times nграммGGNnnNnnграммGG Были разложено в Холецкой формевиде L D L T .AAAL D LTLDLTLDL^T Как …

19 linear-algebra 

Для решения больших линейных систем с использованием итерационных методов часто представляет интерес введение предобусловливания, например, вместо решения вместо решения M - 1 ( A x = b ) , где M здесь используется для предобусловливания системы влево. Как правило, мы должны иметь это M - 1 ≈ A - 1 …

19 linear-algebra  iterative-method  preconditioning 

Читая некоторую литературу по решениям для PDE, я натолкнулся на термин « псевдо-время» . Похоже, это общий термин, однако мне не удалось найти хорошее определение или вводную статью для него. Следовательно: что такое псевдо-временное изменение и как оно обычно используется?

18 pde  time-integration 

В большинстве (всех?) Реализаций быстрого мультипольного метода (FMM) октоды используются для декомпозиции соответствующей области. Теоретически, октреи предоставляют простую объемную границу, которая полезна для доказательства O (n) времени выполнения FMM. Помимо этого теоретического обоснования, есть ли преимущества использования Octree по сравнению с другими структурами дерева или дерева данных? Определение списка взаимодействия …

18 algorithms 

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

18 optimization  cuda 

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

18 pde  wavelet 

Я хотел бы использовать 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 

Я пытаюсь реализовать следующую функцию с плавающей запятой двойной точности с низкой относительной ошибкой: l o g s u m (x, y) = журнал( опыт( х ) + опыт( у) )LограммsUм(Икс,Y)знак равножурнал⁡(ехр⁡(Икс)+ехр⁡(Y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Это широко используется в статистических приложениях для добавления вероятностей или плотностей вероятностей, представленных в …

18 floating-point  stability  numerics 

Хотите улучшить этот пост? Предоставьте подробные ответы на этот вопрос, включая цитаты и объяснение того, почему ваш ответ правильный. Ответы без достаточной детализации могут быть отредактированы или удалены. Адаптивное уточнение сетки (AMR) является распространенным методом решения проблемы широко варьирующихся пространственных масштабов в численном решении PDE. Какие библиотеки общего назначения существуют …

18 pde  libraries  parallel-computing  adaptive-mesh-refinement 

Каковы лучшие (масштабируемость и эффективность) алгоритмы для генерации неструктурированных четырехмерных сеток в 2D? Где я могу найти хороший неструктурированный четырехъядерный генератор сетки? (предпочтительно с открытым исходным кодом)

18 mesh-generation  computational-geometry 

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

18 python  hpc  c++  fortran  simulation 

Физические величины, такие как давление, плотность, энергия, температура и концентрация, всегда должны быть положительными, но численные методы иногда вычисляют отрицательные значения во время процесса решения. Это не хорошо, потому что уравнения будут вычислять сложные или бесконечные значения (как правило, сбой кода). Какие численные методы могут быть использованы для гарантии того, …

18 pde  fluid-dynamics  hyperbolic-pde