Я большой сторонник разработки, основанной на тестировании, в научных вычислениях. Его полезность на практике просто ошеломляет и действительно облегчает классические проблемы, известные разработчикам кода. Однако при тестировании научных кодов, которые не встречаются в общем программировании, есть...
В учебных целях мне понадобится непрерывная функция от одной переменной, которую «трудно» аппроксимировать полиномами, то есть, для степенного ряда достаточно сильных степеней, чтобы «соответствовать» этой функции. Я намерен показать своим ученикам «пределы» того, что может быть достигнуто с...
Больше всего в MPI мне не нравится иметь дело с типами данных (то есть картами данных / масками), потому что они не очень хорошо подходят для объектно-ориентированного C ++. boost::mpiоднако поддерживает только MPI 1.1 со своего веб-сайта: boost :: mpi - это дружественный к C ++ интерфейс...
Я делаю задачу классификации текста с помощью R и получаю матрицу терминов документа размером 22490 на 120 000 (только 4 миллиона ненулевых записей, менее 1% записей). Теперь я хочу уменьшить размерность, используя PCA (анализ основных компонентов). К сожалению, R не может обработать эту огромную...
Я начал использовать python в качестве языка программирования для выполнения всех моих заданий в CFD. У меня очень мало опыта в программировании. Я из области машиностроения и продолжаю получать высшее образование в области аэрокосмической техники. иногда вычислительный аспект CFD становится более...
Мне любопытно узнать, какие хорошие численные алгоритмы существуют для оценки обобщенной гипергеометрической функции (или ряда), определенной как pFq(a1,…,ap;b1,…,bq;z)=∑k=0∞(a1)k⋯(ap)k(b1)k⋯(bq)kzkk!pFq(a1,…,ap;b1,…,bq;z)=∑k=0∞(a1)k⋯(ap)k(b1)k⋯(bq)kzkk!{}_pF_q(a_1,\ldots,a_p;b_1,\ldots,b_q;z) =...
При расчете факторизации QR на практике используют отражения Домохозяина для обнуления нижней части матрицы. Я знаю, что для вычисления собственных значений симметричных матриц лучшее, что вы можете сделать с отражениями Хаусхолдера, - это привести его к трехдиагональной форме. Есть ли очевидный...
Для заданной плотной матрицы каков наилучший способ найти его нулевое пространство в пределах некоторого допуска ?A∈Rm×n,m>>n;max(m)≈100000A∈Rm×n,m>>n;max(m)≈100000A \in R^{m \times n}, m >> n; max(m) \approx 100000 ϵϵ\epsilon Исходя из этого, могу ли я сказать, что некоторые столбцы...
предполагать мин АV e C (U)подлежит Uя , дж≤ max { Uя , к, UК , Дж} ,i , j , k = 1 , … , nminAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,n\begin{align*} \min A &\mathrm{vec}(U) \\ &\text{subject to } U_{i,j} \leq \max\{U_{i,k}, U_{k,j}\}, \quad i,j,k = 1, \ldots, n \end{align*} где -...
Поскольку этот вопрос связан с вычислениями, я решил опубликовать здесь. Надеюсь, это будет рассматриваться как уместно. Я только начал работать с атмосферными и океаническими моделями и понял, что мне нужно больше ядер, памяти и дискового пространства, чем у моего нынешнего рабочего стола. В моем...
Я использую Valgrind + Callgrind для профилирования написанного мной решателя. Как говорится в руководстве пользователя Valgrind, я скомпилировал свой код с опциями отладки для компилятора: «Без отладочной информации лучшее, что смогут сделать инструменты Valgrind, - это угадать, к какой функции...
Я хотел бы быть независимым от коммерческого программного обеспечения для моей научной работы. Я нахожу зависимость от коммерческих пакетов, таких как Matlab и ее наборов инструментов, неудовлетворительной, потому что я не знаю, получу ли я доступ к Matlab в будущем, и потому что мне не нравится...
Я узнал о методе конечных элементов (также немного о других численных методах), но я не знаю, каково именно определение этих двух ошибок и различий между
Я читал о вычислительной науке в Википедии, но мое понимание не очень ясно. Включает ли вычислительная наука программирование? Насколько отличается вычислительная наука от вычислительной _ ____ , где пробелом может быть любая дисциплина (материаловедение, инженерия, химия, биология и т. Д.)? (Я...
Мне известно о том, что инвертировать матрицу для решения линейной системы не очень хорошая идея, поскольку она не так точна и эффективна, как непосредственное решение системы или использование разложения LU, Cholesky или QR. Однако я не смог проверить это на практическом примере. Я пробовал этот...
Я хотел бы знать, есть ли быстрый способ вычисления евклидова расстояния двух векторов в октаве. Кажется, что для этого нет специальной функции, так что я должен просто использовать формулу с...
Предположим, у меня есть двумерная сетка, состоящая из непересекающихся треугольников и набора точек . Как лучше всего определить, в каком треугольнике лежит каждая из точек?{Tk}Nk=1{Tk}k=1N\{T_k\}_{k=1}^N{pi}Mi=1⊂∪Nk=1TK{pi}i=1M⊂∪k=1NTK\{p_i\}_{i=1}^M \subset \cup_{k=1}^N T_K Например, на...
Есть ли причины, по которым следует выбирать неявную Рунге Кутту высокого порядка (IMRK), а не BDF? BDF кажется мне намного проще, так как для ступени IMRK требуется q линейных решений за шаг по времени. Стабильность для BDF и IMRK выглядит спорным вопросом. Я не могу найти никаких ресурсов,...
Какие подходы используются на практике для оценки числа условий больших разреженных
Существует система линейных ограничений . Я хочу найти строго положительный вектор который удовлетворяет этим ограничениям. Это означает, что требуется для каждого компонента из . Как я могу использовать LP-решатель, чтобы найти такой строго положительный вектор (или подтвердить, что существует)? Я...