Pertanyaan yang diberi tag «ode»

Saya memiliki beberapa masalah optimisasi global non-cembung yang menantang untuk dipecahkan. Saat ini saya menggunakan MATLAB's Optimization Toolbox (khusus, fmincon()dengan algoritma = 'sqp'), yang cukup efektif . Namun, sebagian besar kode saya menggunakan Python, dan saya ingin melakukan optimasi dengan Python juga. Apakah ada pemecah NLP dengan binding Python yang …

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 

Saat ini saya sedang mencari metode paralel untuk integrasi ODE. Ada banyak literatur baru dan lama di luar sana yang menggambarkan berbagai pendekatan, tetapi saya belum menemukan survei terbaru atau artikel ikhtisar yang menggambarkan topik secara umum. Ada buku karya Burrage [1], tetapi sudah hampir 20 tahun dan karenanya tidak …

39 parallel-computing  ode  numerical-analysis 

Pertanyaan ini dimigrasikan dari Mathematics Stack Exchange karena dapat dijawab di Computational Science Stack Exchange. Bermigrasi 6 tahun yang lalu . Saya akan mulai bekerja di perpustakaan perangkat lunak pemecah ODE numerik , dan saya berjuang dengan cara merumuskan tes untuk implementasi pemecah. Ambisi saya adalah bahwa perpustakaan, akhirnya, akan …

26 ode 

Bagaimana masalah gravitasi n-tubuh diselesaikan secara numerik secara paralel? Apakah tradeoff presisi-kompleksitas mungkin? Bagaimana presisi mempengaruhi kualitas model?

25 algorithms  ode  parallel-computing  complexity  precision 

Dalam komentar ini saya menulis: ... integrator SciPy default, yang saya asumsikan hanya menggunakan metode symplectic. di mana saya merujuk ke SciPy's odeint, yang menggunakan metode "non-stiff (Adams)" atau "stiff (BDF) method". Menurut sumber : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, …

25 ode  time-integration  integration  differential-equations 

Saya hanya ingin tahu mengapa metode tingkat tinggi (yaitu lebih dari 4) metode Runge-Kutta hampir tidak pernah dibahas / dipekerjakan (setidaknya setahu saya). Saya mengerti ini membutuhkan waktu komputasi yang lebih besar per langkah (misalnya RK14 dengan langkah tertanam urutan ke-12 ), tetapi apakah ada kelemahan lain dari menggunakan metode …

17 ode  runge-kutta 

Pertimbangkan IVP untuk sistem ODE , . Paling umum masalah ini dianggap kaku ketika matriks Jacobi memiliki kedua nilai eigen dengan bagian nyata negatif yang sangat besar dan nilai eigen dengan bagian nyata negatif yang sangat kecil (saya menganggap hanya bagian stabil kasus).y ( x 0 ) = y 0 …

17 ode  stiffness 

Saya ingin mensimulasikan perilaku sistem seperti pendulum ganda. Sistem ini adalah manipulator robot 2 derajat kebebasan yang tidak digerakkan dan oleh karena itu, akan berperilaku seperti pendulum ganda yang dipengaruhi oleh gravitasi. Satu-satunya perbedaan utama dengan pendulum ganda adalah bahwa ia terdiri dari dua benda kaku dengan massa dan sifat …

15 matlab  ode  numerical-modelling 

apa metode canggih untuk solusi numerik ODE dengan sisi kanan terputus-putus? Saya sebagian besar tertarik sedikit demi sedikit fungsi sisi kanan, misalnya tanda. Saya mencoba memecahkan persamaan dari tipe berikut: x˙v˙=v={(|Fexternal|−|Ffriction|)sign(Fexternal)0:|Fexternal|&lt;|Ffriction|:otherwisex˙=vv˙={(|Fexternal|−|Fgesekan|)tanda⁡(Fluar):|Fluar|&lt;|Fgesekan|0:jika tidak\begin{align*} \dot x &= v\\ \dot v &= \begin{cases} (|F_\text{external}| - |F_\text{friction}|) \mathop{\rm sign} (F_\text{external}) & :|F_\text{external}| < |F_\text{friction}|\\ …

15 ode 

Saya ingin menerapkan sistem penyelesaian ODE ke GPU, dalam pengaturan yang 'sepele parallelisable'. Misalnya, melakukan analisis sensitivitas dengan 512 set parameter yang berbeda. Idealnya saya ingin melakukan penyelesaian ODE dengan solver pencatat waktu adaptif yang cerdas seperti CVODE, daripada pencatat waktu yang tetap seperti Forward Euler, tetapi menjalankannya pada GPU …

15 ode  gpu 

Dalam praktiknya, runtime penyelesaian IVP ˙ x ( t ) = f ( t , x ( t ) ) secara numerik x˙(t)=f(t,x(t)) for t∈[t0,t1]x˙(t)=f(t,x(t)) for t∈[t0,t1] \dot{x}(t) = f(t, x(t)) \quad \text{ for } t \in [t_0, t_1] x(t0)=x0x(t0)=x0 x(t_0) = x_0 sering didominasi oleh lamanya mengevaluasi sisi kanan …

14 ode  numerical-analysis 

Saya bekerja cukup banyak dalam biologi matematika / epidemiologi, di mana sebagian besar karya ilmu pemodelan / komputasi masih didominasi oleh set ODE, diakui kadang-kadang set cukup rumit dari mereka. Salah satu kelebihan model-model ini adalah model-model itu agak mudah digambarkan dan ditiru. Tabel nilai parameter, dan persamaan itu sendiri …

14 ode  computational-biology 

Hai dan terima kasih telah meluangkan waktu untuk melihat pertanyaan saya. Ini adalah versi terbaru dari pertanyaan saya yang saya posting sebelumnya di physics.stackexchange.com Saat ini saya sedang mempelajari spin exciton 2D Bose-Einstein Condensate dan saya ingin tahu tentang keadaan dasar sistem ini. Metode matematika untuk mencapai keadaan dasar disebut …

12 matlab  ode  quantum-mechanics 

Saya bertanya-tanya: apa algoritma terbaik untuk menyelesaikan Di mana adalah matriks . A tidak bergantung pada waktu secara eksplisit, biasanya jarang tetapi tidak perlu terikat. Nilai eigennya memiliki bagian nyata yang tidak positif. A juga dapat didiagonalisasi tetapi mungkin terlalu besar untuk diagonalisasi penuh agar efisien secara komputasi.dkamudt= A ududt=Au\begin{equation} …

12 linear-algebra  ode 

Saya telah menggunakan GSL sebagai dasar dari banyak simulasi saya, tetapi itu sedikit berlebihan untuk tujuan saya dan mendefinisikan tipe kompleksnya sendiri untuk alasan warisan. Daripada mengkode pemecah ODE Runge-Kutta saya sendiri, yang mungkin tidak akan sangat efisien, apakah ada pemecah ODE open source yang menggunakan tipe kompleks C99 asli?

12 c  ode