Solver linear iteratif mana yang konvergen untuk matriks semidefinit positif?

10

Saya ingin tahu yang mana dari pemecah linear klasik (mis. Gauss-Seidel, Jacobi, SOR) dijamin akan menyatu untuk masalah mana A adalah positif setengah pasti dan tentu saja b i m ( A )SEBUAHx=bSEBUAHbsayam(SEBUAH)

(Pemberitahuan adalah semi pasti dan tidak pasti)SEBUAH

olamundo
sumber
1
Apakah maksud Anda matriks semi-pasti positif?
meawoppl
1
Apa gunanya menyelesaikan sistem linear dengan matriks seperti itu? Jika saya tidak salah, jika matriks semidefinite positif Anda non-singular maka itu pasti positif.
faleichik
1
Ya saya yakin. Saya harus menyegarkan ingatan saya sebagai bukti aktual, tetapi sesuai dengan apa yang Anda katakan - jika penyebut dalam perhitungan adalah nol, itu berarti A P k adalah nol, yang berarti bahwa semua "arah pencarian" di mana A bukan singular telah habis, dan residu yang tersisa tidak dalam rentang A (dan dengan demikian ini adalah solusi "optimal"). Dalam kasus yang sebenarnya b s p a n ( A ) , ini tidak akan terjadi sebagai sisa akan mencapai nol sebelum pertama kalinya A P k = 0αSEBUAHPkbshalSebuahn(SEBUAH)SEBUAHPk=0
olamundo
1
Set . Kemudian A n b I m ( A ) . CG akan konvergen karena x n A x n > 0x0=bSEBUAHnbsayam(SEBUAH)xnSEBUAHxn>0 untuk semua . Dengan kata lain, Anda tidak pernah meninggalkan I m ( A ) yang A pasti-positif. 0xnsayam(SEBUAH)sayam(SEBUAH)SEBUAH
Deathbreath
2
@faleichik: matriks kerapatan tereduksi dalam mekanika kuantum adalah semi-pasti positif dalam banyak kasus.
Deathbreath

Jawaban:

8

Algoritma gradien konjugasi bekerja untuk masalah semidefinite dan menghasilkan solusi norma minimal.

Arnold Neumaier
sumber
Terima kasih. Setiap ide tentang pemecah "kuno", misalnya SOR Gauss-Seidel dll.
olamundo
Mereka hampir tidak pernah digunakan lagi, dan saya tidak tahu bagaimana ini berperilaku.
Arnold Neumaier
Untuk mengklarifikasi: CG pasti tidak bekerja dalam bentuk vanilla untuk matriks semi-pasti; ini dapat bekerja secara teori jika B berada dalam citra A; tetapi ini tidak mungkin berakhir dengan baik dalam praktik numerik. MINRES berbasis krylov yang sangat mirip adalah pilihan yang jauh lebih baik di sini. Juga, ini "kuno" pemecah banyak digunakan dalam pemecah tipe multigrid, untuk menyebutkan satu contoh.
Eelco Hoogendoorn
1

Berikut ini adalah bukti bahwa Gauss-Seidel cocok dengan kebutuhan Anda, mengingat bahwa ada dalam citra AbSEBUAH .

Hal yang sama sangat tidak benar pada Jacobi; yang memalukan karena siapa yang mau repot dengan Gauss-Seidel pada perangkat keras komputer modern? Jika masalah Anda dapat dipecah menjadi blok-blok yang dominan diagonal, Anda beruntung; Anda dapat menerapkan pembaruan Jacobi ke blok-blok tersebut dengan cara Gauss-Seidel tambahan, dan dapatkan yang terbaik dari keduanya untuk jenis masalah semi-pasti ini.

Eelco Hoogendoorn
sumber