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

Переходячи від сильної форми PDE до форми FEM, здається, завжди слід робити це, спочатку зазначивши варіаційну форму. Для цього ви помножите сильну форму на елемент у якомусь (Соболєві) просторі та інтегруєтесь у своєму регіоні. Це я можу прийняти. Що я не розумію, це те, чому також доводиться використовувати формулу Гріна …

24 pde  finite-element  elliptic-pde 

Опис експерименту: В інтерполяції Лагранжа точне рівняння відбирається в NNN точках (поліноміальний порядок N−1N−1N - 1 ) і інтерполюється в 101 бал. Тут NNN змінюється від 2 до 64. Щоразу, коли підготуються графіки помилок L1L1L_1 , L2L2L_2 і L∞L∞L_\infty . Видно , що, коли функція дискретизируются на равноразнесенние точках, помилка …

24 finite-difference  interpolation  error-estimation  polynomials 

Дуже поширена проблема ланцюга Маркова в Монте-Карло включає обчислення ймовірностей, які є сумою великих експоненціальних доданків, еа1+ еа2+ . . .еа1+еа2+... e^{a_1} + e^{a_2} + ... ааaК: = максi( ai)К: =максi(аi)K := \max_{i}(a_{i}) а'= К+ л о г( еа1- К+ еа2- К+ . . . )а'=К+лог(еа1-К+еа2-К+...)a' =K + log\left( e^{a_1 …

24 monte-carlo  floating-point  statistics 

Наскільки я можу сказати, дві великі загальні рамки програмного забезпечення для обчислювальної науки Міністерства енергетики США - це PETSc і Trilinos . На перший погляд вони схожі за межами відмінностей у мові (C проти C ++). Які основні відмінності між двома рамками, і які фактори повинні впливати на вибір одного …

24 petsc  software  trilinos 

Зараз я не виконую жодного паралельного коду, але передбачаю, що в майбутньому запускається паралельний код, використовуючи гібрид OpenMP та MPI. Налагоджувачі були для мене безцінним інструментом при запуску серійних проектів. Чи може хтось порекомендувати паралельний налагоджувач (або декілька налагоджувачів) використовувати для налагодження паралельного програмного забезпечення? Вільне програмне забезпечення було б …

24 software  parallel-computing 

Як я розумію, існує дві основні категорії ітеративних методів розв’язання лінійних систем рівнянь: Стаціонарні методи (Jacobi, Gauss-Seidel, SOR, Multigrid) Методи підпростори Крилова (кон'югатний градієнт, GMRES тощо) Я розумію, що більшість стаціонарних методів працюють за допомогою ітераційного розслаблення (згладжування) режимів Фур'є помилки. Як я розумію, метод кон'югатного градієнта (метод підпростору Крилова) …

24 linear-algebra  convergence  krylov-method  conjugate-gradient 

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

23 software  publications  documentation  reproducibility 

Що є найсучаснішим у наближенні сильно коливальних інтегралів як в одному вимірі, так і у вищих розмірах до довільної точності?

23 quadrature  precision  numerical-analysis 

Обчислювальний кластер Amazon EC2 коштує приблизно від 800 до 1000 доларів (залежно від робочого циклу) на фізичне ядро ​​процесора протягом 3 років. Під час нашого останнього раунду з придбання обладнання моя лабораторія зібрала обладнання з 48 ядер, що дуже схоже на кластери Amazon, приблизно за $ 300 доларів. Я чогось …

23 cloud-computing  hpc 

(xi,yi)(xi,yi)(x_i,y_i)y(x)y(x)y(x)xxxf(x)≡y(x)−(ax+b)f(x)≡y(x)−(ax+b)f(x) \equiv y(x) - (ax + b)x→∞x→∞x \to \inftyf′(x)f′(x)f'(x)f′′(x)f″(x)f''(x)і т. д. Але я не знаю, що таке функціональна форма для , якщо вона навіть має таку, яку можна описати елементарними функціями.f(x)f(x)f(x) Моя мета - отримати найкращу можливу оцінку асимптотичного схилу . Очевидним неочищеним методом є виділення останніх кількох точок даних …

23 algorithms 

Багато чисельних підходів до CFD можуть бути розширені до довільно високого порядку (наприклад, розривні методи Галеркіна, методи ВЕНО, спектральна диференціація тощо). Як слід вибрати відповідний порядок точності для даної проблеми?

23 fluid-dynamics  accuracy 

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

23 software  matlab  languages  octave 

Збирається викладати студентів бакалаврату курсу під назвою Вступ до комп'ютерного програмування. Я трохи розгубився. У галузі обчислювальної фізики вчені використовують C / C ++ або Python або Fortran, CUDA тощо .... Це час для створення їх бази. Що я повинен використовувати? Я знаю, що ви можете вивчити нову мову програмування …

22 python  c++  computational-physics  languages 

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

22 languages 

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

22 petsc  pde  implicit-methods