Hesaplamalı Bilim

gerçek bir simetrik matris olduğunu ve özdeğer ayrışması verildiğini varsayalım . Toplamın özdeğeri ne olacağını görmek kolaydır nerede (bkz bir sayıl sabittir bu soruyu ). Biz genel durumda herhangi bir sonuç çıkarmayı CanV Λ V T A + c I c A + D DAAAVΛVTVΛVTV \Lambda V^TA+cIA+cIA + cIcccA+DA+DA + …

11 linear-algebra 

Belirli bir modele göre doğruluk karşılaştırması için şu anda Matlab'ı kullanarak 1d ve 2d görüntüler üretiyorum. Yöntemlerimi standart Gaussian .wfn modeliyle karşılaştırmam gerekiyor ve bunu bir molekülün yanı sıra laplacianın yoğunluğunu analiz ederek yapacağım. .Wfn yaklaşımına (muhtemelen kontur haritası aracılığıyla) yaklaşımımdaki farklılıkların 2-D görüntülerini üretmekle ilgileniyorum ve aynı zamanda iki …

11 matlab  computational-chemistry  publications 

Tekrarlanabilir bilimsel sonuçları göz önünde bulundurarak kodun uzun ömürlü olmasını sağlamak için yayınlanmış en iyi uygulamaların bir listesi var mı? (ör. açık kaynak, dokümantasyon uygulamaları, bağımlılık seçme, bir dil seçme, sanal makineler, vb.). Tipik bilimsel kodun veya diğer yazılımların yarı ömrünü tahmin etmeye çalışan herhangi bir çalışmayı (veya eksik, örnekler …

11 software  publications  reproducibility 

Yeterince basit bir soru: bir çeşit Monte Carlo yönteminin uygun olduğuna karar verdiği göz önüne alındığında, çok boyutlu bir integral yapmak, psödondom sayılarını kullanan düzenli bir MC entegrasyonunun, bir quasirandom sekansı kullanan bir yarı Monte Carlo entegrasyonuna göre herhangi bir avantajı vardır. ? Öyleyse, bu avantajın devreye gireceği durumları nasıl …

11 monte-carlo 

Matematiksel Yazılımda Makine İşlemleri Hesaplama Derneği (ACM TOMS) dergisi , yazılım uygulamalarını içeren sayısal algoritmalar hakkında birçok makale yayınlamaktadır. Yayın politikalarına göre , bir algoritma kağıdının sunulması, söz konusu makalede açıklanan algoritmanın uygulanması için kaynak kodunu içerir. Bu kaynak kodu, ACM Yazılımının Telif Hakkı ve Lisans Anlaşmasına tabidir . Kendi …

11 software  documentation 

Aşağıdaki doğrusal sistemi tekrarlı bir yöntemle verimli bir şekilde çözme umudu var mı? A ∈ Rn × n, x ∈ Rn, b ∈ Rn, n > 10 ile6A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 A x = bAx=bAx=b …

11 linear-solver  dense-matrix  linear-system 

Matris toplama, çarpma vb. Gibi basit görevleri yerine getirmek için Paylaşılan Bellek makineleri (hem C hem de FORTRAN) için OpenMP kodlama deneyimim var (Sadece LAPACK ile nasıl rekabet ettiğini görmek için). OpenMP'yi belgelere bakmaya gerek kalmadan basit görevleri yerine getirecek kadar biliyorum. Son zamanlarda, projelerim için Python'a geçtim ve Python …

11 parallel-computing  python 

İgraph, networkX, R veya Python'da uygulanan iki taraflı grafikler (2 modlu ağlar) için topluluk algılaması için herhangi bir algoritma var mı? Özellikle, toplulukların tespitini iki moddan sadece birinde kısıtlayabilecek böyle bir uygulama var mı?

11 python  graph-theory 

Aşağıdaki matris denklemi in Σ - verilen B ve C matrisleri için - çalışmamda bir kovaryans matrisinin karakterizasyonu olarak görünüyor. Bu denklemin, özellikle sürekli zaman kontrol teorisinde, Lyapunov denklemi olarak bilindiğini ve çözülmesi için bu lineer denklemin özel doğasını kullanan çeşitli iyi bilinen algoritmaların bulunduğunu öğrendim . BΣ+ΣBT+C=0BΣ+ΣBT+C=0B\Sigma + \Sigma …

11 matrix  matrix-equations 

A'nın genel bir seyrek matris olduğunu varsayalım ve özdeğerleri hesaplamak istiyorum. Özdeğerlerin çokluğunu nasıl tespit edeceğimi bilmiyorum. Bildiğim kadarıyla, özel bir durum için, eşlik eden matris yöntemiyle polinom köklerini bulmak, köklerin çokluğunu tespit etmek için RRQR uygulayabiliriz.

11 linear-algebra 

Petsc4py vs c / c ++ / fortran ne kadar yavaş? Ben önemli ölçüde yürütülmekte olan koda bağlıdır, ama bir matris-vektör ürün gibi basit bir şey hakkında ne fark?

11 petsc  python 

Bu tür bir denklemi çözmeye çalışıyorum: ( - ∂2∂x2- f( x ) ) ψ(x)=λψ(x)(-∂2∂x2-f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Burada 0'daf( x )f(x)f(x) basit bir kutba sahiptir , en küçük N özdeğerleri ve özvektörleri için. Sınır koşulları şunlardır: ψ ( 0 ) = 0 ve …

11 finite-difference  ode  accuracy 

Çok noktalı alan verisi kümem var, her nokta veri kümesi yapılandırılmamış ağın tek bir hücresi ile ilgilidir. Amaç, verileri doğrudan veya dolaylı olarak hücre merkezine en doğru şekilde enterpole etmektir. Ters Mesafe Ağırlıklı enterpolasyon kullanırsam, kaynak ile hedef (hücre merkezi) arasındaki mesafenin çok küçük olması durumunda, bir kayan nokta istisnasıyla …

11 interpolation  unstructured-mesh 

(Umarım bu soru bu siteye uyar; değilse özürlerimi kabul et). Belirli bir simülasyon çalıştırdım ve bir zaman serisi y (t), t = 0, 1, ... 20 aldım. Bazı işlevleri denedikten sonra şunu buldum: y(t) =~ 1 / (A t + B) Burada A ve B, R ^ 2> 0,99 ile …

11 regression  publications 

sistemi verildiğinde burada , Jacobi yinelemesinin bir çözücü olarak kullanılması durumunda, b sıfırdan farklıysa yöntemin birleşmeyeceğini okudum A'nın sıfır uzayındaki bileşen . Peki, b'nin A'nın sıfır boşluğunu kapsayan sıfır olmayan bir bileşeni olması koşuluyla , Jacobi yönteminin yakınsak olmadığı resmi olarak nasıl ifade edilebilir ? Bunun matematiksel olarak nasıl resmileştirilebileceğini …

11 linear-algebra  optimization  matrix  iterative-method  linear-solver