Запитання з тегом «algorithms»

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

16 algorithms  python  fluid-dynamics 

Я дуже вражений серійними характеристиками багаторівневих попередніх кондиціонерів ILU , особливо для гетерогенних Helmholtz , але я здивований тим, що не в змозі знайти жодних реалізацій з відкритим кодом. Зокрема, ILUPACK робить бінарні файли вільно доступними для вчених, але не видається, що вони випускають свій вихідний код. Чи справді так, …

15 pde  algorithms  linear-solver  libraries  preconditioning 

Я вирішую систему з двох зв'язаних PDE в двох просторових вимірах і в часі, обчислювально. Оскільки оцінки функцій дорогі, я хотів би використовувати багатоступінчастий метод (ініціалізований за допомогою Runge-Kutta 4-5). Метод Адамса-Башфорта, що використовує п'ять попередніх оцінок функції, має глобальну помилку (це випадок, коли у статті Вікіпедії, на яку посилається …

14 pde  algorithms  finite-element  error-estimation 

У мене є набір даних, що складається з мільйонів точок даних у 3D. Для обчислення, який я роблю, мені потрібно обчислити сусід (пошук діапазону) до кожної точки даних у радіусі, спробувати помістити функцію, обчислити помилку для придатності, повторити це для наступної точки даних тощо. Мій код працює належним чином, але …

14 algorithms  data-structures  nearest-neighbors  spatial-data 

При застосуванні класичної формули кута між двома векторами: α=arccosv1⋅v2∥v1∥∥v2∥α=arccos⁡v1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} виявляється, що при дуже малих / гострих кутах втрачається точність, і результат не є точним. Як пояснено у цій відповіді про переповнення стека , одним із варіантів є використання арктангенту: α=arctan2(∥v1×v2∥,v1⋅v2)α=arctan⁡2(‖v1×v2‖,v1⋅v2)\alpha = \arctan2 \left(\|\mathbf{v_1} \times …

14 algorithms  precision  numerical-limitations 

Багато алгоритмів, що використовуються в наукових обчисленнях, мають іншу притаманну структуру, ніж алгоритми, які зазвичай розглядаються в менш інтенсивних математичних формах програмної інженерії. Зокрема, окремі математичні алгоритми, як правило, дуже складні, часто включають сотні чи тисячі рядків коду, але, тим не менш, не містять стану (тобто не діють на складній …

13 algorithms 

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

12 linear-algebra  algorithms  matrix 

Чи існує алгоритм, який перераховує графіки, які відповідають деякому тестуванню точок Делоне в 3D? Якщо так, чи існує ефективна параметризація геометрій, яка відповідає будь-якому "графіку Делоне"? Я хочу систематично перерахувати всі стабільні геометрії молекул визначеного складу без будь-якого апріорного знання про зв'язування тощо. EDIT: Нехай - сукупність графіків з N …

12 algorithms  computational-chemistry  delaunay-triangulation 

Я намагаюся чисельно перевернути таке інтегральне перетворення: F(y) = ∫∞0удосвід[ - 12( у2+ х2) ] Я0( х у) f( х )d xF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Отже, для заданого мені потрібно наблизити f ( x ) де:Ж( у)F(y)F(y)f( х )f(x)f(x) і F ( y ) - реальні і …

11 algorithms  numerical-analysis  approximation  inverse-problem 

Тож теорема про декомпозицію Холеського стверджує, що будь-яка реальна симетрична позитивно-визначена матриця має розклад Холеського де L - нижня трикутна матриця.MMMM=LL⊤M=LL⊤M= LL^\topLLL З огляду на MMM , ми вже знаємо , є швидкі алгоритми для розрахунку його Чолеска фактора LLL . Тепер, припустимо, мені дали прямокутну матрицю A \ \ …

11 linear-algebra  algorithms 

Привіт спільноті scicomp, Я працював у галузі графічних алгоритмів, використовуючи рамки, такі як NetworkX (Python), JUNG та YFiles (Java). Зараз я входжу в область паралельних і високопродуктивних обчислень. Для нового проекту я шукаю бібліотеку графіків C ++ із такими можливостями: має інтуїтивний інтерфейс, який дозволяє розробляти алгоритми підтримує динамічні операції: …

11 algorithms  parallel-computing  libraries  graph-theory 

Я роблю діагоналізацію великої розрідженої матриці Ланцоса (~ 2 мільйони елементів). Майже всі етапи алгоритму Ланкоса виконуються паралельно на графічному процесорі, за винятком діагоналізації матриці Ланцоса для перевірки конвергенції. Для цього я використовував алгоритм TQLI з Numerical Recipes. Чи існують методи пошуку власної системи тридіагональної матриці, які є паралельними або …

11 linear-algebra  algorithms  eigensystem 

Це проста проблема перерахування симетрії. Я даю тут повний досвід, але знання квантової хімії не потрібні. Інтеграл двох частинок є: І має такі 4 симетрії: мене є функція, яка обчислює інтеграл і зберігає їх у 1D масиві , індексується так:⟨ я J | до л ⟩ = ∫ г | …

11 algorithms  performance 

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

11 algorithms  complexity  mpi 

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

11 algorithms  computational-geometry