«algorithms» etiketlenmiş sorular

CFD'de tüm ödevlerimi yapmak için programlama dili olarak python kullanmaya başladım. Programlama konusunda çok az deneyimim var. Makine mühendisliği geçmişindeyim ve Havacılık ve Uzay mühendisliğinde yüksek öğrenim görüyorum. bazen CFD'nin hesaplama yönü denklemleri manipüle etmekten veya matematiği yapmaktan daha sıkıcı olur. Programımızı daha hızlı çalıştırmayı sağlayan genel yönergeler nelerdir? Şeyleri …

16 algorithms  python  fluid-dynamics 

Özellikle heterojen Helmholtz için çok düzeyli ters tabanlı ILU ön koşullandırıcılarının seri performansından çok etkilendim , ancak açık kaynak uygulamaları bulamadığım için şaşırdım. Özellikle, ILUPACK ikili dosyaları akademisyenler için serbestçe kullanılabilir hale getirir, ancak kaynak kodlarını serbest bıraktıkları görülmemektedir. Açık olan hiç kimse uygulamalarını sağlamamış mıydı?

15 pde  algorithms  linear-solver  libraries  preconditioning 

İki birleşmiş PDE'nin bir sistemini iki uzamsal boyutta ve zamanda hesaplamalı olarak çözüyorum. Fonksiyon değerlendirmeleri pahalı olduğundan, çok aşamalı bir yöntem kullanmak istiyorum (Runge-Kutta 4-5 kullanılarak başlatıldı). Önceki beş işlev değerlendirmesini kullanan Adams-Bashforth yönteminde küresel bir hatası vardır (bu, aşağıda atıfta bulunulan Wikipedia makalesinde s = 5 olduğunda ) ve …

14 pde  algorithms  finite-element  error-estimation 

3D olarak milyonlarca veri noktasına çalışan bir veri setim var. Yaptığım hesaplama için, bir yarıçaptaki her veri noktasına komşu (aralık arama) hesaplamam, bir fonksiyona uymaya çalışmam, uyum için hatayı hesaplamam, bir sonraki veri noktası için bunu tekrarlamam vb. Kodum düzgün çalışıyor, ancak veri noktası başına yaklaşık 1 saniye çalıştırmak çok …

14 algorithms  data-structures  nearest-neighbors  spatial-data 

İki vektör arasındaki açı için klasik formülü uygularken: α=arccosv1⋅v2∥v1∥∥v2∥α=arccos⁡v1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} çok küçük / akut açılar için hassasiyet kaybı olduğu ve sonucun doğru olmadığı bulunmuştur. Bu Yığın Taşması cevabında açıklandığı gibi , çözümlerden biri arktanjant kullanmaktır: α=arctan2(∥v1×v2∥,v1⋅v2)α=arctan⁡2(‖v1×v2‖,v1⋅v2)\alpha = \arctan2 \left(\|\mathbf{v_1} \times \mathbf{v_2}\|, \mathbf{v_1} \cdot \mathbf{v_2} \right) …

14 algorithms  precision  numerical-limitations 

Bilimsel hesaplamada kullanılan algoritmaların çoğu, daha az matematik ağırlıklı yazılım mühendisliği formlarında ele alınan algoritmalardan farklı bir doğal yapıya sahiptir. Özellikle, bireysel matematiksel algoritmalar oldukça karmaşık olma eğilimindedir, genellikle yüzlerce veya binlerce satırlık satır içerir, ancak yine de hiçbir durum içermez (yani karmaşık bir veri yapısına etki etmez) ve genellikle …

13 algorithms 

Sayısal kararlılık fedakarlık w / o büyük seyrek matrisler üzerinde matris işlemleri gerçekleştiren bir kütüphane arıyorum. Matrisler 1000+ x 1000+ ve matris değerleri 0 ile 1000 arasında olacaktır. İndeks hesap algoritmasını gerçekleştireceğim, böylece matrisin seri olarak (seyrek) satır vektörlerini oluşturacağım. Her sırayı geliştirirken, doğrusal bağımsızlığı test etmem gerekecek. Matrisimi istenen …

12 linear-algebra  algorithms  matrix 

