«lapack» etiketlenmiş sorular

Üzerinde çalıştığım bir proje için C ++ linear cebir kütüphanelerini inceliyorum. Hala anlayamadığım bir şey, BLAS ve LAPACK'in diğer lineer cebir kütüphanelerine bağlantısı. İçinden bakıldığında lineer cebir kütüphaneleri bu yazıda ben o ilginç bulundu: bazı kütüphaneler BLAS ve LAPACK’tan bağımsızdır bazı BLAS ve LAPACK gerektirir bazılarında BLAS ve LAPACK için …

26 linear-algebra  software  lapack  blas 

Uzun yazı için özür dilerim ama ilk adımda alakalı olduğunu düşündüğüm her şeyi dahil etmek istedim. İstediğim Yoğun Matrisler için Krylov Altuzay Yöntemlerinin paralel bir versiyonunu kullanıyorum .Temelde GMRES, QMR ve CG. (Profillemeden sonra) DGEMV rutinimin acıklı olduğunu fark ettim. Bu yüzden izole ederek buna konsantre olmaya karar verdim. 12 …

15 parallel-computing  matrix  lapack  blas 

Üzerinde çalıştığım bir yazılım projesinde, bazı hesaplamalar yoğun düşük dereceli matrisler için çok daha kolaydır. Bazı problem örnekleri yoğun düşük dereceli matrisleri içerir, ancak bana faktörlerden ziyade tam olarak verilir, bu yüzden düşük dereceli yapıdan yararlanmak istiyorsam dereceyi kontrol etmem ve matrisi faktör yapmalıyım . Söz konusu matrisler tipik olarak …

13 linear-algebra  matrix  lapack  rank  matrix-factorization 

ATLAS'ın derlendiği makine için kendini optimize edebildiğini biliyorum ve bu nedenle maksimum faydalar kaynaktan derlenerek bulundu. LAPACK'i kaynaktan derlemenin herhangi bir faydası var mı? Önceden oluşturulmuş paketi kurmak çok daha kolay olurdu.

13 performance  lapack 

Hem bir matrisin sırasını hem de bir matris denklem sisteminin çözümünü hesaplamak için farklı yöntemler kullanıyorum. Linalg.svd işleviyle karşılaştım. Bunu, Gauss Eliminasyonu ile sistemi çözme çabamla karşılaştırdığımda, hem daha hızlı hem de daha kesin görünüyor. Bunun nasıl mümkün olduğunu anlamaya çalışıyorum. Bildiğim kadarıyla, linalg.svd işlevi, matrisimin özdeğerlerini hesaplamak için bir …

13 linear-algebra  python  lapack 

C ++ için MATLAB mevcut bir kodu taşıma ve çözmek için bir lineer sistem bilgisi am (yerine daha tipik bir şekilde daha A x = b )x A = bxbir=bxA=bA x = bbirx=bAx=b matrisi yoğun ve genel formda olmakla birlikte 1000x1000'den daha büyük değildir. MATLAB'da çözüm, işlev veya eğik çizgi …

12 linear-solver  c++  lapack 

Şu anda ucuz bir matrisi için iyi bir rütbe hesaplamaya çalışıyorum . Bu nedenle QR ayrışmasını kullanarak bir sütunlu döner hesaplıyorumAAA [Q,R,E]=qr(A) Matlab. derecesini kullanarak tahmin ediyorumAAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Bu iyi çalışır ve R'nin tüm diyagonal girişleri üzerinde bir çizim şöyle görünür: RRR Giriş matrisi her …

12 matlab  lapack 

Örneğin nVidia, 7-14x hızlanma vaat eden CUBLAS'a sahiptir. Saf olarak, bu nVidia'nın GPU kartlarından herhangi birinin teorik verimine yakın değildir. GPU'larda lineer cebiri hızlandırmanın zorlukları nelerdir ve halihazırda daha hızlı lineer cebir yönlendirmeleri var mı?

11 linear-algebra  lapack  blas  gpu 

Bazı kuantum mekaniği simülasyonları için fortran 95'e sıkıştım. Dürüst olmak gerekirse, Octave tarafından şımartıldım, bu yüzden verilen matris üssünü aldım. N × n boyutunda (küçük, ) bir çarpıklık -Hermitian matrisi verildiğinde, bu sorunu çözmek için LAPACK'i kullanmanın en etkili yolu nedir? LAPACK95 sarmalayıcısını kullanmıyorum, sadece LAPACK'e doğrudan çağrı yapıyorum.n≤36n≤36n\leq 36n×nn×nn\times …

11 fortran  lapack 

Her gerçek matrisi , dikey simil U dönüşümü kullanılarak gerçek Schur T = U T A U formuna indirgenebilir . Burada T matrisi , ana diyagonal üzerinde 1 x 1 veya 2 x 2 blok ile yarı üçgen formdadır. Gerçek bir özdeğerine her 1 1 ile blok karşılık A ve …

10 eigenvalues  eigensystem  lapack  numerics  scalapack 

PETSc'yi harici bir BLAS / LAPACK kütüphanesi ile derlemek seyrek matrislerdeki performansı önemli ölçüde etkiler mi, yoksa sadece bu kütüphaneleri yoğun matris matematiği için mi kullanır?

10 petsc  sparse-matrix  lapack  blas 

Geçenlerde aynı çizgide eğri-Hermit matrisleri için bir soru sordum . Bu sorunun başarısından esinlenerek, kafamı birkaç saat boyunca bir duvara vurduktan sonra, gerçek asimetrik matrislerin üstel matrisine bakıyorum. Özdeğerleri ve özvektörleri bulma yolu oldukça kıvrımlı görünüyor ve korkuyorum ki kayboldum. Arka plan: Bir süre önce teorik fizik SE hakkında bu …

10 linear-algebra  fortran  lapack 

Hesaplama biliminde yeniyim ve entegrasyon, enterpolasyon, c ++ üzerinde RK4, Numerov vb.Gibi temel yöntemleri öğrendim, ancak son zamanlarda profesör benden matrislerle ilgili problemleri çözmek için LAPACK'i nasıl kullanacağımı öğrenmemi istedi. Örneğin karmaşık bir matrisin özdeğerlerini bulmak gibi. Hiç üçüncü taraf kütüphaneleri kullanmadım ve neredeyse her zaman kendi işlevlerimi yazıyorum. Birkaç …

10 lapack 

LAPACK'in QR rutini Q'yu Hanehalkı reflektörleri olarak depolar. Yansıma ile , böylece sonucun ilk öğesi olur , bu yüzden saklanması gerekmez. Ve gerekli ölçek faktörlerini içeren ayrı bir vektörü saklar . Yani bir reflektör matrisi şöyledir:vvv1 /v11/v11/v_1111ττ\tau'H= Ben- τvvT,'H=ben-τvvT,H=I-\tau v v^T, burada normalleştirilmez. Ders kitaplarında reflektör matrisivvv 'H= Ben- 2 …

9 linear-algebra  matrix  lapack 

Sıfırdan ve tamamen Python3 (utanmaz fiş:) bir kontrol sistemi araç kutusu yazıyorum harold. Geçmişteki araştırmamdan, Riccati çözücüsü hakkında her zaman care.mteknik / alakasız nedenlerden dolayı şikayetlerim var . Bu yüzden kendi rutinlerimi yazıyorum. Bir yol bulamadığım bir şey, en azından kadar iyi bir yüksek performanslı dengeleme algoritması elde etmektir balance.m. …

9 linear-algebra  matlab  lapack  scipy  matrix-equations