
376 Designing Scientific Applications on GPUs
to the diagonal of the local R-matrices. The local R-matrices, together with
the input R-matrix, <
I
, computed on the previous sector, are used to com-
pute the current output R-matrix, <
O
. At the end of a sector evaluation, the
output R-matrix becomes the input R-matrix for the next evaluation.
Algorithm 19: PROP algorithm
1 for all scattering energies do
2 for all sectors do
3 Read amplitude arrays;
4 Read correction data;
5 Construct local R-matrices;
6 From <
I
and local R-matrices, compute <
O
;
7 <
O
becomes <
I
for the next sector;
8 end
9 Compute physical R-Matrix;
10 end
On the first sector, there is no input R-matrix ...