The gradient recovery method proposed by Zhu and Zienkiewicz for one-d
imensional problems and extended to two dimensions by Silvester and Om
eragic is generalized to three-dimensional solutions based on rectangu
lar prism (brick) elements. The extension is not obvious so its detail
s are presented, and the method compared with conventional local smoot
hing and direct differentiation. Illustrative examples are given, with
an extensive experimental study of error. The method is computational
ly cheap and provides better accuracy than conventional local smoothin
g, but its accuracy is position dependent.