Hesaplamalı Bilim

Bilimsel sorunları çözmek için bilgisayar kullanan bilim adamlarının soru-cevapları

2
FEM'de, sertlik matrisi neden pozitiftir?
FEM sınıflarında, genellikle sertlik matrisinin pozitif kesin olduğu kabul edilir, ancak nedenini anlayamıyorum. Birisi açıklama yapabilir mi? Örneğin, Poisson problemini düşünebiliriz: sertlik matrisi: −∇2u=f,−∇2u=f, -\nabla^2 u = f,Kij=∫Ω∇φi⋅∇φjdΩ,Kij=∫Ω∇φi⋅∇φjdΩ,K_{ij} = \int_\Omega\nabla\varphi_i\cdot\nabla\varphi_j\, d\Omega, ki bu simetrik ve pozitif tanımlı. Simetri bariz bir özelliktir, ancak pozitif kesinlik benim için çok açık değildir.


5
Log-log uzayında integral
Genel olarak, log-log uzayında çok daha düzgün ve daha iyi davranan işlevlerle çalışıyorum --- bu yüzden enterpolasyon / ekstrapolasyon vb. Gerçekleştiriyorum ve bu çok iyi çalışıyor. Bu sayısal işlevleri günlük kaydı alanına entegre etmenin bir yolu var mı? yani scipy.integrate.cumtrapzbazı st bulmak için bir kümülatif integral (örneğin python, kullanım ) …

1
Ters matrisin hesaplanması için en hızlı algoritma ve pozitif belirli simetrik matrisler için belirleyicisi nedir?
Pozitif belirli bir simetrik matris verildiğinde, ters matrisi ve determinantını hesaplamak için en hızlı algoritma nedir? İlgilendiğim problemler için matris boyutu 30 veya daha az. Yüksek doğruluk ve hız gerçekten gereklidir. (milyonlarca matris gerçekleştirilir) Her hesaplamada, ters matrisin sadece bir elemanı gereklidir. Teşekkürler!

3
Ağsız yöntemlerin ciddi kullanımına örnekler verebilir misiniz?
Hareketli En Küçük Kareler fonksiyonlarına dayanan Elementsiz Galerkin gibi ağsız yöntemleri kullanan bilimsel kodları ve ticari paketleri duymak istiyorum. "Ciddi" ile kastedilen, örneğin FEM tarafından çözülen boyutlarla karşılaştırılabilir sorunları çözmek için kullanılabilir. Kuruluşundan bu yana on beş yıldan fazla bir süre geçti. Onları geliştiren insanlar onları çok umut vaat ediyorlardı. …

3
N-boyutlu Delaunay Tesselation Yazılım Kütüphaneleri
N-boyutlu uzayda (N> = 2) düzensiz aralıklı bilinen bir dizi nokta / düğüm var ve bu noktaların Delaunay üçgenleme oluşturmak ve karşılık gelen öğeleri döndürmek için bir yol istiyorum. ND Delaunay üçgenleme yapacak mevcut ağ kütüphaneleri var mı? (Bunu yapıyorum, çünkü kafesli elemanları uzayda herhangi bir noktada doğrusal enterpolasyon için …

2
RBF çekirdek matrisleri koşulsuz olma eğilimindedir mi?
Bir çekirdek tabanlı makine öğrenme algoritması (KLPP) uygulamak için RBF çekirdek işlevini kullanıyorum, sonuçta ortaya çıkan çekirdek matrisi son derece kötü koşullu olarak gösterilmiştir. L2 normunun koşul sayısıKKK K(i,j)=exp(−(xi−xj)2σ2m)K(i,j)=exp⁡(−(xi−xj)2σm2)K(i,j)= \exp\left({\frac{-(x_{i}-x_{j})^2}{ \sigma_{m}^2}}\right)1017−10641017−106410^{17}-10^{64} İyi şartlandırmanın bir yolu var mı? Parametre ayarlanması gerekiyor, ama tam olarak nasıl olduğunu bilmiyorum.σσ \sigma Teşekkürler!

1
Çok elektronlu zamana bağlı Schrödinger denklemini sayısal olarak çözmek neden zordur
İnsanların çok elektronlu bir sistemle başa çıkmak için genellikle Tek Aktif Elektron (SAE) yaklaşımını kullandıkları ve sorunu tek bir elektron problemine dönüştürdüğü görülmektedir. Örneğin, bir helyum atomunun lazer alanlarıyla etkileşimi sorununu sayısal olarak çözerken, insanlar yaklaşık olarak yaklaşık olarak bir potansiyel potansiyeli olan elektron-elektron etkisini içerir ve esasen bir elektron …

2
Fortran'da çift kesinlik değerleri nasıl ayarlanır?
Son zamanlarda, FORTRAN95 ile garip bir sorunla karşılaştım. X ve Y değişkenlerini aşağıdaki gibi başlattım: X=1.0 Y=0.1 Daha sonra bunları bir araya getirip sonucu yazdırıyorum: 1.10000000149012 Değişkenler incelendikten sonra, 0.1 tam hassasiyetle çift hassasiyetle temsil edilmemiş gibi görünüyor. Bundan kaçınmamın bir yolu var mı?



3
LAPACK'in kullanmasının nedeni nedir
LAPACK'in QR rutini Q'yu Hanehalkı reflektörleri olarak depolar. Yansıma ile , böylece sonucun ilk öğesi olur , bu yüzden saklanması gerekmez. Ve gerekli ölçek faktörlerini içeren ayrı bir vektörü saklar . Yani bir reflektör matrisi şöyledir:vvv1 /v11/v11/v_1111ττ\tau'H= Ben- τvvT,'H=ben-τvvT,H=I-\tau v v^T, burada normalleştirilmez. Ders kitaplarında reflektör matrisivvv 'H= Ben- 2 …

1
Normal Runge – Kutta yöntemlerinin SDE'lere genelleştirilemeyeceği kolayca anlaşılabilir bir argüman mı?
Stokastik diferansiyel denklemlerin (SDE) çözülmesine naif bir yaklaşım: düzenli çok adımlı Runge – Kutta yöntemini kullanır, Altta yatan Wiener sürecinin yeterince ince bir takdir yetkisini kullanmak, Runge – Kutta yönteminin her adımını Euler – Maruyama'ya benzetir. Şimdi, bu birden fazla seviyede başarısız oluyor ve nedenini anlıyorum. Ancak şimdi Runge-Kutta yöntemleri …

4
Küçük kare matrisler için en hızlı lineer sistem çözme (10x10)
Bazen küçük matrisler olarak adlandırılan küçük matrisler (10x10) için lineer sistem çözmenin cehennemini optimize etmekle çok ilgileniyorum . Bunun için hazır bir çözüm var mı? Matrisin nonsüler olduğu varsayılabilir. Bu çözücü, Intel CPU üzerinde mikrosaniye cinsinden 100000 defadan fazla çalıştırılmalıdır. Bilgisayar oyunlarında kullanılan optimizasyon seviyesinden bahsediyorum. Meclise ve mimariye özgü …

2
Üçgen örgünün düzensizliğini ölçmek için yaygın olarak kullanılan metrikler
Düz bir düzlemde üçgen bir ağınız olduğunu varsayalım. Bu, nihayetinde mekanikteki bazı problemleri çözmek için çizilmiştir. Eşkenar üçgenler ağı, köşeler ve sentroidler arasındaki mesafeler her yerde aynı olduğu için en iyisidir. Bu enterpolasyonları ve degradelerin hesaplanmasını kolay ve doğru bir görev haline getirir. Bununla birlikte, kısıtlamalar ve koşullar nedeniyle, tüm …

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.