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

Какая библиотека является самой быстрой для выполнения триангуляции множеств Делона с миллионами, если 3D-точки? Доступны ли также версии GPU? С другой стороны, наличие вороной тесселяции того же набора точек поможет (с точки зрения производительности) получить триангуляцию Делоне?

26 computational-geometry  delaunay-triangulation  voronoi-diagrams 

Что такое простой алгоритм для вычисления SVD матрицы?2 × 22×22 \times 2 В идеале мне нужен численно устойчивый алгоритм, но я бы хотел увидеть как простые, так и не очень простые реализации. Код C принят. Любые ссылки на документы или код?

26 linear-algebra  matrix  reference-request  svd  matrix-factorization 

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

26 data-storage 

Разработчики программного обеспечения имеют возможность выбрать подходящую лицензию в соответствии с целью (целями) работы. Может ли кто-нибудь дать некоторые рекомендации / опыт, какую лицензию выбрать для программного обеспечения? Каковы плюсы / минусы «раздачи» всей закодированной работы в виде открытых исходных кодов? Как бороться с промышленными игроками, которые хотели бы получить …

26 software 

Что может пойти не так, если использовать методы Крылова из KSP ( пакет линейного решателя PETSc ) с предварительными именами для решения разреженной линейной системы, например, полученной путем дискретизации и линеаризации уравнений в частных производных? Какие шаги я могу предпринять, чтобы определить, что идет не так для моей проблемы? Какие …

26 linear-algebra  petsc  preconditioning  krylov-method  iterative-method 

Я ищу инструмент для визуализации очень больших графиков направленных ссылок. В настоящее время у меня есть ~ 2 миллиона узлов с ~ 10 миллионами ребер. Я пробовал несколько разных вещей, но большинству требуются часы, чтобы даже сделать графы узлов 100k Что я пробовал: я провел день с Gephi, но для …

25 visualization 

Имея реальную функцию реальных переменных, существует ли программное обеспечение, которое может автоматически генерировать числовой точный код для вычисления функции по всем входам на машине, оборудованной арифметикой IEEE 754? Например, если реальная функция для оценки была: Программное обеспечение будет учитывать катастрофическое аннулирование и, возможно, поиск таблиц для определенных наборов входных данных, …

25 software  floating-point  accuracy 

Мне нужно численно оценить интеграл ниже: ∫∞0с я н с'(xr)rE(r)−−−−√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr где , и . Здесь - модифицированная функция Бесселя второго рода. В моем конкретном случае у меня , и .x∈R+λ,κ,ν>0Kλ=0,00313κ=0,00825ν=0,33E(r)=r4(λκ2+r2−−−−−−√)−ν−5/2K−ν−5/2(λκ2+r2−−−−−−√)E(r)=r4(λκ2+r2)−ν−5/2K−ν−5/2(λκ2+r2)E(r) = r^4 (\lambda\sqrt{\kappa^2+r^2})^{-\nu-5/2} K_{-\nu-5/2}(\lambda\sqrt{\kappa^2+r^2})Икс∈R+x∈R+x \in \mathbb{R}_+λ,κ,ν> 0λ,κ,ν>0\lambda, \kappa, \nu >0КKKλ = 0,00313λ=0.00313\lambda = 0.00313κ = 0,00825κ=0.00825\kappa …

25 matlab  quadrature  special-functions 

Смежный вопрос: Состояние Mac OS в научных вычислениях и HPC Значительное количество программных пакетов в вычислительной науке написано на Фортране, и Фортран не уходит. Компилятор Фортрана также необходим для сборки других программных пакетов (одним из примечательных примеров является SciPy ). Тем не менее, Mac OS X не включает в себя …

25 software  fortran 

Какие соображения следует учитывать при выборе между BFGS и сопряженным градиентом для оптимизации? Функция, которую я пытаюсь согласовать с этими переменными, является экспоненциальной функцией; тем не менее, фактическая целевая функция включает в себя, среди прочего, интеграцию и очень дорогая, если это вообще помогает.

25 optimization  conjugate-gradient 

Как гравитационная проблема n-тела может быть решена численно параллельно? Возможен ли компромисс между точностью и сложностью? Как точность влияет на качество модели?

25 algorithms  ode  parallel-computing  complexity  precision 

В этом комментарии я написал: ... интегратор SciPy по умолчанию, который, как я предполагаю, использует только симплектические методы. в котором я ссылаюсь на SciPy's odeint, который использует либо «нежесткий (Адамс) метод», либо «жесткий (BDF) метод». По словам источника : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, …

25 ode  time-integration  integration  differential-equations 

Я не понимаю разного поведения уравнения адвекции-диффузии, когда применяю разные граничные условия. Моя мотивация - симуляция реальной физической величины (плотности частиц) в условиях диффузии и адвекции. Плотность частиц должна сохраняться внутри, если только она не вытекает из краев. По этой логике, если я приведу в исполнение граничные условия Неймана, концы …

25 pde  boundary-conditions  advection  diffusion  crank-nicolson 

Вообще говоря, я слышал, что числовые аналитики высказывают мнение, что «Конечно, с математической точки зрения, время это просто другое измерение, но все же, время является особенным» Как это обосновать? В чем смысл времени особенного для вычислительной науки? Кроме того, почему мы так часто предпочитаем использовать конечные различия (приводящие к «временному …

24 pde  discretization 

При переходе от сильной формы PDE к форме FEM кажется, что всегда следует делать это, сначала указав вариационную форму. Для этого вы умножаете сильную форму на элемент в некотором (соболевском) пространстве и интегрируете по своему региону. Это я могу принять. Я не понимаю, почему нужно использовать формулу Грина (один или …

24 pde  finite-element  elliptic-pde