A brief survey is given to show that harmonic averages enter in a natural way in the numerical solution of various variable coefficient problems, such as in elliptic and transport equations, also of singular perturbation types. Local Green’s functions used as test functions in the Petrov-Galerkin finite element method combined with harmonic averages can be very efficient and are related to exact difference schemes.