Hesaplamalı Bilim

Benim sorum, bu referansta tarif edildiği gibi QMC yöntemlerinden gözlemlenebilirlerin çıkarılmasıyla ilgilidir . Path Integral Monte Carlo gibi çeşitli QMC yöntemlerinin resmi türevini anlıyorum. Ancak, günün sonunda hala bu tekniklerin etkili bir şekilde nasıl kullanılacağı konusunda kafam karıştı. Kuantum MC yöntemlerinin türetilmesinin temel fikri, Trotter yaklaşımı yoluyla, bir kuantum sisteminin …

10 monte-carlo  quantum-mechanics 

Parareal, PITA ve PFASST algoritmalarının tümü , zamana bağlı sorunların zaman içindeki çözümünü paralel hale getirmek için etki alanı genelinde kullanılan tekniklerdir. Bu yöntemlerin ardındaki yol gösterici ilkeler nelerdir? Aralarındaki ana farklar nelerdir? Birinin diğerine dayandığını söyleyebilir miyim? Nasıl? Başvuruları ne olacak? "Hangisi daha iyi?" Sorusuna cevap verilmeyeceğini biliyorum, ancak …

10 parallel-computing  time-integration 

Açık FEM ile örtük FEM arasındaki fark nedir? Buradaki mesaja göre , tek fark örtük veya açık zaman entegrasyonunun kullanılmasıdır. Okuduğum bir kitaptan hatırladığım gibi, örtük FEM, kitlenin düğümlere götürülmediği yerdir. Açık ve kapalı FEM'in kesin tanımları nelerdir?

10 finite-element  implicit-methods  explicit-methods 

Modelin indirgenmesi için, matrisi olan bir matrisin en büyük tekil değerlerini - say 20 ile ilişkili olan sol tekil vektörleri hesaplamak istiyorum , burada ve . Ne yazık ki, matrisim herhangi bir yapı olmadan yoğun olacak. N ≈ 10 6 k ≈ 10 3 AA ∈ RN-, kA∈RN,kA \in \mathbb …

10 linear-algebra  python  reference-request  numpy  svd 

IMPLICIT formunda (zaman ayrıklaştırma için) sonlu hacim yöntemi ile bir dizi kısmi diferansiyel denklemler (PDE) çözen bazı kod için arıyorum. Sonuç olarak, ADI / TDMA şeması tarafından işlenen, x, y, z yönlerinde üç boyutlu bir denklem sistemi vardır. PDE'lerin CUDA ile örtük çözümü hakkında hiçbir şey bulamıyorum. ADI / TDMA …

10 parallel-computing  implicit-methods  cuda 

Bilimsel bilişimde yeniyim. Bir Octree bir Fortran (tercihen f90) uygulaması arıyorum. Benim sorunum, herhangi bir kutuda bazı N parçacıklarından (veya Laplace denklemi gibi bir denklemi çözmek için bir integral denklem yöntemine takılabilen yoğunluk değerini bildiğim kaynaklardan) bulunmayana kadar alanımı bölen bir Octree gerektirir. Bulabildiğim tek şey bu C ++ uygulamasıydı …

10 fortran  integral-equations 

1D ısı denklemini çözmek için Crank-Nicolson sonlu fark şemasını kullanıyorum. Isı denkleminin maksimum / minimum prensibinin (yani maksimum / minimumun başlangıç ​​koşulunda veya sınırlarda gerçekleştiğini) ayrık çözüm için de geçerli olup olmadığını merak ediyorum. Bu muhtemelen Crank-Nicolson'un istikrarlı ve yakınsak bir şema olduğu anlamına gelir. Ama öyle görünüyor ki bunu …

10 linear-algebra  pde  finite-difference  crank-nicolson 

