def METHOD_NAME():
    print ('>>> test_2d_2')
    p_1 = 2
    p_2 = 2
    grid_1 = linspace(0., 1., 3)
    grid_2 = linspace(0., 1., 5)
    V1 = SplineSpace(p_1, grid=grid_1)
    V2 = SplineSpace(p_2, grid=grid_2)
    domain_decomposition = DomainDecomposition([V1.ncells, V2.ncells], [False, False])
    V = TensorFemSpace(domain_decomposition, V1, V2)
    print (V)
    F = FemField(V)