We present a number of test cases and meshes that were designed as a benchmark for numerical schemes dedicated to the approximation of threedimensional anisotropic and heterogeneous diffusion problems. These numerical schemes may be applied to general, possibly non conforming, meshes composed of tetrahedra, hexahedra and quite distorted general polyhedra. A number of methods were tested among which conforming finite element methods, discontinuous Galerkin finite element methods, cell-centered finite volume methods, discrete duality finite volume methods, mimetic finite difference methods, mixed finite element methods, and gradient schemes. We summarize the results presented by the participants to the benchmark, which range from the number of unknowns, the approximation errors of the solution and its gradient, to the minimum and maximum values and energy. We also compare the performance of several iterative or direct linear solvers for the resolution of the linear systems issued from the presented schemes.