- int random = (rand() >> 10) % num_trans;
- double new_x = affine[random].a00 * current_x + affine[random].a01 * current_y + affine[random].b0 ;
- double new_y = affine[random].a10 * current_x + affine[random].a11 * current_y + affine[random].b1 ;
-
+ int random = rand( );
+ double new_x;
+ double new_y;
+ random = (((random >> 10)) & 2) + (((random >> 20) ) & 1);
+ new_x = affine[random].a00 * current_x + affine[random].a01 * current_y + affine[random].b0 ;
+ new_y = affine[random].a10 * current_x + affine[random].a11 * current_y + affine[random].b1 ;
+