Ditampilkan (Yousef Saad, metode berulang untuk sistem linier yang jarang , hal. 260) yangc o n d( A′A ) ≈ c o n d( A )2cond(A′A)≈cond(A)2cond(A'A) \approx cond(A)^2 Apakah ini juga berlaku untuk ?A A′AA′AA' Dalam kasus adalah dengan , saya amati bahwaN × M N ≪ M c o n d ( A ′ A ) ≫ c o n d ( A A...
Dalam proyek saya, saya harus menyelesaikan beberapa matriks tridiagonal pada setiap langkah waktu, jadi sangat penting untuk memiliki pemecah yang baik untuk mereka. Saya melakukan implementasi saya sendiri, hanya cara klasik untuk melakukannya dijelaskan di Wikipedia. Saya kemudian mencoba...
Diberikan matriks jarang generik dengan m << n (koreksi: ) elemen tidak nol (biasanya ). adalah generik dalam arti bahwa ia tidak memiliki sifat spesifik (misalnya kepastian positif), dan tidak ada struktur (misalnya pita) yang diasumsikan.A∈Rn×nA∈Rn×nA \in \mathbb{R}^{n\times n}m≪n2m≪n2m \ll...
Saya menggunakan Intel MKL untuk BLAS dan menggunakan Intel MKL Link Line Advisor untuk bantuan dengan opsi baris perintah. Penasihat menyediakan opsi untuk Static, Dynamic, dan Single Dynamic Library. Apa arti dari istilah-istilah ini? Ada banyak posting di SO mengenai hal ini tetapi semuanya...
Saya punya sistem persamaan diferensial biasa - 7 persamaan, dan ~ 30 parameter yang mengatur perilaku mereka sebagai bagian dari model matematika penularan penyakit. Aku suka mencari negara yang stabil bagi mereka persamaan Mengubah dx/dt = rest of the equationke 0 = equationuntuk masing-masing...
Pertanyaan ini lebih lanjut tentang cara mengatasi masalah secara numerik. Dalam sebuah proyek kecil saya ingin mensimulasikan gerakan coorbital Janus dan Epimetheus. Ini pada dasarnya adalah masalah tiga tubuh. Saya memilih Saturnus untuk diperbaiki pada titik asal, misalkan r1r1r_1 dan r2r2r_2...
Saya membuat pengamatan aneh sambil menghitung solusi untuk persamaan difusi reaksi 1D sederhana: ∂∂ta = ∂2∂x2a - a b∂∂ta=∂2∂x2a−ab\frac{\partial}{\partial t}a=\frac{\partial^2}{\partial x^2}a-ab ∂∂tb = - a b∂∂tb=−ab\frac{\partial}{\partial t}b=-ab ∂∂tc = a∂∂tc=a\frac{\partial}{\partial...
Saya tahu bahwa untuk memecahkan simetris eigen masalah , kita dapat menggunakan Sylvester Inersia Hukum, yaitu jumlah nilai eigen dari A kurang dari satu sama dengan jumlah entri negatif dari D mana matriks diagonal D berasal dari LDL faktorisasi A - sebuah I = L D L T . Kemudian, dengan metode...
Misalkan VVV adalah ruang vektor berdimensi-terbatas dengan norma ∥⋅∥‖⋅‖\|\cdot\| dan biarkan F:V→RF:V→RF : V \rightarrow \mathbb R menjadi fungsional linier terbatas. Itu hanya diberikan sebagai kotak hitam. Saya ingin memperkirakan norma FFF (dari atas dan bawah). Karena FFF adalah kotak hitam,...
Dalam makalah Hierarchical Conforming Finite Element Methods untuk Persamaan Biharmonic , P. Oswald mengklaim elemen tipe Clough-Tocher memiliki kontinuitas sambil menjadi polinomial kubik pada setiap segitiga. Dia tidak memberikan satu set fungsi dasar eksplisit hanya derajat standar kebebasan...
Saya sedang melakukan penelitian tentang struktur dalam komplemen Schur dan menemukan fenomena yang menarik: Misalkan A berasal dari laplacian 5-pt. Jika saya menggunakan metode diseksi bersarang dan metode multifrontal untuk menghitung faktorisasi LU dan kemudian memeriksa blok komplemen schur...
Mendaki bukit tampaknya menjadi alat yang sangat ampuh untuk optimasi. Namun, cara menghasilkan "tetangga" dari solusi selalu membingungkan saya. Misalnya, saya mengoptimalkan solusi . Di sini berada dalam kisaran (0, 0,1) , x_2 dalam kisaran (0, 100) , x_3 dalam kisaran (0, 1000000) . Apa cara...
Saya bertanya-tanya apa perbedaan dan hubungan antara "metode pencarian" dan "metode optimasi"? Apalagi saat menyelesaikan masalah optimasi? Saya menekankan konteks penyelesaian masalah pengoptimalan, karena saya kira metode pencarian tidak hanya untuk menyelesaikan masalah pengoptimalan, tetapi...
Misalkan saya memiliki masalah adveksi berkala 1D berikut: ∂u∂t+c∂u∂x=0∂u∂t+c∂u∂x=0\frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0dalamΩ=[0,1]Ω=[0,1]\Omega=[0,1] u(0,t)=u(1,t)u(0,t)=u(1,t)u(0,t)=u(1,t) u(x,0)=g(x)u(x,0)=g(x)u(x,0)=g(x) manag(x)g(x)g(x)memiliki diskontinuitas...
Saya terjun ke dunia analisis elemen hingga yang menakjubkan dan ingin menyelesaikan masalah termo-mekanis yang besar (hanya termal mekanis, tanpa umpan balik).→→\rightarrow Untuk masalah mekanis, saya sudah memahami jawaban Geoff , bahwa saya harus menggunakan solver iteratif karena ukuran mesh...
Saya mencoba untuk mengoptimalkan distributor aliran di tangki sehingga kecepatan dan distribusi suhu di setiap penampang relatif seragam. Ada banyak parameter yang dapat saya sesuaikan dengan keseragaman lintas bagian maksimum, seperti jumlah pipa saluran masuk, posisi, orientasi, dan arahnya....
Misalkan saya memiliki sistem linier besar, jarang: . Sekarang, saya tidak memiliki A - 1 karena A terlalu besar untuk faktor atau segala bentuk penguraian A , tetapi berasumsi bahwa saya memiliki solusi x 0 yang ditemukan dengan penyelesaian berulang.Ax0=b0Ax0=b0A\textbf{x}_0=\textbf{b}_0A-...
Saya ingin memprediksi runtimes untuk operasi aljabar linier padat pada arsitektur tertentu menggunakan perpustakaan tertentu. Saya ingin mempelajari model yang mendekati fungsinya Fo hlm: :Fop::F_{op} \;::\; ukuran input runtime→→ \rightarrow untuk operasi seperti matrik-gandakan, tambahkan...
Seperti yang dicontohkan oleh jawaban Jed Brown untuk Biaya pencarian versus perhitungan , menggunakan operasi floating point yang di-vectorized vs non-vectorized menghasilkan kode yang jauh lebih cepat. Banyak kompiler modern mengklaim bahwa mereka dapat melakukan vektorisasi otomatis. Bagaimana...
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...