[anishtain4]

Assume you have a scalarField or vectorField named fi(K), and you want to compute its average value in a spherical shell with thickness of one, say:

sigma( magSqr(fi(K)) )
i-0.5<k<i+0.5
k=sqrt(sqr(kx)+sqr(ky)+sqr(kz))

I know you can make a loop or use forAll but this way's gonna check every label once for every shell and that is pretty time consuming. Is there any better way?

[gschaider]

Originally Posted by anishtain4

Assume you have a scalarField or vectorField named fi(K), and you want to compute its average value in a spherical shell with thickness of one, say:

sigma( magSqr(fi(K)) )
i-0.5<k<i+0.5
k=sqrt(sqr(kx)+sqr(ky)+sqr(kz))

I know you can make a loop or use forAll but this way's gonna check every label once for every shell and that is pretty time consuming. Is there any better way?

That is a basic programming: "memory usage buys performance" Create an array for the sums, one for every shell. Loop once over the cells. Find out to which shell the cell belongs. Increase that array. Of course you'll need a second array .... but that is left to the reader as an exercise