14 #ifndef ESYS_LSMCUBICBLOCKITERATOR_HPP
15 #define ESYS_LSMCUBICBLOCKITERATOR_HPP
35 setDimRepeat(
Vec3L(6,3,3));
38 offsetMatrix(0,0,0) = 0.0;
39 offsetMatrix(0,0,1) = 0.0;
41 offsetMatrix(0,0,3) = 0.0;
42 offsetMatrix(0,0,4) = 0.0;
47 offsetMatrix(0,1,2) = 0.0;
50 offsetMatrix(0,1,5) = 0.0;
52 offsetMatrix(0,2,0) = 0.0;
53 offsetMatrix(0,2,1) = 0.0;
55 offsetMatrix(0,2,3) = 0.0;
56 offsetMatrix(0,2,4) = 0.0;
61 offsetMatrix(0,3,2) = 0.0;
64 offsetMatrix(0,3,5) = 0.0;
66 offsetMatrix(0,4,0) = 0.0;
67 offsetMatrix(0,4,1) = 0.0;
69 offsetMatrix(0,4,3) = 0.0;
70 offsetMatrix(0,4,4) = 0.0;
75 offsetMatrix(0,5,2) = 0.0;
78 offsetMatrix(0,5,5) = 0.0;
80 offsetMatrix(1,0,0) = 0.0;
81 offsetMatrix(1,0,1) = 0.0;
82 offsetMatrix(1,0,2) = 0.0;
84 offsetMatrix(1,1,0) = 2.0*SQRT_1_OVER_3*
getRadius();
85 offsetMatrix(1,1,1) = 2.0*SQRT_1_OVER_3*
getRadius();
86 offsetMatrix(1,1,2) = 2.0*SQRT_1_OVER_3*
getRadius();
88 offsetMatrix(1,2,0) = SQRT_1_OVER_3*
getRadius();
89 offsetMatrix(1,2,1) = SQRT_1_OVER_3*
getRadius();
90 offsetMatrix(1,2,2) = SQRT_1_OVER_3*
getRadius();
92 offsetMatrix(2,0,0) = 0.0;
93 offsetMatrix(2,0,1) = 0.0;
94 offsetMatrix(2,0,2) = 0.0;
96 offsetMatrix(2,1,0) = 0.0;
97 offsetMatrix(2,1,1) = 0.0;
98 offsetMatrix(2,1,2) = 0.0;
100 offsetMatrix(2,2,0) = 0.0;
101 offsetMatrix(2,2,1) = 0.0;
102 offsetMatrix(2,2,2) = 0.0;
104 setOffsetMatrix(offsetMatrix);