Misalkan matriks berikut diberikan [ 0.500 - 0,333 - 0,167 - 0.500 ,667 - 0,167 - 0.500 - 0,333 0,833 ] dengan transposnya A T . Produk A T A = G menghasilkan [ 0,750 - 0,334 - 0,417 - 0,334 0,667 - 0,333 - 0,417 - 0,333 0,750 ]
(σ2(x)u′′(x))′′=f(x),0⩽x⩽1(σ2(x)u″(x))″=f(x),0⩽x⩽1 \left( \sigma^{2}(x) u ''(x) \right)'' = f(x), \;\;\; 0 \leqslant x \leqslant 1 u(0)=u(1)=0u(0)=u(1)=0u(0) = u(1) = 0u′′(0)=u′′(1)=0u″(0)=u″(1)=0u''(0) = u''(1) = 0σ(x)⩾σ0>0σ(x)⩾σ0>0\sigma(x) \geqslant \sigma_{0} > 0Au=fAu=fAu =...
Harap cantumkan paket Python (petsc4py, dll ...) dan pemecah langsung yang jarang didukungnya. Tolong, satu (komunitas-wiki) jawaban per
(Saya harap pertanyaan ini cocok dengan situs ini; jika tidak, terima permintaan maaf saya). Saya menjalankan simulasi tertentu, dan mendapat deret waktu y (t), t = 0, 1, ... 20. Setelah mencoba beberapa fungsi, saya menemukan bahwa: y(t) =~ 1 / (A t + B) Di mana A dan B adalah koefisien yang...
Misalkan A adalah matriks umum yang jarang, dan saya ingin menghitung nilai eigen. Saya tidak tahu cara mendeteksi multiplisitas untuk nilai eigen. Sejauh yang saya tahu, untuk kasus khusus, menemukan akar polinomial dengan metode matriks pengiring, kita dapat menerapkan RRQR untuk mendeteksi...
Dengan sistem mana A ∈ R n × n , saya membaca bahwa, dalam kasus iterasi Jacobi digunakan sebagai pemecah, metode ini tidak akan konvergen jika b memiliki komponen bukan nol dalam ruang nol A . Jadi, bagaimana seseorang dapat secara resmi menyatakan bahwa, asalkan b memiliki komponen non-nol yang...
Seberapa lambatkah petsc4py vs c / c ++ / fortran? Saya menyadari itu akan sangat tergantung pada kode yang dieksekusi, tetapi bagaimana dengan sesuatu yang sederhana seperti produk
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 adalah masalah enumerasi simetri sederhana. Saya memberikan latar belakang penuh di sini, tetapi tidak ada pengetahuan tentang kimia kuantum diperlukan. Kedua partikel terpisahkan adalah: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) ψ k ( x ) ψ l ( x ' )⟨ I j | k l ⟩⟨ij|kl⟩\langle ij|kl\rangle...
Saya memiliki set data bidang multipoint, setiap set data titik terkait dengan sel tunggal dari mesh yang tidak terstruktur. Tujuannya adalah untuk menginterpolasi data ke pusat sel, secara langsung atau tidak langsung, dengan cara yang paling akurat. Jika saya menggunakan interpolasi Tertimbang...
Saya mencoba memecahkan persamaan tipe: ( - ∂2∂x2- f( x ) ) ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Di mana f( x )f(x)f(x) memiliki kutub sederhana pada 000 , untuk nilai eigen NNN dan vektor eigen terkecil ....
Apakah ada algoritma untuk deteksi komunitas untuk grafik bipartit (jaringan 2-mode) yang diterapkan dalam igraph, networkX, R atau Python dll? Secara khusus, adakah implementasi seperti itu di mana seseorang akan dapat membatasi deteksi komunitas hanya pada salah satu dari dua
Persamaan matriks berikut dalam untuk matriks dan diberikan muncul dalam karya saya sebagai karakterisasi dari matriks kovarians. Saya telah belajar bahwa persamaan ini dikenal, khususnya dalam teori kontrol waktu kontinu, sebagai persamaan Lyapunov , dan bahwa ada berbagai algoritma terkenal...
Saya memiliki pengalaman dalam pengkodean OpenMP untuk mesin Memori Bersama (baik dalam C dan FORTRAN) untuk melakukan tugas-tugas sederhana seperti penambahan matriks, perkalian dll. (Hanya untuk melihat bagaimana itu bersaing dengan LAPACK). Saya tahu OpenMP cukup untuk melakukan tugas-tugas...
Jurnal Asosiasi untuk Komputasi Mesin Transaksi pada Perangkat Lunak Matematika (ACM TOMS) menerbitkan banyak artikel tentang algoritma numerik yang mencakup implementasi perangkat lunak. Menurut kebijakan editorial mereka , pengajuan makalah algoritma mencakup kode sumber untuk implementasi...
Adakah harapan dalam menyelesaikan sistem linear berikut secara efisien dengan metode berulang? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n >
Misalkan adalah matriks simetris nyata dan dekomposisi nilai eigennya V Λ V T diberikan. Sangat mudah untuk melihat apa yang terjadi dengan nilai eigen dari jumlah A + c I di mana c adalah konstanta skalar (lihat pertanyaan ini ). Bisakah kita menarik kesimpulan dalam kasus umum A + D di mana D...
Saya menghasilkan gambar 1d dan 2d menggunakan Matlab sekarang untuk perbandingan akurasi terhadap model yang diberikan. Saya perlu membandingkan metode saya dengan model Gaussian .wfn standar dan saya akan melakukannya dengan menganalisis kepadatan molekul serta laplacian. Saya tertarik untuk...
Sebuah pertanyaan yang cukup sederhana: untuk melakukan integral multidimensi, mengingat bahwa seseorang telah memutuskan bahwa semacam metode Monte Carlo sesuai, apakah ada keuntungan bahwa integrasi MC biasa menggunakan nomor pseudorandom memiliki lebih dari integrasi kuasi-Monte Carlo...
nVidia, misalnya, memiliki CUBLAS, yang menjanjikan peningkatan 7-14x. Naifnya, ini sama sekali tidak mendekati throughput teoritis kartu GPU nVidia. Apa saja tantangan dalam mempercepat aljabar linier pada GPU, dan apakah sudah ada rute aljabar linier yang lebih