Saya ingin memodifikasi matriks transisi persegi yang padat di tempat dengan mengubah urutan beberapa baris dan kolomnya, menggunakan perpustakaan numpy python. Secara matematis ini berhubungan dengan pra-mengalikan matriks dengan matriks permutasi P dan mengalikannya dengan P ^ -1 = P ^ T, tetapi...
Bagaimana skala array Python / Numpy dengan dimensi array yang meningkat? Ini didasarkan pada beberapa perilaku yang saya perhatikan ketika membandingkan kode Python untuk pertanyaan ini: Cara mengekspresikan ekspresi rumit ini menggunakan irisan numpy Masalahnya sebagian besar melibatkan...
Saya sedang mengembangkan beberapa kode yang lebih besar untuk melakukan perhitungan nilai eigen dari matriks jarang besar, dalam konteks fisika komputasi. Saya menguji rutinitas saya terhadap osilator harmonik sederhana dalam satu dimensi, karena nilai eigen dikenal secara analitis. Melakukannya...
Saya ingin menerapkan ekspresi berikut dalam Python: mana x dan y adalah array numpy dengan ukuran n , dan k adalah array numpy ukuran n × n . Ukuran n mungkin hingga sekitar 10.000, dan fungsinya adalah bagian dari loop dalam yang akan dievaluasi berkali-kali, jadi kecepatan itu penting.xsaya=...
Saya telah menerapkan pemecah Euler mundur dalam python 3 (menggunakan numpy). Untuk kenyamanan saya sendiri dan sebagai latihan, saya juga menulis fungsi kecil yang menghitung perkiraan perbedaan hingga dari gradien sehingga saya tidak selalu harus menentukan Jacobian secara analitis (jika itu...
Saya memecahkan persamaan diferensial yang perlu membalik matriks kuadrat padat. Pembalikan matriks ini menghabiskan sebagian besar waktu perhitungan saya, jadi saya bertanya-tanya apakah saya menggunakan algoritma tercepat yang tersedia. Pilihan saya saat ini adalah numpy.linalg.inv . Dari angka...
Untuk reduksi model, saya ingin menghitung vektor singular kiri terkait dengan - katakanlah 20 - nilai singular terbesar dari matriks , di mana N ≈ 10 6 dan k ≈ 10 3 . Sayangnya, matriks A saya akan padat tanpa struktur apa pun.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx...
Saya seorang ahli fisika yang mencoba memodelkan karakteristik tegangan-arus dari persimpangan superkonduktor-superkonduktor. The Persamaan untuk model ini adalah: saya( V) = 1e Rn - n∫∞- ∞| E|[ E2- Δ21]1 / 2| E+ e V|[ ( E+ e V)2- Δ22]1 / 2[ f( E) - f( E+ e V) ]d