#include <stdio.h>
#include <math.h>

//to create table...
//int minx[1000],miny[1000];
//for(i=0;i<1000;i++)    //fill min arrays with maximum value
//{ minx[i]=32767;
// miny[i]=32767;
//}

//for(x=0;x<32767;x++)  //cycle through all point in first octant
//  for(y=0;y<=x;y++)
//    {
//     fx=x;
//     fy=y;
//     mag=trunc(sqrt(pow(fx,2)+pow(fy,2))); //calculate the magnitude in float
//     if(mag>32699.0) //if it is bigger than 32700
//     {     //find the smalles value
//      if(x<minx[i])minx[i]=x; //this becomes coslimit
//      if(y<miny[i])miny[i]=y; //and this sinlimit
//     }
//    }

int coslimit[] = {
32700,
32700,
32698,
32697,
32694,
32692,
32688,
32685,
32681,
32676,
32671,
32665,
32659,
32652,
32645,
32637,
32629,
32621,
32611,
32602,
32592,
32581,
32570,
32558,
32546,
32533,
32520,
32510,
32496,
32481,
32466,
32451,
32435,
32419,
32402,
32385,
32367,
32348,
32329,
32310,
32290,
32270,
32256,
32235,
32213,
32191,
32169,
32146,
32122,
32098,
32073,
32048,
32023,
32000,
31982,
31955,
31928,
31900,
31872,
31843,
31814,
31784,
31754,
31739,
31708,
31677,
31645,
31613,
31580,
31547,
31513,
31488,
31465,
31430,
31395,
31359,
31322,
31286,
31248,
31232,
31197,
31159,
31120,
31081,
31041,
31001,
30976,
30949,
30908,
30866,
30824,
30781,
30738,
30720,
30685,
30640,
30596,
30551,
30505,
30464,
30452,
30405,
30358,
30311,
30263,
30214,
30208,
30161,
30112,
30062,
30011,
29960,
29952,
29908,
29856,
29804,
29751,
29698,
29696,
29646,
29593,
29538,
29483,
29440,
29433,
29377,
29321,
29265,
29208,
29184,
29159,
29101,
29043,
28984,
28928,
28928,
28878,
28818,
28757,
28696,
28672,
28652,
28590,
28528,
28466,
28416,
28416,
28360,
28297,
28233,
28168,
28160,
28129,
28064,
27998,
27932,
27904,
27896,
27829,
27762,
27695,
27648,
27648,
27594,
27526,
27457,
27392,
27392,
27358,
27288,
27218,
27147,
27136,
27121,
27050,
26979,
26906,
26880,
26880,
26813,
26740,
26666,
26624,
26624,
26576,
26501,
26427,
26368,
26368,
26340,
26264,
26188,
26112,
26112,
26105,
26028,
25951,
25873,
25856,
25856,
25794,
25716,
25637,
25600,
25600,
25562,
25483,
25403,
25344,
25344,
25333,
25253,
25172,
25090,
25088,
25088,
25026,
24944,
24861,
24832,
24832,
24802,
24720,
24636,
24576,
24576,
24576,
24499,
24415,
24330,
24320,
24320,
24283,
24198,
24112,
24064,
24064,
24064,
23986,
23899,
23812,
23808,
23808,
23779,
23692,
23604,
23552,
23552,
23552,
23490,
23402,
23313,
23296,
23296,
23295,
23206,
23170,
23170
	};

