A linearity-preserving diamond scheme with extended least square interpolation for the heterogeneous and anisotropic diffusion problems on polyhedral meshes