2021-10-15 12:56:07 +01:00
|
|
|
import numpy as np
|
|
|
|
|
import pyerrors as pe
|
|
|
|
|
import pytest
|
|
|
|
|
|
2021-10-15 13:05:00 +01:00
|
|
|
np.random.seed(0)
|
|
|
|
|
|
2021-10-15 14:13:06 +01:00
|
|
|
|
2021-10-15 13:05:00 +01:00
|
|
|
def test_root_linear():
|
2021-10-15 12:56:07 +01:00
|
|
|
|
|
|
|
|
def root_function(x, d):
|
|
|
|
|
return x - d
|
|
|
|
|
|
|
|
|
|
value = np.random.normal(0, 100)
|
|
|
|
|
my_obs = pe.pseudo_Obs(value, 0.1, 't')
|
|
|
|
|
my_root = pe.roots.find_root(my_obs, root_function)
|
|
|
|
|
|
|
|
|
|
assert np.isclose(my_root.value, value)
|
2021-10-15 13:05:00 +01:00
|
|
|
difference = my_obs - my_root
|
2021-10-15 15:01:53 +01:00
|
|
|
assert np.allclose(0.0, difference.deltas['t'])
|
