int main(void)
{
int k=0;
double* sln;// A pointer for receiving the ballistic solution.
double bc=0.5; // The ballistic coefficient for the projectile.
double v=1200; // Intial velocity, in ft/s
double sh=1.6; // The Sight height over bore, in inches.
double angle=0; // The shooting angle (uphill / downhill), in degrees.
double zero=100; // The zero range of the rifle, in yards.
double windspeed=0; // The wind speed in miles per hour.
double windangle=0; // The wind angle (0=headwind, 90=right to left, 180=tailwind, 270/-90=left to right)
// If we wish to use the weather correction features, we need to
// Correct the BC for any weather conditions. If we want standard conditions,
// then we can just leave this commented out.
// bc=AtmCorrect(bc, 0, 29.59, 100,.7);
// First find the angle of the bore relative to the sighting system.
// We call this the "zero angle", since it is the angle required to
// achieve a zero at a particular yardage. This value isn't very useful
// to us, but is required for making a full ballistic solution.
// It is left here to allow for zero-ing at altitudes (bc) different from the
// final solution, or to allow for zero's other than 0" (ex: 3" high at 100 yds)
double zeroangle=0;
zeroangle=ZeroAngle(G1,bc,v,1.6,zero,0);
// Now we have everything needed to generate a full solution.
// So we do. The solution is stored in the pointer "sln" passed as the last argument.
// k has the number of yards the solution is valid for, also the number of rows in the solution.
k=SolveAll(G1,bc,v,sh,angle,zeroangle,windspeed,windangle,&sln);
// Now print a simple chart of X / Y trajectory spaced at 10yd increments
int s=0;
for (s=0;s<=100;s++)
{
printf("\nX: %.0f Y: %.2f",GetRange(sln,s*10), GetPath(sln,s*10));
}
return 0;
}
int trajectory (void)
{
char buffer[80];
double dc,pw,iv,zr,sh,sa,wv,wa,at,tp,bp,rh,za,tv,cl=0.0,vcl,hcl,vd;
double vmoa,vmils,hmoa,hmils;
int md,mr,df,yi,x;
double* sln=NULL;
int spbr[5];
int scope;
printf("\n\n\t * Trajectory Calculator * \n");
while(1)
{
printf("\nEnter Ballistic Coefficient : ");
sgets(buffer,80);
dc = atof(buffer);
if (dc > 0 && dc < 2) break;
printf("Ballistic Coefficient must > 0 and < 2!\n");
}
while(1)
{
printf("\n");
printf("G0 Roundball\n");
printf("G1 Ingalls (by far the most popular)\n");
printf("G2 (Aberdeen J projectile)\n");
printf("G5 (short 7.5deg boat-tail, 6.19 calibers long tangent ogive)\n");
printf("G6 (flatbase, 6 calibers long secant ogive)\n");
printf("G7 (long 7.5deg boat-tail, 10 calibers tangent ogive)\n");
printf("G8 (flatbase, 10 calibers long secant ogive)\n");
printf("\nEnter Drag Model 0,1,2,5,6,7 or 8 : ");
sgets(buffer,80);
df = atoi(buffer);
if (df == 0) break;
if (df == 1) break;
if (df == 2) break;
if (df == 5) break;
if (df == 6) break;
if (df == 7) break;
if (df == 8) break;
printf("Drag Model must be 0,1,2,5,6,7 or 8!\n");
}
while(1)
{
printf("\nEnter Bullet Weight (gr) : ");
sgets(buffer,80);
pw = atof(buffer);
if (pw >= 1) break;
printf("Bullet Weight must be >= 1 grain!\n");
}
while(1)
{
printf("\nEnter Muzzle Velocity (ft/s) : ");
sgets(buffer,80);
iv = atof(buffer);
if (iv > 0) break;
printf("Muzzle Velocity must be > 0!\n");
}
while(1)
{
printf("\nEnter Zero Range (yds) : ");
sgets(buffer,80);
zr = atof(buffer);
if (zr >= 1) break;
printf("Zero Range must be >= 1 yard!\n");
}
while(1)
{
printf("\nEnter Sight Height Above Bore (in) : ");
sgets(buffer,80);
sh = atof(buffer);
if (sh > 0) break;
printf("Sight Height must be > 0!\n");
}
while(1)
{
printf("\nEnter Shooting Angle (deg) ");
printf("\n\t+90\n");
printf("\t 0 \n");
printf("\t-90\n :");
sgets(buffer,80);
sa = atof(buffer);
if (sa >= -90 && sa <= 90) break;
printf("Shooting Angle must be between -90 and +90!\n");
}
while(1)
{
printf("\nEnter Wind Velocity (mph) : ");
sgets(buffer,80);
wv = atof(buffer);
if (wv >= 0) break;
printf("Wind Velocity cannot be < 0!\n");
}
while(1)
{
printf("\nEnter Wind Angle is Coming From (0-360 deg)");
printf("\n\t 0\n");
printf("270\t\t90\n");
printf("\t180\n :");
sgets(buffer,80);
wa = atof(buffer);
if (wa >= 0 && wa <= 360) break;
printf("Wind Angle must be between 0 and 360 deg!\n");
}
printf("\nScope (Y/N) : ");
sgets(buffer,80);
scope = 0;
if (buffer[0] == 'Y' || buffer[0] == 'y')
{
while(1)
{
printf("\nEnter # Clicks to Move Scope 1 inch @ 100 yrds : ");
sgets(buffer,80);
cl = atof(buffer);
if (cl >= 0) break;
printf("# Clicks must not be < 0!\n");
}
scope = 1;
}
while(1)
{
printf("\nEnter Vitals Diameter (in): ");
sgets(buffer,80);
vd = atof(buffer);
if (vd >= 1) break;
printf("Vitals Diameter must be at least 1 inch!\n");
}
printf("\nCorrect for Atmospheric Conditions (Y/N) : ");
sgets(buffer,80);
if (buffer[0] == 'Y' || buffer[0] == 'y')
{
while(1)
{
printf("\nEnter Altitude (ft) : ");
sgets(buffer,80);
at = atof(buffer);
if (at >= 0) break;
printf("Altitude must be >= 0!\n");
}
while(1)
{
printf("\nEnter Temperature (F) : ");
sgets(buffer,80);
tp = atof(buffer);
if (tp > -100) break;
printf("Temperature must be > -100 deg!\n");
}
while(1)
{
printf("\nEnter Barometric Pressure (in Hg) : ");
sgets(buffer,80);
bp = atof(buffer);
if (bp > 20 && bp < 35) break;
printf("Barometric Pressure must be > 20 and < 35!\n");
}
while(1)
{
printf("\nEnter Relative Humidity (%%) : ");
sgets(buffer,80);
rh = atof(buffer)/100;
if (rh >= 0 && rh <= 100) break;
printf("Relative Humidity must be between 0 and 100%%!\n");
}
dc = AtmCorrect(dc, at, bp, tp, rh);
}
za = ZeroAngle(df, dc, iv, sh, zr, 0.0);
md = SolveAll(df, dc, iv, sh, sa, za, wv, wa, &sln);
while(1)
{
printf("\nEnter Maximum Range (yds) up to %d yds for Table : ",md);
sgets(buffer,80);
mr = atoi(buffer);
if (mr > 0 && mr < md) break;
printf("Maximum Table Range must be > 0 and < %d yds!\n",md);
}
while(1)
{
printf("\nEnter Yardage Increment (yds) for Table : ");
sgets(buffer,80);
yi = atoi(buffer);
if (yi >= 1 && yi < md) break;
printf("Yardage Increment must be between 1 and %d yds!\n",md);
}
pbr(df, dc, iv, sh, vd, zr, spbr);
printf("\n\n\t * Ballistic Information * \n");
printf("\nFinal Ballistic Coefficient: %3.3f\n", dc);
printf("Ballistic Model: G%d\n",df);
printf("Muzzle Velocity: %3.0f fps\n",iv);
printf("Muzzle Energy: %3.0f ft/lb\n",(pw * iv * iv)/450436);
printf("Absolute Maximum Range: %d yds\n",md);
printf("Zero Range: %3.0f yds\n",zr);
printf("Shooting Angle: %3.1f deg\n",sa);
if (spbr[4]/100.0 >= 0)
printf("\n\nPBR Sight in for %3.0f yds: %3.1f in high\n",zr,fabs(spbr[4]/100.0));
else
printf("\n\nPBR Sight in for %3.0f yds: %3.1f in low\n",zr,fabs(spbr[4]/100.0));
printf("\tPBR Vital Area Diameter: %3.1f in\n",vd);
printf("\tPBR Near Zero: %d yds\n",spbr[0]);
printf("\tPBR Far Zero: %d yds\n", spbr[1]);
printf("\tMinimum PBR: %d yds\n",spbr[2]);
printf("\tMaximum PBR: %d yds\n",spbr[3]);
printf("\n\nWind Speed %3.1f mph from %3.1f deg\n",wv,wa);
printf("\n\n\t * Flight Characteristics *\n");
printf("\n\nDIST\tBD\tWind\tV\tE\tT\n");
printf("(yds)\t(in)\t(in)\t(ft/s)\t(ft/lb)\t(sec)\n");
for (x=yi;x<=mr;x+=yi)
{
printf("%d\t",x);
printf("%3.1f\t",GetPath(sln,x));
printf("%3.1f\t",GetWindage(sln,x));
tv = GetVelocity(sln,x);
printf("%3.0f\t",tv);
printf("%3.0f\t",(pw * tv * tv)/450436);
printf("%3.2f\n",GetTime(sln,x));
}
printf("\n\n\n");
if (scope)
{
printf("\t * Scope Adjustments *\n");
printf("\n\nDIST\tVMOA\tVCLKS\tVMILS\tHMOA\tHCLKS\tHMILS\n");
for (x=yi;x<=mr;x+=yi)
{
vcl=GetMOA(sln,x)*cl;
hcl=GetWindageMOA(sln,x)*cl;
vmoa = GetMOA(sln,x);
vmils = GetMOA(sln,x)*0.2908;
hmoa = GetWindageMOA(sln,x);
hmils = GetWindageMOA(sln,x)*0.2908;
printf("%d\t",x);
if (vmoa > 0)
printf("U%3.1f\t",vmoa);
else
printf("D%3.1f\t",fabs(vmoa));
if (vcl > 0)
printf("U%3.0f\t",vcl);
else
printf("D%3.0f\t",fabs(vcl));
if (vmils > 0)
printf("U%3.1f\t",vmils);
else
printf("D%3.1f\t",fabs(vmils));
if (hmoa > 0)
printf("R%3.1f\t",hmoa);
else
printf("L%3.1f\t",fabs(hmoa));
if (hcl > 0)
printf("R%3.0f\t",hcl);
else
printf("L%3.0f\t",fabs(hcl));
if (hmils > 0)
printf("R%3.1f\t",hmils);
else
printf("L%3.1f\t",fabs(hmils));
printf("\n");
}
}
return 0;
}
static void TweakMovement (usercmd_t *cmd)
{
vec3_t fw, rt;
vec3_t movevec;
vec3_t v1, ang;
float a;
float delta;
vec3_t intentions;
vec3_t intentions_a;
vec3_t cmd_ang;
qbool onground;
static qbool locked = false;
static double startlock = -1;
static float olddelta;
static qbool olddelta_valid = false;
static short oldmove[2];
if (!cl_easymove.value || cl.waterlevel)
return;
if (!cmd->forwardmove && !cmd->sidemove)
{
oldmove[0] = oldmove[1] = 0;
return;
}
onground = (cl.onground && !((cmd->buttons & BUTTON_JUMP) && !em_jump_held));
VectorCopy (cmd->angles, cmd_ang);
cmd_ang[PITCH] = 0;
cmd_ang[ROLL] = 0;
AngleVectors (cmd_ang, fw, rt, NULL);
VectorScale (fw, cmd->forwardmove, intentions);
VectorMA (intentions, cmd->sidemove, rt, intentions);
vectoangles (intentions, intentions_a);
VectorCopy (cl.simvel, v1);
v1[2] = 0;
if (VectorLength(v1) < 30)
return;
vectoangles (v1, ang);
delta = AngleMod(intentions_a[YAW] - ang[YAW]);
if (!locked && startlock == -1)
startlock = curtime;
if ((cmd->forwardmove != oldmove[0] || cmd->sidemove != oldmove[1])
&& (!onground || (oldmove[0] || oldmove[1]))) {
locked = false;
olddelta_valid = false;
startlock = curtime;
}
oldmove[0] = cmd->forwardmove;
oldmove[1] = cmd->sidemove;
if (fabs(delta) < 0.5 || (olddelta_valid && (
(olddelta >= 0 && delta <= 0) || (olddelta <= 0 && delta >= 0)
|| fabs(delta) > fabs(olddelta) + 0.1)) ) {
locked = true;
}
olddelta = delta;
olddelta_valid = true;
if (curtime - startlock > (onground ? 0.1 : 0.3))
locked = true;
if (fabs(delta) > 95) {
locked = true;
return;
}
if (!locked || em_touch_wall)
return;
if (onground)
{
if (fabs(delta) < 0.5)
return;
a = BestAngle (v1, intentions, onground);
}
else
{
a = intentions_a[YAW] + (delta > 0 ? 1 : -1)*ZeroAngle(VectorLength(v1));
}
VectorClear (ang);
ang[YAW] = a;
AngleVectors (ang, movevec, NULL, NULL);
cmd->forwardmove = DotProduct (movevec, fw) * 500;
cmd->sidemove = DotProduct (movevec, rt) * 500;
}
