«sparse-matrix» etiketlenmiş sorular

Depolama, montaj, işlemler ve yalnızca sıfır olmayan elemanların depolandığı seyrek matrislerle başa çıkmanın diğer yönleri ile ilgili sorular. Doğrudan seyrek matrisler ile değil, ancak seyreklik kullanmanın diğer yolları [seyrek-operatör] ile etiketlenmelidir.

17
Python için yüksek kaliteli bir doğrusal olmayan programlama çözücü var mı?
Çözülmesi gereken birkaç dışbükey olmayan küresel optimizasyon problemim var. Şu anda , oldukça etkili olan MATLAB's Optimizasyon Araç Kutusu'nu (özellikle, fmincon()algoritma = ile 'sqp') kullanıyorum . Ancak, kodumun çoğu Python'da ve optimizasyonu Python'da da yapmak isterim. Rekabet edebilecek Python bağlarına sahip bir NLP çözücü var mı fmincon()? O olmalı Doğrusal …

4
Seyrek bir lineer sistem çözücü seçerken hangi yönergeleri izlemeliyim?
Seyrek doğrusal sistemler uygulamalarda artan sıklık ile ortaya çıkıyor. Birinin bu sistemleri çözmek için seçebileceği çok fazla rutini var. En üst düzeyde, doğrudan (örneğin, seyrek Gauss ortadan kaldırılması veya Cholesky ayrışımı, özel sipariş algoritmaları ve çok yönlü yöntemler) ve yinelemeli (örneğin GMRES, (iki) eşlenik gradyan) yöntemleri arasında bir havza vardır. …

7
Genel bir matrisin en büyük özdeğerini hesaplamanın en hızlı yolu nedir?
EDIT: Herhangi bir özdeğerin büyüklüğü bir veya daha fazla olup olmadığını test ediyorum. Büyük, seyrek, simetrik olmayan bir matrisin en büyük mutlak özdeğerini bulmam gerekiyor. eigen()Tüm özdeğerleri bulmak için EISPACK veya LAPACK'den QR algo kullanan R işlevini kullanıyorum ve sonra abs()mutlak değerleri elde etmek için kullanıyorum . Ancak, daha hızlı …

3
Çözme
ve matrislerine sahibim . seyrek ve bir ile çok (birkaç milyon mertebesinde olabilir.), Büyük bir bir ile uzun matris oldukça küçük ( ) ve her sütunun kutunun sadece tek bir sahiptir geri kalan kısım ile girişi öyle ki 's, . çok büyük, bu yüzden tersine çevirmek gerçekten zor ve gibi …

3

5
Güzel bir yazılım tasarımı için% 20 performans cezası
Nesneye yönelik programlamayı en iyi şekilde kullanmayı kendime öğretmenin bir yolu olarak seyrek matris hesaplamaları için küçük bir kütüphane yazıyorum. Parçaların (seyrek matrisler ve bağlantı yapılarını tanımlayan grafikler) çok gevşek bir şekilde birleştiği güzel bir nesne modeline sahip olmak için çok çalıştım. Kendi görüşüme göre, kod bunun için çok daha …


4
Minimum bant genişliği bantlı bir matris üretmek için değişkenleri nasıl yeniden sıralayabilirim?
2B Poisson denklemini sonlu farklarla çözmeye çalışıyorum. Bu süreçte, her denklemde sadece değişkenli seyrek bir matris elde ediyorum . Örneğin, değişkenler U olsaydı, o zaman ayrıklaştırma şöyle olur:555UUU Ui−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jUi−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jU_{i-1,j} + U_{i+1,j} -4U_{i,j} + U_{i,j-1} + U_{i,j+1} = f_{i,j} Bu sistemi yinelemeli bir yöntemle çözebileceğimi biliyorum, ama değişkenleri uygun şekilde sipariş …

1
SciPy eigsh () neden harmonik osilatör durumunda hatalı özdeğerler üretir?
Hesaplamalı fizik bağlamında, büyük seyrek matrislerin özdeğer hesaplamaları gerçekleştirmek için bazı büyük kod geliştiriyorum. Özdeğerler analitik olarak iyi bilindiği için rutinlerimi basit harmonik osilatöre karşı bir boyutta test ediyorum. Bunu yaparak ve kendi rutinlerimi SciPy'nin dahili çözücüleriyle karşılaştırarak, aşağıdaki çizimde gösterilen garipliğe rastladım. Burada ilk 100 sayısal olarak hesaplanmış özdeğer …

2
Çok büyük ve çok seyrek bir bitişiklik matrisinin tüm özdeğerlerini hesaplayın
Her biri yaklaşık n ~ 100000 düğümü olan iki grafiğim var. Her iki grafikte de, her bir düğüm tam olarak 3 düğüme bağlanır, böylece bitişiklik matrisi simetrik ve çok seyrek olur. Zor olan nokta, bitişiklik matrisinin tüm özdeğerlerine ihtiyacım var ama özvektörlere ihtiyacım yok. Doğru olmak için, bu benim hayatımda …

3
Thomas algoritması, simetrik bir çapraz baskın seyrek tridiagonal doğrusal sistemi çözmenin en hızlı yolu mu
Thomas algoritmasının algoritmik karmaşıklık açısından (LAPACK vb. Hem Thomas algoritmasının hem de multigrid'in karmaşıklığı olduğunu biliyorum , ama belki multigrid için sabit faktör daha az? Bana multigrid daha hızlı olabilir gibi görünmüyor ama ben olumlu değilim.O(n)O(n)O(n) Not: Matrislerin çok büyük olduğu durumu düşünüyorum. Doğrudan veya yinelemeli yöntemler kabul edilebilir.

3
Birçok sağ taraf için seyrek doğrusal çözücü
Birçok sağ tarafı (300 ila 1000) ile aynı seyrek doğrusal sistemi (300x300 ila 1000x1000) çözmem gerekiyor. Bu ilk soruna ek olarak, farklı sistemleri de çözmek istiyorum, ancak aynı sıfır olmayan elemanlarla (sadece farklı değerler), sabit seyreklik desenine sahip birçok seyrek sistem. Matrislerim belirsiz. Çarpanlara ayırma ve başlatma işlemlerinin performansı önemli …



4
Dörtlü çift aritmetik seyrek matris paketi var mı?
Bazı koşulsuz büyük seyrek doğrusal denklemler sistemi üzerinde çalışıyorum. Bunları çözmek için çift-çift aritmetik veya dört-çift aritmetik kullanmak istiyorum. Nakata, Maho tarafından geliştirilen ve dörtlü çift aritmetik altında sayısal lineer cebirsel hesaplamalar yapabilen MPACK adında bir paket olduğunu biliyorum. Bununla birlikte, seyrek matris için değil, yoğun matris için tasarlanmıştır. Dörtlü …

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.