14 #ifndef ESYS_LSMHEXAGBLOCKITERATOR_HPP
15 #define ESYS_LSMHEXAGBLOCKITERATOR_HPP
35 setDimRepeat(
Vec3L(2,2,2));
38 offsetMatrix(0,0,0) = 0.0;
40 offsetMatrix(0,0,2) = 0.0;
43 offsetMatrix(0,1,1) = 0.0;
44 offsetMatrix(0,1,2) = 0.0;
46 offsetMatrix(0,2,0) = 0.0;
47 offsetMatrix(0,2,1) = 0.0;
48 offsetMatrix(0,2,2) = 0.0;
50 offsetMatrix(1,0,0) = 0.0;
51 offsetMatrix(1,0,1) = 0.0;
52 offsetMatrix(1,0,2) = 0.0;
54 offsetMatrix(1,1,0) = SQRT_1_OVER_3*
getRadius();
55 offsetMatrix(1,1,1) = SQRT_1_OVER_3*
getRadius();
56 offsetMatrix(1,1,2) = SQRT_1_OVER_3*
getRadius();
58 offsetMatrix(1,2,0) = 0.0;
59 offsetMatrix(1,2,1) = 0.0;
60 offsetMatrix(1,2,2) = 0.0;
62 offsetMatrix(2,0,0) = 0.0;
63 offsetMatrix(2,0,1) = 0.0;
64 offsetMatrix(2,0,2) = 0.0;
66 offsetMatrix(2,1,0) = 0.0;
67 offsetMatrix(2,1,1) = 0.0;
68 offsetMatrix(2,1,2) = 0.0;
70 offsetMatrix(2,2,0) = 0.0;
71 offsetMatrix(2,2,1) = 0.0;
72 offsetMatrix(2,2,2) = 0.0;
74 setOffsetMatrix(offsetMatrix);