Exemple #1
0
 def test_compare_pushforward_pullback_derivatives2(self):
     x = UTPS(numpy.array([1.,0.,0.,0.]), P = 1)
     y = UTPS(numpy.array([2.,1.,0.,0.]), P = 1)
     z = div(x,y)
     
     zbar = UTPS([1.,0.,0., 0.], P = 1)
     xbar = UTPS(numpy.zeros(4), P = 1)
     ybar = UTPS(numpy.zeros(4), P = 1)
     
     epb_div(x,y,z, zbar, xbar, ybar)
     
     facs = numpy.array([1.,1.,2.,6.])
     
     assert_array_almost_equal((z.data*facs)[1:], (ybar.data*facs)[:-1])
Exemple #2
0
 def test_div_epb(self):
     x = UTPS(numpy.array([1.,2.,3.]), P = 1)
     y = UTPS(numpy.array([5.,7.,11.]), P = 1)
     z = div(x,y)
     
     zbar = UTPS(numpy.random.rand(3), P = 1)
     xbar = UTPS(numpy.zeros(3), P = 1)
     ybar = UTPS(numpy.zeros(3), P = 1)
     
     xbar2 = zbar/y
     ybar2 = - zbar*z/y
     epb_div(x,y,z, zbar, xbar, ybar)
     
     assert_array_almost_equal(xbar2.data, xbar.data)
     assert_array_almost_equal(ybar2.data, ybar.data)