def test_Rectangle_non_Function(self):
for order in range(2, 11):
coupler = SpeckleyToRipley(2,
(2 * getMPISizeWorld() * order, order),
order=order,
lengths=[3. * getMPISizeWorld(), 2.])
self.badInterpolations(coupler.getSpeckley(), coupler.getRipley())
def test_Brick_non_Function(self):
for order in range(2, 11):
#values here are arbitrary, just has to be Bricks
coupler = SpeckleyToRipley(
3, (2 * getMPISizeWorld() * order, order, order),
order=order,
lengths=[3. * getMPISizeWorld(), 2., 2.])
self.badInterpolations(coupler.getSpeckley(), coupler.getRipley())
def test_Brick_non_Function(self):
for order in range(2, 11):
# values here are arbitrary, just has to be Bricks
coupler = SpeckleyToRipley(
3,
(2 * getMPISizeWorld() * order, order, order),
order=order,
lengths=[3.0 * getMPISizeWorld(), 2.0, 2.0],
)
self.badInterpolations(coupler.getSpeckley(), coupler.getRipley())
def test_Brick(self):
elements = [3 * getMPISizeWorld(), 2, 1]
lengths = [3.0 * getMPISizeWorld(), 4.0, 5.0]
labels = ["x", "y", "z"]
error_message = "Brick {0}-split of order {1} failed to interpolate " + "correctly, variance of {2}"
for i in range(len(labels)):
e = tuple(elements[i:] + elements[:i])
l = tuple(lengths[i:] + lengths[:i])
for order in range(2, 11):
ele = tuple([j * order for j in e])
coupler = SpeckleyToRipley(3, ele, order=order, lengths=l)
a = self.calculateVariance(coupler.getSpeckley(), coupler.getRipley())
self.assertLess(Lsup(a), 1e-10, error_message.format(labels[i], order, Lsup(a)))
def test_creation(self):
with self.assertRaises(RuntimeError): #invalid order
coupler = SpeckleyToRipley(2,
pointsPerDim=(4, 2 * getMPISizeWorld()),
order=1)
with self.assertRaises(RuntimeError): #invalid order
coupler = SpeckleyToRipley(2,
pointsPerDim=(40,
20 * getMPISizeWorld()),
order=11)
with self.assertRaises(ValueError): #invalid dimension
coupler = SpeckleyToRipley(1,
pointsPerDim=(4, 2 * getMPISizeWorld()))
with self.assertRaises(ValueError): #incalid dimension
coupler = SpeckleyToRipley(4,
pointsPerDim=(4, 2 * getMPISizeWorld()))
def test_Brick(self):
elements = [3 * getMPISizeWorld(), 2, 1]
lengths = [3. * getMPISizeWorld(), 4., 5.]
labels = ["x", "y", "z"]
error_message = "Brick {0}-split of order {1} failed to interpolate " +\
"correctly, variance of {2}"
for i in range(len(labels)):
e = tuple(elements[i:] + elements[:i])
l = tuple(lengths[i:] + lengths[:i])
for order in range(2, 11):
ele = tuple([j * order for j in e])
coupler = SpeckleyToRipley(3, ele, order=order, lengths=l)
a = self.calculateVariance(coupler.getSpeckley(),
coupler.getRipley())
self.assertLess(
Lsup(a), 1e-10,
error_message.format(labels[i], order, Lsup(a)))
def test_Rectangle(self):
for order in range(2, 11):
coupler = SpeckleyToRipley(2,
(2 * getMPISizeWorld() * order, order),
order=order,
lengths=[3. * getMPISizeWorld(), 2.])
d = self.calculateVariance(coupler.getSpeckley(),
coupler.getRipley())
res = Lsup(d)
self.assertLess(res, 1e-10,
"Rectangles of order %d failed to interpolate correctly"%order +\
", variance of %g"%res)
coupler = SpeckleyToRipley(2,
(order, 2 * getMPISizeWorld() * order),
order=order,
lengths=[2., 3. * getMPISizeWorld()])
a = self.calculateVariance(coupler.getSpeckley(),
coupler.getRipley())
self.assertLess(Lsup(a), 1e-10,
"Rectangles of order %d failed to interpolate correctly"%order +\
", variance of %g"%Lsup(a))
def test_Rectangle(self):
for order in range(2, 11):
coupler = SpeckleyToRipley(
2, (2 * getMPISizeWorld() * order, order), order=order, lengths=[3.0 * getMPISizeWorld(), 2.0]
)
d = self.calculateVariance(coupler.getSpeckley(), coupler.getRipley())
res = Lsup(d)
self.assertLess(
res, 1e-10, "Rectangles of order %d failed to interpolate correctly" % order + ", variance of %g" % res
)
coupler = SpeckleyToRipley(
2, (order, 2 * getMPISizeWorld() * order), order=order, lengths=[2.0, 3.0 * getMPISizeWorld()]
)
a = self.calculateVariance(coupler.getSpeckley(), coupler.getRipley())
self.assertLess(
Lsup(a),
1e-10,
"Rectangles of order %d failed to interpolate correctly" % order + ", variance of %g" % Lsup(a),
)
def test_Rectangle_non_Function(self):
for order in range(2, 11):
coupler = SpeckleyToRipley(
2, (2 * getMPISizeWorld() * order, order), order=order, lengths=[3.0 * getMPISizeWorld(), 2.0]
)
self.badInterpolations(coupler.getSpeckley(), coupler.getRipley())
