Запитання з тегом «numpy»

Я хочу змінити щільну квадратну матрицю переходу на місці, змінивши порядок кількох її рядків і стовпців, використовуючи бібліотеку numpy python. Математично це відповідає попередньому множенню матриці на перестановкову матрицю P і післяпомноження на P ^ -1 = P ^ T, але це не є обчислювально обґрунтованим рішенням. Зараз я обмінюю …

27 linear-algebra  python  numpy 

Як масштабувати масиви Python / Numpy зі збільшенням розмірів масиву? Це ґрунтується на деякій поведінці, яку я помітив під час тестування коду Python для цього питання: Як висловити цей складний вираз за допомогою nummy-фрагментів Проблема в основному полягала в індексації для заповнення масиву. Я виявив, що переваги використання (не дуже …

21 python  performance  numpy 

Я розробляю якийсь більший код для виконання обчислень власних значень величезних розріджених матриць в контексті обчислювальної фізики. Я перевіряю свої підпрограми проти простого гармонічного осцилятора в одному вимірі, оскільки власні значення добре відомі аналітично. Роблячи це і порівнюючи власну процедуру з вбудованими рішеннями SciPy, я зіткнувся з диваком, відображеним на …

15 python  sparse-matrix  numpy  eigenvalues 

Я хочу реалізувати наступний вираз у Python: де x і y - нумерові масиви розміром n , а k - масивний масив розміру n × n . Розмір n може становити приблизно 10000, а функція є частиною внутрішнього циклу, який буде оцінюватися багато разів, тому важлива швидкість.хi= ∑j = 1i …

14 python  numpy 

Я реалізував зворотний рішень-ейлер у python 3 (використовуючи numpy). Для власної зручності і як вправи я також написав невелику функцію, яка обчислює наближення градієнта кінцевою різницею, так що мені не завжди доведеться аналітично визначати якобіанські (якщо це навіть можливо!). Використовуючи описи, наведені в Ascher і Petzold 1998 , я написав …

13 pde  stability  numpy  scipy  newton-method 

Я розв'язую диференціальні рівняння, які потребують інвертування щільних квадратних матриць. Ця інверсія матриці займає більшу частину мого часу на обчислення, тому мені було цікаво, чи використовую я найшвидший доступний алгоритм. Мій поточний вибір - numpy.linalg.inv . З моїх числових цифр я бачу, що вона масштабується як де n - кількість …

11 numpy  dense-matrix  inverse 

Для скорочення моделі я хочу обчислити ліві сингулярні вектори, пов'язані з - скажімо, 20 - найбільшими сингулярними значеннями матриці , де N ≈ 10 6 та k ≈ 10 3 . На жаль, моя матриця A буде щільною без будь-якої структури.A ∈ RN, кА∈RN,кA \in \mathbb R^{N,k}N≈ 106N≈106N\approx 10^6k ≈ …

10 linear-algebra  python  reference-request  numpy  svd 

Я фізик, який намагається моделювати характеристики струму напруги з’єднання надпровідник-надпровідник. Рівняння для цієї моделі: Я( V) =1еRn - n∫∞- ∞| Е|[Е2-Δ21]1 / 2| Е+ e V|[ ( Е+ e V)2-Δ22]1 / 2[ ф( Є) - ф( Є+ e V) ]d EЯ(V)=1еRн-н∫-∞∞|Е|[Е2-Δ12]1/2|Е+еV|[(Е+еV)2-Δ22]1/2[f(Е)-f(Е+еV)]гЕ\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} …

9 python  numpy  integral-equations