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

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 

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 

Burada kimse Python aracılığıyla yeni nesil (örn. K20) GPU'larla çift duyarlıklı bilimsel hesaplama kullandı mı? Bu teknolojinin hızla geliştiğini biliyorum, ancak şu anda bunu yapmanın en iyi yolu nedir? GPU, popüler bilimsel Python kütüphaneleri için numpy ve scipy kapsamı dışında ve ben theano kullanmak istedim ama GPU için sadece float32 …

14 python  gpu  precision 

İki vektör arasındaki açı için klasik formülü uygularken: α=arccosv1⋅v2∥v1∥∥v2∥α=arccos⁡v1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} çok küçük / akut açılar için hassasiyet kaybı olduğu ve sonucun doğru olmadığı bulunmuştur. Bu Yığın Taşması cevabında açıklandığı gibi , çözümlerden biri arktanjant kullanmaktır: α=arctan2(∥v1×v2∥,v1⋅v2)α=arctan⁡2(‖v1×v2‖,v1⋅v2)\alpha = \arctan2 \left(\|\mathbf{v_1} \times \mathbf{v_2}\|, \mathbf{v_1} \cdot \mathbf{v_2} \right) …

14 algorithms  precision  numerical-limitations 

Buradaki cevaba göre , büyük koşul sayısı (doğrusal sistem çözümü için) kayan nokta çözümünde garantili doğru basamak sayısını azaltır. Psödospektral yöntemlerde daha yüksek mertebeden farklılaşma matrisleri tipik olarak çok koşulsuzdur. Öyleyse neden hala çok doğru yöntemler oldukları? Koşulsuz matrislerden gelen düşük hassasiyetin sadece garantili bir değer olduğunu anlıyorum, ancak yine …

13 spectral-method  precision  condition-number 

Tek hassasiyetli kayan nokta sayıları hafızanın yarısını kaplar ve modern makinelerde (GPU'larda bile) çift hassasiyetle karşılaştırıldığında neredeyse onlarla iki kat daha hızlı işlem yapılabilir. Sadece bulduğum birçok FDTD kodu tek duyarlıklı aritmetik ve depolama kullanır. Büyük ölçekli seyrek denklem sistemlerini çözmek için tek bir hassasiyetin ne zaman kabul edilebilir olduğuna …

13 matrix  linear-solver  precision 

Sayısal türevleri hesaplamak istendiğinde, Bengt Fornberg tarafından burada (ve burada rapor edilir ) sunulan yöntem çok uygundur (hem kesin hem de uygulanması basittir). 1988'deki orijinal bildiri tarihi olarak, bugün daha iyi bir alternatif olup olmadığını bilmek istiyorum (ya da neredeyse basit ve daha kesin)?

11 finite-difference  precision 

C ile yazılmış büyük bir model (~ 5000 satır) var. Hiçbir yerde rastgele sayı üretimi ile seri bir programdır. FFTW işlevlerini FFTW kullanan işlevler için kullanır - FFTW uygulamasının ayrıntılarını bilmiyorum, ancak buradaki işlevlerin de deterministik olduğunu varsayıyorum (hata yapıyorsam beni düzelt). Anlayamadığım sorun, aynı makinede (aynı derleyici, aynı kütüphaneler) …

10 floating-point  precision  computer-arithmetic 

Bazı yoğun, koşulsuz matrisleri köşegenleştirmeye çalışıyorum. Makine hassasiyetinde sonuçlar yanlıştır (negatif özdeğerleri döndürmek, özvektörlerde beklenen simetriler yoktur). Rasgele hassasiyetten yararlanmak için Mathematica'nın Eigensystem [] işlevine geçtim, ancak hesaplamalar son derece yavaş. Ben çok sayıda çözüme açığım. Koşulsuz sorunlara uygun paketler / algoritmalar var mı? Ön koşullandırma konusunda uzman değilim, bu …

10 linear-algebra  parallel-computing  eigensystem  precision  dense-matrix 

Ben ikinci dereceden alanları kullanmanız gereken bir proje var Özellikle ile formun sayıları .a+b−3−−−√a+b−3a + b \sqrt{-3}a,b∈Qa,b∈Qa,b \in \mathbb{Q} Örneğin, Eisenstein tamsayılarındaki asal sayılar şunlardır : Adaçayı kullanmak istemiyorum. Dahil etmek için kendi veri türümü yazmak istiyorum numpy. PARI faydalı olacaktır - ancak Python ile uyumlu değildir. Bu nesneler için …

9 linear-algebra  precision 

Çok farklı kaynaklardan ve araştırmalarla yapılan görüşmelerden, sayısal kısmi diferansiyel denklemlerde yüksek hassasiyetli hesaplamalar için artan bir talep olduğu izlenimine sahibim. Burada yüksek hassasiyet, standart 64bit çift hassasiyetten daha fazla hassasiyet anlamına gelir. Bu konunun son teknolojisini merak ediyorum. Karşılaştırma yoluyla, sayısal PDE'de, özellikle çok çekirdekli yöntemler, büyük ölçekli paralelleştirme …

9 reference-request  precision