Hesaplamalı Bilim

Varsayalım yinelemeli olarak tanımlanan bir matris denklemim A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Daha sonra A [1] denklemi, sürekli yeniden hesaplamayı önleyen bazı yüksek verimli yöntemlerin bulunduğu sürekli bir fraksiyona benzer (Bazı örnekler için "Sayısal Tarifler" kısmına bakın). Bununla birlikte, b [n] ve a [n] katsayılarının matris olmasına izin …

18 algorithms 

Masaüstüme bir Nvidia GT660 grafik kartı taktım ve bazı mücadelelerden sonra R ile arayüz kurmayı başardım. GPU'ları, özellikle gputools'u kullanan birkaç R paketiyle oynuyordum ve bazı temel işlemleri yapmak için GPU ve CPU'mun aldığı zamanı karşılaştırıyordum: ters matrisler (CPU daha hızlı) qr ayrışması (CPU daha hızlı) büyük korelasyon matrisleri (CPU …

18 r  gpu 

ODE sistemi için bir IVP düşünün y′=f(x,y)y′=f(x,y)y'=f(x,y) , y(x0)=y0y(x0)=y0y(x_0)=y_0 . Bu iş en yaygın sorun olarak kabul edilir katı zaman Jacobi matris ∂f∂y( x0, y0)∂f∂y(x0,y0)\frac{\partial f}{\partial y}(x_0,y_0)sahipher ikiçok küçük negatif gerçek kısmı ile çok büyük negatif gerçek parçanın ve özdeğerler (sadece sabit bir durumda dikkate) ile öz. Öte yandan, sadece …

18 ode  stiffness 

Sonunda yeni keşfedilen moleküler varlığınızda atomların mekansal olarak nasıl düzenlendiğini bulmayı başardınız. Örneğin, spektroskopik yollarla, bir grup atom koordinatlarına, atom tiplerine, bağ uzunluklarına, bağ tiplerine ve molekülünüz için bir şeye sahip olursunuz. Artık molekülünüzün nokta grubunu (simetri grubu) belirlemekle ilgileniyorsunuz. Metan (gibi basit moleküller için TdTdT_d ) ya da benzen …

18 computational-chemistry 

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 

Akışta bulunan biyolojik maddelerle akışkan akışını simüle etmek için bir kod geliştiriyorum. Bu, bazı ilave biyolojik modellere bağlı standart Navier-Stokes denklemlerini içerir. Birçok parametre / sabit vardır. Büyük hesaplamaları işlemek için yazılı işlevlerim var, ama yaşadığım bir sorun bu hesaplamaların bağlı olduğu sabitlerin / parametrelerin çok sayıda. Bir işleve 10-20 …

17 c++ 

İki farklı BLAS uygulaması göz önüne alındığında, tam olarak aynı kayan nokta hesaplamaları yapmasını ve aynı sonuçları döndürmesini bekleyebilir miyiz? Veya örneğin skaler bir ürünü ve bir ( x 1 y 1 + x 2 y) şeklinde hesaplayabilir mi? 2 ) + ( x 3 y 3 + x 4((x1y1+x2y2)+x3y3)+x4y4((x1y1+x2y2)+x3y3)+x4y4 …

17 floating-point  blas 

Fortran'ın hala yaygın olarak kullanıldığı, havacılık mühendisliği alanıyla ilgilenen bir Makine Mühendisliği öğrencisiyim. Öğrenmek için zamanımı hangi Fortran sürümüne yatırmalıyım?

17 fortran 

Sadece yüksek dereceli (yani 4'ten büyük) Runge – Kutta yöntemlerinin neden neredeyse hiç tartışılmadığını (en azından bilgime göre) kullanıldığını merak ettim. Adım başına daha büyük hesaplama süresi gerektirdiğini anlıyorum (örneğin 12. sıradaki gömülü adımla RK14 ), ancak daha yüksek dereceli Runge-Kutta yöntemlerini kullanmanın başka olumsuzlukları var mı (örn. Kararlılık sorunları)? …

17 ode  runge-kutta 

Nesneye yönelik programlamayı en iyi şekilde kullanmayı kendime öğretmenin bir yolu olarak seyrek matris hesaplamaları için küçük bir kütüphane yazıyorum. Parçaların (seyrek matrisler ve bağlantı yapılarını tanımlayan grafikler) çok gevşek bir şekilde birleştiği güzel bir nesne modeline sahip olmak için çok çalıştım. Kendi görüşüme göre, kod bunun için çok daha …

17 linear-algebra  sparse-matrix  oop 

Süreksiz galerkin yönteminde çözümleri temsil etmek için iki genel yaklaşım vardır: düğüm ve kip. Modal : Çözümler, bir dizi polinom kümesiyle çarpılan modal katsayıların toplamı ile temsil edilir, örn. burada genellikle dik polinomlardır , örneğin Legendre. Bunun bir avantajı, ortogonal polinomların diyagonal bir kütle matrisi üretmesidir.u ( x , t …

17 fluid-dynamics  discontinuous-galerkin 

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 

Daha sonra Paraview (.vtu veya vtk dosyası) tarafından okunacak bir hesaplama, bazı veriler yazmak gerekiyor. Dosya boyutu söz konusu olduğunda ASCII formatını mı yoksa Binary formatını mı seçmeliyim?

17 visualization  paraview  vtk 

Kafes oluşturma yazılımı arıyorum özgür ve açık kaynaklıdır, alan adı belirtimi için aklı başında bir komut dosyası arabirimi sağlar, karmaşık geometriler için çalışır, 2D ve 3D kafes üretebilir, Hangi seçeneklerim var?

17 mesh-generation 

Eşlenik gradyanın GMRES yönteminden çok daha iyi çalıştığı durumlarla ilgileniyorum. Genel olarak, CG birçok SPD vakasında tercih edilir (simetrik-pozitif-kesin) çünkü daha az depolama gerektirir ve CG için yakınsama oranına teorik olarak bağlı olması GMRES'in iki katıdır. Bu oranların gerçekte gözlemlendiği herhangi bir sorun var mı? Aynı sayıda spmv (seyrek matris-vektör …

17 linear-solver  conjugate-gradient  gmres