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

40
Какое программное обеспечение и рабочий процесс рекомендуется для публикации научных данных и графиков?

Какое программное обеспечение обеспечивает хороший рабочий процесс от простого построения нескольких точек данных до создания графики уровня публикации с подробными стилями, математической версткой и «профессиональным качеством»? Это немного связано с вопросом Дэвида (« Какие атрибуты делают фигуру...

40
Как связать код с публикациями

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

39
Почему деление намного сложнее, чем другие арифметические операции?

Недавно я столкнулся со случаем, когда мне понадобилось целочисленное деление на чипе, в котором его не было (ARM Cortex-A8). Пытаясь исследовать, почему это должно быть, я обнаружил, что в общем случае деление занимает гораздо больше циклов, чем сложение, вычитание или умножение практически любой...

39
Каково состояние параллельных методов ODE?

В настоящее время я изучаю параллельные методы для интеграции ODE. Существует много новой и старой литературы, описывающей широкий спектр подходов, но я не нашел каких-либо недавних обзоров или обзорных статей, описывающих эту тему в целом. Есть книга Burrage [1], но ей почти 20 лет, и,...

37
Является ли хорошей идеей использовать vector <vector <double >> для формирования матричного класса для высокопроизводительного научного вычислительного кода?

Является ли хорошей идеей использовать vector<vector<double>>(используя std) для формирования матричного класса для высокопроизводительного научного вычислительного кода? Если ответ нет. Зачем?...

36
Места для публикации документов, которые подчеркивают программное обеспечение

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

36
Символьные программные пакеты для выражений Matrix?

Мы знаем, что симметрична и положительно определена. Мы знаем, что B ортогонально:AA\mathbf AВB\mathbf B Вопрос: является ли симметричным и положительно определенным? Ответ: да.B ⋅ A ⋅ B⊤B⋅A⋅B⊤\mathbf B \cdot\mathbf A \cdot\mathbf B^\top Вопрос: Мог ли компьютер сказать нам это? Ответ: возможно....

36
Можно ли сделать карьеру в SciComp без участия в исследованиях в области вооружений?

Я на международной конференции (ICIAM2019) о численных методах и удивлен распространенностью приложений, непосредственно связанных с исследованиями в области вооружений. Примеры: Один из лауреатов премии рассказывает о математической проблеме восстановления / обнаружения движущихся объектов с...

36
Как оператор обратной косой черты в MATLAB решает

Я сравнивал несколько своих кодов с «стандартными» кодами MATLAB. Я удивлен результатами. Я запустил пример кода (разреженная матрица) n = 5000; a = diag(rand(n,1)); b = rand(n,1); disp('For a\b'); tic;a\b;toc; disp('For LU'); tic;LULU;toc; disp('For Conj Grad'); tic;conjgrad(a,b,1e-8);toc;...

36
В какие парадигмы программирования я должен инвестировать, если я хочу, чтобы мой код работал на машинах Petascale в будущем?

Из обзора топ-500 довольно ясно видно, что отрасль имеет тенденцию к экспоненциальному увеличению числа процессорных ядер . Все крупнейшие суперкомпьютеры используют MPI для связи между узлами, хотя не наблюдается явной тенденции к параллелизму на узле, при этом самый простой (но не обязательно...

35
Распараллеливание цикла for в Python

Есть ли в Python инструменты, которые похожи на парфор Matlab? Я нашел эту ветку , но ей четыре года. Я подумал, что, может быть, кто-то здесь может иметь более свежий опыт Вот пример типа вещи, которую я хотел бы распараллелить: X = np.random.normal(size=(10, 3)) F = np.zeros((10, )) for i in...

35
Какие атрибуты делают фигуру «профессиональным качеством»?

Я слышал, как люди говорят, что сюжеты, созданные ORIGIN, выглядят полированными и «профессиональными», тогда как сюжеты, производимые Mathematica, - нет. Тем не менее, большинство программ создания сюжетов вполне настраиваемы, и вполне понятно, что с правильными настройками для таких вещей, как...

35
Математические библиотеки для OpenCL?

Я ищу информацию от всех, кто пытался использовать OpenCL в своем научном коде. Кто-нибудь пробовал (недавно) ViennaCL ? Если да, то как это сравнить с острием ? Как насчет OCLTools ? Это соответствует обещанию? Если это так, будет ли возможным начать писать математические ядра в...

34
CUDA против OpenCL по состоянию на конец 2013 года

Как CUDA и OpenCL сравниваются друг с другом в конце 2013 года с точки зрения программиста? Моя группа думает о том, чтобы попытаться использовать вычисления на GPU. Будем ли мы значительно ограничивать себя, выбирая оборудование, которое поддерживает только OpenCL, но не CUDA? Чтобы быть более...

34
Как мне убедиться, что результаты моих симуляций и результаты в моей работе всегда синхронизированы?

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

33
Странные колебания при решении уравнения переноса методом конечных разностей с полностью замкнутыми граничными условиями Неймана (отражение на границах)

Я пытаюсь решить уравнение адвекции, но странное колебание появляется в решении, когда волна отражается от границ. Если кто-то видел этот артефакт раньше, мне было бы интересно узнать причину и как ее избежать! Это анимированный GIF-файл, открываемый в отдельном окне для просмотра анимации (он...

31
Различия в производительности между ATLAS и MKL?

ATLAS - это бесплатная замена BLAS / LAPACK, которая настраивается на машину при компиляции. MKL - коммерческая библиотека, поставляемая Intel. Сравнимы ли эти две библиотеки с точки зрения производительности или MKL имеет преимущество в некоторых задачах? Если да, то...

31
Современные ресурсы для изучения FEM

Мне нужно начать использовать методы конечных элементов. Я о том, чтобы начать читать численные решения дифференциальных уравнений методом конечных элементов с помощью Клаас Джонсона , но это от 1987. Два вопроса: 1) Какие новые полезные ресурсы / учебники / электронные книги / конспекты лекций на...

31
Как выбрать метод решения линейных уравнений

Насколько мне известно, есть 4 способа решения системы линейных уравнений (поправьте меня, если их больше): Если системная матрица является квадратной матрицей полного ранга, вы можете использовать правило Крамера; Вычислить обратную или псевдообратную матрицу системы; Используйте методы матричного...