Saya mengembangkan kode untuk mensimulasikan aliran fluida dengan zat biologis yang ada dalam aliran. Ini melibatkan persamaan standar Navier-Stokes yang digabungkan ke beberapa model biologis tambahan. Ada banyak parameter / konstanta. Saya memiliki fungsi tertulis untuk menangani perhitungan...
Bahasa C ++ menyediakan pemrograman generik dan metaprogramming melalui template. Teknik-teknik ini telah menemukan jalan mereka ke banyak paket komputasi ilmiah skala besar (misalnya, MPQC , LAMMPS , CGAL , Trilinos ). Tetapi apa yang sebenarnya mereka sumbangkan dalam komputasi ilmiah dalam nilai...
Saya tahu bahwa metode Nitsche adalah metode yang sangat menarik karena memungkinkan untuk memperhitungkan kondisi batas tipe Dirichlet atau kontak dengan kondisi batas gesekan dengan cara yang lemah tanpa menggunakan pengganda Lagrange. Dan keuntungannya, yaitu mengubah kondisi batas Dirichlet...
Saya perlu menulis beberapa data dari perhitungan, yang akan dibaca nanti oleh Paraview (file .vtu atau vtk). Ketika datang ke ukuran file, haruskah saya memilih format ASCII atau format
Misalkan saya memiliki beberapa fungsi fff dan saya ingin mencari sedemikian rupa sehingga . Saya mungkin menggunakan metode Newton-Raphson. Tapi ini mengharuskan saya tahu fungsi turunan . Ekspresi analitik untuk mungkin tidak tersedia. Sebagai contoh, dapat didefinisikan oleh sepotong kode...
Saya memiliki fungsi 2-D non-cembung terbatas yang ingin saya temukan minimum. Fungsinya cukup lancar. Mengevaluasi itu mahal. Kesalahan yang dapat diterima adalah sekitar 3% dari domain fungsi di setiap sumbu. Saya mencoba menjalankan implementasi algoritma DIRECT di perpustakaan NLOPT, tetapi...
Saya ingin tahu apakah ada cara cepat untuk menghitung jarak Euclidean dari dua vektor dalam Oktaf. Sepertinya tidak ada fungsi khusus untuk itu, jadi haruskah saya menggunakan rumus saja
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...
Banyak metode numerik untuk PDE hiperbolik didasarkan pada penggunaan pemecah Riemann. Pemecah seperti itu penting untuk menangkap gelombang kejut secara akurat. Ada berbagai jenis solver yang tersedia untuk sistem yang paling banyak dipelajari (misalnya, solver yang tepat, soloe Roe, solver HLL)....
Metode numerik untuk memecahkan PDE (atau ODE) terbagi dalam dua kategori besar: metode eksplisit dan implisit. Metode implisit memungkinkan timesteps stabil yang lebih besar tetapi membutuhkan lebih banyak pekerjaan per langkah. Untuk PDE hiperbolik, kebijaksanaan umum adalah bahwa metode implisit...
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 =
Diberikan PDE evolusi kamut= A u + B ukamut=SEBUAHkamu+Bkamuu_t = Au + Bu di mana adalah (mungkin nonlinear) operator diferensial yang tidak bolak-balik, pendekatan numerik yang umum adalah untuk bergantian antara penyelesaianA , BSEBUAH,BA,B kamut= A ukamut=SEBUAHkamuu_t = Au dan kamut= B u...
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 di bidang yang tidak selalu melakukan banyak pekerjaan HPC, dan ketika itu benar-benar terjadi, sering kali ini merupakan hasil dari peneliti dari bidang lain yang mengeksplorasi aplikasi baru dengan metode dan sejenisnya. Terutama apa artinya ini adalah bahwa hal itu tidak pernah benar-benar...
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
Pertimbangkan Ax=bAx=bAx=b dengan AAA hampir singular yang berarti ada nilai eigen λ0λ0\lambda_0 dari AAA yang sangat kecil. Kriteria berhenti yang biasa dari metode berulang didasarkan pada residual rn:=b−Axnrn:=b−Axnr_n:=b-Ax_n dan menganggap iterasi dapat berhenti ketika
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...
Ini mungkin pertanyaan tingkat siswa tetapi saya tidak bisa membuatnya jelas untuk diri saya sendiri. Mengapa lebih akurat menggunakan kisi-kisi tidak seragam dalam metode numerik? Saya berpikir dalam konteks beberapa metode beda hingga untuk PDE dari bentuk . Dan anggap saya tertarik pada solusi...
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....