Многие численные подходы к CFD могут быть расширены до произвольно высокого порядка (например, разрывные методы Галеркина, методы WENO, спектральное дифференцирование и т. Д.). Как мне выбрать подходящий порядок точности для данной...
Воспроизводимость становится все более важной в научных исследованиях в области вычислительной техники. (Например, см. Эту статью Роджера Пенга в журнале « Наука» ; мне известны и другие подобные статьи и веб-сайты.) Однако мне не ясно, сколько информации я должен включить в журнальную статью (или...
Каково современное состояние приближения сильно колебательных интегралов как в одном измерении, так и в более высоких измерениях с произвольной
У меня есть ряд точек данных которые, как я ожидаю, (приблизительно) будут следовать функции которая асимптотически пересекается с линией при большом . По существу, приближается к нулю при , и то же самое, вероятно, можно сказать обо всех производных , и т. д. Но я не знаю, что такое функциональная...
Этот вопрос был перенесен из Cross Validated, потому что на него можно получить ответ на бирже вычислительных стеков. Мигрировал 7 лет назад . Я хочу изучать программирование на MATLAB, чтобы я мог проводить некоторые исследования / анализ самостоятельно, а также, чтобы я мог изучать /...
Amazon EC2 Compute Cluster стоит около $ 800- $ 1000 ( в зависимости от рабочего цикла) на физическое ядро процессора в течение 3 -х лет. В нашем последнем раунде приобретения оборудования моя лаборатория взяла на вооружение оборудование на 48 ядер, очень похожее на кластеры Amazon, по цене около $...
Как мы знаем, методы ортогональных преобразований (повороты Гивенса и отражения Хаусхолдера) для систем линейных уравнений более дороги, чем устранение по Гауссу, но теоретически обладают более хорошими свойствами устойчивости в том смысле, что они не изменяют число условий системы. Хотя я знаю...
Я использую пакет нелинейных решателей PETSc SNES для решения системы нелинейных уравнений, полученных путем дискретизации уравнения в частных производных. Как я могу определить, почему решатель не сходится и что я могу сделать, чтобы успешно решить мои...
Дана возрастающая конечная последовательность положительных чисел . Какой из следующих двух алгоритмов лучше для вычисления суммы чисел?Z1, z2, . , , , , ZNz1,z2,.....znz_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end Или: s=0; for \ i=1:n s=s + z_{n-i+1} ; end По моему мнению, было бы...
Собираюсь преподавать студентам бакалавриата курс под названием «Введение в компьютерное программирование». Я немного смущен. В вычислительной физике ученые используют C / C ++ или Python или Fortran, CUDA и т. Д. ..... пришло время создать свою базу. Что я должен использовать? Я знаю, что вы...
Каков предпочтительный и эффективный подход для интерполяции многомерных данных? Вещи, о которых я беспокоюсь: производительность и память для построения, единичная / пакетная оценка размеры обработки от 1 до 6 линейный или более высокого порядка возможность получения градиентов (если не линейная)...
В настоящее время этот вопрос не очень подходит для нашего формата вопросов и ответов. Мы ожидаем, что ответы будут подтверждены фактами, ссылками или опытом, но этот вопрос, скорее всего, вызовет дебаты, споры, опрос или расширенное обсуждение. Если вы считаете, что этот вопрос можно улучшить и,...
Парадигма программирования OpenCL обещает быть открытым стандартом для разнородных вычислений. Должны ли мы инвестировать наше время в разработку программного обеспечения на основе OpenCL? За и...
У меня есть двумерная функция , значения которой я хотел бы попробовать. Функция очень дорога в вычислении и имеет сложную форму, поэтому мне нужно найти способ получить наибольшую информацию о ее форме, используя наименьшее количество точек выборки.е( х , у)е(Икс,Y)f(x,y) Какие хорошие методы...
У меня есть матрицы A и G . A является разреженным и имеет размер n × n с очень большим n (может быть порядка нескольких миллионов). G является матрицей высотой n × m с довольно небольшим m ( 1 < m < 1000 ), и в каждом столбце может быть только один 1 запись с остальным 0 «с, таким образом,...
Что было бы хорошей конечной разностной дискретизацией для следующего уравнения: ∂ρ∂t+∇⋅(ρu)=0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? Мы можем взять 1D случай: ∂ρ∂t+ddx(ρu)=0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho...
Я ищу книги или статьи, или посты в блогах, или любые опубликованные материалы в целом, которые конкретно касаются использования современных функций C ++ (семантика перемещения, STL, итераторы, ленивая оценка и т. Д.) В научных вычислениях. Можете ли вы предложить любой? Я думаю, что эти новые...
Я ищу алгоритмы для рисования стандартных 2d-графиков для функций, которые могут иметь или не иметь особенности. Цель состоит в том, чтобы написать «Мини-CAS», поэтому у меня нет априорных знаний о типах функций, которые хотят видеть пользователи. Эта проблема очень старая, поэтому я предполагаю,...
Я начну с моего личного опыта в нашей лаборатории. В 9 и 10 днях раньше мы были довольно агрессивны с оптимизацией, компилируя с -O3 и специфичными для процессора флагами (например, -xW -xSSE4.2). Но начиная с ifort 11, мы начали замечать: 1. некоторые несоответствия в результатах (поскольку...
Я пытаюсь решить ограниченную задачу оптимизации, в которой я знаю границы некоторых переменных (в частности, рамочное ограничение). ArgминUе( У , х )argminuf(u,x) \arg \min_u f(u,x) при условии a ≤ d ( u , x ) ≤ bс ( и , х ) = 0c(u,x)=0 c(u,x) = 0 a ≤ d( и , х ) ≤ бa≤d(u,x)≤b a \le d(u,x) \le b...