Ilmu Komputasi

Yang merupakan perpustakaan tercepat untuk melakukan triangulasi set delaunay dengan jutaan jika poin 3D? Apakah ada juga versi GPU yang tersedia? Dari sisi lain, memiliki voronoi tessellation dari set poin yang sama, akan membantu (dalam hal kinerja) untuk mendapatkan triangulasi delaunay?

26 computational-geometry  delaunay-triangulation  voronoi-diagrams 

Apakah ada alat yang dirancang khusus untuk mengompresi data ilmiah floating point? Jika suatu fungsi halus, jelas ada banyak korelasi antara angka-angka yang mewakili fungsi itu, sehingga data harus dikompres dengan baik. Zipping / gzipping data floating point biner tidak kompres dengan baik. Saya bertanya-tanya apakah ada metode yang khusus …

26 data-storage 

Apa itu algoritma sederhana untuk menghitung SVD matriks?2 × 22×22 \times 2 Idealnya, saya ingin algoritma yang kuat secara numerik, tetapi saya ingin melihat implementasi yang sederhana dan yang tidak terlalu sederhana. Kode C diterima. Adakah referensi untuk makalah atau kode?

26 linear-algebra  matrix  reference-request  svd  matrix-factorization 

Pengembang perangkat lunak memiliki pilihan untuk memilih lisensi yang sesuai sesuai dengan tujuan pekerjaan. Adakah yang bisa memberikan rekomendasi / pengalaman tentang lisensi mana yang harus dipilih untuk perangkat lunak? Apa pro / kontra dari "memberikan" semua pekerjaan yang dikodekan sebagai kode sumber terbuka? Bagaimana cara berurusan dengan pemain industri …

26 software 

Apa yang bisa salah ketika menggunakan metode Krylov yang dikondisikan sebelumnya dari KSP ( paket solver linier PETSc ) untuk menyelesaikan sistem linier yang jarang seperti yang diperoleh dengan mengkritisi dan membuat linierisasi persamaan diferensial parsial? Langkah apa yang bisa saya ambil untuk menentukan apa yang salah untuk masalah saya? …

26 linear-algebra  petsc  preconditioning  krylov-method  iterative-method 

Pertimbangan apa yang harus saya buat ketika memilih antara BFGS dan konjugasi gradien untuk optimasi? Fungsi yang saya coba cocokkan dengan variabel-variabel ini adalah fungsi eksponensial; Namun, fungsi obyektif yang sebenarnya melibatkan integrasi, antara lain, dan sangat mahal jika itu membantu sama sekali.

25 optimization  conjugate-gradient 

Bagaimana masalah gravitasi n-tubuh diselesaikan secara numerik secara paralel? Apakah tradeoff presisi-kompleksitas mungkin? Bagaimana presisi mempengaruhi kualitas model?

25 algorithms  ode  parallel-computing  complexity  precision 

Saya perlu mengevaluasi secara integral angka-angka di bawah ini: ∫∞0sinc′(xr)rE(r)−−−−√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr di mana , dan . Di sini adalah fungsi Bessel yang dimodifikasi dari jenis kedua. Dalam kasus khusus saya, saya memiliki , dan .E(r)=r4(λκ2+r2−−−−−−√)−ν−5/2K−ν−5/2(λκ2+r2−−−−−−√)E(r)=r4(λκ2+r2)−ν−5/2K−ν−5/2(λκ2+r2)E(r) = r^4 (\lambda\sqrt{\kappa^2+r^2})^{-\nu-5/2} K_{-\nu-5/2}(\lambda\sqrt{\kappa^2+r^2})x∈R+x∈R+x \in \mathbb{R}_+λ,κ,ν>0λ,κ,ν>0\lambda, \kappa, \nu >0KKKλ=0.00313λ=0.00313\lambda = 0.00313κ=0.00825κ=0.00825\kappa = …

25 matlab  quadrature  special-functions 

Diberikan fungsi nyata dari variabel nyata, apakah ada perangkat lunak yang tersedia yang dapat secara otomatis menghasilkan kode numerik-akurat untuk menghitung fungsi atas semua input pada mesin yang dilengkapi dengan aritmatika IEEE 754? Misalnya, jika fungsi sebenarnya yang akan dievaluasi adalah: Perangkat lunak akan mempertimbangkan pembatalan dahsyat dan mungkin pencarian …

25 software  floating-point  accuracy 

Saya mencari alat untuk memvisualisasikan grafik tautan searah yang sangat besar. Saat ini saya memiliki ~ 2 juta node dengan ~ 10 juta tepi. Saya telah mencoba beberapa hal yang berbeda, tetapi kebanyakan membutuhkan waktu berjam-jam untuk bahkan melakukan grafik simpul 100k Apa yang saya coba: Saya menghabiskan satu hari …

25 visualization 

Pertanyaan terkait: Keadaan Mac OS dalam Scientific Computing dan HPC Sejumlah besar paket perangkat lunak dalam ilmu komputasi ditulis dalam Fortran, dan Fortran tidak akan hilang. Kompiler Fortran juga diperlukan untuk membuat paket perangkat lunak lain (salah satu contohnya adalah SciPy ). Namun, Mac OS X tidak termasuk kompiler Fortran. …

25 software  fortran 

Dalam komentar ini saya menulis: ... integrator SciPy default, yang saya asumsikan hanya menggunakan metode symplectic. di mana saya merujuk ke SciPy's odeint, yang menggunakan metode "non-stiff (Adams)" atau "stiff (BDF) method". Menurut sumber : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, …

25 ode  time-integration  integration  differential-equations 

Saya tidak mengerti perilaku berbeda dari persamaan advection-difusion ketika saya menerapkan kondisi batas yang berbeda. Motivasi saya adalah simulasi kuantitas fisik nyata (kepadatan partikel) di bawah difusi dan adveksi. Kepadatan partikel harus dilestarikan di interior kecuali mengalir keluar dari tepi. Dengan logika ini, jika saya menerapkan syarat batas Neumann pada …

25 pde  boundary-conditions  advection  diffusion  crank-nicolson 

Masalah yang sangat umum di Markov Chain Monte Carlo melibatkan probabilitas komputasi yang merupakan jumlah dari istilah eksponensial besar, ea1+ea2+...ea1+ea2+... e^{a_1} + e^{a_2} + ... dimana komponen dapat berkisar dari sangat kecil sampai yang sangat besar. Pendekatan saya adalah memperhitungkan faktor eksponensial terbesar sehingga:aaaK:=maxi(ai)K:=maxi(ai)K := \max_{i}(a_{i}) a′=K+log(ea1−K+ea2−K+...)a′=K+log(ea1−K+ea2−K+...)a' =K + log\left( …

24 monte-carlo  floating-point  statistics 

Ketika beralih dari bentuk kuat PDE ke bentuk FEM, tampaknya orang harus selalu melakukan ini dengan terlebih dahulu menyatakan bentuk variasional. Untuk melakukan ini, Anda mengalikan formulir kuat dengan elemen di beberapa (Sobolev) ruang dan mengintegrasikan di wilayah Anda. Ini bisa saya terima. Yang tidak saya mengerti adalah mengapa seseorang …

24 pde  finite-element  elliptic-pde