Hesaplamalı Bilim

Bilimsel sorunları çözmek için bilgisayar kullanan bilim adamlarının soru-cevapları

1
Strang bölünmesinin optimal kullanımı (reaksiyon difüzyon denklemi için)
Basit bir 1D reaksiyon difüzyon denklemine çözümü hesaplarken garip bir gözlem yaptım: ∂∂ta=∂2∂x2a−ab∂∂ta=∂2∂x2a−ab\frac{\partial}{\partial t}a=\frac{\partial^2}{\partial x^2}a-ab ∂∂tb=−ab∂∂tb=−ab\frac{\partial}{\partial t}b=-ab ∂∂tc=a∂∂tc=a\frac{\partial}{\partial t}c = a Başlangıç değeri sabit (bir ) ve sadece üzerinde integrali ilgilendiğim mesafede ile ( ). Amacı, ve denklem sadece bu integrali değerlendirmektir.bbbb(0,x)=b0b(0,x)=b0b(0,x)=b_0aaa000111∫10a(t,x)dt∫01a(t,x)dt\int_0^1a(t,x)dtccc∂∂tc=a∂∂tc=a\frac{\partial}{\partial t}c = a Difüzyon ve reaksiyon arasındaki bağlantı …

4
A'A ve AA 'formülasyonlarının durum sayısı
(Yousef Saad, seyrek doğrusal sistemler için yinelemeli yöntemler , s. 260)c o n d(bir'Bir ) ≈ c O , n d( A)2cond(A′A)≈cond(A)2cond(A'A) \approx cond(A)^2 Bu doğru mu birbir'AA′AA' de? Durumunda birAA dır-dir N-× MN×MN\times M ile N≪MN≪MN \ll M, Gözlemliyorum cond(A′A)≫cond(AA′)cond(A′A)≫cond(AA′)cond(A'A) \gg cond(AA') Bu, formülasyon açısından AA′AA′AA' bu durumda tercih …

1
Polinomların geri kalanını hesaplamak için sayısal olarak kararlı algoritmalar
Let ve . hesaplamak için asimptotik olarak hızlı ve sayısal olarak kararlı algoritmalar . Amaçlanan uygulamalarda her iki f, g , çift kesinlikli kayar nokta katsayılı yoğun polinomlardır. Ancak, şimdilik, uygulamadan ziyade algoritmalarla ilgileniyorum. Sayısal polinomların GCD'sinin hesaplanması için algoritmalar için referanslar da takdir edilmektedir.f,g∈R[x]f,g∈R[x]f, g \in \mathbb{R}[x]degf>deggdeg⁡f>deg⁡g\deg f > …




4
Cebirsel Multigrid Kodu
Cebirsel Multigrid Yöntemlerinin (AMG) uygulanması hakkında daha fazla ayrıntı anlamak istiyorum. AMG için enterpolasyon, kaba ızgara operatörü ve kaba ızgara seçiminin tüm ayrıntılarını açıklayan "A Multigrid Eğitimi" ni okuyorum. Ancak, sanırım, oynamak ve kod okumak gibi bir şey yok. Bu yüzden, herhangi birinin klasik bir AMG "örnek kodu" tanıyıp tanımadığını, …

1
Hızlı Fourier Dönüşümü kullanarak karışık sınır koşullarında bir 2D poisson problemini çözmek için hangi fourier serisine ihtiyaç vardır?
Sınır koşullarının hepsi tek tip olduğunda poisson problemini çözmek için hızlı bir fourier dönüşümünün kullanılabileceğini duydum ... Dirichlet için sinüs serisi, neumann için kosinüs ve her ikisi de periyodik olarak. 2B dikdörtgen bir alan göz önüne alındığında, karşıt iki tarafın periyodik sınır koşulları ve diğer ikisinin dirichlet koşulları olduğunu varsayalım. …

