Ilmu Komputasi

Saya menggunakan kotak alat pde Matlab untuk memecahkan persamaan elips tertentu dalam 2D. Solusi baik-baik saja, meskipun saya perlu memplotnya di sepanjang garis yang diberikan, yaitu untuk memotong sepotong planar dari jala 3D yang mewakili solusi. Saya tidak bisa menemukan cara yang secara cerdas melibatkan fungsi toolbox (yaitu tidak melibatkan …

9 pde  matlab 

Jadi saya sudah berdebat apakah saya harus repot-repot belajar Python. Dari berbicara dengan profesor saya, Matlab tampaknya menjadi bahasa umum yang digunakan dalam matematika terapan / ilmu komputasi sejauh menyangkut akademisi ; sementara di industri , profesor saya (terutama mereka yang telah bekerja di industri) mengatakan bahwa belajar c ++ …

9 matlab  python  c++ 

Untuk pekerjaan penelitian saya, saya telah mengembangkan perpustakaan C ++ kecil yang bertujuan memfasilitasi komunikasi antara kode komputasi C ++ dan Octave / Matlab (ketika yang terakhir digunakan untuk tujuan pasca-pemrosesan). Saya ingin merilis perpustakaan seperti itu di bawah GPL di beberapa situs hosting gratis seperti Sourceforge atau Github. Jujur …

9 matlab  libraries  c++ 

Dari definisi bilangan kondisi tampaknya diperlukan inversi matriks untuk menghitungnya, saya bertanya-tanya apakah untuk matriks kuadrat generik (atau lebih baik jika simetris positif pasti) dimungkinkan untuk mengeksploitasi beberapa dekomposisi matriks untuk menghitung nomor kondisi dalam suatu cara yang lebih cepat.

9 linear-algebra  matrix  condition-number 

Saya mencoba menerapkan metode Broyden-Fletcher-Goldfarb-Shanno untuk menemukan fungsi minimum. Saya perlu dua tebakan awal & x 0 dan perkiraan awal Matriks Goni B 0 . Satu-satunya persyaratan yang saya temukan untuk B 0 adalah bahwa jika Hessian pasti positif simetris, demikian juga B 0 . Melihat wikipedia, saya melihat bahwa …

9 optimization  algorithms 

Saya ingin memecahkan di mana K adalah matriks kekakuan saya. Namun beberapa kendala mungkin hilang dan karena itu beberapa gerakan tubuh yang kaku mungkin masih ada dalam sistem (karena nilai eigen nol). Karena saya menggunakan CG untuk menyelesaikan sistem linear, ini tidak dapat diterima karena kadang-kadang CG tidak konvergen pada …

9 finite-element  iterative-method  conjugate-gradient 

Apakah ada cara yang disukai bagaimana menerapkan evaluasi cepat (perkiraan) dari polinomial interpolasi Chebyshev pada grid yang seragam (diberi nilai fungsi pada node Chebyshev)? Masalah saya adalah interpolasi menjadi lambat ketika tingkat polinomial interpolasi meningkat. Ide-ide berikut muncul di benak saya: Cobalah untuk menyesuaikan teknik FFT (NFFT) yang tidak seragam …

9 interpolation  fourier-analysis  fftw 

Melakukan simulasi numerik langsung turbulensi isotropik dengan metode pseudo-spektral Fourier (Orzag & Patterson, PRL, 1972) menggunakan FFT. Untuk latar belakang metode ini, yang banyak digunakan dalam komunitas turbulensi, Anda dapat melihat kursus ini: http://www.math.ualberta.ca/~bowman/m655/lab3d.pdf Menggunakan disebut aturan untuk dealiasing yang terdiri dari melakukan u ( k , t ) = …

9 fluid-dynamics  simulation 

Saya memiliki masalah di mana saya harus menemukan semua positif (seperti dalam nilai eigen positif) pasang eigen dari matriks nonsimetrik kecil (biasanya lebih kecil dari 60x60). Saya bisa berhenti menghitung ketika nilai eigen lebih kecil dari ambang tertentu. Saya tahu bahwa nilai eigennya nyata. Adakah saran tentang algoritma yang dapat …

9 performance  eigensystem  gpu 

Saya mendapat kesan, dari sumber daya yang sangat berbeda dan pembicaraan dengan para peneliti, bahwa ada peningkatan permintaan untuk perhitungan presisi tinggi dalam persamaan diferensial parsial numerik. Di sini, presisi tinggi berarti lebih presisi daripada hanya presisi ganda 64bit standar. Saya bertanya-tanya tentang seni topik ini. Sebagai perbandingan, ada komunitas …

9 reference-request  precision 

Sebagian besar pemrograman saya adalah kode penelitian satu kali dalam C untuk saya gunakan sendiri. Saya belum pernah mendistribusikan kode apa pun selain dari kolaborator dekat. Saya telah mengembangkan suatu algoritma yang saya terbitkan dalam jurnal ilmiah. Saya ingin memberikan kode sumber dan mungkin kode yang dapat dieksekusi dalam suplemen …

9 linear-solver  libraries  c++  c 

Saya ingin meminimalkan fungsi tujuan yang rumit, dan saya tidak yakin apakah itu cembung. Apakah ada algoritma yang bagus yang mencoba membuktikan bahwa itu bukan cembung? Tentu saja algoritme dapat gagal untuk membuktikan ini, dalam hal ini saya tidak akan tahu apakah itu cembung atau tidak, dan ini OK; Saya …

9 convex-optimization 

Membolak-balik beberapa buku pelajaran, saya perhatikan bahwa masalah awalnya mengurung minimum selama pencarian baris cenderung menjadi renungan (setidaknya dalam teks sarjana saya). Apakah ada teknik yang sudah mapan atau praktik terbaik untuk jenis masalah ini, atau apakah solusi biasanya tergantung pada aplikasi? Adakah yang bisa merekomendasikan beberapa referensi tentang topik …

9 optimization  reference-request 

Answer memberikan perangkat lunak untuk menghitung bersyarat Lyapunov eksponen (CLE) untuk osilator ditambah dalam sinkronisasi kekacauan. Namun, sulit untuk diikuti dan tidak ada output grafis dari plot (dan itu dalam C, lebih kompleks). Adakah yang tahu cara memodifikasi LET Toolbox yang sangat bagus untuk sistem yang tidak terhubung tetapi saya …

9 matlab  ode  simulation  quantum-mechanics 

Saya ingin tahu apakah ada algoritma yang memberikan satu set titik o dan sudut menghitung cembung-lambung jika sudutnya dan diberi α > 0 menghitung amplop yang mengikuti lebih dekat "perimeter".α=0α=0\alpha = 0α>0α>0\alpha > 0 Dan jika ada definisi perimeter non berpotongan satu set poin, dalam hal ini poligon yang dihasilkan …

9 algorithms  computational-geometry