Hesaplamalı Bilim

Gelen Cevabıma a MSE üzerinde soruya 2D Hamilton fizik simülasyonu ile ilgili olarak, daha yüksek bir dereceden kullanarak sürmüştür simplektik entegratörü . Daha sonra, farklı zaman adımlarının yöntemlerin farklı siparişlerle küresel doğruluk üzerindeki etkilerini göstermenin iyi bir fikir olabileceğini düşündüm ve bu yönde bir Python / Pylab betiği yazdım ve …

10 python  time-integration 

Bunu okudum (1) Kötü şartlandırılmış operasyonlar iyi şartlandırılmış operasyonlardan önce yapılmalıdır. Örnek olarak, çarpma işlemi yapılmadığı sürece çıkarma kötü koşullandırıldığı için ( x - y ) zxz−yzxz-yzxz-yz olarak hesaplanmalıdır .(x−y)z(x-y)z(x-y)z Ancak, her iki algoritmanın birinci dereceden hata analizi, sadece üç (*) faktörü ile farklılık gösterdiklerini ortaya çıkarır ve neden bunu …

10 stability  floating-point 

İş listelerinde bundan bahsedemiyorum. Bahsettiğim tamsayı programlama, MIP, karma-tamsayı doğrusal olmayan programlama, LP, dinamik programlama vs., fakat SDP yok. Akademide sanayiden çok daha mı trend? Elektrik güç sistemlerindeki akademisyenlere ve endüstri katılımcılarına olan sınırlı maruziyetimden, SDP'nin bağımsız sistem operatörleri tarafından optimum güç akışı sorunlarına uygulanabilmesi için iyi bir şans olduğunu …

10 convex-optimization  semidefinite-programming 

Hesaplamalı akışkanlar dinamiği (CFD), Navier-Stokes denklemlerini (veya bunların basitleştirilmesini) çözmeye adanmıştır. CFD, okyanus ve atmosferik modellerin bir alt kümesi, gerçekçi uygulamalar için aynı denklemleri sayısal olarak çözer. Genel CFD yaklaşımları ve uygulanan gerçekçi durumlar arasındaki farklar ve değişimler nelerdir?

10 fluid-dynamics  simulation  computational-physics  numerical-modelling 

Kapalı. Bu soru konu dışı . Şu anda cevapları kabul etmiyor. Bu soruyu geliştirmek ister misiniz? Sorunuzu güncelleyin o yüzden -konu üzerinde Hesaplamalı Bilim Stack Exchange için. 5 yıl önce kapalı . den büyük ve den daha küçük bir karmaşıklık derecesi var mı ?O(n)O(n)O(n)O(nlogn)O(nlog⁡n)O(n \log n)

10 algorithms  complexity  efficiency 

Runge-Kutta 8. dereceden yöntemini (89) C ++ ile yazılmış bir gök mekaniği / astrodinamik uygulamasında bir Windows makinesi kullanarak kullanmak istiyorum. Bu yüzden herkes belgelenmiş ve kullanımı ücretsiz iyi bir kütüphane / uygulama biliyor mu merak ediyorum? Beklenilen herhangi bir derleme problemi olmadığı sürece C ile yazılmışsa sorun yoktur. Şimdiye …

10 c++  reference-request  libraries  c  runge-kutta 

FEA'yı birkaç yıldır kullandım, ancak doğru şekilde kullanmak ve doğru kullanmak iki farklı şeydir, güvenlik faktörü her şeyin çözümü değildir. Bu soruya net bir cevabım yoksa, doğru kullanmayacağım hissine sahibim: Doğru sonuçlar elde etmek için elementlerin ideal şekline (sadece Jacobian'a dayanarak) yakın olmaları gerektiğinin farkındayım .. Ama neden? Bir koordinat …

10 finite-element  simulation 

