Sorry typos. That said, I have given up on this so maybe we just drop it.
inline proc vmRot(jk : range, r : int, param size : int, q : [] ?C, u : []
?R)
{
const (j, k) = (jk.first, jk.last);
param block = 0..#size;
var x = [param b in block] u[r + b, j];
for i in j+1 .. k do
{
const gi = g[i];
const c = gi.re;
const s = -gi.im;
for param b in block do
{
const yi = u[r + b, i];
const xi = x[b];
const cxi = c * xi;
const cyi = c * yi;
u[r + b, i - 1] = cxi - s * yi; // FIXED
x[b] = cyi + s * xi; // FIXED
}
}
[param b in block] u[r + b, k] = x[b]; // FIXED
}
Regards - Damian
Pacific Engineering Systems International, 277-279 Broadway, Glebe NSW 2037
Ph:+61-2-8571-0847 .. Fx:+61-2-9692-9623 | unsolicited email not wanted here
Views & opinions here are mine and not those of any past or present employer
_______________________________________________
Chapel-developers mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/chapel-developers
