ในการพยายามผสมความแม่นยำในโปรแกรมง่ายๆ โดยใช้ทั้งของจริงและสองเท่า และใช้รูทีน ddot จาก BLAS ฉันพบผลลัพธ์ที่ไม่ถูกต้องสำหรับชิ้นส่วนที่มีความแม่นยำสองเท่า นี่คือรหัส: program test !! adding this statement narrowed the issue down to ddot being...

ฉันต้องทำโปรแกรมต่อไปนี้ (จาก http://www.netlib.org/scalapack/examples/pblas.tgz ) ทำงานร่วมกับ SGEMM ฉันต้องเปลี่ยนแปลงอะไรบ้างเพื่อให้มันใช้งานได้? ความรู้ของฉันเกี่ยวกับ Fortran ค่อนข้างจำกัด...