Я хотел бы смоделировать ламинарный поток воды от корней к стеблю растения. В самом конце корней трубы варьируются от миллиметра до сантиметра в диаметре и длине. По мере приближения к стеблю корни становятся больше по длине и диаметру. Я хочу создать случайные трехмерные домены, представляющие...
Многие алгоритмы, используемые в научных вычислениях, имеют внутреннюю структуру, отличную от алгоритмов, обычно рассматриваемых в менее интенсивных математических формах разработки программного обеспечения. В частности, отдельные математические алгоритмы, как правило, очень сложны, часто включают...
Есть ли алгоритм, который перечисляет графики, которые соответствуют некоторой тесселяции Делоне точек в 3D? Если да, есть ли эффективная параметризация геометрии, которая соответствует любому «графу Делоне»? Я стремлюсь систематически перечислять все стабильные геометрии молекул определенного...
Я ищу библиотеку, которая выполняет матричные операции над большими разреженными матрицами без ущерба для числовой стабильности. Матрицы будут 1000+ на 1000+, а значения матрицы будут между 0 и 1000. Я буду выполнять алгоритм исчисления индекса, поэтому буду генерировать (разреженные) векторы строк...
Я работаю над проблемой, которую можно распараллелить, используя одну операцию mpi_allgather или одну операцию mpi_scatter и одну операцию mpi_gather. Эти операции вызываются в цикле while, поэтому их можно вызывать много раз. В реализации со схемой MPI_allgather я собираю распределенный вектор по...
Существует ли усеченный алгоритм SVD, который вычисляет сингулярные значения по одному? Моя проблема: я хотел бы вычислить первые kkk сингулярных значений (и сингулярных векторов) большой плотной матрицы MMM , но я не знаю, каково было бы подходящее значение kkk . MMM велико, поэтому из соображений...
Я пытаюсь численно инвертировать следующее интегральное преобразование: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)dxF(y)=∫0∞yexp[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Так что для данного мне нужно приблизить где:f (...
У меня есть вопрос, похожий на тот, который задавался ранее, за исключением 3D, и мне нужен только объем, а не фактическая форма корпуса. Точнее, мне дан небольшой набор точек (скажем, 10-15) в 3D, все из которых, как известно, лежат на выпуклой оболочке набора точек (поэтому все они «имеют...
Вопрос Как бы вы отсортировали облако точек относительно неструктурированной сетки шестигранных ячеек? Каждая ячейка имеет центр и уникальный ярлык для ее представления. В основном есть две точки облаков (исходное облако точек и облако точек центров ячеек), но информация о геометрии ячейки...
Таким образом, теорема разложения Холецкого утверждает, что любая вещественная симметричная положительно-определенная матрица имеет разложение Холецкого M = L L ⊤, где L - нижняя треугольная матрица.MMMM=LL⊤M=LL⊤M= LL^\topLLL Учитывая , мы уже знаем , есть быстрые алгоритмы для расчета его Чолеска...
Привет сообщество scicomp, Я работал в области графовых алгоритмов, используя такие фреймворки, как NetworkX (Python), JUNG и YFiles (Java). Сейчас я вхожу в область параллельных и высокопроизводительных вычислений. Для нового проекта я ищу библиотеку графов C ++ со следующими функциями: имеет...
Я делаю диагонализацию Ланцоша большой разреженной матрицы (~ 2 миллиона элементов). Почти все шаги в алгоритме Ланцкоса выполняются параллельно на графическом процессоре, за исключением диагонализации матрицы Ланцоша для проверки сходимости. Для этого я использовал алгоритм TQLI из Numeric...
Это простая задача перечисления симметрии. Я даю полный обзор, но знание квантовой химии не требуется. Двухчастичный интеграл : И имеет следующие 4 симметрии: У меня есть функция, которая вычисляет интеграл и сохраняет его в одномерном массиве , проиндексированном следующим образом:⟨ я J | к л ⟩ =...
Учитывая положительно определенную симметричную матрицу, каков самый быстрый алгоритм вычисления обратной матрицы и ее определителя? Для задач, которые меня интересуют, размер матрицы равен 30 или меньше. Высокая точность и скорость действительно необходимы. (выполняются миллионы матриц)...
У меня плотная вещественная симметричная квадратная матрица. Размер составляет около 1000x1000. Мне нужно вычислить первый главный компонент и подумать, каким может быть лучший алгоритм для этого. Похоже, что MATLAB использует алгоритмы Арнольди / Ланцоша (для eigs). Но читая о них, я не уверен,...
Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он подходил для обмена стеками вычислительной науки. Закрыто 5 лет назад . Существует ли степень сложности, которая больше, чем и меньше, чем...
Какие ситуации численного анализа становятся более / менее стабильными, имеют более быструю / медленную сходимость или иным образом сильно различаются при работе с функциями комплексной переменной вместо функций реальной переменной?...
Позволять f,g∈R[x]f,g∈R[x]f, g \in \mathbb{R}[x] а также degf>deggdegf>degg\deg f > \deg g, Я ищу асимптотически быстрые и численно устойчивые алгоритмы для вычисленийfmodgfmodgf \bmod g, В предназначенных приложениях обаf,gf,gf, gплотные полиномы с коэффициентами с плавающей точкой двойной...
Я хотел бы знать, если есть алгоритм, который дает множество точек и угол вычисляет выпуклую оболочку, если угол α = 0α=0\alpha = 0 и дали α > 0α>0\alpha > 0 вычисляет конверт, который более близко следует за «периметром». И если есть определение непересекающегося периметра множества точек, в...
Учитывая общую разреженную матрицу A ∈рn × nA∈рN×NA \in \mathbb{R}^{n\times n}с m << n (поправка:м ≪N2м«N2m \ll n^2ненулевые элементы (обычно m ∈ O ( n )м∈О(N)m \in {\cal O}(n)). AAA является общим в том смысле, что не имеет специфических свойств (например, положительной определенности), и не...