def check_float(self):
# Axis -1
xt = xcf
axis = -1
order = 1
a, e, k = lpc2(xt, order, axis)
at, et, kt = lpc_py(xt, order, axis)
assert a.dtype == e.dtype == k.dtype == numpy.float32
assert_array_almost_equal(a, at)
assert_array_almost_equal(e, et)
assert_array_almost_equal(k, kt)
# Axis 0
xt = xcf
axis = 0
order = 1
a, e, k = lpc2(xt, order, axis)
at, et, kt = lpc_py(xt, order, axis)
assert a.dtype == e.dtype == k.dtype == numpy.float32
assert_array_almost_equal(a, at)
assert_array_almost_equal(e, et)
assert_array_almost_equal(k, kt)
def check_double_rank1(self):
# test rank 1
xt = xc[0]
axis = 0
order = 1
a, e, k = lpc2(xt, order, axis)
at, et, kt = lpc_py(xt, order, axis)
assert_array_almost_equal(a, at)
assert_array_almost_equal(e, et)
assert_array_almost_equal(k, kt)
def check_float_rank1(self):
# test rank 1
xt = xcf[0]
axis = 0
order = 1
a, e, k = lpc2(xt, order, axis)
at, et, kt = lpc_py(xt, order, axis)
assert a.dtype == e.dtype == k.dtype == numpy.float32
assert_array_almost_equal(a, at)
assert_array_almost_equal(e, et)
assert_array_almost_equal(k, kt)
def check_double(self):
# Axis -1
xt = xc
axis = -1
order = 1
a, e, k = lpc2(xt, order, axis)
at, et, kt = lpc_py(xt, order, axis)
assert_array_almost_equal(a, at)
assert_array_almost_equal(e, et)
assert_array_almost_equal(k, kt)
# Axis 0
xt = xc
axis = 0
order = 1
a, e, k = lpc2(xt, order, axis)
at, et, kt = lpc_py(xt, order, axis)
assert_array_almost_equal(a, at)
assert_array_almost_equal(e, et)
assert_array_almost_equal(k, kt)