int sinlimit[] = {
0,
127,
254,
381,
508,
635,
762,
889,
1016,
1143,
1270,
1397,
1524,
1651,
1778,
1905,
2032,
2159,
2286,
2413,
2540,
2667,
2794,
2921,
3048,
3175,
3302,
3429,
3528,
3654,
3780,
3906,
4032,
4158,
4284,
4410,
4536,
4662,
4788,
4914,
5040,
5166,
5292,
5377,
5500,
5625,
5750,
5875,
6000,
6125,
6250,
6375,
6500,
6625,
6735,
6820,
6944,
7068,
7192,
7316,
7440,
7564,
7688,
7812,
7872,
7995,
8118,
8241,
8364,
8487,
8610,
8733,
8824,
8906,
9028,
9150,
9272,
9394,
9516,
9638,
9691,
9801,
9922,
10043,
10164,
10285,
10406,
10481,
10560,
10680,
10800,
10920,
11040,
11160,
11209,
11305,
11424,
11543,
11662,
11781,
11887,
11918,
12036,
12154,
12272,
12390,
12508,
12524,
12636,
12753,
12870,
12987,
13104,
13124,
13224,
13340,
13456,
13572,
13688,
13693,
13800,
13915,
14030,
14145,
14236,
14250,
14364,
14478,
14592,
14706,
14753,
14803,
14916,
15029,
15142,
15255,
15249,
15344,
15456,
15568,
15680,
15725,
15762,
15873,
15984,
16095,
16206,
16183,
16280,
16390,
16500,
16610,
16625,
16677,
16786,
16895,
17004,
17051,
17064,
17172,
17280,
17388,
17496,
17463,
17548,
17655,
17762,
17869,
17862,
17914,
18020,
18126,
18232,
18248,
18270,
18375,
18480,
18585,
18690,
18623,
18720,
18824,
18928,
19032,
18987,
19055,
19158,
19261,
19364,
19341,
19380,
19482,
19584,
19686,
19686,
19695,
19796,
19897,
19998,
20099,
20021,
20100,
20200,
20300,
20400,
20347,
20394,
20493,
20592,
20691,
20665,
20678,
20776,
20874,
20972,
21070,
20975,
21049,
21146,
21243,
21340,
21278,
21312,
21408,
21504,
21600,
21696,
21573,
21660,
21755,
21850,
21945,
21861,
21902,
21996,
22090,
22184,
22278,
22142,
22227,
22320,
22413,
22506,
22417,
22448,
22540,
22632,
22724,
22816,
22686,
22750,
22841,
22932,
23023,
22949,
22950,
23040,
23130
	};

int limiter( short *vd, short *vq)
{
 int vdt;
 int vqt;
 int quad=0;
 int i;

 //figure out which quadrant and set temp variables to Quad 0

 vdt=*vd;
 vqt=*vq;

  if (vdt<0)
 {
     quad+=2;
     vdt=-vdt;
 }

 if (vqt<0)
 {
     quad+=1;
     vqt=-vqt;
 }

 if (vdt>vqt)
 {
    if(vdt<23122) return(0);
    i=vdt>>8;
    i=vqt/i;
    if (i>255)i=256;
    if ((vdt<coslimit[i+1])||(vqt<sinlimit[i])) return (0);  //magnituide is acceptablle
    switch (quad)
    {
        case 0:
            *vd=coslimit[i+1];  //too big limit
            *vq=sinlimit[i];
            return (1);
        case 1:
            *vd=coslimit[i+1];  //too big limit
            *vq=-sinlimit[i];
            return (1);
        case 2:
            *vd=-coslimit[i+1];  //too big limit
            *vq=sinlimit[i];
            return (1);
        case 3:
            *vd=-coslimit[i+1];  //too big limit
            *vq=-sinlimit[i];
            return (1);
    }
  }
  else
  {
    if(vqt<23122) return(0);
    i=vqt>>8;
    i=vdt/i;
    if (i>255)i=256;
    if ((vqt<coslimit[i+1])||(vdt<sinlimit[i])) return (0);  //magnituide is acceptablle
    switch (quad)
     {
        case 0:
            *vq=coslimit[i+1];  //to big limit
            *vd=sinlimit[i];
            return (1);
        case 1:
            *vq=-coslimit[i+1];  //to big limit
            *vd=sinlimit[i];
            return (1);
        case 2:
            *vq=coslimit[i+1];  //to big limit
            *vd=-sinlimit[i];
            return (1);
        case 3:
            *vq=-coslimit[i+1];  //to big limit
            *vd=-sinlimit[i];
            return (1);
     }
  }
 return (-1);
 }