Моя система представляет собой симметричную задачу КЭ с множителями Лагранжа (например, несжимаемый поток Стокса): ( АВВTС)(AВTВС)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} где - типичный случай (я даже удостоверился, что уравнения пронумерованы так, чтобы множители Лагранжа появились...
Вот типичное использование указателей на функции в C. Я хотел бы сделать нечто подобное в Fortran. У меня есть некоторые идеи, но я хотел бы знать, есть ли какой-то канонический способ сделать это. Указатели и контексты функций, переданные пользователем, сохраняются, а затем вызываются позже....
Параллельные вычисления часто моделируются с использованием детерминированной локальной скорости вычислений, времени ожидания и пропускной способности сети. На самом деле, они пространственно переменные и недетерминированные. Такие исследования, как Skinner и Kramer (2005), наблюдают...
Насколько медленнее petsc4py против c / c ++ / fortran? Я понимаю, что это будет в значительной степени зависеть от выполняемого кода, но как насчет чего-то простого, такого как матрично-векторный...
Мне дают матрицу Q которая является симметричной, обратимой, положительно определенной и плотной. Мне нужно проверить, если det, где J является матрицей всех единиц.12 × 1212×1212 \times 12QQQdet ( Q ) = det ( 12 I- Q - J)( 1 )йе(Q)знак равнойе(12я-Q-J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ В...
Влияет ли компиляция PETSc с внешней библиотекой BLAS / LAPACK существенным образом на производительность разреженных матриц, или он использует эти библиотеки только для математики с плотной...
Я хотел бы использовать метод 8-го порядка Рунге-Кутты (89) в приложении по небесной механике / астродинамике, написанном на C ++, на машине Windows. Поэтому мне интересно, знает ли кто-нибудь хорошую библиотеку / реализацию, которая документирована и бесплатна для использования? Это нормально,...
Что вычисляется быстрее, или или ? , и - положительные вещественные числа с .abaбa^blogacжурналaс\log_a cc√bсб\sqrt[b]{c}aaabбbcсcb>1б>1b>1 Какие алгоритмы вы будете использовать при сравнении? Каковы их сложности? Например, когда илиc≡abс≡aбc \equiv a^bc≈abс≈aбc \approx a^b Этот вопрос был...
Существует много хорошо известных численных методов решения уравнений типа например, метод деления пополам, метод Ньютона и т. Д.е( х ) = 0 ,x ∈ RN,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, В моем приложении рассчитывается стохастическим методом (результат является средним).е( х )е(Икс)f(x)...
Я ищу заслуживающие доверия ссылки на то, сколько ресурсов суперкомпьютеры тратят на координацию, а не на выполнение реальной работы, связанной с задачами. Ресурсы могут иметь доступную вычислительную мощность, но даже ватты кажутся действительными единицами. Я полагаю, что один из моих профессоров...
Я новичок в пакете PETSc. У меня матрица A размером ~ 4000x4000 в формате матричного рынка, и я хочу, чтобы PETSc решил эту проблему с помощью нескольких процессоров. Я знаю, как решить систему на одном процессоре, но я не знаю, как распределить матрицу и векторы между различными процессорами. Есть...
У меня есть значительная модель (~ 5000 строк), написанная на C. Это последовательная программа, нигде не генерирующая случайные числа. Она использует библиотеку FFTW для функций, использующих FFT - я не знаю деталей реализации FFTW, но я предполагаю, что функции в ней также являются...
В последние годы появилось несколько библиотечно-программных проектов, которые предлагают ту или иную форму параллелизма общей памяти на основе данных общего назначения. Основная идея состоит в том, что вместо написания явно поточного кода программисты реализуют свои алгоритмы как взаимозависимые...
Я очень новичок в программировании GPGPU, поэтому, пожалуйста, прости меня, если вопрос не особенно уместен. Из того, что я понимаю, программирование на GPU - очень сложная часть инженерной работы по сравнению с обычным программированием на CPU. Нужно быть очень осторожным в вопросах расхождения,...
По сути, FEM кажется проблемой, которая в значительной степени «решена». Существует множество мощных сред, таких как Trilinos, PETSc, FEniCS, Libmesh или MOOSE. У них есть одна общая черта: они чрезвычайно «тяжелые». Во-первых, установка обычно очень болезненная. Во-вторых, их интерфейс / API...
Я погружаюсь в увлекательный мир анализа методом конечных элементов и хотел бы решить большую термомеханическую проблему (только термическая механическая, без обратной связи).→→\rightarrow Для механической задачи я уже понял из ответа Джеффа , что мне нужно будет использовать итеративный решатель...
Мне нужно решить Ax = b, но я понимаю, что, даже если оно редкое, сохранение матричных коэффициентов моей задачи займет слишком много памяти. Поэтому сейчас я рассматриваю возможность использования метода без матрицы, потому что одни и те же коэффициенты появляются в матрице много раз, поэтому я...
Большая часть моего программирования - это одноразовые исследовательские коды на Си для моего собственного использования. Я никогда не распространял никакого кода, кроме близких соавторов. Я разработал алгоритм, который я публикую в научном журнале. Я хочу предоставить исходный код и, возможно,...
Как видно из ответа Джеда Брауна « Затраты на поиск по сравнению с вычислениями» , использование векторизованных и не векторизованных операций с плавающей запятой приводит к гораздо более быстрому коду. Многие современные компиляторы утверждают, что они могут выполнять автоматическую векторизацию....