PGF90 (Version 10.1) 01/29/2010 11:33:44 page 6
( 292) do loop=1,nn
( 293) gosa= 0.0
( 294) do k=2,kmax-1
( 295) do j=2,jmax-1
( 296) do i=2,imax-1
( 297) s0=a(I,J,K,1)*p(I+1,J,K) &
( 298) +a(I,J,K,2)*p(I,J+1,K) &
( 299) +a(I,J,K,3)*p(I,J,K+1) &
( 300) +b(I,J,K,1)*(p(I+1,J+1,K)-p(I+1,J-1,K) &
( 301) -p(I-1,J+1,K)+p(I-1,J-1,K)) &
( 302) +b(I,J,K,2)*(p(I,J+1,K+1)-p(I,J-1,K+1) &
( 303) -p(I,J+1,K-1)+p(I,J-1,K-1)) &
( 304) +b(I,J,K,3)*(p(I+1,J,K+1)-p(I-1,J,K+1) &
( 305) -p(I+1,J,K-1)+p(I-1,J,K-1)) &
( 306) +c(I,J,K,1)*p(I-1,J,K) &
( 307) +c(I,J,K,2)*p(I,J-1,K) &
( 308) +c(I,J,K,3)*p(I,J,K-1)+wrk1(I,J,K)
( 309) ss=(s0*a(I,J,K,4)-p(I,J,K))*bnd(I,J,K)
( 310) GOSA=GOSA+SS*SS
( 311) wrk2(I,J,K)=p(I,J,K)+OMEGA *SS
( 312) enddo
( 313) enddo
( 314) enddo
-----------------------
~/Himeno> pgf90 -fast -Minfo=intensity himenoBMTxp.f90
jacobi:
292, Intensity = [symbolic], and not printable, try the -Mpfi -Mpfo options
294, Intensity = [symbolic], and not printable, try the -Mpfi -Mpfo options
295, Intensity = [symbolic], and not printable, try the -Mpfi -Mpfo options
296, Intensity = 1.06
コンパイル時に情報が特定できないループ情報(外側ループの情報)は、1度実行する。
まず、-Mpfi オプションを付してコンパイル&リンクします。
~/Himeno> pgf90 -fast -Minfo=intensity himenoBMTxp.f90 -Mpfi
~/Himeno> ./a.out (実行)
実行後、pgfi.outと言う統計情報ファイルができます。これを元に再度、フィードバック
コンパイル(-Mpfo) を行うと、以下のように外側ループの Intensity が出力されます。
~/Himeno> pgf90 -fast -Minfo=intensity himenoBMTxp.f90 -Mpfo
jacobi:
292, Intensity = 12.37
294, Intensity = 1.06
295, Intensity = 1.06
296, Intensity = 1.06
|