Boundary value problems on large periodic networks arise in many applications such as soil mechanics in geophysics or the analysis of photonic crystals in nanotechnology. As a model example, singularly perturbed elliptic differential equations of second order are addressed. Typically, the length of periodicity is very small compared to the size of the covered region. The overall complexity of the networks raises serious problems on the computational side. The high density of the graph, the huge number of edges and vertices and highly oscillating coefficients necessitate solution schemes, where even a numerical approximation is no longer feasible. Realizing that such a system depends on two spatial scales - global scale (full domain) and loc...