1 | From http://www.cray.com/products/systems/crayt3e/1200/performance.html |
2 | REAL*8 AA(513,513), DD(513,513) |
3 | REAL*8 X (513,513), Y (513,513) |
4 | REAL*8 RX(513,513), RY(513,513) |
5 | DO J = 2,N-1 |
6 | DO I = 2,N-1 |
7 | XX = X(I+1,J)-X(I-1,J) |
8 | YX = Y(I+1,J)-Y(I-1,J) |
9 | XY = X(I,J+1)-X(I,J-1) |
10 | YY = Y(I,J+1)-Y(I,J-1) |
11 | A = 0.25 * (XY*XY+YY*YY) |
12 | B = 0.25* (XX*XX+YX*YX) |
13 | C = 0.125 * (XX*XY+YX*YY) |
14 | Continued on Next Page |