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

Наша рабочая группа выпускает настольное приложение, которое имитирует энергоэффективность здания. Это приложение .NET, и когда пользователь запускает много симуляций, они могут занимать довольно много времени. Моделирование полностью распараллеливаемо, и у нас есть несколько очень значительных вычислительных ресурсов HPC в «офисе». У нас была одна идея - позволить пользователям разгружать модели, …

17 hpc  partitioning 

Я знаю, что метод Ницше является очень привлекательным, поскольку он позволяет учитывать граничные условия типа Дирихле или слабо связываться с граничными условиями трения без использования множителей Лагранжа. И его преимущество, заключающееся в преобразовании граничного условия Дирихле в слабые члены, аналогично граничному условию Неймана, обусловлено тем фактом, что реализация зависит от …

17 finite-element  boundary-conditions  elliptic-pde  nitsche-method 

Я узнал о методе конечных элементов (также немного о других численных методах), но я не знаю, каково именно определение этих двух ошибок и различий между ними?

16 finite-element  numerical-analysis  error-estimation 

В учебных целях мне понадобится непрерывная функция от одной переменной, которую «трудно» аппроксимировать полиномами, то есть, для степенного ряда достаточно сильных степеней, чтобы «соответствовать» этой функции. Я намерен показать своим ученикам «пределы» того, что может быть достигнуто с помощью степенных рядов. Я думал о том, чтобы придумать что-то «шумное», но …

16 interpolation 

Мне известно о том, что инвертировать матрицу для решения линейной системы не очень хорошая идея, поскольку она не так точна и эффективна, как непосредственное решение системы или использование разложения LU, Cholesky или QR. Однако я не смог проверить это на практическом примере. Я пробовал этот код (в MATLAB) M = …

16 matrix  accuracy  inverse 

Поскольку этот вопрос связан с вычислениями, я решил опубликовать здесь. Надеюсь, это будет рассматриваться как уместно. Я только начал работать с атмосферными и океаническими моделями и понял, что мне нужно больше ядер, памяти и дискового пространства, чем у моего нынешнего рабочего стола. В моем учреждении есть высокопроизводительная вычислительная служба, где …

16 hpc 

Я хотел бы знать, есть ли быстрый способ вычисления евклидова расстояния двух векторов в октаве. Кажется, что для этого нет специальной функции, так что я должен просто использовать формулу с sqrt?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Недавно я сравнивал различные нелинейные решатели от scipy и был особенно впечатлен примером Ньютона-Крылова в книге рецептов Scipy, в которой они решают уравнение дифференциального уравнения второго порядка с нелинейным членом реакции примерно в 20 строках кода. Я изменил пример кода для решения нелинейного уравнения Пуассона ( также называемого уравнением Пуассона-Больцмана …

16 python  nonlinear-equations  poisson  scipy 

Я большой сторонник разработки, основанной на тестировании, в научных вычислениях. Его полезность на практике просто ошеломляет и действительно облегчает классические проблемы, известные разработчикам кода. Однако при тестировании научных кодов, которые не встречаются в общем программировании, есть определенные трудности, поэтому тексты TDD не очень полезны в качестве учебных пособий. Например: Как …

16 testing 

Можно получить много понимания из опыта, мне просто интересно, видел ли кто-нибудь что-то подобное раньше. График показывает начальное условие (зеленый) для уравнения адвекции-диффузии, затем решение на итерации 200 (синий) и затем снова на итерации 400 (красный). Решение уравнения адвекции-диффузии взрывается после нескольких итераций. Число Пекле , а условие КЛЛ выполнено, …

16 numerical-analysis  advection-diffusion 

Я хотел бы быть независимым от коммерческого программного обеспечения для моей научной работы. Я нахожу зависимость от коммерческих пакетов, таких как Matlab и ее наборов инструментов, неудовлетворительной, потому что я не знаю, получу ли я доступ к Matlab в будущем, и потому что мне не нравится язык. Поэтому я ищу …

16 python  machine-learning 

При программировании вычислений с плотной матрицей, есть ли причина выбирать макет с основной строкой над макетом с основной колонкой? Я знаю, что в зависимости от макета выбранной матрицы нам нужно написать соответствующий код, чтобы эффективно использовать кэш-память для повышения скорости. Раскладка с основными строками кажется более естественной и простой (по …

16 matrix  fortran 

Больше всего в MPI мне не нравится иметь дело с типами данных (то есть картами данных / масками), потому что они не очень хорошо подходят для объектно-ориентированного C ++. boost::mpiоднако поддерживает только MPI 1.1 со своего веб-сайта: boost :: mpi - это дружественный к C ++ интерфейс стандартного интерфейса передачи …

16 hpc  c++  mpi 

Я читал о вычислительной науке в Википедии, но мое понимание не очень ясно. Включает ли вычислительная наука программирование? Насколько отличается вычислительная наука от вычислительной _ ____ , где пробелом может быть любая дисциплина (материаловедение, инженерия, химия, биология и т. Д.)? (Я буду заниматься вычислительной материаловедением.)

16 computational-chemistry  education 

Вероятно, это вопрос студенческого уровня, но я не могу сделать его понятным для себя. Почему более точно использовать неоднородные сетки в численных методах? Я думаю в контексте некоторого метода конечных разностей для PDE вида . И предположим, что меня интересует решение в точке x ∗ . Итак, я могу видеть, …

16 pde  finite-difference