«numpy» etiketlenmiş sorular

Python's numpy kitaplığını kullanarak, satırlarının ve sütunlarının birçoğunun sırasını değiştirerek, yoğun bir kare geçiş matrisini yerinde değiştirmek istiyorum. Matematiksel olarak bu, matrisin permütasyon matrisi P ile önceden çarpılması ve P ^ -1 = P ^ T ile çarpılması ve sonradan hesaplanması açısından makul bir çözüm değildir. Şu anda elle satırları …

27 linear-algebra  python  numpy 

Python / Numpy dizileri artan dizi boyutlarıyla nasıl ölçeklenir? Bu, bu soru için Python kodunu karşılaştırırken fark ettiğim bazı davranışlara dayanıyor: Bu karmaşık ifadenin en iyi dilimler kullanılarak nasıl ifade edileceği Sorun çoğunlukla bir diziyi doldurmak için dizin oluşturmayı içeriyordu. Python döngüsünde (çok iyi olmayan) Cython ve Numpy sürümlerini kullanmanın …

21 python  performance  numpy 

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 …

15 python  sparse-matrix  numpy  eigenvalues 

Python'da aşağıdaki ifadeyi uygulamak istiyorum: burada x ve y n boyutunda nümerik dizilerdir ve k nümerik bir dizidir n × n boyutunda . N boyutu yaklaşık 10000 olabilir ve fonksiyon birçok kez değerlendirilecek bir iç döngünün parçasıdır, bu nedenle hız önemlidir.xben= ∑j = 1i - 1ki - j , jbiri …

14 python  numpy 

Python 3'e (numpy kullanarak) geriye doğru bir Euler çözücü uyguladım. Kendi rahatlığım için ve bir egzersiz olarak, degradenin sonlu fark yaklaşımını hesaplayan küçük bir işlev yazdım, böylece her zaman Jacobian'ı analitik olarak belirlemek zorunda değilim (mümkünse!). Ascher ve Petzold 1998'de verilen açıklamaları kullanarak, verilen x noktasında gradyanı belirleyen bu fonksiyonu …

13 pde  stability  numpy  scipy  newton-method 

Yoğun kare matrisleri ters çevirmeyi gerektiren diferansiyel denklemleri çözüyorum. Bu matris inversiyonu hesaplama süremin çoğunu tüketiyor, bu yüzden mevcut en hızlı algoritmayı kullanıp kullanmadığımı merak ediyordum. Şu anki seçimim numpy.linalg.inv . Nümeriklerimden , n'nin satır sayısı olduğu olarak ölçeklendiğini görüyorum, bu yüzden yöntem Gauss eliminasyonu gibi görünüyor.Ö ( n3)Ö(n3)O(n^3) Wikipedia'ya …

11 numpy  dense-matrix  inverse 

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 

Ben süperiletken-süperiletken birleşme noktasının akım-gerilim özelliklerini modellemeye çalışan bir fizikçiyim. Denklem bu model için geçerli: I(V)=1eRn−n∫∞−∞|E|[E2−Δ21]1/2|E+eV|[(E+eV)2−Δ22]1/2[f(E)−f(E+eV)]dEI(V)=1eRn−n∫−∞∞|E|[E2−Δ12]1/2|E+eV|[(E+eV)2−Δ22]1/2[f(E)−f(E+eV)]dE\begin{align} I(V) = \frac{1}{eR_{\mathrm{n-n}}}\int_{-\infty}^{\infty}\frac{|E|}{[E^{2} - \Delta_{1}^{2}]^{1/2}}\frac{|E + eV|}{[(E + eV)^{2} - \Delta_{2}^{2}]^{1/2}}[f(E) - f(E + eV)]\,\mathrm{d}E \end{align} Güncel (IIIveya Ikodda) değerler, bu integral verilen voltajlar için değerlendirilerek hesaplanır (VVVveya vkodda). Bunu Python'da denedim. Kod aşağıda …

9 python  numpy  integral-equations