битовые операции

for (k = 0; k<16; ++k)                            
    {
        d = zt >> 15;
std::cout << d<<  "\t"; 
        t = xt; 
std::cout << xt<< "\t"; 
std::cout << yt<< "\t"; 
std::cout << zt<< "\n";
 
        xt -=(((yt >> k) ^ d) - d);
        yt += (((t >> k) ^ d) - d);
        zt -= ((cordic_tab[k] ^ d) - d);

int cordic_tab[] = {
    0x00003243, 0x00001DAC, 0x00000FAD, 0x000007F5, 0x000003FE, 0x000001FF, 0x000000FF, 0x0000007F, 0x0000003F, 0x0000001F, 0x0000000F, 0x00000007, 0x00000003, 0x00000001, 0x00000000, 0x00000000
};
 
int k, d, t, xt=9949, yt=0, zt=13653;
 
for (k = 0; k<16; ++k)                            
    {
        d = zt >> 15;
std::cout << d<<  "\t"; 
        t = xt; 
std::cout << xt<< "\t"; 
std::cout << yt<< "\t"; 
std::cout << zt<< "\n";
 
        xt -=(((yt >> k) ^ d) - d);
        yt += (((t >> k) ^ d) - d);
        zt -= ((cordic_tab[k] ^ d) - d);

#include <iostream>
#include <cmath>
 
#define pi 3
#define cordic_1K 9949
#define half_pi 0x00006487
#define MUL 16384.000000
#define CORDIC_NTAB 16
int cordic_tab[] = {
    0x00003243, 0x00001DAC, 0x00000FAD, 0x000007F5, 0x000003FE, 0x000001FF, 0x000000FF, 0x0000007F, 0x0000003F, 0x0000001F, 0x0000000F, 0x00000007, 0x00000003, 0x00000001, 0x00000000, 0x00000000
};
 
void cordic16_core_rotation(int *x, int *y, int *z, int n)//26DD ; 0 ;число
{
    int k, d, t,ro;
    int xt = *x, yt = *y, zt = *z; //xt=26DD, yt=0, zt=число
    for (k = 0; k<n; ++k)   
    {
        d = zt >> 15;
        t = xt; 
 
        xt -=(((yt >> k) ^ d) - d);
        yt += (((t >> k) ^ d) - d);
        zt -= ((cordic_tab[k] ^ d) - d);
    }
    *x = xt, *y = yt, *z = zt;
}
 
void cordic16_sincos(int alpha, int *s, int *c) //(число*3)/180*16384 ; указатель 1 ; указатель 2
{
    int z = alpha;                          //z = (число*3)/180*16384
    *s = 0, *c = cordic_1K;                 //d1=0  d2=26DD
    cordic16_core_rotation(c, s, &z, 16);   //26DD ; 0 ; (число*3)/180*16384 , 16 
}
 
int main()
{
    int d1, d2, angle1;
    double angle;
    std::cout << "Najmite 1 esli v gladisah 2 - v radianah\n";
    std::cin >> angle1;
    std::cout << "Vvedite ugol\n";
    std::cin >> angle;
    if (angle1 == 2) {
        cordic16_sincos(angle * MUL, &d1, &d2);         
        std::cout << "cos = " << (double)d2 / MUL << std::endl;
        std::cout << "sin = " << (double)d1 / MUL << std::endl;
    }
    if (angle1 == 1) {
        cordic16_sincos(angle*pi / 180 * MUL, &d1, &d2);      //(число*3)/180*16384 ; указатель 1 ; указатель 2
        std::cout << "cos = " << (double)d2 / MUL << std::endl;
        std::cout << "sin = " << (double)d1 / MUL << std::endl;
    }
    
    return 0;
}

for (k = 0; k<n; ++k)                     
    {
        d = zt >> 15;
std::cout << d<<  "\t"; 
        t = xt; 
std::cout << xt<< "\t"; 
std::cout << yt<< "\t"; 
std::cout << zt<< "\n";
 
        xt -=(((yt >> k) ^ d) - d);
        yt += (((t >> k) ^ d) - d);
        zt -= ((cordic_tab[k] ^ d) - d);
    }