Недавно я столкнулся со странной проблемой с FORTRAN95. Я инициализировал переменные X и Y следующим образом: X=1.0 Y=0.1 Позже я складываю их вместе и печатаю результат: 1.10000000149012 После изучения переменных кажется, что 0,1 не представляется в двойной точности с полной точностью. Есть ли...
После некоторых шагов байесовского обновления у меня осталось апостериорное распределение в виде смеси нормальных распределений,Таким образом, параметр \ theta взят из распределения, PDF которого задан как взвешенная смесь нормальных PDF, и не является суммой нормальных RV. Я хотел бы нарисовать...
У меня плотная вещественная симметричная квадратная матрица. Размер составляет около 1000x1000. Мне нужно вычислить первый главный компонент и подумать, каким может быть лучший алгоритм для этого. Похоже, что MATLAB использует алгоритмы Арнольди / Ланцоша (для eigs). Но читая о них, я не уверен,...
Я решаю множество PDE численно, но прикладная математика не является моей областью. Я не понял, какие журналы по прикладной математике я должен читать, чтобы не отставать от последних событий в этой области. Какие хорошие журналы нужно почитать, чтобы не отставать от последних достижений в...
У меня есть ODE: U'= - 1000 U + ев я п ( т )U'знак равно-1000U+sяN(T)u'=-1000u+sin(t) и ( 0 ) = - 11000001U(0)знак равно-11000001u(0)=-\frac{1}{1000001} Я знаю, что этот конкретный ODE жесток, аналитически. Я также знаю, что если мы используем явный (прямой) шаговый метод времени (Эйлер,...
Я решаю проблему мультимасштаба, используя Гетерогенный Метод Мультимасштаба (HMM) . По сути, моя конкретная процедура использует следующий итерационный процесс: Решите много локальных матричных систем. Вычислите значение интереса из решений локальных систем. Соберите глобальную матричную систему...
В методе самосогласованного поля Хартри-Фока для решения не зависящего от времени электронного уравнения Шредингера мы стремимся минимизировать энергию основного состояния, , системы электронов во внешнем поле относительно выбора спина орбитали, .E0E0E_{0}{χi}{χi}\{\chi_{i}\} Мы делаем это путем...
Недавно я задал вопрос в том же духе для косоэрмитовых матриц. Вдохновленный успехом этого вопроса, и после того, как пару часов ударился головой о стену, я смотрю на экспоненциальную матрицу вещественных асимметричных матриц. Путь к поиску собственных значений и собственных векторов кажется...
Я уже разработал рабочее решение метода конечных элементов для решения задач теплопередачи, используя GPU и OpenCL, используя метод сопряженных градиентов. Основным недостатком этого метода является высокий спрос на память. Более того, в случае видеокарт память часто очень ограничена. Я вижу два...
Я начал изучать OpenFOAM из учебника Cavity, который размещен на сайте . При экспериментировании с различными числами Рейнольдса в разделе «2.1.8.2 Выполнение кода» в учебном пособии говорится, что необходимо перезапустить решатель, поскольку «разумно увеличить время решения». Но когда я сделал...
Чем (обобщенное) геометрическое программирование отличается от общего выпуклого программирования? Геометрическая программа может быть преобразована в выпуклую программу и обычно решается методом внутренней точки. Но в чем преимущество прямой формулировки задачи как выпуклой программы и ее решения...
Я хочу знать , какие из классических линейных решателей (например , Гаусс-Зейделя, Jacobi, SOR) гарантированно сходятся для задачи , где положительно полу определена и, конечноA b ∈ i m ( A )A x = bAx=bAx=bAAAb ∈ i m ( A )b∈im(A)b \in im(A) (Примечание является полуопределенным и не...
У меня есть проблема, похожая по формулировке на этот пост, с несколькими заметными отличиями: Какие простые методы существуют для адаптивной выборки 2D-функции? Как в этом посте: У меня есть и оценка этой функции несколько дороже, чтобы вычислитье( х , у)f(x,y)f(x,y) В отличие от этого поста: Меня...
Как человек, которому приходится преподавать курсы по вычислительной науке, я сталкиваюсь с извечным вопросом: как я могу оценить способность студентов изучать предмет, который зависит от приложений, которые трудно проверить с помощью «стандартных» методов тестирования ( письменные или устные...
Я знаю, что мы можем использовать методы математического анализа, чтобы доказать, есть ли у IVP или BVP решение, оно уникально и постоянно зависит от граничных / начальных значений. Для некоторых PDE, особенно нелинейных pde, очень трудно, если не невозможно, доказать правильность. Есть ли...
как следует из названия, я пытаюсь вычислить интеграл от компактно поддерживаемой функции (квинтический полином Вендланда) на треугольнике. Обратите внимание, что центр функции находится где-то в трехмерном пространстве. Я интегрирую эту функцию в произвольный, но маленький треугольник ( ). В...
Какие рекомендуемые способы выполнения нелинейных наименьших квадратов, мин , с коробкой ограничений л O J < = р J < = ч я J ? Мне кажется (вбегают дураки), что можно сделать квадратные ограничения квадратичными и минимизировать ∑ i e r r i ( p ) 2 + C ∗ ∑ j t u b ( p j , l...
Предположим, я знал, что последовательность случайных чисел была сгенерирована линейным конгруэнтным генератором. Это, xn+1=(aXn+c)modmxn+1=(aXn+c)modmx_{n+1}=(aX_n+c) \bmod m Если мне дают весь период (или, по крайней мере, большую смежную подпоследовательность этого), как я могу восстановить...
Пусть - выпуклая полигонально ограниченная липшицева область в , пусть .R 2 f ∈ L 2 ( Ω )ΩΩ\OmegaR2R2\mathbb R^2f∈L2(Ω)f∈L2(Ω)f \in L^2(\Omega) Тогда решение задачи Дирихле в , в имеет единственное решение в и корректно, т. Е. Для некоторой константы имеем .Δu=fΔu=f\Delta u =...
Мне интересно, как найти собственные значения некоторой разреженной матрицы в заданном интервале [a, b] итерационным методом. По моему личному пониманию, более очевидно использовать подпространственный метод Крылова для нахождения крайних собственных значений, а не...