A cache optimized multilevel algorithm to solve variable coefficient elliptic boundary value problems on adaptively refined structured meshes is described here. The algorithm is optimized to exploit the cache memory subsystem. Numerical results are given demonstrating the efficiency of the cache optimization.