def test_Probe_vectorfunctionspace_2D(VF2_self):
u0 = interpolate(Expression(('x[0]', 'x[1]'), degree=1), VF2_self)
x = array([0.5, 0.75])
p = Probe(x, VF2_self)
p(u0)
p(u0)
assert round(p[0][0] - 0.5, 7) == 0
assert round(p[1][1] - 0.75, 7) == 0
def test_Probe_functionspace_3D(V3_self):
u0 = interpolate(Expression('x[0]', degree=1), V3_self)
x = array([0.25, 0.5, 0.5])
p = Probe(x, V3_self)
p(u0)
p(u0)
assert round(p[0][0] - 0.25, 7) == 0
assert round(p[1][0] - 0.25, 7) == 0