อัลกอริทึม 1/t ที่ใช้ความหนาแน่นสองสถานะสำหรับประมาณปริพันธ์หลายมิติ

ปัญหาที่งานวิจัยแก้ไข

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

ผู้ได้รับประโยชน์:

• นักฟิสิกส์เชิงคำนวณ
• นักวิจัยด้านการอินทิเกรตเชิงตัวเลข
• นักวิทยาศาสตร์ที่วิเคราะห์ฟังก์ชันหลายมิติซึ่งมียอดแหลม
• ผู้พัฒนาการจำลองที่ต้องการลดต้นทุนมอนติคาร์โล

วิธีการ

หลังสุ่มตัวอย่างเบื้องต้นจำนวนเล็กน้อย วิธีนี้เลือกเกณฑ์ y_s ที่เข้าถึงได้ในบริเวณที่มีความหนาแน่นสูง แล้วแบ่งปริภูมิ y เป็นเพียงสองสถานะ คือ y < y_s และ y >= y_s จากนั้นใช้การเดินสุ่ม 1/t ตามปกติระหว่างสองสถานะ พร้อมปรับค่าเฉลี่ยภายในสถานะ จึงไม่ต้องแบ่งช่วงคงที่ละเอียดหรือกำหนด y_min และ y_max งานวิจัยทดสอบปริพันธ์เกาส์เซียนหลายตัวแปรตั้งแต่หนึ่งถึงหกมิติ ปริพันธ์วงแหวนเกาส์เซียนตั้งแต่สองถึงหกมิติ และปริพันธ์แผนภาพไฟน์แมนแบบ setting-sun ด้วยการจำลองอิสระ 50 ครั้งและ N_final = 10^8

ผลการศึกษา

วิธีสองสถานะให้ค่าที่แม่นยำในการทดสอบเกาส์เซียนทุกกรณีถึงหกมิติ โดยได้ 0.9995586 สำหรับเกาส์เซียน 6D และ 0.9992425 สำหรับวงแหวนเกาส์เซียน 6D ขณะที่วิธี 1/t แบบหลายช่วงเยี่ยมไม่ครบทุกสถานะในการทดสอบเกาส์เซียน 5D และ 6D ภายใน 10^8 การทดลอง สำหรับตัวอย่าง setting-sun ค่าประมาณ 1.594486 x 10^9 ใกล้เคียงกับวิธีควอดราเจอร์แบบปรับตัวที่ 1.596160 x 10^9 แต่มอนติคาร์โลพื้นฐานให้เพียง 4.041376 x 10^7 ความคลาดเคลื่อนลดลงตามอัตรา N^(-1/2) โดยไม่เกิดการอิ่มตัว