Я хотел бы выполнить PCA для набора данных, состоящего из приблизительно 40 000 выборок, каждая из которых отображает около 10 000 функций. Использование функции Matlab princomp последовательно занимает более получаса, после чего я завершаю процесс. Я хотел бы найти реализацию / алгоритм, который...
CVXOPT: http://abel.ee.ucla.edu/cvxopt/index.html OpenOpt: http://openopt.org/Welcome Какая связь между ними? Каковы их преимущества / недостатки, соответственно? Кстати, есть ли еще какая-нибудь высококачественная библиотека выпуклой оптимизации общего назначения для Python / C ++, на которую...
Существует ли опубликованный список лучших практик, обеспечивающих долговечность кода, с учетом воспроизводимых научных результатов? (например, открытый исходный код, практика документирования, выбор зависимостей, выбор языка, виртуальные машины и т. д.). Знайте о каких-либо исследованиях (или об...
Большинство математических библиотек имеют несколько версий логарифмических функций. Большую часть времени мы предполагаем, что они идеальны, но на самом деле многие из них просто предлагают определенное количество цифр точности. Для некоторых функций существуют численно более стабильные варианты....
Привет сообщество scicomp, Я работал в области графовых алгоритмов, используя такие фреймворки, как NetworkX (Python), JUNG и YFiles (Java). Сейчас я вхожу в область параллельных и высокопроизводительных вычислений. Для нового проекта я ищу библиотеку графов C ++ со следующими функциями: имеет...
Таким образом, теорема разложения Холецкого утверждает, что любая вещественная симметричная положительно-определенная матрица имеет разложение Холецкого M = L L ⊤, где L - нижняя треугольная матрица.MMMM=LL⊤M=LL⊤M= LL^\topLLL Учитывая , мы уже знаем , есть быстрые алгоритмы для расчета его Чолеска...
Я решаю для огромной разреженной положительно определенной матрицы используя метод сопряженного градиента (CG). Можно вычислить детерминант А, используя информацию, полученную в ходе решения?AA x = bAИксзнак...
Параллельные вычисления часто моделируются с использованием детерминированной локальной скорости вычислений, времени ожидания и пропускной способности сети. На самом деле, они пространственно переменные и недетерминированные. Такие исследования, как Skinner и Kramer (2005), наблюдают...
Хотите улучшить этот пост? Предоставьте подробные ответы на этот вопрос, включая цитаты и объяснение того, почему ваш ответ правильный. Ответы без достаточной детализации могут быть отредактированы или удалены. Я совершенно не знаком с понятием вычислительной науки и ищу хорошую отправную точку. Я...
Скажем, у меня есть функция которую я хочу интегрировать по тетраэдру T ⊂ R 3 . Если бы f было произвольным, квадратура Гаусса была бы хорошим решением, но я случайно узнал, что f - гармоника. Насколько можно ускорить квадратуру Гаусса, используя эту информацию?f:R3→Rf:R3→Rf : \mathbf{R}^3 \to...
Теоретически, как время для вычисления теории функционала плотности (DFT) зависит от количества электронов? Меня интересуют «типичные» реализации DFT, такие как VASP, ABINIT и т. Д., А не O (N)...
У меня есть неоднородная линейная система A x = bAИксзнак равноб Ax=b где - вещественная матрица с . Гарантируется, что нулевое пространство в A имеет нулевую размерность, поэтому уравнение имеет единственный обратный x = A - 1 b . Поскольку результат попадает в правую часть ODE, которую я намерен...
Учитывая неизвестную функцию , мы можем оценить ее значение в любой точке ее области, но у нас нет ее выражения. Другими словами, f для нас как черный ящик.f:Rd→Rf:Rd→Rf:\mathbb R^d \to \mathbb Rfff Как называется проблема поиска минимизатора ? Какие существуют методы?fff Как называется задача...
Как MATLAB, например, вычисляет SVD данной матрицы? Я предполагаю, что ответ, вероятно, включает в себя вычисление собственных векторов и собственных значений A*A'. Если это так, я также хотел бы знать, как он вычисляет...
У меня есть проблемы определяется уравнением Пуассона в двух измерениях −∇2u=f(x,y),inΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y), \; in \; \Omega у меня есть измерения двух компонентов градиента∂u/∂x∂u/∂x\partial{u}/\partial{x} и∂u/∂y∂u/∂y\partial{u}/\partial{y} вдоль некоторой части...
На работе я по сути функционирую как независимый консультант. Для руководства и клиентов мне нужно оценить, сколько времени потребуется на разработку программного обеспечения в рамках моих исследований в области вычислительной техники. Тем не менее, мои оценки времени, как правило, выключены. Я...
У меня есть набор данных 3-мерных точек, для которых я хотел бы построить сетку, используя python. Все программное обеспечение, которое я видел, требует, чтобы вы предоставили края. Есть ли в Python программа, которая принимает в качестве входных данных набор точек в 3D и выводит сетку? Если...
Начиная с уравнения адвекции в форме сохранения. UT= ( a ( x ) u )Иксut=(a(x)u)x u_t = (a(x)u)_x где - скорость, зависящая от пространства, а u - концентрация вида, которая сохраняется.а ( х )a(x)a(x)Uuu Дискретизация потока (где поток , определенный на краях ячеек между точками сетки) дает, u t =...
В настоящее время я разрабатываю метод декомпозиции области для решения задачи рассеяния. По сути, я решаю систему BVP Гельмгольца итеративно. Я дискретизирую уравнения, используя метод конечных элементов по треугольным или тетраэдрическим сеткам. Я разрабатываю код для моей кандидатской...
Мне дают матрицу 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 В...