HPF 1 Example: Irregular Mesh (2)
REAL x(nnode), flux(nedge)
INTEGER permute_node(nnode), permute_edge(nedge)
!HPF$ DISTRIBUTE x(BLOCK)
!HPF$ DISTRIBUTE flux(BLOCK)
!HPF$ ALIGN iedge(i,*) WITH flux(i)
!HPF$ ALIGN permute_edge(i) WITH flux(i)
!HPF$ ALIGN permute_node(i) WITH x(i)
CALL renumber_nodes( iedge, permute_node )
FORALL (i=1:nedge) iedge(i,:) = permute_node(iedge(i,:))
permute_edge = GRADE_UP( iedge(:,1) )
FORALL (i=1:nedge) iedge(i,:) = iedge(permute_edge(i),:)
DO WHILE (err > tol)
flux=(x(iedge(1:nedge,1))-x(iedge(1:nedge,2)))/2
err = SUM( flux*flux ) / nedge
x=SUM_SCATTER(-flux(1:nedge),x,iedge(1:nedge,1))
x=SUM_SCATTER( flux(1:nedge),x,iedge(1:nedge,2))