Nonlinear multigrid and scalar auxiliary variable solvers of the Cahn-Hilliard equation (MATLAB, Julia, Python)