В моей программе PhD по вычислительной технике мы работаем почти исключительно на C ++ и Fortran. Кажется, некоторые профессора предпочитают одного другому. Мне интересно, какой из них «лучше» или один лучше другого в определенных...
Является ли хорошей идеей использовать vector<vector<double>>(используя std) для формирования матричного класса для высокопроизводительного научного вычислительного кода? Если ответ нет. Зачем?...
Проблема В настоящее время я работаю над симуляцией Навье-Стокса с помощью конечных элементов, и я хотел бы изучить влияние различных параметров. Некоторые параметры указываются во входном файле или через параметры командной строки; другие параметры предоставляются в виде флагов в Makefile, поэтому...
Хотите улучшить этот пост? Предоставьте подробные ответы на этот вопрос, включая цитаты и объяснение того, почему ваш ответ правильный. Ответы без достаточной детализации могут быть отредактированы или удалены. Версия стандарта MPI 3.0 формально удалила интерфейс C ++ (ранее он был устаревшим)....
Предположим, что я работаю над научным кодом на C ++. В недавнем обсуждении с коллегой было высказано мнение, что шаблоны выражений могут быть очень плохой вещью, потенциально делая программное обеспечение компилируемым только на определенных версиях gcc. Предположительно, эта проблема затронула...
Я ищу книги или статьи, или посты в блогах, или любые опубликованные материалы в целом, которые конкретно касаются использования современных функций C ++ (семантика перемещения, STL, итераторы, ленивая оценка и т. Д.) В научных вычислениях. Можете ли вы предложить любой? Я думаю, что эти новые...
Собираюсь преподавать студентам бакалавриата курс под названием «Введение в компьютерное программирование». Я немного смущен. В вычислительной физике ученые используют C / C ++ или Python или Fortran, CUDA и т. Д. ..... пришло время создать свою базу. Что я должен использовать? Я знаю, что вы...
Предположим, что у меня есть функция, которая принимает в качестве входных данных несколько значений с плавающей точкой (одинарные или двойные), выполняет некоторые вычисления и выдает выходные значения с плавающей точкой (также одинарные или двойные). Я работаю в основном с MSVC 2008, но также...
Я ищу тензорную библиотеку C ++, которая поддерживает независимый от измерения код. В частности, мне нужно выполнять операции по каждому измерению (до 3), например, вычислять взвешенную сумму. Размеры - это параметр шаблона (и, следовательно, константа времени компиляции). Другое ограничение...
Я хочу разработать параллельное программное обеспечение для научных вычислений с нуля. Я хочу несколько мыслей, на каком языке начать. Программа включает в себя чтение / запись данных в текстовые файлы и параллельное выполнение тяжелых вычислений со многими факторизациями LU и использованием...
Язык C ++ обеспечивает общее программирование и метапрограммирование с помощью шаблонов. Эти методы нашли свое применение во многих крупномасштабных научных вычислительных пакетах (например, MPQC , LAMMPS , CGAL , Trilinos ). Но что они на самом деле внесли в научное вычисление в значении, которое...
Я разрабатываю код для моделирования потока жидкости с биологическими веществами, присутствующими в потоке. Это включает в себя стандартные уравнения Навье-Стокса, связанные с некоторыми дополнительными биологическими моделями. Есть много параметров / констант. Я написал функции для обработки...
Больше всего в MPI мне не нравится иметь дело с типами данных (то есть картами данных / масками), потому что они не очень хорошо подходят для объектно-ориентированного C ++. boost::mpiоднако поддерживает только MPI 1.1 со своего веб-сайта: boost :: mpi - это дружественный к C ++ интерфейс...
Многие из наиболее известных библиотек C ++ в вычислительной науке, такие как Eigen , Trilinos и deal.II, используют стандартный объект библиотеки заголовков шаблонов C ++ std::complex<>для представления сложных чисел с плавающей запятой. В ответе Джека Полсона на вопрос о конструкторах по...
Функция имеет особенность вблизи . Эту особенность можно снять, однако: для следует иметь , поскольку И, следовательно, Однако форма не только не определена в , он также численно нестабилен в окрестности этой точки; чтобы оценить для очень малых численно, можно использовать разложение Тейлора, то...
Я решаю физическую проблему, используя неявную числовую схему. Это приводит меня к решению линейного уравнения с трехдиагональной матрицей. Я закодировал этот алгоритм из Википедии. Интересно, есть ли эффективная библиотека, которая позволяет оптимизировать этот тип уравнения? Важным примечанием...
Я портирую существующий код из MATLAB на C ++ и имею линейную систему для решения (вместо более типичной формы )A x = bх = ЬxA=bxA=bA x = bAx=bAx=b Матрица плотная и общего вида, но не больше 1000x1000. Таким образом, в MATLAB решение находится по функции или через косую чертуAAAmrdivide(b,A)x =...
У меня есть программа Mathematica, которая выполняет некоторые интегралы в 3 или 4 измерениях, используя QuasiMonteCarloметод. Проблема в том, что запуск занимает очень много времени, и некоторые из этих вычислений не могут быть завершены в течение максимального рабочего времени, доступного в нашем...
В статье «Пересмотренные шаблоны выражений: анализ эффективности современных методологий» в SIAM Journal of Scientific Computing упоминается библиотека линейной алгебры «Blaze». Я не слышал об этом раньше, и не могу найти ссылки в Интернете. (Очевидный поиск в Google возвращает вышеуказанную...
Как я могу заменить метод Эйлера на 4-й порядок Рунге-Кутты, чтобы определить движение свободного падения с не постоянной гравитационной величиной (например, свободное падение с высоты 10 000 км над землей)? До сих пор я писал простую интеграцию по методу Эйлера: while() { v += getMagnitude(x) *...