Ilmu Komputasi

Saya menggunakan paket PETSc non -linier SNES untuk menyelesaikan sistem persamaan nonlinier yang diperoleh dengan mendiskritisasi persamaan diferensial parsial. Bagaimana saya bisa menentukan mengapa pemecah tidak konvergen dan apa yang bisa saya lakukan untuk berhasil menyelesaikan persamaan saya?

22 petsc  pde  implicit-methods 

Seperti yang kita ketahui, metode transformasi ortogonal (rotasi Givens dan pantulan rumah tangga) untuk sistem persamaan linear lebih mahal daripada eliminasi Gaussian, tetapi secara teoritis memiliki sifat stabilitas yang lebih baik dalam arti bahwa mereka tidak mengubah nomor kondisi sistem. Meskipun saya tahu hanya satu contoh akademik dari sebuah matriks …

22 linear-algebra  reference-request 

Pergi untuk mengajar siswa tingkat sarjana kursus berjudul Pengantar Pemrograman Komputer. Saya sedikit bingung. Dalam Ilmu Fisika Komputasi, para ilmuwan menggunakan C / C ++ atau Python atau Fortran, CUDA dll ..... inilah saatnya untuk membangun basis mereka. Apa yang harus saya gunakan? Saya tahu Anda dapat mempelajari bahasa pemrograman …

22 python  c++  computational-physics  languages 

Seperti yang ada saat ini, pertanyaan ini tidak cocok untuk format tanya jawab kami. Kami berharap jawaban didukung oleh fakta, referensi, atau keahlian, tetapi pertanyaan ini kemungkinan akan mengumpulkan debat, argumen, polling, atau diskusi panjang. Jika Anda merasa bahwa pertanyaan ini dapat diperbaiki dan mungkin dibuka kembali, kunjungi pusat bantuan …

22 languages 

Saya memiliki fungsi dua dimensi yang nilainya ingin saya sampel. Fungsi ini sangat mahal untuk dikomputasi dan memiliki bentuk yang kompleks, jadi saya perlu mencari cara untuk mendapatkan informasi paling banyak tentang bentuknya menggunakan jumlah titik sampel yang paling sedikit.f( x , y)f(x,y)f(x,y) Metode apa yang bagus untuk melakukan ini? …

22 adaptive-mesh-refinement 

Paradigma pemrograman OpenCL berjanji untuk menjadi standar bebas royalti untuk komputasi heterogen. Haruskah kita menginvestasikan waktu kita dalam mengembangkan perangkat lunak berbasis OpenCL? Pro kontra?

22 hpc  programming-paradigms 

Apa pendekatan yang disukai dan efisien untuk menginterpolasi data multidimensi? Hal yang saya khawatirkan: kinerja dan memori untuk konstruksi, evaluasi tunggal / batch menangani dimensi dari 1 hingga 6 linier atau tingkat tinggi kemampuan untuk mendapatkan gradien (jika tidak linier) grid biasa vs tersebar menggunakan sebagai Fungsi Interpolasi, misalnya untuk …

22 python  performance  interpolation  delaunay-triangulation  extrapolation 

Saya mencari buku atau artikel, atau posting blog, atau materi apa pun yang diterbitkan secara umum, yang membahas secara khusus penggunaan fitur-fitur modern C ++ (memindahkan semantik, STL, iterator, evaluasi malas, dll.) Dalam komputasi ilmiah. Bisakah Anda menyarankan? Saya pikir fitur-fitur baru ini akan membuatnya lebih mudah untuk menulis kode …

22 c++  reference-request 

Saya akan mulai dengan pengalaman pribadi saya di lab kami. Kembali ke ifort 9 dan 10 hari, kami dulu cukup agresif dengan optimasi, mengkompilasi dengan -O3 dan flag spesifik prosesor (-xW -xSSE4.2 misalnya). Tetapi mulai dengan ifort 11, kami mulai memperhatikan: 1. beberapa inkonsistensi dalam hasil (karena semantik tidak dipertahankan) …

21 performance  hpc  compiling  fortran 

Bagaimana skala array Python / Numpy dengan dimensi array yang meningkat? Ini didasarkan pada beberapa perilaku yang saya perhatikan ketika membandingkan kode Python untuk pertanyaan ini: Cara mengekspresikan ekspresi rumit ini menggunakan irisan numpy Masalahnya sebagian besar melibatkan pengindeksan untuk mengisi array. Saya menemukan bahwa keuntungan menggunakan (tidak-sangat-baik) versi Cython …

21 python  performance  numpy 

Saya mencari algoritma untuk menggambar 2d-grafik standar untuk fungsi yang mungkin atau mungkin tidak memiliki singularitas. Tujuannya adalah untuk menulis "Mini-CAS", jadi saya tidak memiliki pengetahuan apriori tentang jenis fungsi, pengguna ingin membuat grafik. Masalah ini sudah sangat tua, jadi saya membayangkan harus ada beberapa algoritma standar dalam literatur. Untuk …

21 algorithms  visualization 

Ada sejumlah perpustakaan yang berbeda di luar sana yang memecahkan sistem persamaan linear yang jarang, namun saya merasa sulit untuk mencari tahu apa perbedaannya. Sejauh yang saya tahu ada tiga paket utama: Trilinos , PETSc , dan Intel MKL . Mereka semua dapat melakukan penyelesaian matriks yang jarang, mereka semua …

21 linear-algebra  linear-solver  libraries 

Rumus perbedaan hingga standar dapat digunakan untuk menghitung turunan secara numerik di bawah ekspektasi bahwa Anda memiliki nilai fungsi pada titik yang diberi spasi secara merata, sehingga adalah konstanta. Bagaimana jika saya memiliki titik spasi yang tidak rata, sehingga sekarang bervariasi dari satu pasang titik yang berdekatan ke yang berikutnya? …

21 finite-difference  discretization 

Saya mencari untuk memecahkan masalah optimasi terbatas di mana saya tahu batas-batas pada beberapa variabel (khususnya kendala kotak). argminuf(u,x)arg⁡minkamuf(kamu,x) \arg \min_u f(u,x) tunduk pada c(u,x)=0c(kamu,x)=0 c(u,x) = 0 a≤d(u,x)≤bSebuah≤d(kamu,x)≤b a \le d(u,x) \le b di mana adalah vektor variabel desain, adalah vektor variabel keadaan, dan adalah kendala kesetaraan (biasanya PDE). …

21 pde  optimization  nonlinear-programming 

Apakah ada metode untuk menyelesaikan sistem linear dari bentuk mana adalah matriks SPD tetap dan adalah matriks diagonal positif?O(n3+n2k)O(n3+n2k)O(n^3+n^2 k)kkk(Di+A)xi=bi(Di+A)xi=bi(D_i + A) x_i = b_iAAADiDiD_i Misalnya, jika setiap adalah skalar, itu sudah cukup untuk menghitung SVD dari . Namun, ini rusak untuk umum karena kurangnya komutatif.DiDiD_iAAADDD Pembaruan : Jawaban sejauh …

21 linear-algebra  algorithms  performance  complexity