44 return (k + j*nn[2] + i*nn[1]*nn[2]);
66 vector cellL = 2 * (box.
max() - cornerCellCentre);
74 NN[i] = label(L[i]/cellL[i] + 0.5);
82 rdeltaByL[i] = NN[i]/(L[i]*L[i]);
88 <<
"calculated number of cells is incorrect"
92 for (i=0; i<NN[0]; i++)
94 scalar k1 = (i - NN[0]/2)*twoPi/L[0];
96 for (label j=0; j<NN[1]; j++)
98 scalar k2 = (j - NN[1]/2)*twoPi/L[1];
100 for (label
k=0;
k<NN[2];
k++)
102 scalar k3 = (
k - NN[2]/2)*twoPi/L[2];
104 (*this)[rep(i, j,
k, NN)] =
vector(k1, k2, k3);
109 Kmax =
mag((*
this)[rep(NN[0]-1, NN[1]-1, NN[2]-1, NN)]);