Обчислювальна наука

Багато моєї власної роботи обертається над тим, щоб покращити масштаб алгоритмів, і один із кращих способів показу паралельного масштабування та / або паралельної ефективності - це побудувати графік продуктивності алгоритму / коду за кількістю ядер, наприклад де -ось представляє кількість ядер, а y -ось деяка метрика, наприклад робота, виконана за …

17 parallel-computing  hpc  scaling  plotting 

Припустимо, у мене є деяка функція і я хочу знайти x таку, що f ( x ) ≈ 0 . Я можу використовувати метод Ньютона-Рафсона. Але для цього потрібно, щоб я знав похідну функцію f ′ ( x ) . Аналітичний вираз дляfffxxxf(x)≈0f(x)≈0f(x)\approx 0f′(x)f′(x)f'(x) може бути недоступним. Наприклад, f може …

17 reference-request  approximation 

Мені було цікаво, чи існує швидкий та ефективний метод заздалегідь знайти кількість ненульових значень для операції розмноження матричної матриці, якщо обидві матриці знаходяться у форматі CSC або CSR. Я знаю, що в пакеті smmp є такий, але мені потрібно щось, що вже реалізовано в C або C ++. Будь-яка допомога …

17 matrix  sparse-matrix 

Паралельне зменшення передбачає, що відповідна операція є асоціативною. Це припущення порушено для додавання чисел з плаваючою комою. Ви можете запитати, чому я переймаюся цим. Ну, це робить результати менш відтворюваними. І стає гірше, коли імітований відпал використовується для оптимізації (або підгонки параметрів) над підпрограмами, що дають такі невідтворювані результати. Які …

17 parallel-computing  reproducibility 

Днями мій інструктор з обчислювальної динаміки рідини був відсутній, і він направив свого кандидата кандидатів наук на заміну. У своїй лекції він, схоже, вказав на кілька недоліків, пов'язаних з різними схемами дискретизації для моделювання потоку рідини: Метод скінченної різниці: важко задовольнити збереження та застосувати неправильну геометрію Метод кінцевого обсягу: він …

17 fluid-dynamics 

У мене є обмежена невипукла 2-D функція, яку я хотів би знайти мінімум. Функція досить плавна. Оцінити це дорого. Прийнятна помилка становить близько 3% домену функції в кожній осі. Я спробував запустити реалізацію алгоритму DIRECT в бібліотеці NLOPT, але це не дало значного покращення в порівнянні з пошуком грубої сили …

17 optimization 

Мова C ++ забезпечує загальне програмування та метапрограмування через шаблони. Ці методи знайшли свій шлях у багато великих наукових обчислювальних пакети (наприклад, MPQC , LAMMPS , CGAL , Trilinos ). Але що вони насправді сприяли науково-обчислювальним значенням, що виходить за рамки загальних, неметальних мов, таких як C або Fortran, з …

17 performance  c++  programming-paradigms 

Коли обчислення з обмеженою пропускною здатністю пам’яті виконуються в середовищах спільної пам’яті (наприклад, з потоком через OpenMP, Pthreads або TBB), існує дилема, як забезпечити правильну розподіл пам’яті по фізичній пам’яті, таким чином, щоб кожен потік в основному отримував доступ до пам’яті "локальна" шина пам'яті. Хоча інтерфейси не є портативними, у …

17 performance  memory-management  multicore 

Я завжди обчислюю речі вручну, але зараз мої товариші стають неприємними і роблять багато повторюваних вправ, що включають просто підключення речей, як вираз вище. Мене особливо цікавить програмне забезпечення з відкритим кодом, таким як Python або R для спрощення таких рівнянь. Я спробував використовувати Wolfram Alpha , але мені не …

17 education  symbolic-computation 

У літературі про числові методи багато дослідницьких робіт складаються з опису нової алгоритмічної варіації з подальшими кількома тестовими проблемами порівняння нового методу з одним або двома існуючими методами. Це ускладнює визначення Як новий метод діятиме на інші цікаві проблеми Як новий метод порівнявся б із іншими існуючими методами Звичайно, обидва …

17 pde  testing 

Я планую використовувати мультирешітку для визначення деяких власних значень та векторів, і я помітив, що PETSc має підтримку на високому рівні для мультисетки. У документації PETSc зазначено, що ця частина PETSc не повинна використовуватися, оскільки вона незабаром буде замінена. Які ще бібліотеки мають високу підтримку для мультирешітки, і приблизно, як …

17 petsc  libraries  multigrid 

Розглянемо IVP для системи ODE y′=f(x,y)y′=f(x,y)y'=f(x,y) , . Найчастіше ця проблема вважається жорсткою, коли матриця Якобі має як власні значення з дуже великою негативною реальною частиною, так і власні значення з дуже малою негативною реальною частиною (я вважаю лише стабільну випадок).y(x0)=y0y(x0)=y0y(x_0)=y_0∂f∂y(x0,y0)∂f∂y(x0,y0)\frac{\partial f}{\partial y}(x_0,y_0) З іншого боку, у випадку лише одного …

17 ode  stiffness 

Розробка наукових алгоритмів - це дуже ітераційний процес, який часто включає зміну безлічі параметрів, які я хочу змінити або як частина мого експериментального дизайну, або як частина роботи алгоритму налаштування. Які стратегії я можу взяти за структурування цих параметрів, щоб я міг легко змінювати їх між ітераціями, щоб я міг …

17 algorithms  testing  programming-paradigms 

Моя дипломна робота полягає у розробці чисельних методів зменшення моделі згоряння. Я керую своїми методами виключно на хімічній частині моделювання горіння, і у мене є безліч прикладів для 0-D моделювання (без потоку). Мені хотілося б запустити симуляції, які мають в них потік, бажано 2-D або 3-D моделювання. Ці симуляції повинні …

17 pde  fluid-dynamics  software  parallel-computing 

Повернувшись до світанку OS X, здавалося, було дуже багато гомін, принаймні в світі Mac (я тоді ніде не був поблизу наукових обчислень) щодо Mac OS як платформи для наукових обчислень та додатків HPC. XGrid вийшов з коробки, Virginia Tech мали свої вигадливі комп'ютерні кластери на основі Mac, Стенфорд робив круті …

17 software  hpc