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

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

10 computational-geometry  monte-carlo  geometry  random-sampling 

Недавно я написал код, который пытался вычислить без использования библиотечных функций. Вчера я просматривал старый код и пытался сделать его как можно быстрее (и исправить). Вот моя попытка:log(x)log(x)log(x) const double ee = exp(1); double series_ln_taylor(double n){ /* n = e^a * b, where a is an non-negative integer */ double …

10 numerical-analysis 

Я не вижу упоминаний об этом в списках вакансий. Я видел упомянутое целочисленное программирование, MIP, смешанно-целочисленное нелинейное программирование, LP, динамическое программирование и т. Д., Но без SDP. Это намного моднее в академии, чем в промышленности? Из-за моего ограниченного контакта с учеными и промышленными участниками в электроэнергетических системах, я думаю, есть …

10 convex-optimization  semidefinite-programming 

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он подходил для обмена стеками вычислительной науки. Закрыто 5 лет назад . Существует ли степень сложности, которая больше, чем и меньше, чем ?O(n)O(n)O(n)O(nlogn)O(nlog⁡n)O(n \log n)

10 algorithms  complexity  efficiency 

Я хотел бы использовать метод 8-го порядка Рунге-Кутты (89) в приложении по небесной механике / астродинамике, написанном на C ++, на машине Windows. Поэтому мне интересно, знает ли кто-нибудь хорошую библиотеку / реализацию, которая документирована и бесплатна для использования? Это нормально, если он написан на C, если нет никаких проблем …

10 c++  reference-request  libraries  c  runge-kutta 

Я использовал FEA в течение нескольких лет, но его использование и правильное использование - это две разные вещи, фактор безопасности не является решением для всего. У меня такое чувство, что я не буду использовать это правильно, если у меня нет четкого ответа на этот вопрос: Я знаю, что элементы должны …

10 finite-element  simulation 

Область вычислительной гидродинамики (CFD) посвящена решению уравнений Навье-Стокса (или некоторому их упрощению). Подмножество моделей CFD, океана и атмосферы численно решают те же уравнения для реалистичных приложений. Каковы различия и компромиссы между общими подходами CFD и применяемыми реалистичными случаями?

10 fluid-dynamics  simulation  computational-physics  numerical-modelling 

Например, библиотеки разреженных матриц C ++, которые я использовал - Eigen и SuiteSparse, похоже, что они не имеют какой-либо функциональности SVD для разреженных матриц. Так что просто любопытно, является ли SVD сложнее, чем QR / LU для разреженной матрицы?

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

Я часто сталкиваюсь с общей пословицей, что методы внутренней точки трудно начать с теплого начала. Есть ли интуитивное объяснение этому совету? Существуют ли ситуации, в которых можно ожидать выгоды от теплого старта в методе внутренней точки? Кто-нибудь может порекомендовать несколько полезных ссылок по теме?

10 optimization  constrained-optimization 

Я пытаюсь понять разницу между этими двумя видеокартами для академических вычислений, особенно для компонента DGEMM. Если мы посмотрим на необработанные статистические данные, оба имеют одинаковый чип GK110, имеют сопоставимые статистические данные практически во всех категориях и, как мне кажется, имеют одинаковую архитектуру ядра. До каких-либо скидок K20X примерно в 4 …

10 performance  gpu  efficiency 

Существует много хорошо известных численных методов решения уравнений типа например, метод деления пополам, метод Ньютона и т. Д.е( х ) = 0 ,x ∈ RN,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, В моем приложении рассчитывается стохастическим методом (результат является средним).е( х )е(Икс)f(x) Существуют ли методы решения численных уравнений, которые …

10 monte-carlo  nonlinear-equations  stochastic  numerics 

Я работаю над некой обусловленной большой разреженной линейной системой уравнений. Я хочу использовать двойную арифметику или квадро-двойную арифметику для их решения. Я знаю, что существует пакет под названием MPACK, разработанный Накатой, Махо, который может выполнять числовые линейные алгебраические вычисления в четырехугольной арифметике. Тем не менее, он предназначен для плотной матрицы, …

10 sparse-matrix  matrix  condition-number 

Я хотел бы знать, как обычно применяются условия Дирихле при использовании метода конечных объемов на неоднородной сетке с центром на ячейках, Моя текущая реализация просто накладывает граничное условие, фиксируя значение первой ячейки, φ1= гD( хL)ϕ1=gD(xL) \phi_1 = g_D(x_L) где - это переменная решения, а - значение граничного условия Дирихле в …

10 boundary-conditions  finite-volume  poisson 

Я решаю нелинейную систему связанных уравнений и вычисляю якобиан дискретной системы. Результат действительно сложный, ниже приведены (только!) Первые 3 столбца матрицы ,3 × 93×93\times 9 (Сложность возникает отчасти потому, что численная схема требует экспоненциального подбора для устойчивости.) У меня довольно общий вопрос о реализации числовых кодов с использованием якобианов. Я …

10 symbolic-computation  jacobian  design-patterns  sympy 

Я новичок в FE. Мое приложение - оценка финансовых производных, где пространство является пятимерным. Итак, добавляя время, проблема имеет шесть измерений. Я пытался осмотреться (Fenics, escript, deal.II, ...), но, насколько я понимаю, эти программы ограничены 3 + 1 (трехмерное пространство + 1d время). Это правильно? Мой целевой язык Python или …

10 pde  finite-element