«scipy» etiketlenmiş sorular

1
Newton-Krylov ne zaman uygun bir çözücü değildir?
Son zamanlarda scipy'den farklı lineer olmayan çözücüleri karşılaştırdım ve özellikle 20 kod satırında lineer olmayan reaksiyon terimi ile ikinci dereceden diferansiyel denklem denklemini çözdükleri Scipy Yemek Kitabı'ndaki Newton-Krylov örneğinden özellikle etkilendim . Formdaki yarı iletken heterostrüktürler için doğrusal olmayan Poisson denklemini ( Poisson-Boltzmann denklemi de denir, bu notlarda sayfa 17'ye …

1
Sonlu farklılıklara sahip yaklaşık bir Jacobian, Newton yönteminde kararsızlığa neden olabilir mi?
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 …

2
Python'da çok büyük ve seyrek bir bitişiklik matrisinin tüm özdeğerlerini hesaplamanın en hızlı yolu nedir?
Scipy.sparse.linalg.eigsh kullanmaktan çok büyük ve seyrek bir bitişiklik matrisinin tüm özdeğerlerini ve özvektörlerini hesaplamanın daha hızlı bir yolu olup olmadığını anlamaya çalışıyorum. matrisin simetri özellikleri. Bir bitişiklik matrisi de ikiliktir, bu da beni daha hızlı bir yol olduğunu düşündürüyor. Rastgele 1000x1000 seyrek bitişiklik matrisi oluşturdum ve x230 ubuntu 13.04 dizüstü …

1
Matris Dengeleme Algoritması
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. …
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.