2
Çapraz baskın matrislere yinelemeli yöntemlerin güvenli uygulanması
Aşağıdaki doğrusal sistemin verildiğini varsayalım Lx=c,(1)(1)Lx=c,Lx=c,\tag1 nerede LLL pozitif olduğu bilinen Laplacian ağırlıklı semi−semi−semi-tek boyutlu boş uzay ile 1n=(1,…,1)∈Rn1n=(1,…,1)∈Rn1_n=(1,\dots,1)\in\mathbb{R}^nve çeviri varyansı x∈Rnx∈Rnx\in\mathbb{R}^{n}, yani, x+a1nx+a1nx+a1_n (türevi olan fonksiyon değerini değiştirmez) (1)(1)(1)). Tek olumlu girişlerLLL negatif diyagonal girişlerin mutlak değerlerinin bir toplamı olan diyagonalindedir. Alanında oldukça alıntı yapılan bir akademik çalışmada, LLL …

1
Periyodik sınır koşulları ile Schrodinger denklemi
Aşağıdakilerle ilgili birkaç sorum var: 1D Schrodinger denklemi krank nicolson ayrıklaştırma kullanarak sonuçlanan tridiagonal matris ters çevirerek çözmeye çalışıyorum. Benim sorunum şimdi periyodik sınır koşulları ile ilgili bir soruna dönüştü ve bu yüzden Sherman Morrison algoritmasını kullanmak için kodumu değiştirdim. Varsayalım vben tridiagonal matrisinin tersini istediğinizde her zaman adımında benim …

3
Blok yapısız belirsiz sistemler için yinelemeli yöntemler
Belirsiz matris sistemleri, örneğin eyer noktası problemlerinin karışık sonlu elemanlar tarafından ayrıklaştırılmasında ortaya çıkar. Sistem matrisi daha sonra forma konabilir (ABBtC)(ABtBC)\begin{pmatrix} A & B^t \\ B & C\end{pmatrix} nerede AAA negatif (yarı) -definit, CCC pozitif (yarı) kesin ve BBBkeyfi. Elbette, konvansiyona bağlı olarak doğruluk koşullarını kullanabilirsiniz, ancak bu hemen hemen …

1
Bir işlevi tekilliklerle entegre etmek için hangi sayısal kareleme seçilir?
Örneğin, sayısal olarak hesaplanması için istiyorum arasında -norm bir etki olup, sıfır içerir, ben Gauss quadrature denedim ve başarısız, entegrasyon için küresel koordinatlar kullanarak birim topu gerçek normalden çok uzak, bunu yapmak için iyi bir yolu var mı? Bu sorun genellikle yeniden açılan köşeleri olan alanlar için oyuncak problemlerini hesaplayan …

2
'sıfır' ile olası bölünme ile sayısal entegrasyon
Entegre olmaya çalışıyorum ∫10t2 n + 2tecrübe(αr0t) dt∫01t2n+2exp⁡(αr0t)dt\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt ki bu basit bir dönüşüm ∫∞1x2 ntecrübe( - αr0x ) dx∫1∞x2nexp⁡(−αr0x)dx\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx kullanma t =1xt=1xt = \frac1{x}çünkü uygun olmayan integrallerin sayısal olarak tahmin edilmesi zordur. Ancak bu yeni integrali sıfıra yakın değerlendirme problemine yol açar. Aralık sadece uzunluk …

1
Düzenli ızgara üzerinde iç içe diseksiyon
Seyrek doğrusal sistemleri doğrudan çarpanlara ayırma yöntemlerini kullanarak çözerken, kullanılan sıralama stratejisi, faktörlerdeki sıfır olmayan elemanların doldurma faktörünü önemli ölçüde etkiler. Böyle bir sıralama stratejisi iç içe diseksiyondur. Sadece ızgara parametreleri (ilk sıra farkları ile bir M x N kare sonlu fark ızgarası varsayalım) verilen önceden iç içe diseksiyon siparişi …

2
FEM DG yöntemlerini Riemann çözücülerine bağlama
Süreksiz galerkin sonlu eleman çözücülerini Riemann çözücülerle birleştiren iyi kağıtlar veya kodlar var mı? Eliptik ve hiperbolik bağlanma problemlerini keşfetmem gerekiyor, ancak çoğu bölme yöntemi en iyi şekilde geçici. Çok miktarda FEniCS kodum olduğu için, sadece Riemann çözücüsünü onunla birleştirmek istiyorum. Basit bir Karaca çözücü bir başlangıç ​​olsa da, daha …

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.