Меня интересуют случаи, когда градиент сопряжения работает намного лучше, чем метод GMRES. Вообще, CG является предпочтительным выбором во многих случаях SPD (симметрично-положительно определенного), потому что он требует меньшего объема памяти и теоретическая оценка скорости сходимости для CG...
Параллельное сокращение предполагает, что соответствующая операция является ассоциативной. Это предположение нарушается при добавлении чисел с плавающей запятой. Вы можете спросить, почему я забочусь об этом. Что ж, это делает результаты менее воспроизводимыми. И это ухудшается, когда моделируемый...
Разработка научных алгоритмов - это очень итеративный процесс, часто включающий в себя изменение множества параметров, которые я хочу изменить либо как часть моего экспериментального проекта, либо как часть производительности алгоритма настройки. Какие стратегии я могу использовать для...
Я студент машиностроения, интересующийся аэрокосмической техникой, где, как мне сказали, Фортран все еще широко используется. Какую версию Fortran я должен потратить на изучение своего...
Большая часть моей собственной работы вращается вокруг того, чтобы алгоритмы масштабировались лучше, и один из предпочтительных способов показать параллельное масштабирование и / или параллельную эффективность - это построить производительность алгоритма / кода по количеству ядер, например где ось...
На днях мой инструктор по вычислительной гидродинамике отсутствовал, и он послал своего кандидата в доктора наук, чтобы заменить его. В своей лекции он, казалось, указал на несколько недостатков, связанных с различными схемами дискретизации для моделирования потока жидкости: Метод конечных...
У меня есть ограниченная невыпуклая 2-D функция, минимум которой я бы хотел найти. Функция довольно плавная. Оценка это дорого. Допустимая ошибка составляет около 3% от области функции в каждой оси. Я попытался запустить реализацию алгоритма DIRECT в библиотеке NLOPT, но он не дал значительного...
Моя диссертация посвящена разработке численных методов для моделирования снижения горения. Я использую свои методы исключительно в химической модели моделирования горения, и у меня есть много примеров для моделирования 0-D (без потока). Я хотел бы запустить симуляции, в которых есть потоки,...
Я планирую использовать multigrid для вычисления некоторых собственных значений и векторов, и я заметил, что PETSc имеет высокоуровневую поддержку multigrid. Документация PETSc гласит, что эту часть PETSc не следует использовать, поскольку она скоро будет заменена. Какие другие библиотеки имеют...
Предположим, у меня есть некоторая функция ееf и я хочу найти ИксИксx такой, что е( х ) ≈ 0е(Икс)≈0f(x)\approx 0 . Я мог бы использовать метод Ньютона-Рафсона. Но для этого необходимо, чтобы я знал производную функцию е'( х )е'(Икс)f'(x) . Аналитическое выражение для ееf может быть недоступно....
Возвращаясь к заре OS X, казалось, было много шума, по крайней мере, в мире Mac (я в то время не был близок к научным вычислениям) о Mac OS как платформе для научных вычислений и приложений HPC. XGrid вышла из коробки, у Virginia Tech был свой модный компьютерный кластер на базе Mac, Стэнфорд делал...
Какой из них лучше: Фортран или Питон? И я думаю, что в обоих случаях вам нужен Gnuplot, я прав? Я сейчас работаю на машине с Windows. Я хотел бы использовать его для получения численных решений физических задач, включая моделирование по методу Монте-Карло, численное интегрирование и...
Язык C ++ обеспечивает общее программирование и метапрограммирование с помощью шаблонов. Эти методы нашли свое применение во многих крупномасштабных научных вычислительных пакетах (например, MPQC , LAMMPS , CGAL , Trilinos ). Но что они на самом деле внесли в научное вычисление в значении, которое...
Я знаю, что метод Ницше является очень привлекательным, поскольку он позволяет учитывать граничные условия типа Дирихле или слабо связываться с граничными условиями трения без использования множителей Лагранжа. И его преимущество, заключающееся в преобразовании граничного условия Дирихле в слабые...
Я пишу небольшую библиотеку для разреженных матричных вычислений, чтобы научить себя наилучшим образом использовать объектно-ориентированное программирование. Я очень много работал над созданием хорошей объектной модели, в которой части (разреженные матрицы и графики, описывающие их структуру...
Я ищу программное обеспечение для генерации сетки, которое является бесплатным и открытым исходным кодом, предоставляет разумный интерфейс сценариев для спецификации домена, работает для сложных геометрий, может генерировать 2D и 3D сетки, Какие варианты у меня...
Рассмотрим IVP для системы ODE , . Чаще всего эта проблема считается жесткой, когда матрица Якоби имеет как собственные значения с очень большой отрицательной вещественной частью, так и собственные значения с очень маленькой отрицательной вещественной частью (я рассматриваю только стабильную...
В литературе по численным методам многие исследовательские работы состоят из описания нового алгоритмического варианта, за которым следует несколько тестовых задач, сравнивающих новый метод с одним или двумя существующими методами. Это затрудняет определение Как новый метод будет работать на другие...
Мне было интересно, существует ли быстрый и эффективный способ заранее определить количество ненулевых элементов для операции умножения разреженных матриц при условии, что обе матрицы находятся в формате CSC или CSR. Я знаю, что есть один в пакете smmp, но мне нужно что-то, что уже реализовано в C...
Я всегда вычисляю вещи вручную, но теперь мои товарищи становятся противными и делают много повторяющихся упражнений, в том числе просто вставляя вещи, как показано выше. Я особенно заинтересован в программном обеспечении с открытым исходным кодом, таком как Python или R, чтобы упростить эти виды...