«ode» etiketlenmiş sorular

Çözülmesi gereken birkaç dışbükey olmayan küresel optimizasyon problemim var. Şu anda , oldukça etkili olan MATLAB's Optimizasyon Araç Kutusu'nu (özellikle, fmincon()algoritma = ile 'sqp') kullanıyorum . Ancak, kodumun çoğu Python'da ve optimizasyonu Python'da da yapmak isterim. Rekabet edebilecek Python bağlarına sahip bir NLP çözücü var mı fmincon()? O olmalı Doğrusal …

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 

Şu anda ODE entegrasyonu için paralel yöntemler arıyorum. Geniş bir yaklaşım yelpazesini tanımlayan pek çok yeni ve eski literatür var, ancak konuyu genel olarak anlatan herhangi bir son anket veya genel bakış makalesi bulamadım. Burrage [1] adlı kitabın kitabı var, ancak neredeyse 20 yaşında ve bu nedenle parareal algoritması gibi …

39 parallel-computing  ode  numerical-analysis 

Bu soru, Matematiksel Yığın Değişim Borsası'ndan taşınmıştır , çünkü Hesaplamalı Bilim Yığın Değişimi alanında yanıtlanabilir. 6 yıl önce göç etmiş . Sayısal ODE çözücülerinin bir yazılım kütüphanesinde çalışmaya başlamak üzereyim ve çözücü uygulamaları için testlerin nasıl formüle edileceği ile mücadele ediyorum. Benim hırs kütüphane, sonuçta, hem nonstiff ve sert problemleri …

26 ode 

Yerçekimi n-vücut problemi sayısal olarak paralel olarak nasıl çözülebilir? Hassas-karmaşıklık tradeoff mümkün mü? Hassasiyet modelin kalitesini nasıl etkiler?

25 algorithms  ode  parallel-computing  complexity  precision 

Bu yazımda şöyle yazdım: ... varsayılan olarak sadece sempatik yöntemler kullandığını varsaydığım SciPy entegratörü. İçinde odeint"sert olmayan (Adams) bir yöntem" veya "sert (BDF) bir yöntem" kullanan SciPy'ye atıfta bulunuyorum . Kaynağa göre : 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, hmin=0.0, ixpr=0, …

25 ode  time-integration  integration  differential-equations 

ODE sistemi için bir IVP düşünün y′=f(x,y)y′=f(x,y)y'=f(x,y) , y(x0)=y0y(x0)=y0y(x_0)=y_0 . Bu iş en yaygın sorun olarak kabul edilir katı zaman Jacobi matris ∂f∂y( x0, y0)∂f∂y(x0,y0)\frac{\partial f}{\partial y}(x_0,y_0)sahipher ikiçok küçük negatif gerçek kısmı ile çok büyük negatif gerçek parçanın ve özdeğerler (sadece sabit bir durumda dikkate) ile öz. Öte yandan, sadece …

18 ode  stiffness 

