«python» etiketlenmiş sorular

Oktav'da iki vektörün Öklid mesafesini hesaplamanın hızlı bir yolu olup olmadığını bilmek istiyorum. Görünüşe göre bunun için özel bir fonksiyon yok, bu yüzden sadece formülü kullanmalı sqrtmıyım?

18 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 

Hangisi daha iyi: FORTRAN veya Python? Sanırım her iki durumda da Gnuplot'a ihtiyacınız var, değil mi? Şu anda bir Windows makinesinde çalışıyorum. Monte-Carlo simülasyonları, sayısal entegrasyon ve farklılaşma, moleküler dinamikler gibi fizik problemlerine sayısal çözümler elde etmek için kullanmak istiyorum. Hesaplamalı fizik üzerine hem FORTRAN'ı (77 inanıyorum) hem de Python'u …

17 python  fortran 

Son zamanlarda scipy'den farklı lineer olmayan çözücüleri karşılaştırdım ve özellikle 20 kod satırında lineer olmayan reaksiyon terimi ile ikinci dereceden diferansiyel denklem denklemini çözdükleri Scipy Yemek Kitabı'ndaki Newton-Krylov örneğinden özellikle etkilendim . Formdaki yarı iletken heterostrüktürler için doğrusal olmayan Poisson denklemini ( Poisson-Boltzmann denklemi de denir, bu notlarda sayfa 17'ye …

16 python  nonlinear-equations  poisson  scipy 

Bilimsel çalışmalarım için ticari yazılımlardan bağımsız olmak istiyorum. Matlab ve araç kutuları gibi ticari bir paketleri tatmin edici bulmuyorum, çünkü gelecekte Matlab'a erişip erişemeyeceğimi bilmiyorum ve dili sevmiyorum. Bu nedenle, alternatifler arıyorum. Neyse ki, Python'da oldukça akıcıyım (ve dili seviyorum) ve NumPy, SciPy, Matplotlib, Basemap ve NetCDF okuma ve yazma …

16 python  machine-learning 

CFD'de tüm ödevlerimi yapmak için programlama dili olarak python kullanmaya başladım. Programlama konusunda çok az deneyimim var. Makine mühendisliği geçmişindeyim ve Havacılık ve Uzay mühendisliğinde yüksek öğrenim görüyorum. bazen CFD'nin hesaplama yönü denklemleri manipüle etmekten veya matematiği yapmaktan daha sıkıcı olur. Programımızı daha hızlı çalıştırmayı sağlayan genel yönergeler nelerdir? Şeyleri …

16 algorithms  python  fluid-dynamics 

Araştırma grubum, daha sonra analiz edilmesi gereken tek bir yörüngenin parçası olarak gigabaytlarca veri üretebilen moleküler dinamiklere odaklanıyor. İlgilendiğimiz sorunların birçoğu veri kümesindeki korelasyonları içerir, bu da daha ardışık bir yaklaşım kullanmak yerine büyük miktarda veriyi hafızada izlememiz ve analiz etmemiz gerektiği anlamına gelir. Bilmek istediğim şey, büyük veri kümelerinin …

15 python  c  efficiency 

Hesaplamalı fizik bağlamında, büyük seyrek matrislerin özdeğer hesaplamaları gerçekleştirmek için bazı büyük kod geliştiriyorum. Özdeğerler analitik olarak iyi bilindiği için rutinlerimi basit harmonik osilatöre karşı bir boyutta test ediyorum. Bunu yaparak ve kendi rutinlerimi SciPy'nin dahili çözücüleriyle karşılaştırarak, aşağıdaki çizimde gösterilen garipliğe rastladım. Burada ilk 100 sayısal olarak hesaplanmış özdeğer …

15 python  sparse-matrix  numpy  eigenvalues 

Burada kimse Python aracılığıyla yeni nesil (örn. K20) GPU'larla çift duyarlıklı bilimsel hesaplama kullandı mı? Bu teknolojinin hızla geliştiğini biliyorum, ancak şu anda bunu yapmanın en iyi yolu nedir? GPU, popüler bilimsel Python kütüphaneleri için numpy ve scipy kapsamı dışında ve ben theano kullanmak istedim ama GPU için sadece float32 …

14 python  gpu  precision 

Python'da aşağıdaki ifadeyi uygulamak istiyorum: burada x ve y n boyutunda nümerik dizilerdir ve k nümerik bir dizidir n × n boyutunda . N boyutu yaklaşık 10000 olabilir ve fonksiyon birçok kez değerlendirilecek bir iç döngünün parçasıdır, bu nedenle hız önemlidir.xben= ∑j = 1i - 1ki - j , jbiri …

14 python  numpy 

Tekrarlanabilir bir kağıt yazıyorum ve kağıt bir Python komut dosyası (benzer bir MATLAB komut dosyası neredeyse aynı sonuçları üretir) tarafından oluşturulan hesaplama sonuçları vardır. Kağıdın içindeki hesaplamaları koddaki hesaplarla eşleştirebilseler, makalenin okuyucular için daha kolay anlaşılabileceğini hissediyorum. Çalışma soyut bir formalizm öneriyor ve makaledeki örneklerin bu formalizmi okuyucular için (birçoğu …

14 python  software  reproducibility 

Bir işlevi sayısal olarak uygulamakta sorun yaşıyorum. Büyük giriş değerlerinde sonucun çok büyük bir sayının çok küçük bir sayının çarpımı olduğu gerçeğinden muzdariptir. Felaket iptalinin doğru terim olup olmadığından emin değilim, lütfen öyleyse düzeltin. Bir şeyin ters gittiğine dair kanıtlar: 6'nın daha büyük girdileri için salınımları ve 0.0 değerini atamayı …

13 python  numerical-analysis  floating-point 

Matlab'ın sayısal entegrasyonu Scipy ile karşılaştırması konusunda biraz hayal kırıklığı yaşıyorum. Aşağıdaki test kodumda aşağıdaki farklılıkları gözlemliyorum: Matlab'ın versiyonu python eşdeğümden ortalama 24 kat daha hızlı çalışıyor ! Matlab'ın sürümü, integrali uyarmadan hesaplayabilirken, python döner nan+nanj Bahsedilen iki noktaya göre python'da aynı performansı almamı sağlamak için ne yapabilirim? Dokümantasyona göre, …

13 matlab  python  quadrature  numba 

Hem bir matrisin sırasını hem de bir matris denklem sisteminin çözümünü hesaplamak için farklı yöntemler kullanıyorum. Linalg.svd işleviyle karşılaştım. Bunu, Gauss Eliminasyonu ile sistemi çözme çabamla karşılaştırdığımda, hem daha hızlı hem de daha kesin görünüyor. Bunun nasıl mümkün olduğunu anlamaya çalışıyorum. Bildiğim kadarıyla, linalg.svd işlevi, matrisimin özdeğerlerini hesaplamak için bir …

13 linear-algebra  python  lapack 

Scipy.sparse.linalg.eigsh kullanmaktan çok büyük ve seyrek bir bitişiklik matrisinin tüm özdeğerlerini ve özvektörlerini hesaplamanın daha hızlı bir yolu olup olmadığını anlamaya çalışıyorum. matrisin simetri özellikleri. Bir bitişiklik matrisi de ikiliktir, bu da beni daha hızlı bir yol olduğunu düşündürüyor. Rastgele 1000x1000 seyrek bitişiklik matrisi oluşturdum ve x230 ubuntu 13.04 dizüstü …

12 linear-algebra  python  performance  eigensystem  scipy 

Çözmem gerek 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} Ben düşünüyorum onunla çözülebilir olması gereken bir kuadratik sorundur CVXOPT , ama nasıl dışarı çalışamaz.

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