+
+
+
+template<int h>
+void encrypt_ecb_rm(uchar* seq_in, uchar *seq_out, int len,uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, uchar *IV, int debug) {
+
+ uchar X[h];
+ uchar Y[h];
+ uchar fX[h];
+ uchar gY[h];
+ uchar IV1[h];
+ uchar IV2[h];
+ uchar *RM1;
+ uchar *RM2;
+ uchar tmp[h];
+
+
+
+ for(int a=0;a<h;a+=4) {
+ IV1[a]=IV[a];
+ IV1[a+1]=IV[a+1];
+ IV1[a+2]=IV[a+2];
+ IV1[a+3]=IV[a+3];
+ }
+
+
+ for(int a=0;a<h;a+=4) {
+ IV2[a]=IV[h+a];
+ IV2[a+1]=IV[h+a+1];
+ IV2[a+2]=IV[h+a+2];
+ IV2[a+3]=IV[h+a+3];
+
+ }
+
+
+
+
+
+
+ for(int it=0;it<len/2;it++) {
+ int ind1=Pbox[it]*h;
+ int ind2=Pbox[it+len/2]*h;
+
+
+ RM1=&RM[PboxSRM[it]*h];
+ RM2=&RM[h*h+PboxSRM[len/2-it]*h];
+
+
+
+ for(int a=0;a<h;a+=4) {
+ X[a]=seq_in[ind2+a];
+ X[a+1]=seq_in[ind2+a+1];
+ X[a+2]=seq_in[ind2+a+2];
+ X[a+3]=seq_in[ind2+a+3];
+ }
+
+ for(int a=0;a<h;a+=4) {
+ Y[a]=seq_in[ind1+a];
+ Y[a+1]=seq_in[ind1+a+1];
+ Y[a+2]=seq_in[ind1+a+2];
+ Y[a+3]=seq_in[ind1+a+3];
+ }
+
+
+ for(int a=0;a<h;a+=4) {
+ tmp[a]=Sbox1[X[a]^RM1[a]];
+ tmp[a+1]=Sbox1[X[a+1]^RM1[a+1]];
+ tmp[a+2]=Sbox1[X[a+2]^RM1[a+2]];
+ tmp[a+3]=Sbox1[X[a+3]^RM1[a+3]];
+ }
+
+ for(int a=0;a<h;a+=4) {
+ fX[a]=Sbox2[tmp[a]^Y[a]];
+ fX[a+1]=Sbox2[tmp[a+1]^Y[a+1]];
+ fX[a+2]=Sbox2[tmp[a+2]^Y[a+2]];
+ fX[a+3]=Sbox2[tmp[a+3]^Y[a+3]];
+ }
+
+
+ /*for(int a=0;a<h;a+=4) {
+ fX[a]=Sbox2[Sbox1[X[a]^RM1[a]]^Y[a]];
+ fX[a+1]=Sbox2[Sbox1[X[a+1]^RM1[a+1]]^Y[a+1]];
+ fX[a+2]=Sbox2[Sbox1[X[a+2]^RM1[a+2]]^Y[a+2]];
+ fX[a+3]=Sbox2[Sbox1[X[a+3]^RM1[a+3]]^Y[a+3]];
+ }
+ */
+
+
+
+ for(int a=0;a<h;a+=4) {
+ tmp[a]=Sbox2[fX[a]^Y[a]];
+ tmp[a+1]=Sbox2[fX[a+1]^Y[a+1]];
+ tmp[a+2]=Sbox2[fX[a+2]^Y[a+2]];
+ tmp[a+3]=Sbox2[fX[a+3]^Y[a+3]];
+
+ }
+ for(int a=0;a<h;a+=4) {
+ gY[a]=Sbox1[tmp[a]^RM2[a]];
+ gY[a+1]=Sbox1[tmp[a+1]^RM2[a+1]];
+ gY[a+2]=Sbox1[tmp[a+2]^RM2[a+2]];
+ gY[a+3]=Sbox1[tmp[a+3]^RM2[a+3]];
+
+ }
+
+
+ /* for(int a=0;a<h;a+=4) {
+ gY[a]=Sbox1[Sbox2[fX[a]^Y[a]]^RM2[a]];
+ gY[a+1]=Sbox1[Sbox2[fX[a+1]^Y[a+1]]^RM2[a+1]];
+ gY[a+2]=Sbox1[Sbox2[fX[a+2]^Y[a+2]]^RM2[a+2]];
+ gY[a+3]=Sbox1[Sbox2[fX[a+3]^Y[a+3]]^RM2[a+3]];
+
+ } */
+
+
+
+
+
+
+ for(int a=0;a<h;a+=4) {
+ seq_out[ind2+a]=gY[a];
+ seq_out[ind2+a+1]=gY[a+1];
+ seq_out[ind2+a+2]=gY[a+2];
+ seq_out[ind2+a+3]=gY[a+3];
+ }
+
+ for(int a=0;a<h;a+=4) {
+ seq_out[ind1+a]=fX[a];
+ seq_out[ind1+a+1]=fX[a+1];
+ seq_out[ind1+a+2]=fX[a+2];
+ seq_out[ind1+a+3]=fX[a+3];
+ }
+
+
+
+ }
+
+
+
+
+}
+
+
+
+
+
+
+
+
+template<int h>
+void decrypt_ecb_rm(uchar* seq_in, uchar *seq_out, int len, uchar* RM, int *Pbox, int *PboxSRM, uchar *Sbox1, uchar *Sbox2, uchar *Inv_Sbox1, uchar *Inv_Sbox2, uchar *IV, int debug) {
+
+ uchar invfX[h];
+ uchar invgY[h];
+ uchar fX[h];
+ uchar gY[h];
+ uchar IV1[h];
+ uchar IV2[h];
+ uchar *RM1;
+ uchar *RM2;
+ uchar tmp[h];
+
+
+
+ for(int a=0;a<h;a+=4) {
+ IV1[a]=IV[a];
+ IV1[a+1]=IV[a+1];
+ IV1[a+2]=IV[a+2];
+ IV1[a+3]=IV[a+3];
+ }
+
+
+ for(int a=0;a<h;a+=4) {
+ IV2[a]=IV[h+a];
+ IV2[a+1]=IV[h+a+1];
+ IV2[a+2]=IV[h+a+2];
+ IV2[a+3]=IV[h+a+3];
+
+ }
+
+
+
+ for(int it=0;it<len/2;it++) {
+ int ind1=Pbox[it]*h;
+ int ind2=Pbox[it+len/2]*h;
+
+
+ RM1=&RM[PboxSRM[it]*h];
+ RM2=&RM[h*h+PboxSRM[len/2-it]*h];
+
+
+
+ for(int a=0;a<h;a+=4) {
+ gY[a]=seq_in[ind2+a];
+ gY[a+1]=seq_in[ind2+a+1];
+ gY[a+2]=seq_in[ind2+a+2];
+ gY[a+3]=seq_in[ind2+a+3];
+ }
+
+ for(int a=0;a<h;a+=4) {
+ fX[a]=seq_in[ind1+a];
+ fX[a+1]=seq_in[ind1+a+1];
+ fX[a+2]=seq_in[ind1+a+2];
+ fX[a+3]=seq_in[ind1+a+3];
+ }
+
+
+
+
+ for(int a=0;a<h;a+=4) {
+ tmp[a]=Inv_Sbox1[gY[a]]^RM2[a];
+ tmp[a+1]=Inv_Sbox1[gY[a+1]]^RM2[a+1];
+ tmp[a+2]=Inv_Sbox1[gY[a+2]]^RM2[a+2];
+ tmp[a+3]=Inv_Sbox1[gY[a+3]]^RM2[a+3];
+ }
+
+
+ for(int a=0;a<h;a+=4) {
+ invgY[a]=Inv_Sbox2[tmp[a]]^fX[a];
+ invgY[a+1]=Inv_Sbox2[tmp[a+1]]^fX[a+1];
+ invgY[a+2]=Inv_Sbox2[tmp[a+2]]^fX[a+2];
+ invgY[a+3]=Inv_Sbox2[tmp[a+3]]^fX[a+3];
+ }
+
+
+
+ /* for(int a=0;a<h;a+=4) {
+ invgY[a]=Inv_Sbox2[Inv_Sbox1[gY[a]]^RM2[a]]^fX[a];
+ invgY[a+1]=Inv_Sbox2[Inv_Sbox1[gY[a+1]]^RM2[a+1]]^fX[a+1];
+ invgY[a+2]=Inv_Sbox2[Inv_Sbox1[gY[a+2]]^RM2[a+2]]^fX[a+2];
+ invgY[a+3]=Inv_Sbox2[Inv_Sbox1[gY[a+3]]^RM2[a+3]]^fX[a+3];
+ } */
+
+
+
+
+ for(int a=0;a<h;a+=4) {
+ tmp[a]=Inv_Sbox2[fX[a]]^invgY[a];
+ tmp[a+1]=Inv_Sbox2[fX[a+1]]^invgY[a+1];
+ tmp[a+2]=Inv_Sbox2[fX[a+2]]^invgY[a+2];
+ tmp[a+3]=Inv_Sbox2[fX[a+3]]^invgY[a+3];
+
+ }
+
+ for(int a=0;a<h;a+=4) {
+ invfX[a]=Inv_Sbox1[tmp[a]]^RM1[a];
+ invfX[a+1]=Inv_Sbox1[tmp[a+1]]^RM1[a+1];
+ invfX[a+2]=Inv_Sbox1[tmp[a+2]]^RM1[a+2];
+ invfX[a+3]=Inv_Sbox1[tmp[a+3]]^RM1[a+3];
+
+ }
+
+
+ /*
+ for(int a=0;a<h;a+=4) {
+ invfX[a]=Inv_Sbox1[Inv_Sbox2[fX[a]]^invgY[a]]^RM1[a];
+ invfX[a+1]=Inv_Sbox1[Inv_Sbox2[fX[a+1]]^invgY[a+1]]^RM1[a+1];
+ invfX[a+2]=Inv_Sbox1[Inv_Sbox2[fX[a+2]]^invgY[a+2]]^RM1[a+2];
+ invfX[a+3]=Inv_Sbox1[Inv_Sbox2[fX[a+3]]^invgY[a+3]]^RM1[a+3];
+
+ }
+ */
+
+ for(int a=0;a<h;a+=4) {
+ seq_out[ind2+a]=invfX[a];
+ seq_out[ind2+a+1]=invfX[a+1];
+ seq_out[ind2+a+2]=invfX[a+2];
+ seq_out[ind2+a+3]=invfX[a+3];
+ }
+
+ for(int a=0;a<h;a+=4) {
+ seq_out[ind1+a]=invgY[a];
+ seq_out[ind1+a+1]=invgY[a+1];
+ seq_out[ind1+a+2]=invgY[a+2];
+ seq_out[ind1+a+3]=invgY[a+3];
+ }
+
+
+
+ }
+
+
+
+
+}
+
+
+
