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

Алгоритмы Parareal, пита, и PFASST все поголовное-области техники для распараллеливания решения нестационарных задач во времени. Каковы руководящие принципы, лежащие в основе этих методов? Каковы основные различия между ними? Могу ли я сказать, что одно основано на другом? Как? Как насчет их приложений? Я знаю, что не будет ответа на вопрос …

10 parallel-computing  time-integration 

Мой вопрос касается извлечения наблюдаемых из методов QMC, как описано в этой ссылке . Я понимаю формальное происхождение различных методов QMC, таких как Path Integral Monte Carlo. Однако в конце дня я все еще не понимаю, как эффективно использовать эти методы. Основная идея получения квантовых методов МК заключается в дискретизации …

10 monte-carlo  quantum-mechanics 

В чем разница между явным FEM и неявным FEM? Согласно сообщению здесь , кажется, что единственная разница в том, используется ли неявная или явная интеграция времени. Как я помню из одной книги, которую я прочитал, неявная FEM - это то место, где масса не сосредоточена на узлах. Каковы точные определения …

10 finite-element  implicit-methods  explicit-methods 

Для редукции модели я хочу вычислить левые сингулярные векторы, связанные с, скажем, 20-ю самыми большими сингулярными значениями матрицы , где N ≈ 10 6 и k ≈ 10 3 . К сожалению, моя матрица А будет плотной без какой-либо структуры.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx 10^3AAA Если я просто вызываю svdподпрограмму …

10 linear-algebra  python  reference-request  numpy  svd 

Я ищу портировать некоторый код, который разрешает систему уравнений в частных производных (PDE) методом конечных объемов в форме IMPLICIT (для дискретизации по времени). В результате существует трехдиагональная система уравнений в направлениях x, y, z, которая обрабатывается схемой ADI / TDMA. Я не могу найти что-то относительно неявного решения PDE с …

10 parallel-computing  implicit-methods  cuda 

Я новичок в научных вычислениях. Я ищу Fortran (предпочтительно f90) реализацию Octree. Моя проблема требует Octree, которая делит мою область, пока в любом блоке не будет больше, чем несколько N частиц (или источников, где я знаю значение плотности, которое можно включить в метод интегрального уравнения для решения уравнения, подобного уравнению …

10 fortran  integral-equations 

Я использую конечно-разностную схему Кранка-Николсона для решения одномерного уравнения теплопроводности. Мне интересно, если принцип максимума / минимума уравнения теплопроводности (то есть, что максимум / минимум возникает при начальном условии или на границах) также имеет место для дискретного решения. Вероятно, это связано с тем, что Крэнк-Николсон является стабильной и сходящейся схемой. …

10 linear-algebra  pde  finite-difference  crank-nicolson 

Я пытаюсь решить двумерное уравнение Пуассона, используя метод прерывистого Галеркина (DG) и следующую дискретизацию (у меня есть файл png, но мне не разрешено загружать его, извините): Уравнение: ∇ ⋅ ( х ∇ Т) + f= 0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Новые уравнения: Q= κ ∇ T∇ …

10 pde  finite-element  fenics 

В методе готовых решений (MMS) постулируется точное решение, подставляется его в уравнения и вычисляется соответствующий член источника. Затем решение используется для проверки кода. Для несжимаемых уравнений Навье-Стокса MMS легко приводит к (ненулевому) исходному члену в уравнении непрерывности. Но не все коды допускают исходные члены в уравнениях неразрывности, поэтому для этих …

10 navier-stokes  incompressible  verification 

Мы все знакомы со многими вычислительными методами для решения стандартной линейной системы Однако мне любопытно, существуют ли какие-либо «стандартные» вычислительные методы для решения более общей (конечномерной) линейной системы видаAx=b.Ax=b. Ax=b. где, скажем,являются м 1 × п 1 матрица B является м 2 × п 2 матрицы, а L представляет собой …

10 linear-algebra 

Поскольку все знают, что несжимаемый поток не существует в действительности, это предположение введено для упрощения управляющих уравнений. Мы не можем применить это предположение прямо. Обычно число Маха (M <0,3 для несжимаемого потока), изменение плотности (нулевое изменение плотности) и расхождение скорости (равно нулю для несжимаемого потока) являются общими критериями для определения …

10 fluid-dynamics 

C ++ 11 вводит семантику перемещения, которая может, например, улучшить производительность кода в ситуациях, когда C ++ 03 потребуется выполнить конструкцию копирования или назначение копирования. В этой статье сообщается, что при компиляции с C + 11 следующий код ускоряется в 5 раз: vector<vector<int> > V; for(int k = 0; k …

10 finite-element  performance  c++ 

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он подходил для обмена стеками вычислительной науки. Закрыто 2 года назад . Как говорится в заголовке, я хочу иметь возможность остановить код в определенном месте и безопасно остановить его. …

10 matlab 

в Matlab и Linsolve, и Mldivide используются для решения системы линейных уравнений во всех детерминированных, переопределенных и недоопределенных случаях. Читая их документы, мне было интересно, какие различия между ними? Используют ли они почти одинаковые алгоритмы матричной факторизации и триангуляции в трех случаях? Если A имеет свойства в опциях, linsolve работает …

10 matlab  linear-solver 

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

10 fluid-dynamics  simulation