create16bit := proc(C,C2)

read(create8bit);
create8bit(B);

for i to 16 do
   for j to 16 do
      A[i][2*j-1] := B[(i-1)*16 + j];
      A[i][2*j]   := B[(i-1)*16 + j];
   od:
od:

read(shiftm2);


for i to 16 do
for j to 16 do

   shiftm2(A,i,16-(14-i+1+j mod 16),1-j mod 8,8,C,1+8*(j+(i-1)*16-1)*32,sign(1-j));
   shiftm2(a,i,16-(14-i+1+j mod 16),1-j mod 8,8,C2,1+8*(j+(i-1)*16-1)*32,sign(1-j));

od:
od:

output := fopen("apgc16bit",WRITE,TEXT);

for i to 2^11 do

#print("list:",i);
#mystring:=cat("list: ",i);
#print(mystring);
#writeline(output,"Hi");

for j to 32 do
#print(j,C2[(i-1)*32+j],C[(i-1)*32+j]);
#mystring:=cat("item[",i,",",j,"]=",C2[(i-1)*32+j],"=",C[(i-1)*32+j]);
#print(mystring);
writeline(output,C[(i-1)*32+j]);
od:

#j:=1:
#print(j,C2[(i-1)*32+j],C[(i-1)*32+j]);
#j:=32:
#print(j,C2[(i-1)*32+j],C[(i-1)*32+j]);

od:

close(output);

end;