def METHOD_NAME():
    print ('>>> test_2d_3')
    p = 2
    grid_1 = linspace(0., 1., 3)
    grid_2 = linspace(0., 1., 5)
    # ... first component
    V1 = SplineSpace(p-1, grid=grid_1)
    V2 = SplineSpace(p, grid=grid_2)
    domain_decomposition = DomainDecomposition([V1.ncells, V2.ncells], [False, False])
    Vx = TensorFemSpace(domain_decomposition, V1, V2)
    # ...
    # ... second component
    V1 = SplineSpace(p, grid=grid_1)
    V2 = SplineSpace(p-1, grid=grid_2)
    Vy = TensorFemSpace(domain_decomposition, V1, V2)
    # ...
    V = VectorFemSpace(Vx, Vy)
    print (V)