def METHOD_NAME(quads=False):
    mesh = Mesh(unit_square.GenerateMesh(quad=quads))
    for spacename in spaces2d.keys():
        if not quads or (quads and spaces2d[spacename]["quad"]):
            for order in spaces2d[spacename]["order"]:
                space = FESpace(type=spacename,mesh=mesh,order=order)
                for vb in spaces2d[spacename]["vorb"]:
                    for el in space.Elements(vb):
                        assert space.GetFE(el).ndof == len(space.GetDofNrs(el)), [spacename,vb,order]
    return