Herkesin bildiği gibi, sıkıştırılamaz akış gerçekte mevcut değildir, yöneten denklemleri basitleştirmek için getirilen bir varsayımdır. Bu varsayımı doğrudan uygulayamayız. Genellikle Mach sayısı (sıkıştırılamaz akış için M <0.3), yoğunluk değişimi (sıfır yoğunluk değişimi) ve hızın ayrıştırılması (sıkıştırılamaz akış için sıfıra eşittir) akışı sıkıştırılamaz akış olarak tanımlamak için ortak kriterdir. Isı transferi …

10 fluid-dynamics 

Süreksiz Galerkin yöntemini (DG) ve aşağıdaki takdir yetkisini kullanarak 2D Poisson denklemini çözmeye çalışıyorum (bir png dosyası var, ancak yüklememe izin verilmiyor, üzgünüm): Denklem: ∇⋅(κ∇T)+f=0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Yeni denklemler: q=κ∇T∇⋅q=−fq=κ∇T∇⋅q=−fq = \kappa \nabla T\\\nabla \cdot q = -f ve sayısal akıları ile zayıf form …

10 pde  finite-element  fenics 

Üretilen çözeltiler (MMS) yönteminde, kesin bir çözüm önerilmektedir, denklemlerde ikame edilir ve karşılık gelen kaynak terimi hesaplanır. Çözelti daha sonra kod doğrulaması için kullanılır. Sıkıştırılamaz Navier-Stokes denklemleri için MMS, süreklilik denkleminde kolayca (sıfırdan farklı) bir kaynak terimine yol açar. Ancak tüm kodlar süreklilik denklemlerinde kaynak terimlerine izin vermez, bu nedenle …

10 navier-stokes  incompressible  verification 

C ++ 11, örneğin C ++ 03'ün bir kopya oluşturma veya kopyalama ataması gerçekleştirmesi gereken durumlarda kod performansını artırabilecek taşıma semantiği sunar. Bu makalede , aşağıdaki kod C + 11 ile derlendiğinde 5 kat hızlanma yaşadığı bildirilmektedir: vector<vector<int> > V; for(int k = 0; k < 100000; ++k) { vector<int> …

10 finite-element  performance  c++ 

Standart lineer sistemi çözmek için birçok hesaplama yöntemine aşinayız Ax=b.Ax=b. Ax=b. Ancak, formun daha genel (sonlu boyutlu) doğrusal bir sistemini çözmek için herhangi bir "standart" hesaplama yöntemi olup olmadığını merak ediyorum LA=B,LA=B, LA=B, diyelim ki, AAA bir olan m1×n1m1×n1m_1\times n_1 matrisi, BBB bir bir m2×n2m2×n2m_2\times n_2 matris ve LLL alarak …

10 linear-algebra 

Matlab'de, hem linsolve ve mldivide , belirlenen overdetermine tabutlarına örneklerin tümünde, lineer denklem sisteminin çözülmesi için kullanılır. Belgelerini okurken aralarında ne gibi farklılıklar olduğunu merak ediyordum. Üç durumda hemen hemen aynı matris çarpanlarına ayırma ve üçgenleştirme algoritmaları kullanıyorlar mı? A opts'ta özelliklere sahipse, linsolve mldivide göre daha hızlıdır, çünkü linsolve …

10 matlab  linear-solver 

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. 2 yıl önce kapalı . Başlığın dediği gibi, kodu belirli bir yerde durdurabilmek ve güvenli bir şekilde durdurmak istiyorum. Bunu yapmak için bir …

10 matlab 

Bugün diferansiyel formlar teorisi hakkında okurken, bana ikinci dereceden Sonlu Hacim Metodu'nu (FVM) hatırlatması beni etkiledi. Bu şekilde önemsiz olduğunu düşünmek için uğraşıyorum ya da daha derin bir bağlantı var mı. Diferansiyel formlar, bir yüzeydeki sıvı akısı gibi ikinci dereceden FVM'de köklü bazı kavramları genelleştirmeye hizmet eder ve hepimiz FVM'deki …

10 fluid-dynamics  finite-volume