2 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
6 * A PICTURE RECONSTRUCTION PROGRAM *
8 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
10 transm.cpp,v 1.2 2008/08/25 16:11:09 jklukowska Exp
12 SUBROUTINE CALLED BY RTFORT & NOT USEFUL INDEPENDENTLY.
13 APPLIES TRANSFORM SYMMETRY RELATIONS
14 J1: INDEX OF REAL PART OF ( I, K ) COEFFICIENT
15 J2: INDEX OF REAL PART OF ( M-I, N-K ) "
26 void transm(REAL* x, INTEGER j1, INTEGER j2, REAL sn, REAL cn)
37 ba = x[j1 + 1] + x[j2 + 1];
38 bb = x[j1 + 1] - x[j2 + 1];
39 real = cn * ba + sn * ab;
40 amag = sn * ba - cn * ab;
41 x[j1] = (aa + real) * (REAL) 0.5;
42 x[j2] = (aa - real) * (REAL) 0.5;
43 x[j1 + 1] = (amag + bb) * (REAL) 0.5;
44 x[j2 + 1] = (amag - bb) * (REAL) 0.5;