Пусть треугольник , и пусть гладкая функция на .TTTееfTTT Мы можем использовать среднюю точку квадратурной , где является средней точкой .∫еdх ≈ | T| ⋅ф( хM)∫еdИкс≈|T|⋅е(ИксM)\int f dx \approx |T|\cdot f(x_M)ИксMИксMx_MTTT Можете ли вы предоставить мне (ссылку на) формулы высшего порядка на...
Я пытаюсь сопоставить предметы. Учитывая набор из предметов, я могу оценить по шкале от 0 до 100, насколько они похожи друг на друга. Например, если элемент n 1 является молоком, а элемент n 2 также является молоком, то сходство между n 1 и n 2 будет 100%. Если n 3 является содой, то сходство между...
Я пытаюсь написать полную реализацию SVM на Python, и у меня есть несколько проблем с вычислением коэффициентов Лагранжа. Сначала позвольте мне перефразировать то, что я понимаю из алгоритма, чтобы убедиться, что я на правильном пути. Если x1,x2,...,xnx1,x2,...,xnx_1, x_2, ..., x_n - это набор...
Для оптимизации из Википедии : В информатике метаэвристика обозначает вычислительный метод, который оптимизирует проблему путем итеративной попытки улучшить подходящее решение с учетом заданного показателя качества. Метаэвристика делает мало или вообще не делает предположений об оптимизируемой...
Я должен сделать преобразование координат между двумя системами отсчета (осями). Для этого необходимо умножить три матрицы ( ) из-за использования некоторых промежуточных осей. Я думал о двух подходах, чтобы решить эту проблему:3×33×33\times3 Способ № 1 : Выполнение умножения напрямую, то есть...
Умножается ли матричное умножение (как Mat * Mat, так и Mat * Vec) на количество ненулевых элементов или на размер матрицы? Или какая-то комбинация двух. Как насчет формы. Например, у меня есть матрица 100 x 100 с 100 значениями в ней или матрица 1000 x 1000 с 100 значениями в ней. При возведении в...
Влияет ли компиляция PETSc с внешней библиотекой BLAS / LAPACK существенным образом на производительность разреженных матриц, или он использует эти библиотеки только для математики с плотной...
По заданию в колледже я сделал 1d симуляцию. Постановка задачи состояла в том, чтобы решить проблему 1-мерной ударной трубы, используя сжимаемый идеальный газ в качестве рабочей жидкости Для этой задачи я решил систему уравнений Эйлера с помощью решателя Римана Роу. Я хочу знать, чтобы решить...
Сейчас я работаю над методом решения дифференциальных уравнений, который называется базис-сплайн-коллокацией. У меня возникли проблемы с созданием метода построения сплайна произвольного порядка с соотношением с начальным условием B 1 i (x)={ 1Вк + 1я( х ) = х - хяИкск + я- хяВКя+ хk + i + 1- хИксk...
Я пытаюсь диагонализировать некоторые плотные, плохо обусловленные матрицы. В машинной точности результаты являются неточными (возвращая отрицательные собственные значения, собственные векторы не имеют ожидаемой симметрии). Я переключился на функцию Eigensystem [] Mathematica, чтобы использовать...
Я очень новичок в программировании GPGPU, поэтому, пожалуйста, прости меня, если вопрос не особенно уместен. Из того, что я понимаю, программирование на GPU - очень сложная часть инженерной работы по сравнению с обычным программированием на CPU. Нужно быть очень осторожным в вопросах расхождения,...
Я решаю систему реакции-диффузии Тьюринга с помощью следующего кода C ++. Это слишком медленно: для текстуры 128x128 пикселей допустимое количество итераций составляет 200, что приводит к задержке в 2,5 секунды. Мне нужно 400 итераций, чтобы получить интересное изображение, но 5 секунд ожидания -...
Я начинаю свою аспирантуру в области машиностроения и буду работать над проектами в области вычислительной науки. Я заметил , что там была некоторая дискуссия о преимуществах и недостатках реализации собственных алгоритмов здесь . Отличается ли анализ затрат и выгод от реализации ваших собственных...
Существует ли открытая C-реализация для решения квартичных уравнений: ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0 Я думаю о реализации решения Ferrari. В Википедии я читал, что решение устойчиво в вычислительном отношении только для некоторых из возможных комбинаций знаков коэффициентов....
Рассмотрим следующую задачу где форсирующий член может зависеть от (см. Формулировку 1 ниже для формулировки), а также от и его первых производных. Это 1 + 1 мерное волновое уравнение. У нас есть начальные данные, прописанные в .Wuv=FWuv=F W_{uv} = F u,vu,vu,vWWW{u+v=0}{u+v=0}\{u+v = 0\} Меня...
Я пытаюсь оптимизировать приложение MPI с высоко асинхронным шаблоном связи. Каждый ранг имеет список вещей для вычисления и отправляет сообщения по мере необходимости, если входы или выходы находятся на другом ранге. Кроме того, каждый ранг является многопоточным (в настоящее время с одним...
В последние годы появилось несколько библиотечно-программных проектов, которые предлагают ту или иную форму параллелизма общей памяти на основе данных общего назначения. Основная идея состоит в том, что вместо написания явно поточного кода программисты реализуют свои алгоритмы как взаимозависимые...
Каков современный способ реализации специальных функций двойной точности? Мне нужен следующий интеграл: для и , что можно записать в терминах нижней неполной гамма-функции. Вот моя реализация на Фортране и Си: м=0,1,2,. , , t>0Fм( т ) = ∫10U2 ме- т у2dу = γ( м + 12, т )2 тм + 12Fм(T)знак...
У меня есть числовая функция, f(x, y)возвращающая двойное число с плавающей запятой, которая реализует некоторую формулу, и я хочу проверить, является ли она корректной по отношению к аналитическим выражениям для всех комбинаций параметров xи yкоторая мне интересна. Как правильно сравнивать...
Я заинтересован в максимизации функции , где .θ ∈ R pе( θ )f(θ)f(\mathbf \theta)θ ∈ Rпθ∈Rp\theta \in \mathbb R^p Проблема в том, что я не знаю аналитической формы функции или ее производных. Единственное, что я могу сделать, это оценить функцию по точкам, подключив значение и получить оценку NOISY...