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

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

10 software  education 

Я вирішую систему реакції дифузії Тьюрінга з наступним кодом С ++. Це занадто повільно: для текстури розміром 128x128 пікселів прийнятна кількість ітерацій становить 200 - це призводить до затримки в 2,5 секунди. Мені потрібно 400 ітерацій, щоб отримати цікаве зображення - але 5 секунд очікування - це занадто багато. Також …

10 pde  stiffness 

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

10 c++  b-spline 

Я намагаюсь діагоналізувати деякі щільні, погано обумовлені матриці. У машинній точності результати є неточними (повертаючи негативні власні значення, власні вектори не мають очікуваної симетрії). Я перейшов на функцію Eigensystem [] Mathematica, щоб скористатися довільною точністю, але обчислення надзвичайно повільні. Я відкритий для будь-якої кількості рішень. Чи є пакети / алгоритми, …

10 linear-algebra  parallel-computing  eigensystem  precision  dense-matrix 

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

10 parallel-computing  petsc  gpu  cuda 

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

10 fluid-dynamics 

Чи компіляція PETSc із зовнішньою бібліотекою BLAS / LAPACK суттєво впливає на продуктивність на розріджених матрицях, чи використовує лише ці бібліотеки для математики щільної матриці?

10 petsc  sparse-matrix  lapack  blas 

Чи множиться матричне множення (і Mat * Mat, і Mat * Vec) з кількістю ненульових значень, або з розміром матриці? Або якесь поєднання двох. Що з формою Наприклад, у мене є матриця 100 x 100 із 100 значеннями, або матриця 1000 x 1000 зі 100 значеннями. Під час складання цих …

10 linear-algebra  performance  sparse-matrix 

Я маю зробити перетворення координат між двома опорними системами (осями). Для цього потрібно помножити три матриці ( ) через деякі проміжні осі. Я подумав про два підходи для вирішення цього питання:3 × 33×33\times3 Спосіб №1 : Створення множення безпосередньо, тобто vf= R1 R2 R3 vivf=R1 R2 R3 viv_f = R_1\ …

10 matrix  performance  matlab 

Я намагаюся відповідати предметам. З огляду на набір елементів, я можу класифікувати за шкалою від 0 до 100, наскільки вони схожі один на одного. Наприклад, якщо пункт n 1 - молоко, а пункт n 2 - також молоко, то подібність між n 1 і n 2 складе 100%. Якщо n …

10 regression  data-analysis 

Нехай трикутник , і нехай гладка функція на .ТТTfffТТT Ми можемо використовувати середню точку квадратурной , де є середньою точкою .∫fгx ≈ | Т| ⋅f( хМ)∫fгх≈|Т|⋅f(хМ)\int f dx \approx |T|\cdot f(x_M)хМхМx_MТТT Чи можете ви надати мені (посилання на) формули вищого порядку в симплексі?

10 quadrature 

Для оптимізації з Вікіпедії : В інформатиці метагеврист позначає обчислювальний метод, який оптимізує проблему, ітеративно намагаючись вдосконалити рішення кандидата щодо даного показника якості. Метагевристика робить мало припущень щодо оптимізації проблеми та може шукати дуже великі простори кандидатських рішень. Однак метагевристика не гарантує, що оптимальне рішення ніколи не знайдеться. Багато метагевристики …

10 optimization  heuristics 

Мені цікаво, як знайти власне значення деякої розрідженої матриці в заданому інтервалі [a, b] ітераційним методом. На моє особисте розуміння, більш очевидно використовувати метод підпростору Крилова для пошуку крайніх власних значень, а не внутрішніх.

10 linear-algebra 

Я намагаюся написати повну реалізацію SVM в Python, і у мене є кілька проблем, які обчислюють коефіцієнти Лагранжа. Спочатку дозвольте перефразувати те, що я розумію з алгоритму, щоб переконатися, що я на вірному шляху. Якщо - це набір даних, а - мітка класу , тоy i ∈ { - 1 …

10 optimization  machine-learning  support-vector-machines  quadratic-programming 

Припустимо, я знав, що послідовність випадкових чисел генерується лінійним вродженим генератором. Це є, xn+1=(aXn+c)modmxn+1=(aXn+c)modmx_{n+1}=(aX_n+c) \bmod m Якщо мені дають весь період (або принаймні велику суміжну послідовність його), як я можу реконструювати параметри a,c,ma,c,ma,c,m та x0x0x_0 які створили цю послідовність? Я шукаю загальний метод, який зможе визначити початкові параметри, якщо відомий …

10 random-number-generation  inverse-problem