     matrix<1, 0,   0, // matrix-shear_y_to_z 
            0, 1, 0.5,
            0, 0, 1,
            0, 0, 0>
