«floating-point» 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 

Google’a gerçekten zor olan basit bir sorum var (kanuni yanı sıra, Her Bir Bilgisayar Bilim Adamının Kayan Nokta Aritmetik makalesi hakkında bilmesi gerekenler ). Gibi fonksiyonlar ne zaman log1pya da expm1yerine kullanılabilir logve exp? Ne zaman kullanılmamalıdır? Bu işlevlerin farklı uygulamaları kullanımları açısından nasıl farklılık gösterir?

30 floating-point 

Gerçek değişkenlerin gerçek bir işlevi göz önüne alındığında, IEEE 754 aritmetiği ile donatılmış bir makinedeki tüm girdiler üzerindeki işlevi hesaplamak için otomatik olarak sayısal olarak doğru kod üretebilen bir yazılım var mı? Örneğin, değerlendirilecek asıl işlev şuysa: Yazılım, hesaplamalı doğruluk kaybını önlemek için felaket iptalini ve bazı girdi kümeleri için …

25 software  floating-point  accuracy 

Markov Zincirinde çok yaygın bir sorun olan Monte Carlo, büyük üstel terimlerin toplamı olan hesaplama olasılıklarını içermektedir. ebir1+ ebir2+ . . .ea1+ea2+... e^{a_1} + e^{a_2} + ... Bir kutudaki bileşenler çok küçükten çok geniş yelpazeye kadar değişebilir. Benim yaklaşımım, en büyük üstel terimi etmektir, böylece:biraaK: = maksben( aben)K:=maxi(ai)K := \max_{i}(a_{i}) …

24 monte-carlo  floating-point  statistics 

Pozitif sayılar artan bir sonlu dizisi Verilen . Sayıların toplamını hesaplamak için aşağıdaki iki algoritmadan hangisi daha iyidir?z1, z2, .. . . . znz1,z2,.....znz_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end Veya: s=0; for \ i=1:n s=s + z_{n-i+1} ; end Bence sayıları en büyükten en küçük …

22 algorithms  floating-point 

Giriş olarak birkaç kayan nokta değeri (tek veya çift) alan, bazı hesaplama yapan ve çıktı kayan nokta değerleri (ayrıca tek veya çift) üreten bir işleve sahip olduğumu varsayalım. Öncelikle MSVC 2008 ile çalışıyorum, aynı zamanda MinGW / GCC ile çalışmayı da planlıyorum. C ++ ile programlama yapıyorum. Sonuçlarda ne kadar …

20 floating-point  c++ 

Düşük göreli hata ile çift hassasiyetli kayan nokta aşağıdaki işlevi uygulamaya çalışıyorum : l o g s u m (x,y) = günlük( exp( x ) + tecrübe( y) )lÖgsum(x,y)=günlük⁡(tecrübe⁡(x)+tecrübe⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Bu, günlük uygulamalarda temsil edilen olasılıkları veya olasılık yoğunluklarını eklemek için istatistiksel uygulamalarda yaygın olarak kullanılır. Tabii …

18 floating-point  stability  numerics 

İki farklı BLAS uygulaması göz önüne alındığında, tam olarak aynı kayan nokta hesaplamaları yapmasını ve aynı sonuçları döndürmesini bekleyebilir miyiz? Veya örneğin skaler bir ürünü ve bir ( x 1 y 1 + x 2 y) şeklinde hesaplayabilir mi? 2 ) + ( x 3 y 3 + x 4((x1y1+x2y2)+x3y3)+x4y4((x1y1+x2y2)+x3y3)+x4y4 …

17 floating-point  blas 

GROMACS moleküler dinamiği (MD) kodunu çalıştırdım24 Intel Xeon CPU içeren düğümlerden oluşan bir Ubuntu Linux kümesinde . Özel ilgi alanım kayan nokta aritmetik hassasiyetine biraz duyarlı olduğu için, çift hassasiyetin daha yüksek hesaplama maliyetine rağmen GROMACS'ı tek hassasiyet yerine çift hassasiyetle çalıştırmak zorunda kaldım. Küme üzerinde, GROMACS'ı çift hassasiyetle derledim. …

17 hpc  floating-point  molecular-dynamics  gpu 

Bilgisayarımın yapabileceği FLOP'ların (Kayan Nokta İşlemleri) teorik sayısını belirlemek istiyorum. Birisi bana bu konuda yardımcı olabilir mi? (Aralarındaki fark hakkında fikir edinmek için bilgisayarımı bazı süper bilgisayarlarla karşılaştırmak istiyorum)

15 floating-point 

Izin vermek x, yiki kayan nokta sayı olmak. Ortalamalarını hesaplamanın doğru yolu nedir? Saf yol , çok büyük ve (x+y)/2aşırı olduğunda taşmalara neden olabilir . Belki daha iyi olduğunu düşünüyorum , ama iki çarpma (belki de verimsiz) içerir ve yeterince iyi olup olmadığından emin değilim. Daha iyi bir yol var …

15 floating-point 

Bir işlevi sayısal olarak uygulamakta sorun yaşıyorum. Büyük giriş değerlerinde sonucun çok büyük bir sayının çok küçük bir sayının çarpımı olduğu gerçeğinden muzdariptir. Felaket iptalinin doğru terim olup olmadığından emin değilim, lütfen öyleyse düzeltin. Bir şeyin ters gittiğine dair kanıtlar: 6'nın daha büyük girdileri için salınımları ve 0.0 değerini atamayı …

13 python  numerical-analysis  floating-point 

Allen ve Tildesley'nin Sıvılar Bilgisayar Simülasyonu kitabını okuyorum. 71. sayfadan itibaren yazarlar, Newton'un hareket denklemlerini moleküler dinamik (MD) simülasyonlarına entegre etmek için kullanılan çeşitli algoritmaları tartışıyorlar. 78. sayfadan itibaren, yazarlar belki de MD'deki kanonik entegrasyon algoritması olan Verlet algoritmasını tartışıyorlar. Belirtiyorlar: Belki de hareket denklemlerini entegre etmek için en yaygın …

13 floating-point  computer-arithmetic 

Basit bir fonksiyondaki FLOP sayısını değerlendirirken, genellikle temel aritmetik operatörleri ifade eden ifadeden aşağı inilebilir. Bununla birlikte, bölünmeyi bile içeren matematiksel ifadeler söz konusu olduğunda, bunu yapamaz ve sadece toplama ve çarpma ile fonksiyonlardan FLOP sayımları ile karşılaştırmayı bekleyebiliriz. İşlem bir kütüphanede uygulandığında durum daha da kötüdür. Bu nedenle, özel …

13 performance  complexity  floating-point 

Çoğu matematik kütüphanesinde logaritma fonksiyonlarının bazı versiyonları vardır. Çoğu zaman onların mükemmel olduklarını varsayıyoruz, ancak aslında birçoğu sadece belirli sayıda hassas basamak sunuyor. Bazı işlevler için sayısal olarak daha kararlı varyantlar vardır. Örneğin, Fortran, R, Java ve C'nin her ikisi de Math.log1phesaplama için log(1.0+x)(küçük x değerleri için daha yüksek hassasiyet …

11 floating-point  accuracy