這是我的代碼 program load_flow
!!!!!!!region starts
!$acc region
p1=0;p=0
do i=2,n
do j=1,n
p1=(V(i)*V(j)*Yabs(i,j)*cos(angle(i,j)-d(i)+d(j)))
p=p1+p
end do
P0(i-1)=p
p
我試圖使用共享內存來緩存OpenACC的東西。 基本上就是我工作的是一個矩陣乘法，和我有什麼是這樣的： typedef float ff;
// Multiplies two square row-major matrices a and b, puts the result in c.
void mmul(const restrict ff* a,
const restric