Hesaplamalı Bilim

Söyleyebileceğim kadarıyla, iki büyük jenerik ABD Enerji Bakanlığı hesaplamalı bilim yazılımı çerçevesi PETSc ve Trilinos'tur . İlk bakışta, dil farklılıklarının ötesinde benzer görünüyorlar (C'ye karşı C ++). İki çerçeve arasındaki temel farklar nelerdir ve birini seçmenin diğerine göre hangisini etkilemesi gerekir? (Kurumsal önyargıları ve mevcut altyapıyı görmezden gelin.)

24 petsc  software  trilinos 

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 

C ++ dilinde bilimsel bir kod üzerinde çalıştığımı varsayalım. Bir meslektaşı ile yakın zamanda yapılan bir tartışmada, ifade şablonlarının gerçekten kötü bir şey olabileceği, potansiyel olarak yazılımı yalnızca belirli gcc sürümlerinde derlenebilir hale getirdiği iddia edildi. Güya bu problem, Downfall'ın bu parodisinin altyazılarında olduğu gibi, birkaç bilimsel kodu etkiledi . …

24 performance  c++  automatic-differentiation  interval-arithmetic 

Deneme açıklaması: Lagrange enterpolasyonunda, tam denklem noktalarından (polinom sırası ) örneklenir ve 101 noktadan enterpolasyon yapılır. Burada 2 ile 64 arasında değişmektedir. Her seferinde , ve hata grafiği hazırlanır. İşlev eşit aralıklı noktalardan örneklendiğinde, hatanın başlangıçta düştüğü ( yaklaşık 15'ten az olduğu durumlarda gerçekleşir) ve daha sonra hatanın daha da …

24 finite-difference  interpolation  error-estimation  polynomials 

Şu anda herhangi bir paralel kod çalıştırmıyorum, ancak gelecekte OpenMP ve MPI hibridlerini kullanarak paralel kod çalıştırmayı tahmin ediyorum. Hata ayıklayıcılar seri projeler yürütürken benim için paha biçilmez araçlardı. Paralel yazılımda hata ayıklamak için kullanmak isteyen herhangi biri paralel hata ayıklayıcı (veya çoklu hata ayıklayıcı) önerebilir mi? Özgür yazılım tercih …

24 software  parallel-computing 

Anladığım kadarıyla, doğrusal denklem sistemlerini çözmek için iki ana yinelemeli yöntem kategorisi vardır: Durağan Yöntemler (Jacobi, Gauss-Seidel, SOR, Multigrid) Krylov Alt uzay yöntemleri (Konjuge Gradyan, GMRES, vb.) Durağan yöntemlerin çoğunun hatanın Fourier modlarını yinelemeli olarak gevşeterek (düzgünleştirerek) çalıştığını biliyorum. Bunu anlamak gibi Eşlenik gradyan yöntemi ile (Krylov alt uzay yöntemi) …

24 linear-algebra  convergence  krylov-method  conjugate-gradient 

Bilgisayar bilimi araştırmalarında tekrarlanabilirlik giderek daha önemli hale geldi. (Örneğin, bkz Roger Peng bu makaleyi Science , ben de bu tür başka makale ve web sitelerinin farkındayım.) Ancak, bu marka için bir dergi makale ile içermelidir ne kadar bilgi bana net değil (veya çevrimiçi) hesaplamalı araştırmam tekrar üretilebilir (fikri mülkiyet …

23 software  publications  documentation  reproducibility 

Hakkında bir Amazon EC2 hesaplama küme maliyetleri $ 800- $ 1000 3 yıl boyunca fiziksel işlemci çekirdek başına (görev döngüsüne bağlı olarak). Son donanım edinme turumuzda laboratuvarım, Amazon'un kümelerine benzeyen 48 çekirdekli donanım çekirdeğini yaklaşık 300 dolarlık bir çekirdek için aldı. Burada bir şey mi eksik? Moleküler dinamik simülasyonları gibi …

23 cloud-computing  hpc 

Veri noktaları, bir dizi bilgisi I (yaklaşık) bir fonksiyonu takip edeceğini büyük bir satırına asimptotlar . Temel olarak, sıfıra olarak yaklaşır ve , türevlerinde de aynı şey söylenebilir. , vb. Ama f (x) için işlevsel biçimin ne olduğunu bilmiyorum , eğer temel işlevler açısından tanımlanabilecek bile olsa.(xi,yi)(xi,yi)(x_i,y_i)y(x)y(x)y(x)xxxf(x)≡y(x)−(ax+b)f(x)≡y(x)−(ax+b)f(x) \equiv y(x) - …

23 algorithms 

CFD'ye birçok sayısal yaklaşım keyfi olarak yüksek dereceli olarak genişletilebilir (örneğin, süreksiz Galerkin yöntemleri, WENO yöntemleri, spektral farklılık, vb.). Belirli bir problem için uygun bir doğruluk sırasını nasıl seçmeliyim?

23 fluid-dynamics  accuracy 

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 

Hesaplamalı Bilim Yığın Borsası'nda yanıtlanabileceği için , bu soru Çapraz Doğrulama'dan taşındı . 7 yıl önce göç etti . MATLAB programlamasını öğrenmek istiyorum, böylece kendi başıma bazı araştırma / analizler yapabilirim, ayrıca çevrimiçi bulduğum bazı MATLAB komut dosyalarını inceleyebilir / değiştirebilirim. Ancak sorun şu ki, MATLAB'ı karşılayamam. GNU Octave, duyduğum …

23 software  matlab  languages  octave 

Şu anda olduğu gibi, bu soru bizim soru-cevap formatımız için uygun değil. Yanıtların gerçekler, referanslar veya uzmanlık tarafından desteklenmesini bekliyoruz, ancak bu soru muhtemelen tartışma, tartışma, oylama veya genişletilmiş tartışmalar gerektirecektir. Bu sorunun iyileştirilebileceğini ve yeniden açılabileceğini düşünüyorsanız, rehberlik için yardım merkezini ziyaret edin . 7 yıl önce kapandı . …

22 languages 

Lisans düzeyindeki öğrencilere, Bilgisayar Programlamaya Giriş dersi veren bir ders verecek. Kafam biraz karıştı. Hesaplamalı Fizikte bilim adamları C / C ++ veya Python veya Fortran, CUDA vb. Kullanıyorlar. Ne kullanmalıyım? Hayatında istediğin zaman yeni bir programlama dili öğrenebileceğini biliyorum ama bu daha sonra tüm temel programlama kavramlarını ve OOP …

22 python  c++  computational-physics  languages 

Aşağıdaki denklem için iyi sonlu farkların ayrıklaştırılması ne olacaktır: ∂ρ∂t+∇⋅(ρu)=0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? 1D davayı alabiliriz: ∂ρ∂t+ ddx( ρ u ) = 0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho u\right)=0 Nedense bulabildiğim tüm planlar Lagrangian koordinatlarındaki formülasyon içindir. Şu an için bu şema ile karşılaştım ( j …

22 finite-difference  fluid-dynamics