Ilmu Komputasi

Untuk beberapa domain cembung sederhana dalam 2D, kami memiliki beberapa u ( x ) yang memenuhi persamaan berikut: - d i v ( A ∇ u ) + c u n = f dengan Dirichlet tertentu dan / atau kondisi batas Neumann. Setahu saya, menerapkan metode Newton dalam ruang elemen …

11 pde  finite-element 

Saya ingin memvisualisasikan getaran molekuler yang bukan mode normal. Saya ingin menyajikan representasi statis gerakan vektor, dan saya ingin beberapa fleksibilitas dalam gaya vektor (ukuran, warna dll). Saya juga tertarik untuk menghasilkan video getaran. Apa sajakah sumber daya yang baik untuk menampilkan getaran molekuler? Preferensi saya adalah untuk alat open …

11 computational-chemistry 

Berikut ini adalah penggunaan khas pointer fungsi di C. Saya ingin melakukan sesuatu yang serupa di Fortran. Saya punya beberapa ide, tetapi saya ingin tahu apakah ada cara kanonik untuk melakukannya. Pointer fungsi dan konteks yang dilewatkan oleh pengguna disimpan, lalu dipanggil nanti. typedef PetscErrorCode (*TSIFunction)(TS,PetscReal,Vec,Vec,Vec,void*); PetscErrorCode TSSetIFunction(TS ts,Vec res,TSIFunction …

11 petsc  fortran 

Adakah yang bisa membantu saya menemukan buku tentang solusi numerik (beda hingga dan metode Crank-Nicolson) dari persamaan Poisson dan difusi termasuk contoh-contoh pada geometri tidak beraturan, seperti domain yang terdiri dari area antara persegi panjang dan lingkaran (terutama buku atau tautan) pada contoh kode MATLAB dalam kasus ini)?

11 pde  finite-difference 

Bersamaan dengan penelitian saya pada Metode Ruang Krylov, saya memiliki pilihan untuk mengeksplorasi matematika di belakang HPC selangkah lebih maju atau teori komputasi (perangkat keras, OS, kompiler dll.). Saat ini, saya tahu baik cukup untuk hanya mendapatkan oleh. Sebagai contoh, saya tahu bagaimana menurunkan persamaan untuk CG dan dasar-dasar metode …

11 parallel-computing  education 

Ketika memecahkan PDE tergantung waktu menggunakan metode elemen hingga, misalnya mengatakan persamaan panas, jika kita menggunakan stepping waktu eksplisit maka kita harus menyelesaikan sistem linear karena matriks massa. Misalnya jika kita tetap dengan contoh persamaan panas, ∂kamu∂t= c ∇2kamu∂kamu∂t=c∇2kamu\frac{\partial{u}}{\partial{t}} = c\nabla{}^{2}u kemudian gunakan forward Euler yang kita dapatkan M.( kamun …

11 finite-element  time-integration  navier-stokes 

Saya membaca di suatu tempat bahwa Pandas pertama kali dikembangkan untuk dunia keuangan, setidaknya tidak terutama untuk ilmu alam (fisika, biologi, dan sebagainya), jadi apakah ada analisis data serupa paket Python yang lebih “berorientasi ilmu alam”? Saya baru saja mulai menggunakan Pandas dan saya telah menemukan dua masalah yang tidak …

11 python  data-analysis  data-management 

Apakah ada algoritma SVD terpotong yang menghitung nilai singular satu per satu? Masalah saya: Saya ingin menghitung nilai singular pertama (dan vektor singular) dari matriks padat , tapi saya tidak tahu apa nilai tepat. itu besar, jadi untuk alasan efisiensi, saya lebih suka tidak mengevaluasi SVD lengkap hanya untuk memotong …

11 linear-algebra  algorithms  svd 

Saya punya pertanyaan yang mirip dengan yang diajukan sebelumnya kecuali dalam 3D, dan saya hanya perlu volume, bukan bentuk lambung yang sebenarnya. Lebih tepatnya, saya diberi satu set kecil poin (katakanlah, 10-15) dalam 3D, yang semuanya diketahui terletak pada cembung lambung set point (sehingga mereka semua "penting" dan mendefinisikan lambung). …

11 algorithms  reference-request  computational-geometry 

Untuk sebuah proyek, saya harus mengimplementasikan kedua metode ini dan membandingkan kinerjanya pada fungsi yang berbeda. Sepertinya metode gradien konjugasi dimaksudkan untuk menyelesaikan sistem persamaan linear untuk Ax=bAx=b A\mathbf{x} = \mathbf{b} Di mana adalah matriks n-by-n yang simetris, pasti-positif dan nyata.AAA Di sisi lain, ketika saya membaca tentang gradient descent, …

11 optimization  conjugate-gradient 

Saya mengimplementasikan makalah " Transportasi Massal Optimal untuk Registrasi dan Warping ", tujuan saya adalah untuk membuatnya online karena saya tidak dapat menemukan kode transportasi massal euler online dan ini akan menarik setidaknya untuk komunitas riset dalam pemrosesan gambar. Makalah ini dapat diringkas sebagai berikut: - temukan peta awal uuu …

11 pde  finite-difference  matlab  fluid-dynamics  image-processing 

Ada algoritma untuk memecahkan fungsi Odes pegangan , di manay∈Rn. Tetapi dalam banyak sistem fisik, persamaan diferensial adalah otonom, jadidydydt= f( y, t )dydt=f(y,t)\frac{dy}{dt} = f(y, t)y∈Rny∈Rny \in \mathbb R^n,y∈Rn, dengant yangtertinggal. Dengan asumsi yang disederhanakan ini, perbaikan apa yang dapat dilihat dalam metode numerik yang ada? Misalnya, jikan=1, masalahnya …

10 ode  numerics 

Saya menemukan bahwa banyak pemrograman sains komputasi saya memiliki persyaratan pengujian yang tidak dicakup oleh kerangka uji standar: Pengujian waktu komputasi Untuk memastikan bahwa algoritma tidak menjadi lebih lambat. Saya bisa melakukan sesuatu seperti assureSmallerEqual(RuntimeWrapper(algorithm),53)tetapi saya ingin ambang 53 detik dikurangi terus menerus karena saya bekerja pada algoritma, yaitu sesuatu …

10 testing 

Setiap matriks nyata dapat direduksi menjadi bentuk Schur nyata T = U T A U menggunakan transformasi U yang mirip ortogonal . Di sini matriks T adalah bentuk kuasi-segitiga dengan 1 oleh 1 atau 2 oleh 2 blok pada diagonal utama. Setiap 1 oleh 1 blok berkorespondensi dengan nilai eigen …

10 eigenvalues  eigensystem  lapack  numerics  scalapack 

Dalam penelitian tentang masalah invers, adalah umum untuk membangun set data sintetik dari set parameter yang diketahui dan kemudian menguji apakah teknik inversi dapat merekonstruksi parameter tersebut. Dengan demikian, penting untuk menambahkan tingkat kebisingan acak yang sesuai ke data sintetis. Lebih jauh, jika metode yang digunakan untuk menghitung data sintetis …

10 reference-request  inverse-problem