for (m=0;m<4;m++){
       for (n=0;n<4;n++){
           x1[m][n]=c1[n]*gc1[0]+c2[n]*gc2[0]+s1[n]*gs1[0]+s2[n]*gs2[0];
           x2[m][n]=c1[n]*gc1[1]+c2[n]*gc2[1]+s1[n]*gs1[1]+s2[n]*gs2[1];
       }
   }
   for (m=0;m<4;m++){
       for (n=0;n<4;n++){
           if (x1[m][n]>=4){
               x1[m][n]=x1[m][n]-4;
           }
           if (x2[m][n]>=4){
               x2[m][n]=x1[m][n]-4;
           }
       }
   }