Örneğin, kullandığım C ++ seyrek matris kütüphaneleri - Eigen ve SuiteSparse, seyrek matris için herhangi bir SVD işlevine sahip görünmüyorlar. Çok meraklı, seyrek matris için SVD QR / LU'dan daha mı zor?

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

Genellikle iç nokta yöntemlerinin ısınmanın zor olduğu genel atasözü ile karşılaşıyorum. Bu tavsiyenin arkasında sezgisel bir açıklama var mı? Bir iç nokta yönteminde sıcak başlangıçtan fayda bekleyebileceğiniz durumlar var mı? Konuyla ilgili bazı yararlı referanslar önerebilir misiniz?

10 optimization  constrained-optimization 

FE ile yeni başlayan biriyim. Uygulamam, uzayın beş boyutlu olduğu finansal türevlerin fiyatlandırılması. Yani, zaman ekleyerek, sorunun altı boyutu vardır. Etrafa bakmaya çalıştım (Fenics, escript, deal.II, ...), ama benim anlayışım bu yazılımların 3 + 1 (3d uzay + 1d zaman) ile sınırlı olmasıdır. Bu doğru mu? Hedeflediğim dil Python veya …

10 pde  finite-element 

türündeki türündeki denklemleri çözmek için iyi bilinen birçok sayısal yöntem vardır örneğin ikiye ayırma yöntemi, Newton yöntemi vb.f( x ) = 0 ,x ∈ Rn,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, Uygulamamda stokastik bir yöntemle hesaplanmaktadır (sonuç ortalamadır).f( x )f(x)f(x) Bu durumu iyi idare eden sayısal denklem çözme yöntemleri …

10 monte-carlo  nonlinear-equations  stochastic  numerics 

Bazı koşulsuz büyük seyrek doğrusal denklemler sistemi üzerinde çalışıyorum. Bunları çözmek için çift-çift aritmetik veya dört-çift aritmetik kullanmak istiyorum. Nakata, Maho tarafından geliştirilen ve dörtlü çift aritmetik altında sayısal lineer cebirsel hesaplamalar yapabilen MPACK adında bir paket olduğunu biliyorum. Bununla birlikte, seyrek matris için değil, yoğun matris için tasarlanmıştır. Dörtlü …

10 sparse-matrix  matrix  condition-number 

Doğrusal olmayan bir çiftli denklem sistemini çözüyorum ve bağımsız sistemin Jacobian'ını hesapladım. Sonuç gerçekten karmaşıktır, aşağıda (sadece!) 3 × 9 matrisin ilk 3 sütunu ,3×93×93\times 9 (Karmaşıklık kısmen ortaya çıkar, çünkü sayısal şema kararlılık için üstel uyum gerektirir.) Jacobians kullanarak sayısal kodların uygulanması hakkında oldukça genel bir sorum var. Devam …

10 symbolic-computation  jacobian  design-patterns  sympy 

Ben özellikle DGEMM bileşeni için, akademik bilgi işlem için bu iki grafik kartı arasındaki farkı anlamaya çalışıyorum. Ham istatistiklere bakarsak, her ikisi de aynı GK110 yongasına sahip, neredeyse her kategoride karşılaştırılabilir istatistiklere sahip ve inanıyorum ki, aynı çekirdek mimariye sahip. Herhangi bir indirimden önce K20X bir Titan'ın maliyetinin yaklaşık 4 …

10 performance  gpu  efficiency 

Hücre merkezli tek tip olmayan bir ızgarada sonlu hacim yöntemi kullanılırken Dirichlet koşullarının normal olarak nasıl uygulandığını bilmek istiyorum, Mevcut uygulamam, ilk hücrenin değerini sabitlediğim sınır koşulunu getiriyor, φ1= gD( xL)φ1=gD(xL) \phi_1 = g_D(x_L) burada , çözüm değişkenidir ve , alan Dirichlet sınır koşulu değeridir ( NB ). Ancak bu …

10 boundary-conditions  finite-volume  poisson