def METHOD_NAME():
    domains_with_quadratic_coordinates = []
    for D in get_domains():
        V = VectorElement("CG", D.ufl_cell(), 2)
        E = Mesh(V)
        domains_with_quadratic_coordinates.append(E)
    return domains_with_quadratic_coordinates