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())