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