Irregular Mesh: Sequential Program

REAL x(nnode), y(nnode), flux

INTEGER iedge(nedge,2)

err = tol * 1e6

DO WHILE (err > tol)
DO i = 1, nedge

flux = (y(iedge(i,1))-y(iedge(i,2))) / 2

x(iedge(i,1)) = x(iedge(i,1)) - flux
x(iedge(i,2)) = x(iedge(i,2)) + flux
err = err + flux(i)*flux(i)
END DO
err = err / nedge
DO j = 1, nnode
y(i) = x(i)

END DO

END DO

Previous slide

Next slide

Back to first slide

View graphic version