def test_GantryDominant2D_adjacency_graph(self):
nodes = [[1, 2, 3], [-1, 0, 1]]
lattice = Lattice(nodes)
ratio = 1.0
sampler = GantryDominant2D(lattice, ratio=ratio)
graph_true = np.array([
[0, 1, 0, 0, 1, 0, 0, 0, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
[0, 1, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
[0, 0, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 0, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
],
dtype=bool)
graph_test = sampler._adjacency_graph()
self.assertTrue(np.all(graph_true == graph_test))
ratio = 2.0
sampler = GantryDominant2D(lattice, ratio=ratio)
graph_true = np.array([
[0, 1, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
[0, 1, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
[0, 1, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
],
dtype=bool)
graph_test = sampler._adjacency_graph()
self.assertTrue(np.all(graph_true == graph_test))
nodes = [[1, 2, 3], [-0.5, 0, 0.5]]
lattice = Lattice(nodes)
ratio = 1.0
sampler = GantryDominant2D(lattice, ratio=ratio)
graph_true = np.array([
[0, 1, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
[0, 1, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
[0, 1, 0, 0, 1, 0, 0, 1, 0],
[0, 0, 1, 0, 0, 1, 0, 0, 1],
[0, 0, 0, 0, 0, 0, 0, 0, 0],
],
dtype=bool)
graph_test = sampler._adjacency_graph()
self.assertTrue(np.all(graph_true == graph_test))
def test_Lattice3D_print(self):
lattice = Lattice(self.nodes3D)
self.assertTrue(hasattr(lattice, '__repr__'))
self.assertTrue(lattice.__repr__())
self.assertTrue(hasattr(lattice, '__str__'))
self.assertTrue(lattice.__str__())
def test_Lattice3D__getitem__(self):
lattice = Lattice(self.nodes3D)
nodes = lattice.nodes
nnodes_dim = lattice.nnodes_dim
nx_st, nx_end = 2, 4
ny_st, ny_end = 0, 3
nz_st, nz_end = 2, 3
lattice_z = lattice[:, :, nz_st]
lattice_zslice = lattice[:, :, nz_st:nz_end]
lattice_xyzslice = lattice[nx_st:nx_end, ny_st:ny_end, nz_st:nz_end]
self.assertEqual(lattice_z.nnodes_dim, (nnodes_dim[0], nnodes_dim[1]))
self.assertEqual(lattice_z, Lattice([nodes[0], nodes[1]]))
self.assertEqual(lattice_zslice.nnodes_dim,
(nnodes_dim[0], nnodes_dim[1], 1))
self.assertEqual(lattice_zslice,
Lattice([nodes[0], nodes[1], nodes[2][nz_st:nz_end]]))
self.assertEqual(
lattice_xyzslice,
Lattice([
nodes[0][nx_st:nx_end], nodes[1][ny_st:ny_end],
nodes[2][nz_st:nz_end]
]))
def test_Lattice2D__eq__(self):
lattice = Lattice(self.nodes2D)
lattice_is = lattice
lattice_eq = Lattice(self.nodes2D)
lattice_non_eq = Lattice([self.x_non_eq, self.y])
lattice_short = Lattice([self.x_short, self.y])
self.assertTrue(lattice is lattice_is)
self.assertTrue(lattice == lattice_eq)
self.assertEqual(lattice, lattice_eq)
self.assertFalse(lattice is lattice_eq)
self.assertFalse(lattice == lattice_non_eq)
self.assertFalse(lattice == lattice_short)
def test_Trajectory_to_latticemap(self):
nodes = [[-2, -1, 0, 1, 2], [-3, 1, 5, 9]]
lattice = Lattice(nodes)
pts = [(-2, 1), (-1, -4), (0, 3), (1, 4), (2, 8.8)]
values = np.array([
0,
1,
0,
0,
0,
1,
0,
1,
0,
0,
0,
0,
0,
1,
0,
0,
0,
0,
0,
1,
],
dtype=bool)
map_test = Trajectory(pts).to_latticemap(lattice)
map_true = LatticeMap(lattice, values)
self.assertEqual(map_true, map_test)
def test_GantryDominant2D__call__conditioning(self):
# Computational lattice
nodes = [
[-4, -3, -2, -1, 0, 1, 2, 3, 4],
[0, 1, 2, 3, 4, 5, 6],
]
lattice = Lattice(nodes)
# Sampler
ratio = 2.0
sampler = GantryDominant2D(lattice=lattice, ratio=ratio)
# Sample Trajectory
trajectory = sampler(conditions=[(0, 3)])
# Check length
self.assertTrue(len(trajectory), len(nodes[0]))
# Check conditioing
self.assertEqual(trajectory[4], (0, 3))
# Check max trajectory spreading
pt_last = trajectory[0]
self.assertEqual(pt_last[0], nodes[0][0])
for pt_next, nd_next in zip(trajectory[1:], nodes[0][1:]):
self.assertEqual(pt_next[0], nd_next)
max_spread = abs(pt_next[0] - pt_last[0]) * ratio
max_spread += NP_ATOL + NP_RTOL * max_spread
self.assertTrue(abs(pt_next[1] - pt_last[1]) <= max_spread)
pt_last = pt_next
def test_GantryDominant2D_gen(self):
nodes = [[-1.5, 1.5, 5, 8, 9], [1, 2, 3, 4, 5, 6]]
lattice = Lattice(nodes)
sampler = GantryDominant2D(lattice, ratio=2.0)
self.assertTrue(isinstance(sampler, Sampler))
self.assertTrue(isinstance(sampler, GantryDominant2D))
def test_Lattice2D_ndim(self):
ndim = 2
lattice = Lattice(self.nodes2D)
self.assertTrue(hasattr(lattice, 'ndim'))
self.assertEqual(lattice.ndim, ndim)
self.assertTrue(isinstance(lattice.ndim, int))
def test_LatticeMap2D_slice(self):
nodes = [[-1.5, 1.5, 5, 8, 9], [1, 2, 3, 4, 5, 6]]
lattice = Lattice(nodes)
values = [
0.51,
0.52,
0.53,
0.54,
0.55,
0.61,
0.62,
0.63,
0.64,
0.65,
0.71,
0.72,
0.73,
0.74,
0.75,
0.71,
0.72,
0.73,
0.74,
0.75,
0.61,
0.62,
0.63,
0.64,
0.65,
0.51,
0.52,
0.53,
0.54,
0.55,
]
latticemap = LatticeMap(lattice, values=values)
lattice_slice_0 = Lattice([[1, 2, 3, 4, 5, 6]])
values_slice_0 = [0.52, 0.62, 0.72, 0.72, 0.62, 0.52]
latticemap_slice_0 = LatticeMap(lattice_slice_0, values_slice_0)
lattice_slice_1 = Lattice([[-1.5, 1.5, 5, 8, 9]])
values_slice_1 = [0.71, 0.72, 0.73, 0.74, 0.75]
latticemap_slice_1 = LatticeMap(lattice_slice_1, values_slice_1)
self.assertEqual(latticemap_slice_0, latticemap.slice(0, 1))
self.assertEqual(latticemap_slice_1, latticemap.slice(1, 3))
def __init__(self, lattice):
# Also accept objects that can be used to generate a lattice
if not isinstance(lattice, Lattice):
lattice = Lattice(lattice)
self._lattice = lattice
# Make default prior maps (all equal probability, all allowed)
ones = np.ones(lattice.nnodes_dim)
self._prior_prob = LatticeMap(lattice, ones, dtype=int)
self._prior_cond = LatticeMap(lattice, ones, dtype=bool)
def test_LatticeMap3D_ndim(self):
ndim = 3
lattice = Lattice(self.nodes3D)
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
self.assertTrue(hasattr(latticemap, 'ndim'))
self.assertEqual(latticemap.ndim, ndim)
self.assertTrue(isinstance(latticemap.ndim, int))
def test_LatticeMap2D_print(self):
lattice = Lattice(self.nodes2D)
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
self.assertTrue(hasattr(latticemap, '__repr__'))
self.assertTrue(latticemap.__repr__())
self.assertTrue(hasattr(latticemap, '__str__'))
self.assertTrue(latticemap.__str__())
def test_Lattice2D_nnodes(self):
ndim = 2
nnodes_dim = (6, 4)
nnodes = 24
lattice = Lattice(self.nodes2D)
self.assertTrue(isinstance(lattice.nnodes_dim, tuple))
self.assertEqual(len(lattice.nnodes_dim), ndim)
self.assertEqual(lattice.nnodes_dim, nnodes_dim)
self.assertTrue(isinstance(lattice.nnodes, int))
self.assertEqual(lattice.nnodes, nnodes)
def test_LatticeMap3D_print(self):
lattice = Lattice(self.nodes3D)
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
# Force the implementation of __repr__() and __str__()
self.assertTrue(hasattr(latticemap, '__repr__'))
self.assertTrue(latticemap.__repr__())
self.assertTrue(hasattr(latticemap, '__str__'))
self.assertTrue(latticemap.__str__())
def test_Lattice3D_nnodes(self):
ndim = 3
nnodes_dim = (6, 4, 3)
nnodes = 72
lattice = Lattice(self.nodes3D)
self.assertTrue(isinstance(lattice.nnodes_dim, tuple))
self.assertEqual(len(lattice.nnodes_dim), ndim)
self.assertEqual(lattice.nnodes_dim, nnodes_dim)
self.assertTrue(isinstance(lattice.nnodes, int))
self.assertEqual(lattice.nnodes, nnodes)
def test_Lattice2D__getitem__(self):
lattice = Lattice(self.nodes2D)
nodes = lattice.nodes
nnodes_dim = lattice.nnodes_dim
nx_st, nx_end = 2, 4
ny_st, ny_end = 0, 3
lattice_x = lattice[nx_st, :]
lattice_xslice = lattice[nx_st:nx_end, :]
lattice_y = lattice[:, ny_st]
lattice_yslice = lattice[:, ny_st:ny_end]
lattice_xyslice = lattice[nx_st:nx_end, ny_st:ny_end]
self.assertEqual(lattice_x.nnodes_dim, (nnodes_dim[1], ))
self.assertEqual(lattice_x, Lattice([nodes[1]]))
self.assertEqual(lattice_xslice.nnodes_dim, (2, nnodes_dim[1]))
self.assertEqual(lattice_xslice,
Lattice([nodes[0][nx_st:nx_end], nodes[1]]))
self.assertEqual(lattice_y.nnodes_dim, (nnodes_dim[0], ))
self.assertEqual(lattice_y, Lattice([nodes[0]]))
self.assertEqual(lattice_yslice.nnodes_dim, (nnodes_dim[0], 3))
self.assertEqual(lattice_yslice,
Lattice([nodes[0], nodes[1][ny_st:ny_end]]))
self.assertEqual(
lattice_xyslice,
Lattice([nodes[0][nx_st:nx_end], nodes[1][ny_st:ny_end]]))
def test_LatticeMap3D__add__(self):
lattice = Lattice(self.nodes3D)
nodes_short = [self.x_short, self.y, self.z]
lattice_short = Lattice(nodes_short)
num = -25.89
values_left = np.random.randn(lattice.nnodes)
values_right = np.random.randn(lattice.nnodes)
values_short = np.random.randn(lattice_short.nnodes)
latticemap_left = LatticeMap(lattice, values_left)
latticemap_right = LatticeMap(lattice, values_right)
latticemap_short = LatticeMap(lattice_short, values_short)
latticmap_sum = LatticeMap(lattice, values_left + values_right)
latticmap_sum_num = LatticeMap(lattice, values_left + num)
self.assertEqual(latticmap_sum, latticemap_left + latticemap_right)
self.assertTrue(latticmap_sum == latticemap_left + latticemap_right)
self.assertEqual(latticmap_sum_num, latticemap_left + num)
self.assertTrue(latticmap_sum_num == latticemap_left + num)
with self.assertRaises(ValueError):
latticemap_left + latticemap_short
with self.assertRaises(TypeError):
latticemap_left + 'foobar'
def test_LatticeMap3D__eq__(self):
lattice = Lattice(self.nodes3D)
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
latticemap_is = latticemap
latticemap_eq = LatticeMap(lattice, values)
latticemap_non_eq = LatticeMap(lattice, values + 1)
self.assertTrue(latticemap is latticemap_is)
self.assertTrue(latticemap == latticemap_eq)
self.assertEqual(latticemap, latticemap_eq)
self.assertFalse(latticemap is latticemap_eq)
self.assertFalse(latticemap == latticemap_non_eq)
def test_LatticeMap2D__getitem__(self):
lattice = Lattice(self.nodes2D)
nodes = lattice.nodes
nnodes_dim = lattice.nnodes_dim
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
nx_st, nx_end = 2, 4
ny_st, ny_end = 0, 3
ind_x = np.arange(nx_st, nx_end)
ind_y = np.arange(ny_st, ny_end)
latticemap_x = latticemap[nx_st, :]
lin_ind_x = np.array(
[nx_st + iy * nnodes_dim[0] for iy in range(nnodes_dim[1])])
latticemap_xslice = latticemap[nx_st:nx_end, :]
lin_ind_xslice = np.array([
ix + iy * nnodes_dim[0] for iy in range(nnodes_dim[1])
for ix in ind_x
])
latticemap_y = latticemap[:, ny_st]
lin_ind_y = np.array(
[ix + ny_st * nnodes_dim[0] for ix in range(nnodes_dim[0])])
latticemap_yslice = latticemap[:, ny_st:ny_end]
lin_ind_yslice = np.array([
ix + iy * nnodes_dim[0] for iy in ind_y
for ix in range(nnodes_dim[0])
])
latticemap_xyslice = latticemap[nx_st:nx_end, ny_st:ny_end]
lin_ind_xyslice = np.array(
[ix + iy * nnodes_dim[0] for iy in ind_y for ix in ind_x])
self.assertTrue(latticemap_x == LatticeMap(lattice=lattice[nx_st, :],
values=values[lin_ind_x]))
self.assertTrue(latticemap_xslice == LatticeMap(
lattice=lattice[nx_st:nx_end, :], values=values[lin_ind_xslice]))
self.assertTrue(latticemap_y == LatticeMap(lattice=lattice[:, ny_st],
values=values[lin_ind_y]))
self.assertTrue(latticemap_yslice == LatticeMap(
lattice=lattice[:, ny_st:ny_end], values=values[lin_ind_yslice]))
self.assertTrue(
latticemap_xyslice == LatticeMap(lattice=lattice[nx_st:nx_end,
ny_st:ny_end],
values=values[lin_ind_xyslice]))
def test_LatticeMap3D__mul__(self):
lattice = Lattice(self.nodes3D)
num = -25.89
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
latticmap_mul = LatticeMap(lattice, values * num)
latticemap_squared = LatticeMap(lattice, values**2)
self.assertEqual(latticmap_mul, latticemap * num)
self.assertTrue(latticmap_mul == latticemap * num)
self.assertEqual(latticmap_mul, num * latticemap)
self.assertTrue(latticmap_mul == num * latticemap)
self.assertEqual(latticemap * latticemap, latticemap_squared)
with self.assertRaises(TypeError):
latticemap * 'foobar'
def test_LatticeMap3D__getitem__(self):
lattice = Lattice(self.nodes3D)
nodes = lattice.nodes
nnodes_dim = lattice.nnodes_dim
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
nx_st, nx_end = 2, 4
ny_st, ny_end = 0, 3
nz_st, nz_end = 2, 3
ind_x = np.arange(nx_st, nx_end)
ind_y = np.arange(ny_st, ny_end)
ind_z = np.arange(nz_st, nz_end)
latticemap_z = latticemap[:, :, nz_st]
lin_ind_z = np.array([
ix + iy * nnodes_dim[0] + nz_st * nnodes_dim[0] * nnodes_dim[1]
for iy in range(nnodes_dim[1]) for ix in range(nnodes_dim[0])
])
latticemap_zslice = latticemap[:, :, nz_st:nz_end]
lin_ind_zslice = np.array([
ix + iy * nnodes_dim[0] + iz * nnodes_dim[0] * nnodes_dim[1]
for iz in ind_z for iy in range(nnodes_dim[1])
for ix in range(nnodes_dim[0])
])
latticemap_xyzslice = latticemap[nx_st:nx_end, ny_st:ny_end,
nz_st:nz_end]
lin_ind_xyzslice = np.array([
ix + iy * nnodes_dim[0] + iz * nnodes_dim[0] * nnodes_dim[1]
for iz in ind_z for iy in ind_y for ix in ind_x
])
self.assertTrue(latticemap_z == LatticeMap(lattice=lattice[:, :,
nz_st],
values=values[lin_ind_z]))
self.assertTrue(latticemap_zslice == LatticeMap(
lattice=lattice[:, :,
nz_st:nz_end], values=values[lin_ind_zslice]))
self.assertTrue(
latticemap_xyzslice == LatticeMap(lattice=lattice[nx_st:nx_end,
ny_st:ny_end,
nz_st:nz_end],
values=values[lin_ind_xyzslice]))
def test_LatticeMap2D_gen(self):
lattice = Lattice(self.nodes2D)
values = np.random.randn(lattice.nnodes)
latticemap = LatticeMap(lattice, values)
self.assertTrue(isinstance(latticemap, LatticeMap))
self.assertTrue(np.all(values == latticemap.values))
with self.assertRaises(ValueError):
LatticeMap(lattice, values[:-1])
with self.assertRaises(ValueError):
LatticeMap(lattice, [])
latticemap = LatticeMap(self.nodes2D, values)
self.assertTrue(isinstance(latticemap, LatticeMap))
self.assertTrue(np.all(values == latticemap.values))
with self.assertRaises(ValueError):
LatticeMap(self.nodes2D, values[:-1])
with self.assertRaises(ValueError):
LatticeMap(self.nodes2D, [])
def test_LatticeMap2D_make_latticemap_from_txt(self):
nodes = [[-1.5, 1.5, 5, 8, 9], [1, 2, 3, 4, 5, 6]]
lattice = Lattice(nodes)
values = [
0.5,
0.5,
0.5,
0.5,
0.5,
0.6,
0.6,
0.6,
0.6,
0.6,
0.7,
0.7,
0.7,
0.7,
0.7,
0.7,
0.7,
0.7,
0.7,
0.7,
0.6,
0.6,
0.6,
0.6,
0.6,
0.5,
0.5,
0.5,
0.5,
0.5,
]
latticemap_true = LatticeMap(lattice, values)
file = PATH_TESTFILES + 'latticemap2d_simple.txt'
latticemap_test = LatticeMap.from_txt(file)
self.assertTrue(isinstance(latticemap_test, LatticeMap))
self.assertEqual(latticemap_true, latticemap_test)
def test_Lattice3D_indices_from_point(self):
lattice = Lattice(self.nodes3D)
components = (2, -0.5, 1000)
ind_point = lattice.indices_from_point(components)
self.assertEqual(ind_point, (1, 1, 2))
components = (1100, -0.5, -100)
ind_point = lattice.indices_from_point(components)
self.assertEqual(ind_point, (5, 1, 0))
components = (2, -0.499, -11001.55)
ind_point = lattice.indices_from_point(components)
self.assertEqual(ind_point, (1, 2, 0))
components = (2, -0.5)
with self.assertRaises(ValueError):
lattice.indices_from_point(components)
def test_Lattice2D_indices(self):
lattice = Lattice(self.nodes2D)
components = (2, -0.5)
ind_point = lattice.indices_from_point(components)
self.assertEqual(ind_point, (1, 1))
components = (1100, -0.5)
ind_point = lattice.indices_from_point(components)
self.assertEqual(ind_point, (5, 1))
components = (2, -0.499)
ind_point = lattice.indices_from_point(components)
self.assertEqual(ind_point, (1, 2))
components = (2, -0.5, -1)
with self.assertRaises(ValueError):
lattice.indices_from_point(components)
def _cond_map_from_str(self, file):
"""Loads a conditioning map from a saved file.
!!! ATTENTION !!!
This map invertes the values. AMS is used to indicate blocked nodes
by `1` and allowed nodes by `0`. However, for the sampling it is more
convenient to use it the other way around. Therefore, the values are
inverted withing this function.
"""
# Read the conditioning file
nodes, values = readfile_nodes_values(file)
# Invert the values (c.f. docstring)
values = utl.ams_val_map_to_bool_map(values)
# Check consistency of the lattices
if Lattice(nodes) != self._lattice:
raise ValueError(msg.err3001(file, self._lattice))
return LatticeMap(nodes, values)
def test_GantryDominant2D_set_prior_cond_with_str(self):
nodes = [[1, 2, 3], [4, 5, 6, 7, 8, 9, 10]]
lattice = Lattice(nodes)
conditions = [PATH_TESTFILES + 'condmap2d_simple.txt']
values_true = np.array([
1,
1,
1,
1,
1,
1,
1,
1,
1,
1,
1,
0,
1,
0,
0,
0,
0,
0,
0,
0,
0,
],
dtype=bool)
prior_cond_true = LatticeMap(nodes, values_true)
ratio = 1.0
sampler = GantryDominant2D(lattice, ratio=ratio)
sampler.set_prior_cond(conditions)
self.assertEqual(sampler._prior_cond, prior_cond_true)
def test_LatticeMap3D_make_latticemap_from_txt(self):
nodes = [[-1.5, 1.5], [5, 8, 9], [-2, 3]]
lattice = Lattice(nodes)
values = [
0.5,
0.5,
0.8,
0.8,
0.1,
0.1,
0.6,
0.6,
0.9,
0.9,
0.2,
0.2,
]
latticemap_true = LatticeMap(lattice, values)
file = PATH_TESTFILES + 'latticemap3d_simple.txt'
latticemap_test = LatticeMap.from_txt(file)
self.assertTrue(isinstance(latticemap_test, LatticeMap))
self.assertEqual(latticemap_true, latticemap_test)
def test_GantryDominant2D__cond_map_from_point(self):
nodes = [[-2, -1, 0, 1, 2], [-2, 0, 2]]
lattice = Lattice(nodes)
conditions = [(0, 0)]
values_true = np.array([
1,
1,
0,
1,
1,
1,
1,
1,
1,
1,
1,
1,
0,
1,
1,
],
dtype=bool)
prior_cond_true = LatticeMap(nodes, values_true)
ratio = 2.0
sampler = GantryDominant2D(lattice, ratio=ratio)
sampler.set_prior_cond(conditions)
self.assertEqual(sampler._prior_cond, prior_cond_true)
conditions = [(-2, -2)]
values_true = np.array([
1,
1,
1,
1,
1,
0,
1,
1,
1,
1,
0,
0,
1,
1,
1,
],
dtype=bool)
prior_cond_true = LatticeMap(nodes, values_true)
ratio = 2.0
sampler = GantryDominant2D(lattice, ratio=ratio)
sampler.set_prior_cond(conditions)
self.assertEqual(sampler._prior_cond, prior_cond_true)
conditions = [(1, -2)]
values_true = np.array([
1,
1,
1,
1,
1,
1,
1,
1,
0,
1,
1,
1,
1,
0,
1,
],
dtype=bool)
prior_cond_true = LatticeMap(nodes, values_true)
ratio = 4.0
sampler = GantryDominant2D(lattice, ratio=ratio)
sampler.set_prior_cond(conditions)
self.assertEqual(sampler._prior_cond, prior_cond_true)
def test_Lattice2D_gen(self):
lattice = Lattice(self.nodes2D)
self.assertTrue(isinstance(lattice, Lattice))
with self.assertRaises(ValueError):
Lattice([[1, 2, 3], [1, -1, 2]])