วิธีการหนึ่งที่สามารถทำให้ขนานหลายวิธีสำหรับการแก้ระบบสมการเชิงเส้น?


11

ตามที่ฉันเข้าใจวิธี multigrid แก้ปัญหาระบบเชิงเส้นโดยการแก้ไขปัญหาแบบเดียวกันของ coarser (โดยการกำจัดข้อผิดพลาดความถี่ต่ำ) จากนั้นฉายกลับไปที่กริดที่ดีเพื่อลดข้อผิดพลาดความถี่สูง สำหรับระบบที่มีขนาดใหญ่ฉันสามารถดูว่าวิธีการวนซ้ำสามารถนำมาใช้ในการขนานในแต่ละระดับ วิธีนี้มีขนาดที่ดีในแบบคู่ขนานหรือไม่? มีแหล่งที่มาของการเกิดขึ้นพร้อมกันในอัลกอริทึมที่หนึ่งสามารถใช้ประโยชน์ในแบบคู่ขนาน?

คำตอบ:


14

Multigrid ทางเรขาคณิตแบบขนานนั้นง่ายต่อการติดตั้งบนกริดที่มีโครงสร้าง พีชคณิตแบบหลายจุดและไม่มีโครงสร้างนั้นเป็นเทคนิคมากกว่าดูคำตอบสำหรับลิงค์ไปยังการนำไปใช้

ในวิธีการคูณ (เช่น -004) สามารถคำนวณได้ครั้งละหนึ่งระดับเท่านั้น เนื่องจากจำนวนของระดับคือบันทึก Nที่Nคือจำนวนขององศาอิสระและเป็นปัจจัยอนุภาค (มักจะเกี่ยวกับ2 dหรือ3 dในdมิติ) ระยะลอการิทึมนี้ไม่ได้ถอดออกได้ วิธีการเติมจะเสียสละปัจจัยบางอย่าง แต่สามารถคำนวณได้ทุกระดับพร้อมกันซึ่งจะช่วยลดปัจจัยลอการิทึมเพื่อบันทึก2 log c NVเข้าสู่ระบบยังไม่มีข้อความยังไม่มีข้อความ2d3ddเข้าสู่ระบบ2เข้าสู่ระบบยังไม่มีข้อความ. ฉันยังไม่ได้เห็นการสาธิตเกี่ยวกับฮาร์ดแวร์จริงที่การทำงานพร้อมกันที่เพิ่มขึ้นช่วยให้ค่าคงที่ที่ด้อยลงและความทนทานของวิธีการเติมเพิ่มลดลง

O(ยังไม่มีข้อความ/P)

ในทางปฏิบัติกริดแบบหยาบจะมีค่าถึงขีดจำกัดความสามารถในการปรับขนาดได้อย่างรวดเร็ว (ซึ่งนอกเหนือจากการเพิ่มกระบวนการอื่น ๆ จะเพิ่มเวลาทำงาน) ดังนั้นจึงควรอยู่ในตัวสื่อสาร MPI ที่เล็กกว่าเดิม สิ่งนี้จะเพิ่มความซับซ้อนเล็กน้อยให้กับการนำไปใช้งาน สำหรับปัญหาที่ระดับหยาบมีโครงสร้างมากเกินไปที่จะดำเนินการหยาบต่อไปการแก้ปัญหาระดับหยาบอาจกลายเป็นปัญหาคอขวด

สำหรับการทดสอบวิธีการหลายมัลติพอร์ตแบบขนานผมแนะนำให้ใช้ไลบรารี่อย่างPETScซึ่งช่วยให้คุณสามารถรันอัลกอริทึมที่แตกต่างกันด้วยรหัสผู้ใช้น้อยมาก


ลิงก์ Adams (2001) ใช้งานไม่ได้อีกต่อไป ผมเชื่อว่าบทความที่คุณหมายถึงเป็นหนึ่งในนี้: ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=1592790&tag=1 "อัลกอริทึม Gauss-Seidel แบบกระจายไม่มีโครงสร้างสำหรับ Multigrid Smoothers" แจ้งให้เราทราบหากฉันผิด
nukeguy
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.