Mengikuti dari pertanyaan saya sebelumnya, saya mencoba menerapkan kondisi batas pada jaring volume terbatas yang tidak seragam ini, Saya ingin menerapkan kondisi batas tipe Robin ke lhs domain ( x = xL.)x=xL.)x=x_L) , sehingga, σL.= ( dkamux+ a u ) ∣∣∣x=
Saya seorang penganjur besar pengembangan yang digerakkan oleh tes dalam komputasi ilmiah. Utilitas dalam praktiknya hanya mengejutkan, dan benar-benar mengurangi masalah klasik yang diketahui pengembang kode. Namun, ada kesulitan yang melekat dalam pengujian kode ilmiah yang tidak ditemukan dalam...
Apakah ada alasan mengapa seseorang harus memilih Runge Kutta (IMRK) implisit tingkat tinggi dari waktu ke waktu melangkah? BDF tampaknya lebih mudah bagi saya karena stage IMRK membutuhkan linear solves per step waktu. Stabilitas untuk BDF dan IMRK tampaknya menjadi titik perdebatan. Saya tidak...
Untuk tujuan pengajaran, saya membutuhkan fungsi kontinu dari satu variabel yang "sulit" untuk diperkirakan dengan polinomial, yaitu seseorang akan membutuhkan kekuatan yang sangat tinggi dalam rangkaian daya untuk "menyesuaikan" fungsi ini dengan baik. Saya bermaksud menunjukkan kepada siswa saya...
Hal yang saya paling tidak suka tentang MPI adalah berurusan dengan tipe data (yaitu peta data / topeng) karena mereka tidak cocok dengan C + berorientasi objek. boost::mpiNamun, hanya mendukung MPI 1.1 dari situs web mereka: boost :: mpi adalah antarmuka ramah C ++ ke Antarmuka Lulus Pesan...
Saya melakukan tugas klasifikasi teks dengan R, dan saya mendapatkan matriks istilah dokumen dengan ukuran 22490 sebanyak 120.000 (hanya 4 juta entri bukan nol, entri kurang dari 1%). Sekarang saya ingin mengurangi dimensionalitas dengan memanfaatkan PCA (Principal Component Analysis). Sayangnya, R...
Saya sudah mulai menggunakan python sebagai bahasa pemrograman untuk melakukan semua tugas saya di CFD. Saya memiliki sedikit pengalaman dalam pemrograman. Saya dari latar belakang teknik mesin dan sedang mengejar pendidikan tinggi di bidang teknik Aerospace. terkadang aspek komputasi CFD menjadi...
Baru-baru ini saya telah membandingkan berbagai pemecah non-linear dari scipy dan sangat terkesan dengan contoh Newton-Krylov dalam Scipy Cookbook di mana mereka memecahkan persamaan persamaan diferensial orde dua dengan istilah reaksi non-linear dalam sekitar 20 baris kode. Saya memodifikasi kode...
Saya ingin tahu apa algoritma numerik yang baik yang ada untuk evaluasi fungsi hypergeometric umum (atau seri), didefinisikan sebagai halFq( a1, ... , ahal; b1, ... , bq; z) = ∑k = 0∞( a1)k⋯ ( ahal)k( b1)k⋯ ( bq)kzkk
Saat menghitung faktorisasi QR dalam praktiknya, seseorang menggunakan pantulan Householder untuk menghilangkan bagian bawah dari sebuah matriks. Saya tahu bahwa untuk menghitung nilai eigen dari matriks simetris, yang terbaik yang dapat Anda lakukan dengan pantulan Householder adalah membuatnya...
Diberi matriks padat apa cara terbaik untuk menemukan basis ruang-nol dalam beberapa toleransi ?A∈Rm×n,m>>n;max(m)≈100000A∈Rm×n,m>>n;max(m)≈100000A \in R^{m \times n}, m >> n; max(m) \approx 100000 ϵϵ\epsilon Berdasarkan pada dasar itu dapatkah saya mengatakan bahwa col tertentu...
Seharusnya minAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i , j , k =1 , … , nminAvec(U)subject to Ui,j≤maks{Usaya,k,Uk,j},saya,j,k=1,...,n\begin{align*} \min A &\mathrm{vec}(U) \\ &\text{subject to } U_{i,j} \leq \max\{U_{i,k}, U_{k,j}\}, \quad i,j,k = 1, \ldots, n \end{align*} di mana UUU adalah...
Saya menggunakan Valgrind + Callgrind ke profil pemecah yang saya tulis. Seperti yang dinyatakan oleh manual pengguna Valgrind, saya telah mengkompilasi kode saya dengan opsi debugging untuk kompiler: "Tanpa info debug, alat Valgrind terbaik yang dapat Anda lakukan adalah menebak fungsi kode...
Saya ingin mandiri dari perangkat lunak komersial untuk karya ilmiah saya. Saya menemukan ketergantungan suatu paket komersial seperti Matlab dan kotak alatnya tidak memuaskan, karena saya tidak tahu apakah saya akan memiliki akses ke Matlab di masa depan, dan karena saya tidak suka bahasanya....
Saya telah belajar tentang Metode Elemen Hingga (juga sedikit tentang metode numerik lainnya) tetapi saya tidak tahu apa sebenarnya definisi dari dua kesalahan dan perbedaan di antara
Saya membaca tentang ilmu komputasi di Wikipedia, tetapi pemahaman saya tidak begitu jelas. Apakah ilmu komputasi melibatkan pemrograman? Seberapa berbeda ilmu komputasi dengan komputasi _ ____ , di mana yang kosong bisa berupa disiplin apa saja (ilmu material, teknik, kimia, biologi, dan...
Misalkan saya memiliki 2D jala terdiri dari non-overlapping segitiga , dan satu set poin { p i } M i = 1 ⊂ ∪ N k = 1 T K . Apa cara terbaik untuk menentukan di mana segitiga masing-masing titik terletak?{ Tk}Nk =1{Tk}k=1N\{T_k\}_{k=1}^N{ psaya}M.i = 1⊂ ∪Nk =
Saya sadar bahwa membalikkan matriks untuk menyelesaikan sistem linear bukanlah ide yang baik, karena itu tidak seakurat dan seefisien langsung memecahkan sistem atau menggunakan LU, Cholesky atau penguraian QR. Namun, saya belum dapat memeriksa ini dengan contoh praktis. Saya telah mencoba kode...
Pendekatan mana yang digunakan dalam praktik untuk memperkirakan jumlah kondisi dari matriks jarang
Ada sistem kendala linier Ax≤bAx≤b{\bf Ax} \leq {\bf b} . Saya ingin menemukan vektor yang benar-benar positif x>0x>0{\bf x} > 0 yang memenuhi batasan-batasan ini. Itu berarti, xi>0xi>0x_i > 0 diperlukan untuk setiap komponen xixix_i dari xx{\bf x} . Bagaimana saya bisa menggunakan LP...