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

У мене є декілька складних проблем, що стосуються невипуклої глобальної оптимізації. В даний час я використовую панель інструментів оптимізації MATLAB (конкретно, fmincon()з алгоритмом = 'sqp'), що є досить ефективним . Однак більша частина мого коду знаходиться в Python, і я хотів би зробити оптимізацію і в Python. Чи є розв'язувач …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

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

25 software  floating-point  accuracy 

Багато чисельних підходів до CFD можуть бути розширені до довільно високого порядку (наприклад, розривні методи Галеркіна, методи ВЕНО, спектральна диференціація тощо). Як слід вибрати відповідний порядок точності для даної проблеми?

23 fluid-dynamics  accuracy 

Мені відомо про те, що інвертування матриці для вирішення лінійної системи не є хорошою ідеєю, оскільки вона не настільки точна і настільки ефективна, як безпосередньо розв’язання системи або з використанням розкладання LU, Cholesky або QR. Однак я не зміг перевірити це на практичному прикладі. Я спробував цей код (у MATLAB) …

16 matrix  accuracy  inverse 

Більшість математичних бібліотек мають ряд версій функцій логарифму. Більшу частину часу ми вважаємо їх ідеальними, але насправді їх досить багато просто пропонують певну кількість цифр точності. Для деяких функцій існують чисельніші стабільні варіанти. Наприклад, у Fortran, R, Java і C є Math.log1pобчислення log(1.0+x)(що забезпечує більш високу точність для малих значень …

11 floating-point  accuracy 

Я намагаюся вирішити рівняння типу: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Якщо має простий полюс при 0 , для найменших N власних значень та власних векторів. Граничними умовами є: ψ ( 0 ) = 0 і ψ ( R ) = 0 , і я …

11 finite-difference  ode  accuracy 

Який сучасний спосіб здійснення спеціальних функцій подвійної точності? Мені потрібен такий інтеграл: для і , які можна записати через нижню неповну гамма-функцію. Ось моя реалізація Fortran та C: м=0,1,2,. . . t>0Fm(t)=∫10u2me−tu2du=γ(m+12,t)2tm+12Fm(t)=∫01u2me−tu2du=γ(m+12,t)2tm+12 F_m(t) = \int_0^1 u^{2m} e^{-tu^2} d u = {\gamma(m+{1\over 2}, t)\over 2 t^{m+{1\over 2}}} m=0,1,2,...m=0,1,2,...m=0, 1, 2, ...t>0t>0t>0 …

10 efficiency  accuracy  special-functions 

Чи можна отримати точність другого порядку для шестигранних кінцевих елементів із меншою кількістю 8 точок Гаусса без введення нефізичних режимів? Одинарна центральна точка Гаусса вводить нефізичний режим стрижки, а стандартне симетричне розташування 8 точок Гаусса дороге порівняно з тетраедральними розсудами. Редагувати : хтось просив рівнянь. Мене цікавлять рівняння - це …

10 finite-element  accuracy 

Я намагаюся інтегруватися ∫10т2 n + 2досвід(αr0т) dт∫01т2н+2досвід⁡(αr0т)гт\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt що є простою трансформацією ∫∞1х2 ндосвід( - αr0x ) dх∫1∞х2ндосвід⁡(-αr0х)гх\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx використовуючи тому що важко наблизити неправильні інтеграли. Це, однак, призводить до проблеми оцінювання нового інтегралу біля нуля. Буде дуже легко отримати належну кількість квадратурних вузлів, бачачи, що …

9 quadrature  accuracy