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

Яка найшвидша бібліотека для виконання триангуляції множин з мільйонами, якщо 3D-очки? Чи доступні також версії GPU? З іншого боку, маючи вороной тесселяцію одного і того ж набору точок, допомогло б (з точки зору продуктивності) отримати зворотну тріангуляцію?

26 computational-geometry  delaunay-triangulation  voronoi-diagrams 

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

26 data-storage 

Який простий алгоритм для обчислення SVD матриць ?2×22×22 \times 2 В ідеалі мені хотілося б чисельно надійний алгоритм, але я хотів би бачити як прості, так і не дуже прості реалізації. C код прийнято. Будь-які посилання на папери чи код?

26 linear-algebra  matrix  reference-request  svd  matrix-factorization 

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

26 software 

Що може піти не так, якщо використовувати попередньо налаштовані методи Крилова з KSP ( пакет лінійних розв'язків PETSc ) для розв’язання розрідженої лінійної системи, такої, яку отримують шляхом дискретизації та лінеаризації часткових диференціальних рівнянь? Які кроки я можу зробити, щоб визначити, що не так у моїй проблемі? Які зміни я …

26 linear-algebra  petsc  preconditioning  krylov-method  iterative-method 

Мені потрібно чисельно оцінити інтеграл нижче: ∫∞0s i n c'( x r ) r E( r )----√гr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr де , і . Тут - модифікована функція Бесселя другого роду. У моєму конкретному випадку я , і .x∈R+λ,κ,ν>0Kλ=0,00313κ=0,00825ν=0,33E(r)=r4(λκ2+r2−−−−−−√)−ν−5/2K−ν−5/2(λκ2+r2−−−−−−√)E(r)=r4(λκ2+r2)−ν−5/2K−ν−5/2(λκ2+r2)E(r) = r^4 (\lambda\sqrt{\kappa^2+r^2})^{-\nu-5/2} K_{-\nu-5/2}(\lambda\sqrt{\kappa^2+r^2})x∈R+x∈R+x \in \mathbb{R}_+λ,κ,ν>0λ,κ,ν>0\lambda, \kappa, \nu >0KKKλ = …

25 matlab  quadrature  special-functions 

Пов'язане питання: Стан Mac OS у науковому обчисленні та HPC Значна кількість програмних пакетів з обчислювальної науки написана у Fortran, і Fortran не йде. Компілятор Fortran також необхідний для створення інших програмних пакетів (одним із помітних прикладів є SciPy ). Однак Mac OS X не включає компілятор Fortran. Як слід …

25 software  fortran 

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

25 optimization  conjugate-gradient 

Враховуючи реальну функцію реальних змінних, чи доступне програмне забезпечення, яке може автоматично генерувати числовий точний код для обчислення функції на всіх входах на машині, оснащеній арифметикою IEEE 754? Наприклад, якщо реальна функція, яку слід оцінювати, була: Програмне забезпечення розглядає катастрофічне скасування та, можливо, пошук таблиць виходів для певних наборів входів, …

25 software  floating-point  accuracy 

Як можна гравітаційну задачу n-тіла вирішити чисельно паралельно? Чи можливий компроміс точності та складності? Як точність впливає на якість моделі?

25 algorithms  ode  parallel-computing  complexity  precision 

Я шукаю інструмент для візуалізації дуже великих графіків спрямованих посилань. На даний момент у мене є ~ 2 млн. Вузлів з ~ 10 млн. Країв. Я спробував декілька різних речей, але більшості потрібні години, щоб зробити навіть графіки на 100k вузлів Що я спробував: я провів день з gephi, але …

25 visualization 

У цьому коментарі я написав: ... інтегратор SciPy за замовчуванням, для якого я припускаю, що він використовує лише симплектичні методи. в якому я замислююся про SciPy odeint, який використовує або "не жорсткий (метод Адамса)", або "жорсткий (BDF) метод". За даними джерела : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, …

25 ode  time-integration  integration  differential-equations 

Я не розумію різної поведінки рівняння адвекції-дифузії, коли я застосовую різні граничні умови. Моя мотивація - це моделювання реальної фізичної величини (щільності частинок) при дифузії та адвекції. Щільність частинок повинна зберігатися у внутрішніх приміщеннях, якщо вона не витікає з країв. За цією логікою, якщо я виконую умови кордону Неймана, кінці …

25 pde  boundary-conditions  advection  diffusion  crank-nicolson 

Взагалі кажучи, я чув, як чисельні аналітики висловлюють цю думку «Звичайно, з математичної точки зору, час це просто інший вимір, але все ж, час є особливим» Як це виправдати? У якому сенсі час особливий для обчислювальної науки? Більше того, чому ми так часто вважаємо за краще використовувати кінцеві відмінності (що …

24 pde  discretization 

Припустимо, я працюю над науковим кодом на C ++. У недавній дискусії з колегою було стверджено, що шаблони виразів можуть бути дійсно поганою справою, що потенційно робить програмне забезпечення компільованим лише для певних версій gcc. Нібито ця проблема торкнулася декількох наукових кодексів, про які йдеться у підзаголовках цієї пародії на …

24 performance  c++  automatic-differentiation  interval-arithmetic