Hesaplamalı Bilim

Ben çok soru özellikle uygun değildir eğer öyleyse beni affet programlama GPGPU yeni. Anladığım kadarıyla GPU programlama, olağan CPU programlamasına kıyasla çok karmaşık bir mühendislik çalışmasıdır. Diverjans sorunları, döşeme, sabitlenmiş bellek tahsisi ve ana bilgisayar aygıt iletişimi / aygıt hesaplama çakışması konusunda çok dikkatli olmak gerekir. Biraz araştırma yaptıktan sonra …

10 parallel-computing  petsc  gpu  cuda 

Zorlama teriminin ( formülasyon için aşağıdaki Düzenleme 1'e bakınız) ve ve bunun ilk türevlerine bağlı olabileceği aşağıdaki problemi düşünün . Bu 1 + 1 boyutlu dalga denklemidir. reçete edilen ilk verilerimiz var .Wuv=FWuv=F W_{uv} = F W { u + v = 0 }u,vu,vu,vWWW{u+v=0}{u+v=0}\{u+v = 0\} aralığında bağımlılık alanı içindeki …

10 finite-difference  reference-request  hyperbolic-pde 

Üniversitede ödev olarak 1d simülasyon yaptım. Sorun ifadesi, çalışma sıvısı olarak sıkıştırılabilir ideal gazı içeren 1d şok tüpü problemini çözmekti. Bu problem için, Roe Riemann çözücüsünü kullanarak Eulers denklemleri sistemini çözdüm. Bilmek istiyorum, Euler denklemlerini 2 veya 3 boyutta çözmek için nereden başlamalıyım? Test problemi hangisi, önce düşünmeliyim? (Lütfen ticari …

10 fluid-dynamics 

Şu anda temel-spline kollokasyonu adı verilen bir diferansiyel denklem çözme yöntemi üzerinde çalışıyorum. Ne ile sorun yaşıyorum ile ilişkisi ile keyfi bir sipariş spline oluşturmak için bir yöntem inşa etmektir ile başlangıç ​​koşulu ve bu sorunla başlamakta bile sorun yaşıyorum, çünkü özyinelemeli "üst" veya "alt" ile başlayabilir ve genel bir …

10 c++  b-spline 

Bazı yoğun, koşulsuz matrisleri köşegenleştirmeye çalışıyorum. Makine hassasiyetinde sonuçlar yanlıştır (negatif özdeğerleri döndürmek, özvektörlerde beklenen simetriler yoktur). Rasgele hassasiyetten yararlanmak için Mathematica'nın Eigensystem [] işlevine geçtim, ancak hesaplamalar son derece yavaş. Ben çok sayıda çözüme açığım. Koşulsuz sorunlara uygun paketler / algoritmalar var mı? Ön koşullandırma konusunda uzman değilim, bu …

10 linear-algebra  parallel-computing  eigensystem  precision  dense-matrix 

Turing'in reaksiyon difüzyon sistemini aşağıdaki C ++ koduyla çözüyorum. Çok yavaş: 128x128 piksel doku için kabul edilebilir yineleme sayısı 200'dür - bu 2,5 saniyelik gecikme ile sonuçlanır. İlginç bir görüntü elde etmek için 400 yinelemeye ihtiyacım var - ancak 5 saniye beklemek çok fazla. Ayrıca, dokunun boyutu aslında 512x512 olmalıdır …

10 pde  stiffness 

PETSc'yi harici bir BLAS / LAPACK kütüphanesi ile derlemek seyrek matrislerdeki performansı önemli ölçüde etkiler mi, yoksa sadece bu kütüphaneleri yoğun matris matematiği için mi kullanır?

10 petsc  sparse-matrix  lapack  blas 

Matris çarpımı (hem Mat * Mat hem de Mat * Vec) sıfır olmayan sayılarla veya matrisin boyutuyla ölçeklendiriliyor mu? Veya ikisinin bir kombinasyonu. Peki ya şekil. Örneğin, içinde 100 değeri olan bir 100 x 100 matrisi veya içinde 100 değeri olan bir 1000 x 1000 matrisi var. Bu matrisleri karelerken …

10 linear-algebra  performance  sparse-matrix 

İki referans sistemi (eksen) arasında koordinat dönüşümü yapmam gerekiyor. Bunun için, bazı ara eksenlerin kullanılmasından dolayı üç matrisin ( ) çarpılması gerekir. Bunu çözmek için iki yaklaşım düşündüm:3×33×33\times3 Yöntem # 1 : doğrudan çarpma yapmak olduğunu vf=R1 R2 R3 vivf=R1 R2 R3 viv_f = R_1\ R_2\ R_3\ v_i Yöntem # …

10 matrix  performance  matlab 

Eşyaları eşleştirmeye çalışıyorum. Bir dizi öğesi göz önüne alındığında , 0 ile 100 arasında bir skalaya benzerlik gösterebilirim. Örneğin, n 1 maddesi süt ve n 2 maddesi de sütse, n 1 ve n 2 arasındaki benzerlik % 100 olacaktır. Eğer n, 3 soda, sonra arasındaki benzerlik , n , 1 …

10 regression  data-analysis 

Let üçgen olmak ve izin f üzerinde düzgün bir fonksiyonu T .TTTfffTTT Orta nokta kareselini kullanabiliriz ; burada x M , T'nin orta noktasıdır .∫fdx ≈ | T| ⋅f( xM)∫fdx≈|T|⋅f(xM)\int f dx \approx |T|\cdot f(x_M)xMxMx_MTTT Bana bir simpleks üzerinde daha yüksek dereceli formüller sağlayabilir misiniz?

10 quadrature 

Python'da tam bir SVM uygulaması yazmaya çalışıyorum ve Lagrange katsayılarını hesaplamak için birkaç sorunum var. Önce doğru yolda olduğumdan emin olmak için algoritmadan anladığımı tekrar yazmama izin verin. Eğer isimli bir veri seti ve sınıf etikettir sonrax1,x2,...,xnx1,x2,...,xnx_1, x_2, ..., x_nx i ∀ i , y i ( w T x …

10 optimization  machine-learning  support-vector-machines  quadratic-programming 

Optimizasyon için Wikipedia'dan : Bilgisayar biliminde, meta-sezgisel bir kaliteyi belirli bir kalite ölçüsüne göre aday bir çözümü geliştirmeye çalışarak bir problemi optimize eden bir hesaplama yöntemi belirler . Meta-sezgiler, optimize edilen problem hakkında az ya da hiç varsayımda bulunmaz ve çok geniş aday çözüm alanlarını arayabilir. Ancak, meta-sezgisel tarama şimdiye …

10 optimization  heuristics 

Bazı seyrek matrisin özdeğerlerini tekrarlanan yöntemle [a, b] belirli aralıklarla nasıl bulacağımı merak ediyorum. Kişisel anlayışımla, Krylov altuzay yöntemini, iç özdeklerden ziyade aşırı özdeğerleri bulmak için kullanmak daha açıktır.

10 linear-algebra 

Farzedin ki, doğrusal bir konjügasyon jeneratörü tarafından rastgele bir sayı dizisi üretildi. Yani, xn+1=(aXn+c)modmxn+1=(aXn+c)modmx_{n+1}=(aX_n+c) \bmod m Tüm dönemi (veya en azından bitişik olarak büyük bir ardışık) verirsem , bu diziyi üreten ve x 0 parametrelerini nasıl yeniden yapılandırabilirim ? Sahte rasgele sayı üreteci biliniyorsa, başlangıç ​​parametrelerini belirlemek mümkün olacak genel …

10 random-number-generation  inverse-problem