62£®INDEPENDENTµÄµäÐÍʹÓà |
SUBROUTINE update(x, nbr)
¡¡¡¡¡¡REAL x(:), tmp
¡¡¡¡¡¡INTEGER nbr(2,:), i
¡¡¡¡¡¡¡¡¡¡DO i = 1, UBOUND(nbr,2)
¡¡¡¡¡¡¡¡¡¡¡¡¡¡tmp = flux( x(nbr(1,i)), x(nbr(2,i)))
¡¡¡¡¡¡¡¡¡¡¡¡¡¡x(nbr(1,i)) = x(nbr(1,i)) + tmp
¡¡¡¡¡¡¡¡¡¡¡¡¡¡x(nbr(2,i)) = x(nbr(2,i)) - tmp
¡¡¡¡¡¡¡¡¡¡END DO
¡¡¡¡¡¡END DO
Copyright: NPACT |