«scipy» 태그된 질문

최근에 나는 scipy의 다른 비선형 솔버를 비교해 왔으며 특히 Scipy Cookbook 의 Newton-Krylov 예제 에서 약 20 줄의 코드로 비선형 반응 항으로 2 차 미분 방정식 방정식을 풀었습니다 . 예제 코드를 수정하여 반도체 이종 구조에 대한 비선형 Poisson 방정식 ( Poisson-Boltzmann 방정식이라고도 함, 이 노트의 17 페이지 참조)을 해결합니다. 디2ϕ디엑스2− …

16 python  nonlinear-equations  poisson  scipy 

파이썬 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 

scipy.sparse.linalg.eigsh를 사용하는 것보다 매우 크고 희소 한 인접 행렬의 모든 고유 값과 고유 벡터를 계산하는 더 빠른 방법이 있는지 알아 내려고 노력 중입니다. 행렬의 대칭 속성. 인접 행렬도 바이너리이므로 더 빠른 방법이 있다고 생각합니다. 임의의 1000x1000 희소 인접 행렬을 만들고 x230 우분투 13.04 랩톱의 여러 방법을 비교했습니다. scipy.sparse.linalg.eigs : 0.65 …

12 linear-algebra  python  performance  eigensystem  scipy 

제어 시스템 도구 상자를 처음부터 순수하게 Python3 (shameless plug :)으로 작성했습니다 harold. 과거의 연구에서 나는 care.m기술적 / 관련이없는 이유로 Riccati 솔버 에 대해 항상 불만을 제기했습니다 . 따라서 나는 내 자신의 루틴 세트를 작성해 왔습니다. 내가 찾을 수없는 한 가지 방법은 적어도만큼 좋은 고성능 밸런싱 알고리즘을 얻는 것 balance.m입니다. 언급하기 …

9 linear-algebra  matlab  lapack  scipy  matrix-equations