54 scalar l0(
K.sizeOfBox()[0]*(scalar(
K.nn()[0])/(scalar(
K.nn()[0])-1.0)));
62 y /= scalar(
K.size());
77 const label tnp =
Ek.size();
78 const label NoSubintervals = label
80 pow(scalar(tnp), 1.0/vector::dim)*
pow(1.0/vector::dim, 0.5) - 0.5
87 scalar kmax =
K.max()*
pow(1.0/vector::dim,0.5);
88 scalar delta_k = kmax/(NoSubintervals);
92 k1D[a] = (a + 1)*delta_k;
99 scalar kmag =
mag(
K[l]);
101 for (label a=0; a<NoSubintervals; a++)
105 kmag <= ((a + 1)*delta_k + delta_k/2.0)
106 && kmag > ((a + 1)*delta_k - delta_k/2.0)
109 scalar dist = delta_k/2.0 -
mag((a + 1)*delta_k - kmag);
127 for (label a=0; a<NoSubintervals; a++)
131 Ek1D[a] /= EWeight[a];
135 return graph(
"E(k)",
"k",
"E(k)", k1D, Ek1D);