«python» 태그된 질문

FORTRAN 또는 Python 중 어느 것이 더 낫습니까? 그리고 두 경우 모두 Gnuplot이 필요하다고 생각합니다. 현재 Windows 컴퓨터에서 작업하고 있습니다. Monte-Carlo 시뮬레이션, 수치 적분 및 미분, 분자 역학 등 물리 문제에 대한 수치 솔루션을 얻는 데 사용하고 싶습니다. FORTRAN (믿습니다)과 Python을 모두 소개하는 계산 물리학 과정을 봤습니다. 하나를 시작한 다음 …

17 python  fortran 

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

16 python  nonlinear-equations  poisson  scipy 

옥타브에서 두 벡터의 유클리드 거리를 계산하는 빠른 방법이 있는지 알고 싶습니다. 특별한 기능이없는 것 같습니다. 따라서 수식을 사용해야 sqrt합니까?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

과학적인 작업을 위해 상용 소프트웨어와 독립하고 싶습니다. Matlab과 같은 상용 패키지 및 도구 상자와 같은 상용 패키지가 만족스럽지 못하다는 것을 알게되었습니다. 왜냐하면 나중에 Matlab에 액세스 할 수 있는지 알지 못하기 때문에 언어를 좋아하지 않기 때문입니다. 따라서 대안을 찾고 있습니다. 다행스럽게도 저는 파이썬에 능통하고 언어를 좋아하며 NumPy, SciPy, Matplotlib, Basemap 및 …

16 python  machine-learning 

CFD에서 모든 과제를 수행하기 위해 프로그래밍 언어로 파이썬을 사용하기 시작했습니다. 프로그래밍 경험이 거의 없습니다. 저는 기계 공학 분야 출신이며 항공 우주 공학 분야의 고등 교육을 추구하고 있습니다. 때때로 CFD의 컴퓨팅 측면은 방정식을 조작하거나 수학을 수행하는 것보다 지루합니다. 프로그램을보다 빠르게 실행시키는 일반적인 지침은 무엇입니까? 병렬 처리를 수행하는 요령은 무엇입니까? 더 빨리 …

16 algorithms  python  fluid-dynamics 

저의 연구 그룹은 분자 역학에 중점을 두었습니다. 분자 역학은 분명히 단일 궤도의 일부로 기가 바이트의 데이터를 생성하여 분석해야합니다. 우리가 우려하는 몇 가지 문제는 데이터 세트의 상관 관계와 관련이 있습니다. 즉, 순차적 인 접근 방식을 사용하지 않고 메모리에서 많은 양의 데이터를 추적하고 분석해야합니다. 내가 알고 싶은 것은 큰 데이터 세트의 I …

15 python  c  efficiency 

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

15 python  sparse-matrix  numpy  eigenvalues 

여기 누구든지 파이썬을 통해 차세대 (예 : K20) GPU와 함께 배정도 과학 컴퓨팅을 사용 했습니까? 이 기술이 빠르게 발전하고 있음을 알고 있지만 현재이를 수행하는 가장 좋은 방법은 무엇입니까? GPU는 인기있는 과학 Python 라이브러리 numpy 및 scipy의 범위를 벗어 났으며 theano를 사용하고 싶었지만 GPU에는 float32 정밀도 만 사용하는 것 같습니다. 나는 …

14 python  gpu  precision 

파이썬에서 여기서 와 는 크기의 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 

나는 재현 가능한 논문을 작성하고 있는데,이 논문은 파이썬 스크립트에 의해 생성 된 계산 결과를 가지고 있습니다 (유사한 MATLAB 스크립트는 거의 동일한 결과를 생성합니다). 나는 독자들이 논문의 계산과 코드의 계산을 일치시킬 수 있다면 독자가 이해하기 쉽다고 생각합니다. 이 연구는 추상적 인 형식주의를 제안하고, 논문의 예는이 형식주의를 독자들 (많은 엔지니어가 될 것)에게보다 …

14 python  software  reproducibility 

matlab이 숫자 통합과 Scipy를 처리하는 방식에 약간의 좌절이 있습니다. 아래 테스트 코드에서 다음과 같은 차이점이 있습니다. Matlab의 버전은 내 파이썬 버전 보다 평균 24 배 빠릅니다 ! Matlab의 버전은 경고없이 정수를 계산할 수 있으며 파이썬은 반환합니다 nan+nanj 언급 된 두 가지 점과 관련하여 파이썬에서 동일한 성능을 얻으려면 어떻게해야합니까? 문서에 따르면 …

13 matlab  python  quadrature  numba 

함수를 숫자로 구현하는 데 문제가 있습니다. 큰 입력 값에서 결과는 매우 많은 수의 매우 작은 수라는 사실로 인해 어려움을 겪습니다. 치명적인 취소가 올바른 용어인지 확실하지 않으므로 수정하십시오. 무언가 잘못되었다는 증거 : 6의 더 큰 입력에 대해 진동과 0.0의 어설 션을 어떻게 피할 수 있습니까? 내 기능은 다음과 같습니다. import numpy …

13 python  numerical-analysis  floating-point 

행렬의 순위와 방정식의 행렬 시스템의 솔루션을 모두 계산하기 위해 다른 방법을 사용했습니다. linalg.svd 기능을 발견했습니다. 이것을 가우시안 제거 (Gaussian Elimination)로 시스템을 해결하려는 나의 노력과 비교하면 더 빠르고 정확한 것으로 보입니다. 이것이 어떻게 가능한지 이해하려고합니다. 내가 아는 한 linalg.svd 함수는 QR 알고리즘을 사용하여 행렬의 고유 값을 계산합니다. 나는 이것이 수학적으로 어떻게 …

13 linear-algebra  python  lapack 

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

12 linear-algebra  python  performance  eigensystem  scipy 

나는 해결해야한다 s.t.minx∥Ax−b∥22,∑ixi=1,xi≥0,∀i.minx‖Ax−b‖22,s.t.∑ixi=1,xi≥0,∀i.\begin{alignat}{1} & \min_{x}\|Ax - b\|^2_{2}, \\ \mathrm{s.t.} & \quad\sum_{i}x_{i} = 1, \\ & \quad x_{i} \geq 0, \quad \forall{i}. \end{alignat} 내가 생각하는 그것으로 풀 수 있어야한다 이차 문제 CVXOPT ,하지만 난 방법을 작동하지 않을 수 있습니다.

12 optimization  python  convex-optimization  least-squares  quadratic-programming