52 scalar l0(K.
sizeOfBox()[0]*(scalar(K.
nn()[0])/(scalar(K.
nn()[0])-1.0)));
60 y /= scalar(K.
size());
75 const label tnp = Ek.
size();
76 const label NoSubintervals = label
78 pow(scalar(tnp), 1.0/vector::dim)*
pow(1.0/vector::dim, 0.5) - 0.5
85 scalar kmax = K.
max()*
pow(1.0/vector::dim,0.5);
86 scalar delta_k = kmax/(NoSubintervals);
90 k1D[a] = (a + 1)*delta_k;
97 scalar kmag =
mag(K[l]);
99 for (label a=0; a<NoSubintervals; a++)
103 kmag <= ((a + 1)*delta_k + delta_k/2.0)
104 && kmag > ((a + 1)*delta_k - delta_k/2.0)
107 scalar
dist = delta_k/2.0 -
mag((a + 1)*delta_k - kmag);
125 for (label a=0; a<NoSubintervals; a++)
129 Ek1D[a] /= EWeight[a];
133 return graph(
"E(k)",
"k",
"E(k)", k1D, Ek1D);