Bazı Delaunay noktalarının 3D olarak mozaiklenmesine karşılık gelen grafikleri numaralandıran bir algoritma var mı? Eğer öyleyse, herhangi bir "Delaunay grafiğine" karşılık gelen geometrilerin etkin bir parametrelenmesi var mı? Önceden bir bağlanma bilgisi olmadan belirli bir kompozisyondaki moleküllerin tüm kararlı geometrilerini sistematik olarak numaralandırmak istiyorum. EDIT: , N köşeli grafik kümesi …

12 algorithms  computational-chemistry  delaunay-triangulation 

Sayısal olarak aşağıdaki integral dönüşümü ters çevirmeye çalışıyorum: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)dxF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Bu nedenle, belirli bir için yaklaşık f ( x ) değerine ihtiyacım var:F(y)F(y)F(y)f(x)f(x)f(x) ve F ( y ) gerçek ve pozitiftirf(x)f(x)f(x)F(y)F(y)F(y)(sürekli olasılık dağılımlarıdır) gerçek ve pozitiftirx,yx,yx,y(büyüklüklerdir) Şu anda bunu yapmak için çok dağınık ve kaba …

11 algorithms  numerical-analysis  approximation  inverse-problem 

Dolayısıyla Cholesky ayrışma teoremi, herhangi bir gerçek simetrik pozitif tanımlı matris , nin daha düşük bir üçgen matris olduğu bir Cholesky ayrışma olduğunu belirtir .MMMM= L L⊤M=LL⊤M= LL^\topLLL verildiğinde , Cholesky faktörünü hesaplamak için hızlı algoritmalar olduğunu zaten biliyoruz .MMMLLL Şimdi bir dikdörtgen matrisi verildiğini ve pozitif kesin olduğunu bildiğimizi …

11 linear-algebra  algorithms 

Merhaba scicomp topluluğu, Böyle olarak çerçeveler kullanılarak grafik algoritmaları alanında çalışmış NetworkX (Python), JUNG ve YFiles (Java). Şimdi paralel ve yüksek performanslı bilgisayar alanına giriyorum. Yeni bir proje için, aşağıdaki özelliklere sahip bir C ++ grafik kütüphanesi arıyorum: algoritma geliştirmeyi sağlayan sezgisel bir arayüze sahiptir dinamik işlemleri destekler: örneğin, rastgele …

11 algorithms  parallel-computing  libraries  graph-theory 

Bu basit bir simetri numaralandırma problemidir. Burada tam bir arka plan veriyorum, ama kuantum kimyası bilgisi gerekli değil. İki parçacık integrali geçerli: ⟨ i j | k l ⟩ = ∫ ψ * ı ( x ) ψ * j ( X ' ) ψ k ( x ) ψ …

11 algorithms  performance 

Büyük bir seyrek matrisin (~ 2 milyon element) Lanczos köşegenleştirmesini yapıyorum. Lanzcos algoritmasındaki adımların neredeyse tamamı, yakınsama olup olmadığını kontrol etmek için Lanczos matrisinin köşegenleştirilmesi dışında GPU'da paralel olarak yapılır. Bunun için Sayısal Tarifler TQLI algoritması kullanıyorum. Üçgen bir matrisin paralel veya kolayca paralelleştirilebilen öz sistemi bulmak için yöntemler var …

11 linear-algebra  algorithms  eigensystem 

Tek bir mpi_allgather işlemi veya bir mpi_scatter ve bir mpi_gather işlemi kullanarak paralelleştirilebilir bir sorun üzerinde çalışıyorum. Bu işlemler bir while döngüsü içinde çağrılır, bu nedenle birçok kez çağrılabilirler. Bir MPI_allgather şeması ile uygulamada, yinelenen matris çözme için tüm süreçlere dağıtılmış bir vektör topluyorum. Diğer uygulamada, dağıtılmış vektörü tek bir …

11 algorithms  complexity  mpi 

Soru Yapılandırılmamış altı yüzlü hücre ağına göre bir puan bulutunu nasıl sıralarsınız? Her hücrenin bir merkezi ve onu temsil eden benzersiz bir etiketi vardır. Temelde iki bulut noktası vardır (orijinal nokta bulutu ve hücre merkezlerinin bir nokta bulutu), ancak hücre geometrisi bilgileri (sınırlama kutusu) kullanımda olabilir, emin değilim. Sonuçlar Bazı …

11 algorithms  computational-geometry