«quadrature» 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 

Aşağıdaki integrali sayısal olarak değerlendirmem gerekiyor: ∫∞0sinc′(xr)rE(r)−−−−√dr∫0∞sbennc'(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr burada , x \ in \ mathbb {R} _ + ve \ lambda, \ kappa, \ nu> 0 . Burada K , ikinci türün değiştirilmiş Bessel işlevidir. Benim özel durumumda \ lambda = 0.00313 , \ kappa = 0.00825 …

25 matlab  quadrature  special-functions 

Hem tek boyutlu hem de daha yüksek boyutlarda yüksek derecede salınımlı integrallerin keyfi kesinliğe yaklaşmasında son teknoloji nedir?

23 quadrature  precision  numerical-analysis 

Kareleme için en sayısal yöntemler integrali kara kutu işlevi olarak görür. Daha fazla bilgimiz varsa ne olur? Özellikle, eğer varsa, integralin ilk birkaç türevini bilmekten ne fayda sağlayabiliriz? Başka hangi bilgiler değerli olabilir? Özellikle türevler için: temel kareleme için hata tahminleri (dikdörtgen / trapzoid / simpson kuralları) yakından ilişkilidir. Belki …

19 quadrature  error-estimation 

Oktav'da iki vektörün Öklid mesafesini hesaplamanın hızlı bir yolu olup olmadığını bilmek istiyorum. Görünüşe göre bunun için özel bir fonksiyon yok, bu yüzden sadece formülü kullanmalı sqrtmıyım?

18 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Anekdotsal olarak duydum ki, kişi sayısal olarak formun bir integralini yapmaya çalışırken ∫∞0f( x ) J0( x )g x∫0∞f(x)J0(x)dx\int_0^\infty f(x) J_0(x)\,\mathrm{d}x ile (vs, örneğin yüksek titreşimli tekil olmayan kendi başına) düz ve iyi huylu, o zaman doğruluğu olarak yeniden için yardımcı olacaktırf( x )f(x)f(x) 1π∫π0∫∞0f( x ) cos( x günahθ …

15 quadrature  special-functions 

Matlab'ın sayısal entegrasyonu Scipy ile karşılaştırması konusunda biraz hayal kırıklığı yaşıyorum. Aşağıdaki test kodumda aşağıdaki farklılıkları gözlemliyorum: Matlab'ın versiyonu python eşdeğümden ortalama 24 kat daha hızlı çalışıyor ! Matlab'ın sürümü, integrali uyarmadan hesaplayabilirken, python döner nan+nanj Bahsedilen iki noktaya göre python'da aynı performansı almamı sağlamak için ne yapabilirim? Dokümantasyona göre, …

13 matlab  python  quadrature  numba 

Bir fonksiyonu olan böyle sonludur ve bu integrali yaklaşık istiyoruz. ∫ R 3 f ( x , y , z ) d Vf( x , y, z)f(x,y,z)f(x,y,z) ∫R,3f(x , y, z) dV∫R3f(x,y,z)dV\int_{R^3} f(x,y,z)dV Kareleme kurallarına ve integrallerin monte carlo yaklaşımlarına aşinayım, ancak bunları sonsuz bir alanda uygulamakta bazı zorluklar görüyorum. …

13 monte-carlo  quadrature 

Bu uzun sorunun sadece ilk iki bölümü önemlidir. Diğerleri sadece örnekleme amaçlıdır. Arka fon Yüksek dereceli kompozit Newton – Cotes, Gauß – Legendre ve Romberg gibi gelişmiş kareleme, esas olarak, fonksiyonun iyi bir şekilde örneklenebileceği ancak analitik olarak entegre edilemediği durumlar için tasarlanmıştır. Bununla birlikte, örnekleme aralığından (örnek için Ek …

12 reference-request  quadrature 

Sayısal kareleme için birkaç yöntem ailesi mevcuttur. Belirli bir integrand sınıfım varsa ideal yöntemi nasıl seçerim? Hem integral (örneğin, pürüzsüz mü? Bu soruların cevapları çeşitli yöntem ailelerini nasıl dışlar veya teşvik eder? Basitlik için sadece tek veya düşük boyutlu integralleri ele alalım. Örneğin Quadpack üzerinde Wikipedia makalesi oldukça genel belirten …

12 quadrature 

Ben sıfıra yakın belirli değerlerde NaNs gösteren zor bir integral ile ilgileniyorum ve şu anda onlarla oldukça kabaca bu olduğunda integrand sıfıra ayarlar bir ISNAN deyimi kullanarak ilgileniyorum. Bu FORTRAN (q1da rutin - q1dax farklı değil) NMS kitaplığı ve C (QAGS rutin kullanarak) C GSL kitaplığı ile denedim. Özellikle integral …

12 quadrature 

Ben 3D bir 4-düğüm eleman üzerine bir polinom ifade entegre etmek istiyorum. FEA üzerine birkaç kitap, entegrasyonun keyfi düz 4-elemansız bir eleman üzerinden gerçekleştirildiği durumu kapsamaktadır. Bu durumda olağan prosedür, Jacobi matrisini bulmak ve entegrasyon temelini, daha basit entegrasyon limitlerine [-1; 1] sahip olduğum ve Gauss-Legendre kareleme tekniğinin kolayca kullanıldığı …

12 finite-element  quadrature 

(2 boyutlu) uygun olmayan integralim var ben= ∫birW( x , y)F( x , y)d x d yI=∫AW(x,y)F(x,y)dxdyI=\int_A \frac{W(x,y)}{F(x,y)}\,\mbox{d}x\mbox{d}y burada entegrasyonu alanı x = [ - 1 , 1 ] , y = [ - 1 , 1 ] ' den daha küçüktür ancak F ( x , y ) > …

12 monte-carlo  quadrature 

Gelen karmaşık fonksiyon teorisi uygulamaları bu gelişmiş ders bir egzersiz salınım yapan ayrılmaz bir noktasında ben( λ ) = ∫∞- ∞marul( λ cosx ) günahxxdxI(λ)=∫−∞∞cos⁡(λcos⁡x)sin⁡xxdxI(\lambda)=\int_{-\infty}^{\infty} \cos (\lambda \cos x) \frac{\sin x}{x} d x karmaşık düzlemde eyer noktası yöntemi kullanılarak büyük değerleri için yaklaşık olarak hesaplanmalıdır .λλ\lambda Oldukça salınımlı yapısı nedeniyle, …

11 quadrature  integration  oscillations 

Ben bir işlev var ki ı entegre etmek istediklerini üzerinden bir dört yüzlü T ⊂ R 3 . Eğer f keyfi idi Gauss nümerik iyi bir çözüm olurdu ama bunu biliyor f harmonik olduğunu. Gauss kareleme bu bilgiler kullanılarak ne kadar hızlandırılabilir?f: R3→ Rf:R,3→R,f : \mathbf{R}^3 \to \mathbf{R}T⊂ R3T⊂R,3T \subset …

11 quadrature