USE HPF_LIBRARY
|
REAL x(nnode), flux(nedge)
|
INTEGER iedge(nedge,2)
|
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 )
|
x( permute_node(:) ) = x
|
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),:)
|
err = tol * 1e6
|
DO WHILE (err > tol)
-
flux=(x(iedge(1:nedge,1))-x(iedge(1:nedge,2)))/2
-
x=SUM_SCATTER(-flux(1:nedge),x,iedge(1:nedge,1))
-
x=SUM_SCATTER( flux(1:nedge),x,iedge(1:nedge,2))
-
err = SUM( flux*flux ) / nedge
|
END DO
|