Sadece yüksek dereceli (yani 4'ten büyük) Runge – Kutta yöntemlerinin neden neredeyse hiç tartışılmadığını (en azından bilgime göre) kullanıldığını merak ettim. Adım başına daha büyük hesaplama süresi gerektirdiğini anlıyorum (örneğin 12. sıradaki gömülü adımla RK14 ), ancak daha yüksek dereceli Runge-Kutta yöntemlerini kullanmanın başka olumsuzlukları var mı (örn. Kararlılık sorunları)? …

17 ode  runge-kutta 

ODE'lerin sistemlerini GPU'lara 'çok paralellikle paralel' bir ortamda çözmek istiyorum. Örneğin, 512 farklı parametre seti ile duyarlılık analizi yapmak. İdeal olarak, ODE çözmeyi Forward Euler gibi sabit bir zaman çizelgesi yerine CVODE gibi akıllı bir uyarlanabilir zaman çizelgesi çözücüsü ile yapmak istiyorum, ancak CPU yerine bir NVIDIA GPU üzerinde çalıştırıyorum. …

16 ode  gpu 

Çift sarkaçlı bir sistemin davranışını simüle etmek istiyorum. Sistem, harekete geçirilmeyen ve bu nedenle çoğunlukla yerçekiminden etkilenen bir çift sarkaç gibi davranacak olan 2 serbestlik dereceli bir robot manipülatördür. Çift sarkaçlı tek temel fark, kütle merkezlerinde kütle ve atalet özelliklerine sahip iki katı cisimden oluşmasıdır. Temel olarak, ode45Matlab altında aşağıdaki …

15 matlab  ode  numerical-modelling 

süreksiz sağ tarafı olan ODE'lerin sayısal çözümü için son teknoloji yöntemler nelerdir? Çoğunlukla parçalı pürüzsüz sağ taraf fonksiyonlarıyla ilgileniyorum, örneğin işaret. Aşağıdaki türdeki denklemi çözmeye çalışıyorum: x˙v˙=v={(|Fexternal|−|Ffriction|)sign(Fexternal)0:|Fexternal|&lt;|Ffriction|:otherwisex˙=vv˙={(|Fexternal|−|Ffriction|)sign⁡(Fexternal):|Fexternal|&lt;|Ffriction|0:otherwise\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}|\\ 0 & : \text{otherwise} \end{cases} …

15 ode 

Uygulamada, bir IVP ˙ x ( t ) = f ( t , x ( t ) ) sayısal olarak çözülme süresi x ( t 0 ) = x 0 genellikle sağ taraftaki (ST) değerlendirme süresi hakimdir f . Bu nedenle, diğer tüm işlemlerin anında (yani hesaplama maliyeti olmadan) olduğunu …

14 ode  numerical-analysis 

Modelleme / hesaplama bilimi çalışmalarının çoğuna hâlâ ODE'lerin hâkim olduğu, bazen de oldukça ayrıntılı kümelerin hâkim olduğu matematiksel biyoloji / epidemiyolojide oldukça çalışıyorum. Bu modellerin artılarından biri, tanımlanması ve çoğaltılması oldukça kolay olmasıdır. Bir parametre değerleri tablosu ve denklemlerin kendileri ve birisine araştırmanızı uygulamak istedikleri şekilde çoğaltmak için ihtiyaç duydukları …

14 ode  computational-biology 

Selam ve soruma bakmak için zaman ayırdığınız için teşekkürler. Bu, daha önce physics.stackexchange.com'da gönderdiğim sorumun güncellenmiş bir versiyonudur. Şu anda 2D eksiton spinor Bose-Einstein Kondensini inceliyorum ve bu sistemin temel durumunu merak ediyorum. Temel duruma ulaşmanın matematiksel yöntemine hayali zaman yöntemi denir . Kuantum mekaniğinde zamanın hayali bir ile değiştirildiği …

12 matlab  ode  quantum-mechanics 

Acaba: ı çözmek için en iyi algoritma nedir Burada gerçek matrisidir. A, açıkça zamana bağlı değildir, genellikle seyrek değildir, ancak zorunlu olarak bantlanmış değildir. Özdeğerlerinin pozitif olmayan gerçek kısımları vardır. A da köşegenleştirilebilir ancak tam bir köşegenleştirmenin hesaplama açısından verimli olması için çok büyük olabilir.dudt= A ududt=Au\begin{equation} \frac{du}{dt} = Au …

12 linear-algebra  ode 

GSL'yi simülasyonlarımın çoğunun temeli olarak kullanıyorum , ancak benim amacım için biraz abartılı ve eski nedenlerden dolayı kendi karmaşık türünü tanımlıyor. Muhtemelen çok verimli olmayacak olan kendi Runge-Kutta ODE çözücümü kodlamak yerine, yerel C99 karmaşık tipini kullanan herhangi bir açık kaynak ODE çözücüsü var mı?

12 c  ode