Hesaplamalı Bilim

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 

Birçok uyarlanabilir FEM kütüphanesi, düğümleri, kenarları, üçgenleri, tetrahedra vb. Ekleme / çıkarma işlemlerini gerçekleştirmek için daha gelişmiş ağ veri yapıları kullanır . Örneğin, en yaygın kitaplık, uyarlanabilir ağ inceltme için oktree veri yapılarını kullanır; statik bir ağda hesaplamalar için kullanılan oktrees genellikle bulamazsınız. Uyarlamalı FEM için lineer cebir tarafında neler …

9 linear-algebra  finite-element  sparse-matrix 

Her zaman kolay paralelleştirmenin DG yöntemlerinin avantajlarından biri olduğunu duydum, ancak bu nedenlerden herhangi birinin neden sürekli Galerkin için geçerli olmadığını gerçekten anlamıyorum.

9 finite-element  parallel-computing  discontinuous-galerkin 

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Yukarıdaki trapezi kullanarak func()sınırları arasındaki 1D …

9 c++  performance 

Büyük bir kübik özdeğer problemim var: (bir0+ λbir1+λ2bir2+λ3bir3) x =0.(bir0+λbir1+λ2bir2+λ3bir3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Bunu lineer bir özdeğer problemine dönüştürerek çözebilirdim ama bir sistemle sonuçlanırdı 32323^2 büyük gibi: ⎡⎣⎢-bir0000ben000ben⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥= λ⎡⎣⎢bir1ben0bir20benbir300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[-bir0000ben000ben][xyz]=λ[bir1bir2bir3ben000ben0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ 0 & 0 & …

9 c++  eigensystem  eigenvalues 

FEM ve XFEM arasındaki temel farklar nelerdir? Ne zaman FEM'in XFEM intead'ini kullanmamalıyız? Başka bir deyişle, yeni bir sorunla karşılaştığımda, bunlardan hangisini kullanmayı nasıl bilebilirim?

9 finite-element  numerical-analysis  adaptive-mesh-refinement  numerics 

Gauss-Seidel yerine SOR yapmaya değip değmeyeceğini söylemek için basit bir kural var mı ? (ve rexation parametresini tahmin etmenin olası yoluωω\omega) Demek istediğim, sadece matrise veya matrisin temsil ettiği belirli bir problem bilgisine bakarak ? Bu soruların cevabını okuyordum: Ardışık aşırı gevşeme (SOR) yöntemini optimize etmek için herhangi bir buluşsal …

9 linear-solver  iterative-method  heuristics 

Bu sorular compcci'de biraz konu dışı olabilir. gerekiyorsa lütfen nereye uyduğunu önerin. Soru, tüm simülasyonların verimli bir şekilde nasıl yönetileceği ile ilgilidir. diyelim ki bir simülasyon belirli bir önerilen değer aralığında tanımlanması gereken 2 parametreyi sabitlemeyi gerektirir. İki parametrenin bir çiftinin ürettiği daha iyi bir sonuç bulmak için (örneğin simülasyon …

9 monte-carlo  data-management 

Kesinlikle dışbükey kısıtsız optimizasyon problemini düşününX_ \ text {opt} ' un benzersiz minimasını göstermesine izin verin ve x_0 , x_ \ text {opt}' a verilen bir başlangıç ​​yaklaşımı olsun. Biz bir vektör arayacak x bir \ epsilon yakın çözüm \ mathcal {O} eğer \ begin {denklem} \ frac {|| x …

9 convex-optimization  complexity 

Hesaplama Bilimi için çok iyi bir posta listesi veya google grupları olup olmadığını merak ediyorum , burada soru sormak ve cevaplamak yerine soruları tartışabiliriz. Aslında, PDE'lerin paralel hesaplama ve sayısal çözümleriyle daha fazla ilgileniyorum. Ama bu bölgedeki insanların ne yaptığını ve nasıl yaptığını bilmiyorum. Bu alandaki yol haritasını kavramak için …

9 parallel-computing  reference-request 

Bazı sınıfsal sayısal problemleri çözecek küçük bir MATLAB paketi yazıyorum. Algoritmanın 3 aşaması vardır ve kullanıcının her aşama için 5 seçeneği vardır. Tüm sorunu işlev ve 3 anahtar durumu (her algoritma aşaması için bir tane) kullanarak uyguladım . İyi çalışıyor ama daha fazla şey yapmayı (5'ten fazla seçenek ve bir …

9 matlab 

Son birkaç gündür bu hatayı hata ayıklamaya çalışıyorum. Poisson denklemini, bilinmeyenlerin hücre merkezlerinde ve hücre yüzlerindeki akılarda tanımlandığı, muntazam olmayan bir sonlu hacim ağında bir adım yük dağılımı (elektrostatik / yarı iletken fiziğinde ortak bir sorun) için çözüyorum. 0=(ϕx)x+ρ(x)0=(ϕx)x+ρ(x) 0 = (\phi_x)_x + \rho(x) ücret profili (kaynak terim) tarafından verilir, …

9 discretization  finite-volume  poisson 

Fortran'da program yapmayı öğrenmenin en iyi yolunun ne olduğunu merak ediyordum. Kitap arıyorum, bana uygun görünen örnekler "Fortran 95/2003 Açıklaması" ve "Bilim Adamları ve Mühendisler İçin Fortran" idi . Bunlar hakkında ne düşündüğünüzü bilmiyorum. Merak ettiğim başka bir şey de hangi derleyiciyi kullanmam gerektiğidir? Bir Windows 7 makinesinde (64 bit) …

9 fortran 

3 çiftli denklem sistemine sonlu farklar yöntemi uyguluyorum. Denklemlerin ikisi birbirine bağlı değildir, ancak üçüncü denklem diğer ikisiyle de eşleşir. Denklemlerin sırasını değiştirerek, den kadar katsayı matrisinin simetrik hale geldiğini fark ettim .( x , y, z)(x,y,z)(x, y, z)( x , z, y)(x,z,y)(x, z, y) Bunu yapmanın bir avantajı var …

9 finite-difference  symmetry 

Hem doğrusal hata düzeltme (LDC) hem de tam yaklaşım şeması (FAS) kullanarak bir V-Cycle multigrid çözücü uyguladım. Benim sorunum şu: LDC kullanarak artık döngü başına ~ 0.03 faktör azalır. FAS uygulaması da doğrusal bir faktörle yakınsar, ancak faktör sadece ~ 0.58'dir. Bu nedenle FAS, döngü sayısının yaklaşık 20 katına ihtiyaç …

9 multigrid