การปรับระดับต่ำส่งผลต่อการรวมกันของวิธี Krylov อย่างไร


14

บอกว่าฉันมีระบบเชิงเส้นx = Bซึ่งลู่อย่างรวดเร็วโดยใช้วิธีการที่เหมาะสม Krylov (เช่น CG หรือ GMRES) สำหรับทุกข ถ้าBเป็นเมทริกซ์ที่มีระดับต่ำrวิธี Krylov เดียวกันบนระบบ( A + B ) x = b จะมาบรรจบกันอย่างรวดเร็ว (โดยมีการวนซ้ำจำนวนมากที่ขึ้นอยู่กับrเท่านั้น)Ax=bbBr(A+B)x=br

ตัวอย่างของระบบดังกล่าวคือความยืดหยุ่นของเมมเบรนที่ผ่านการปรับสภาพล่วงหน้าและการดัดรวมถึงเงื่อนไขความดันอากาศที่ไม่ได้รับการกำหนดเงื่อนไขพร้อมโครงสร้างผลิตภัณฑ์ด้านนอกที่หนาแน่น

โปรดทราบว่าคำถามจะเหมือนกันมีหรือไม่มี preconditioning ตั้งแต่เป็นอันดับRเปลี่ยนแปลงของP QP(A+B)Q=PAQ+PBQrPAQ

คำตอบ:


7

หากพื้นที่ย่อย Krylov ของคุณขึ้นอยู่กับพลังของบรรจบกันจะล่าช้าตามจำนวนการวนซ้ำที่มากที่สุดในระดับของการแก้ไข หากมันขึ้นอยู่กับพลังของA T Aดังนั้นจำนวนนี้จะมากเป็นสองเท่าAATA

ฉันไม่ได้เห็นคำกล่าวดังกล่าวในวรรณคดี แต่เพื่อดูความถูกต้องในกรณีแรกก็เพียงพอที่จะแสดงให้เห็นว่าพื้นที่ Krylov ของเมทริกซ์A + U S V Tซึ่งU , VkA+USVTU,Vมีคอลัมน์อยู่ในพื้นที่ที่สอดคล้องกันโดยไม่มีการแก้ไขระดับต่ำ แต่มี ดัชนีสูงขึ้นตามลําดับk + R ตรงไปตรงมาเพื่อตรวจสอบrk+r


คุณช่วยอธิบายสิ่งที่คุณหมายถึงโดย "อิงจากพลังของ " ได้หรือไม่? ตัวแก้ปัญหา Krylov ได้รับข้อมูลเกี่ยวกับA + Bเท่านั้นไม่ใช่Aโดยตรง AA+BA
Geoffrey Irving

ไม่เป็นไร: สมมุติว่าคุณหมายถึงพลังของเมทริกซ์ที่เป็นปัญหาดังนั้นในกรณีนี้ A+B
Geoffrey Irving

ใช่. วิธีการที่มีเมทริกซ์เป็นพารามิเตอร์และเมทริกซ์นี้มักจะเขียนแทนด้วย A
Arnold Neumaier

บางทีเพื่อความสนใจเพิ่มเติมคุณสามารถเขียนสมการของคุณ (หรือวิธีแก้ปัญหา) ด้วยความต้องการบางอย่างในถึงx = ( E + k = 1 ( A - 1 B ) k ) A - 1 bซึ่งอาจช่วยถ้าBเป็น nilpotent หรือA - 1 Bของบรรทัดฐานขนาดเล็ก หนึ่งยังตระหนักถึงการพึ่งพาการแก้ไขปัญหา Bx=(E+k=1(A1B)k)A1bBA1Bundisturbed
Bastian Ebeling
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.