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

Мені хотілося б знати, як зазвичай застосовуються умови Діріхле при використанні методу кінцевих обсягів на нерівномірній сітці, орієнтованій на комірки, Моя поточна реалізація просто нав'язує граничну умову мого фіксації значення першої комірки, ϕ1= gD( хL)ϕ1=гD(хL) \phi_1 = g_D(x_L) де є змінним розчином і г Д ( х L ) є …

10 boundary-conditions  finite-volume  poisson 

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

10 parallel-computing  time-integration 

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

10 linear-algebra  python  reference-request  numpy  svd 

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

10 parallel-computing  implicit-methods  cuda 

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

10 finite-element  implicit-methods  explicit-methods 

Я новачок у наукових обчисленнях. Я шукаю Фортран (бажано f90) впровадження указу. Моя проблема потребує Октре, який розділяє мій домен, доки в будь-якому полі не буде більше ніж декілька N частинок (або джерел, де я знаю значення щільності, яке можна включити в метод інтегрального рівняння для вирішення рівняння, як рівняння …

10 fortran  integral-equations 

Я намагаюся розв’язати 2D рівняння Пуассона за допомогою методу Переривчастий Галеркін (DG) та наступної дискретизації (у мене є png-файл, але мені не дозволяється завантажувати його, вибачте): Рівняння: ∇ ⋅ ( κ ∇ T) + f= 0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Нові рівняння: q= κ ∇ T∇ …

10 pde  finite-element  fenics 

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

10 linear-algebra  pde  finite-difference  crank-nicolson 

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

10 finite-element  performance  c++ 

У методі виготовлених розчинів (MMS) один постулює точне рішення, підставляє його в рівняннях і обчислює відповідний вихідний термін. Потім рішення використовується для перевірки коду. Для нестислимих рівнянь Нав'є-Стокса, MMS легко призводить до (ненульового) джерела в рівнянні безперервності. Але не всі коди дозволяють вихідні умови в рівняннях безперервності, тому для цих …

10 navier-stokes  incompressible  verification 

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

10 fluid-dynamics 

у matlab і linsolve, і mldivide використовуються для вирішення системи лінійних рівнянь у всіх визначених, переопределених і недоозначених випадках. Читаючи їх документи, мені було цікаво, які відмінності між ними? Чи використовують вони майже однакові алгоритми матричної факторизації та тріангуляризації у трьох випадках? Якщо A має властивості в opts, linsolve швидше, …

10 matlab  linear-solver 

Зачинено. Це питання поза темою . Наразі відповіді не приймаються. Хочете вдосконалити це питання? Оновіть питання, щоб воно було тематичним для обміну стеками обчислювальної науки. Закрито 2 роки тому . Як говорить заголовок, я хочу мати змогу зупинити код на певному місці та безпечно зупинити його. Я не можу знайти …

10 matlab 

Всі ми знайомі з багатьма обчислювальними методами для вирішення стандартної лінійної системи A x = b .Ах=б. Ax=b. Однак мені цікаво, чи існують якісь "стандартні" обчислювальні методи для вирішення більш загальної (кінцевомірної) лінійної системи форми L A = B ,LА=Б, LA=B, де, скажімо,ААA -матрицям1× n1м1×н1m_1\times n_1 ,ББB -матрицям2× n2м2×н2m_2\times n_2 …

10 linear-algebra 

Зараз я працюю над вирішенням дуже великих симетричних (але не позитивних певних) систем, породжених деякими певними алгоритмами. Ці матриці мають приємну розрідженість блоку, яку можна використовувати для паралельного вирішення. Але я не можу вирішити, чи слід використовувати прямий підхід (наприклад, мультифронтальний) або ітеративний (попередньо обумовлений GMRES або MINRES). Всі мої …

10 parallel-computing  linear-solver  sparse-matrix  linear-system  gmres