62£®INDEPENDENTµÄµäÐÍʹÓÃ

BACKWARDFORWARD



 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 BACKWARDFORWARD