«numpy» 태그된 질문

파이썬의 numpy 라이브러리를 사용하여 여러 행과 열의 순서를 변경하여 밀도가 높은 사각형 전이 행렬을 제자리에서 수정하고 싶습니다. 수학적으로 이것은 순열 행렬 P에 의해 행렬을 미리 곱하고 P ^ -1 = P ^ T에 의해 곱셈하는 것에 해당하지만, 이것은 계산적으로 합리적인 해결책은 아닙니다. 지금은 수동으로 행과 열을 교환하고 있지만 numpy에는 멋진 …

27 linear-algebra  python  numpy 

Python / Numpy 배열은 배열 크기가 증가함에 따라 어떻게 확장됩니까? 이것은이 질문에 대한 파이썬 코드를 벤치마킹 할 때 발견 한 행동을 기반으로합니다 : numpy slices를 사용 하여이 복잡한 표현을 표현하는 방법 문제는 주로 배열을 채우는 색인 생성과 관련이 있습니다. 파이썬 루프에 비해 (아주 좋지 않은) Cython 및 Numpy 버전을 사용하는 …

21 python  performance  numpy 

계산 물리학의 맥락에서 거대한 희소 행렬의 고유 값 계산을 수행하는 더 큰 코드를 개발 중입니다. 고유 값은 분석적으로 잘 알려져 있기 때문에 간단한 고조파 발진기에 대해 일차원으로 테스트합니다. 그렇게하고 내 자신의 루틴을 SciPy의 빌트인 솔버와 비교하면서 아래 그림에 표시된 이상한 점을 발견했습니다. 여기에서 처음 100 개의 숫자로 계산 된 고유 …

15 python  sparse-matrix  numpy  eigenvalues 

파이썬에서 여기서 와 는 크기의 numpy 배열 이며 크기의 numpy 배열입니다 . 크기 은 최대 약 10000 일 수 있으며이 함수는 여러 번 평가 될 내부 루프의 일부이므로 속도가 중요합니다.엑스나는= ∑j = 1난 − 1케이i − j , jㅏ나는 − jㅏ제이,엑스나는=∑제이=1나는−1케이나는−제이,제이ㅏ나는−제이ㅏ제이, x_i = \sum_{j=1}^{i-1}k_{i-j,j}a_{i-j}a_j, 엑스엑스x와이와이y엔엔n케이케이kn × n엔×엔n\times n엔엔n 이상적으로는 for …

14 python  numpy 

파이썬 3에서 numpy를 사용하여 Back-Euler 솔버를 구현했습니다. 나 자신의 편의와 운동을 위해 그래디언트의 유한 차분 근사값을 계산하는 작은 함수를 작성하여 항상 야곱을 분석적으로 결정할 필요가 없도록했습니다. Ascher와 Petzold 1998에 제공된 설명을 사용하여 주어진 지점 x에서 기울기를 결정하는이 함수를 작성했습니다. def jacobian(f,x,d=4): '''computes the gradient (Jacobian) at a point for a …

13 pde  stability  numpy  scipy  newton-method 

밀도가 높은 정방 행렬을 반전시키는 데 필요한 미분 방정식을 풀고 있습니다. 이 행렬 반전은 계산 시간을 가장 많이 소비하므로 사용 가능한 가장 빠른 알고리즘을 사용하고 있는지 궁금합니다. 내 현재 선택은 numpy.linalg.inv 입니다. 내 숫자에서 나는 그것이 n ( 행 수 ) 인 로 확장되는 것을 알기 O(n3)O(n3)O(n^3)때문에 방법은 가우시안 제거 …

11 numpy  dense-matrix  inverse 

모형 축소를 위해 행렬 의 20 개 가장 큰 특이 값과 연관된 왼쪽 특이 벡터를 계산하려고합니다 . 여기서 N ≈ 10 6 및 k ≈ 10 3 입니다. 불행히도 내 행렬 A 는 아무런 구조없이 밀도가 높습니다.∈ R엔, kㅏ∈아르 자형엔,케이A \in \mathbb R^{N,k}엔≈ 106엔≈106N\approx 10^6k ≈ 10삼케이≈10삼k\approx 10^3ㅏㅏA 이 크기의 …

10 linear-algebra  python  reference-request  numpy  svd 

저는 초전도체-초전도체 접합부의 전류-전압 특성을 모델링하려고하는 물리학 자입니다. 이 모델 의 방정식 은 다음과 같습니다. 나는( V) =1이자형아르 자형n - n∫∞− ∞| 이자형|[이자형2−Δ21]1 / 2| 이자형+ 전자 V|[ ( E+ 전자 V)2−Δ22]1 / 2[ f( E) − f( E+ 전자 V) ]d 전자I(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 …

9 python  numpy  integral-equations