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

У меня есть PETSc, Matи я хочу оценить его

В сети обсуждается множество схем ФД для уравнения переноса . Например, здесь: http://farside.ph.utexas.edu/teaching/329/lectures/node89.html∂T∂T+ U ∂T∂Икс= 0∂T∂T+U∂T∂Иксзнак равно0\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=0 Но я не видел, чтобы кто-нибудь предлагал "неявную"...

Вычислительные люди науки: Первоначально я разместил этот вопрос на Math Stack Exchange, и кто-то сказал, что я мог бы получить «гораздо лучшие» ответы здесь: Я новичок в численных методах и Matlab. Я пытаюсь оценить следующую сумму двух тройных интегралов (она, очевидно, может быть записана более...

Какие подходы используются на практике для оценки числа условий больших разреженных

Я хотел бы вывести системы решения ODE на графические процессоры в «тривиально распараллеливаемой» обстановке. Например, анализ чувствительности с 512 различными наборами параметров. В идеале я хочу решать ODE с помощью интеллектуального адаптивного временного шага, такого как CVODE, а не с...

Я взял базовое введение в метод конечных элементов, в котором не подчеркивалось глубокое понимание «слабой формулировки». Я понимаю, что с помощью метода Галеркина мы умножаем обе стороны (эллиптического) PDE на тестовую функцию и затем интегрируем (по частям или по теореме дивергенции). Иногда мне...

Я написал программу / библиотеку, которую использовал для получения результатов в статье. (Вот и все , но мой вопрос общий.) У меня есть тесты, которые я регулярно использую ctest(для запуска требуется несколько минут). Чтобы воспроизвести некоторые таблицы или рисунки в статье, мне нужно создать...

Я заинтересован в решении уравнения Пуассона с использованием метода конечных разностей. Я хотел бы лучше понять, как написать матричное уравнение с граничными условиями Неймана. Будет ли кто-то пересмотреть следующее, это правильно? Конечно-разностная матрица Уравнение Пуассона, ∂2ты ( х )∂Икс2=...

Мне нужно вычислить следующий интеграл: Где - матрица (кинетическая и потенциальная энергия одной частицы выражена в базисе), - матрица, зависящая от (многочастичная одночастичная тело функции Грина), а контурный интеграл представляет собой левый полукруг. Подынтегральная функция имеет полюсы на...

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

У меня есть набор данных и я хочу найти параметр m такой, чтобы он минимизировал сумму k ∑ i = 1 | м - х я | , то естьx1,x2,…,xkx1,x2,…,xkx_{1}, x_{2}, \ldots, x_{k}mmm∑i=1k∣∣m−xi∣∣.∑i=1k|m−xi|.\sum_{i=1}^{k}\big|m-x_i\big|....

Насколько мне известно, многосеточные решатели используют итеративные сглаживатели, такие как Якоби, Гаусс-Зайдель и SOR, чтобы смягчить ошибку на различных частотах. Можно ли использовать метод подпространств Крылова (например, сопряженный градиент, GMRES и т. Д.)? Я не думаю, что они...

Я пытаюсь решить двумерное уравнение Пуассона с помощью конечных разностей. В процессе, я получаю разреженную матрицу только с переменными в каждом уравнении. Например, если переменные были U , то дискретизация даст:555UUU...

Я слышал анекдотично, что когда кто-то пытается численно сделать интеграл вида ∫∞0е( х ) J0(х )д х∫0∞е(Икс)J0(Икс)dИкс\int_0^\infty f(x) J_0(x)\,\mathrm{d}x с сглаженным и хорошо себя ведет (например, не сам колеблющийся, неособой и т. д.), то это поможет точность переписать его каке( х )е(Икс)f(x)...

Существует система линейных ограничений . Я хочу найти строго положительный вектор который удовлетворяет этим ограничениям. Это означает, что требуется для каждого компонента из . Как я могу использовать LP-решатель, чтобы найти такой строго положительный вектор (или подтвердить, что существует)? Я...

Я создаю простой симулятор астрономии, который должен использовать ньютоновскую физику для симуляции движения планет в системе (или любых других объектов). Все тела представляют собой круги в евклидовой плоскости, которые обладают такими свойствами, как положение, скорость, масса, радиус и...

Извините за длинный пост, но я хотел включить все, что я думал, было актуально с первого взгляда. Что я хочу Я реализую параллельную версию методов подпространств Крылова для плотных матриц. В основном GMRES, QMR и CG. Я понял (после профилирования), что моя рутина DGEMV была жалкой. Поэтому я...

Мне нужно простое объяснение многосеточного метода или некоторая литература по этому поводу. Я знаком с итерационными методами, включая BiCGStab, CG, GS, Jacobi и предварительные условия, но я новичок в многосеточном методе. Может кто-нибудь объяснить это подробно или хотя бы предоставить явно...

В моделировании полупроводников, как правило, уравнения масштабируются, чтобы они имели нормированные значения. Например, в крайних случаях плотность электронов в полупроводниках может варьироваться более чем на 18 порядков, а электрическое поле может изменяться более чем на 6 (или более) порядков....

Моя исследовательская группа фокусируется на молекулярной динамике, которая, очевидно, может генерировать гигабайты данных как часть одной траектории, которую затем необходимо проанализировать. Некоторые из проблем, с которыми мы сталкиваемся, связаны с корреляциями в наборе данных, что означает,...