Implement function to query if an element has derivatives of order n in direction d

[Martin Sandve Alnæs]

The naive implementation would be just:

def has_derivatives(self, order, direction):
    return self.degree() >= order

but a product element may have different degrees in each spatial direction.