Before revision:  A self-contained Python program implements the scheme , and its optimal performanceis demonstrated .
After revision:  A self-contained Python program implements the scheme . Optimal performance, i.e. work proportional to the number of unknowns, is demonstrated for both kinds of cycles, including convergence nearly to discretization error in a single F-cycle .