def testFreeRefineTrapezoidBoundaries(self): """Tests surplus based refine for Hash-Storage""" from pysgpp import GridStorage, HashGenerator from pysgpp import SurplusRefinementFunctor, HashRefinementBoundaries, DataVector s = GridStorage(2) g = HashGenerator() g.regularWithBoundaries(s, 1, True) d = DataVector(9) d[0] = 0.0 d[1] = 0.0 d[2] = 0.0 d[3] = 0.0 d[4] = 0.0 d[5] = 0.0 d[6] = 0.0 d[7] = 0.0 d[8] = 1.0 f = SurplusRefinementFunctor(d) r = HashRefinementBoundaries() r.free_refine(s, f) self.failUnlessEqual(s.size(), 21)
def testRegularTrapezoidBoundaries3D(self): """Tests 3D grid generation""" from pysgpp import GridStorage, HashGenerator s = GridStorage(3) g = HashGenerator() g.regularWithBoundaries(s, 2, True) self.failUnlessEqual(s.size(), 81)
def testRegularTrapezoidBoundaries2D(self): """Tests 2D grid generation""" from pysgpp import GridStorage, HashGenerator s = GridStorage(2) g = HashGenerator() g.regularWithBoundaries(s, 2, True) self.failUnlessEqual(s.size(), 21) s = GridStorage(2) g.regularWithBoundaries(s, 3, True) self.failUnlessEqual(s.size(), 49) from pysgpp import GridIndex i = GridIndex(2) i.set(0,1,1) i.set(1,1,1) self.failUnless(s.has_key(i)) i.set(1,2,1) self.failUnless(s.has_key(i)) i.set(0,2,3) self.failUnless(s.has_key(i)) i.set(0,3,5) self.failIf(s.has_key(i)) i.set(1,1,1) self.failUnless(s.has_key(i)) i.set(1,0,0) self.failUnless(s.has_key(i)) i.set(0,0,0) self.failUnless(s.has_key(i))
def testRegularTrapezoidBoundaries2D(self): """Tests 2D grid generation""" from pysgpp import GridStorage, HashGenerator s = GridStorage(2) g = HashGenerator() g.regularWithBoundaries(s, 2, True) self.failUnlessEqual(s.size(), 21) s = GridStorage(2) g.regularWithBoundaries(s, 3, True) self.failUnlessEqual(s.size(), 49) from pysgpp import GridIndex i = GridIndex(2) i.set(0, 1, 1) i.set(1, 1, 1) self.failUnless(s.has_key(i)) i.set(1, 2, 1) self.failUnless(s.has_key(i)) i.set(0, 2, 3) self.failUnless(s.has_key(i)) i.set(0, 3, 5) self.failIf(s.has_key(i)) i.set(1, 1, 1) self.failUnless(s.has_key(i)) i.set(1, 0, 0) self.failUnless(s.has_key(i)) i.set(0, 0, 0) self.failUnless(s.has_key(i))
def testFreeRefineTrapezoidBoundaries(self): """Tests surplus based refine for Hash-Storage""" from pysgpp import GridStorage, HashGenerator from pysgpp import SurplusRefinementFunctor, HashRefinementBoundaries, DataVector s = GridStorage(2) g = HashGenerator() g.regularWithBoundaries(s, 2, False) d = DataVector(17) for i in xrange(len(d)): d[i] = 0.0 d[12] = 1.0 f = SurplusRefinementFunctor(d) r = HashRefinementBoundaries() r.free_refine(s, f) self.failUnlessEqual(s.size(), 21)