- affine[1].a00 = 0.25 ; affine[1].a01 = 0.00 ; affine[1].a10 = 0.00 ; affine[1].a11 = 0.25 ;
- affine[1].b0 = 0.5 ; affine[1].b1 = 0.0 ;
+
+ affine[1].a00 = 1/6. ; affine[1].a01 = -1/3.*sin(FGH_PI/3.) ; affine[1].a10 = 1/3.*sin(FGH_PI/3.) ; affine[1].a11 = 1/6. ;
+ affine[1].b0 = 1/3. ; affine[1].b1 = 0.0 ;
+
+ affine[2].a00 = 1/6. ; affine[2].a01 = -1/3.*sin(-FGH_PI/3.) ; affine[2].a10 = 1/3.*sin(-FGH_PI/3.) ; affine[2].a11 = 1/6. ;
+ affine[2].b0 = 0.5 ; affine[2].b1 = sqrt(3)/6. ;
+
+ affine[3].a00 = 1/3. ; affine[3].a01 = 0.00 ; affine[3].a10 = 0.00 ; affine[3].a11 = 1/3. ;
+ affine[3].b0 = 2/3. ; affine[3].b1 = 0.0 ;