1 void washQtoM3(double m
[9], double q
[4]);
5 double clp
[2], xyz
[3], q
[4], len
;
8 for (xi
=0; xi
<samp
; xi
++)
11 q
[1] = ( ((double)(1)-(-1))*((double)((float)xi
)-(-0.5)) / ((double)(samp
-0.5)-(-0.5)) + (-1));
12 q
[2] = ( ((double)(1)-(-1))*((double)((float)0)-(-0.5)) / ((double)(samp
-0.5)-(-0.5)) + (-1));
13 q
[3] = ( ((double)(1)-(-1))*((double)((float)0)-(-0.5)) / ((double)(samp
-0.5)-(-0.5)) + (-1));
14 len
= (sqrt((((q
))[0]*((q
))[0] + ((q
))[1]*((q
))[1] + ((q
))[2]*((q
))[2] + ((q
))[3]*((q
))[3])));
15 ((q
)[0] = (q
)[0]*1.0/len
, (q
)[1] = (q
)[1]*1.0/len
, (q
)[2] = (q
)[2]*1.0/len
, (q
)[3] = (q
)[3]*1.0/len
);