/* ===========================================================;*/
/* == MultiCelestron 999 :  ==================================;*/
/* == as multi points movement symlation Ver.1.02 ============;*/
/* == Logic spreading Version as test ========================;*/
/* == Extended From comet3xx Basic Source File ===============;*/
/* == 1999 Feb 13th- 1999 Mar 10th. ==========================;*/
/* == Copy Right:  Kazuki Kimijima/Kaztomi Miyama ============;*/
/* == Address VYE00525@Nifty.ne.jp & =========================;*/
/* == URL:http://home.interlink.or.jp/~kimijima/ =============;*/
/* == First, For Borland/Imprise TurboC++ for DOS 4.00 =======;*/
/* ===========================================================;*/

#include<graphics.h>;
#include<stdlib.h>;
#include<stdio.h>;
#include<conio.h>;
#include<math.h>;

/* == Functions for Using ====================================;*/
void pixelcircle(double x,double y,int flag,double frame,int flag98)
{
int xmoni,ymoni;
double dotperauscale;
int dotx,doty;
/*int gd=DETECT,gm;
initgraph(&gd,&gm,"")     ;*/
/*This is for VGA*/
xmoni=640;
if(flag98==98){ymoni=400;}else{ymoni=480;}
dotperauscale=(xmoni/2)/frame;
     /*Default, for "frame" AU per Half Xmonitor Scale.*/

dotx=int(xmoni/2+dotperauscale*x);
if (dotx<0 && xmoni<dotx){goto nodrawforover;}
doty=int(ymoni/2-dotperauscale*y);
if (doty<0 && ymoni<doty){goto nodrawforover;}

if (flag==1){putpixel(dotx,doty,4);}
if (flag==2){circle(dotx,doty,2);}
nodrawforover:;
}
/*Long Mathmatical Function Declear*/
double actualrandom01(void)
{return((double) rand()/32767.0);}
double cosrandom(void){return(cos(2*3.141592*actualrandom01()));}
double sinrandom(void){return(sin(2*3.141592*actualrandom01()));}


/*Main Routine*/
void main(int argc, char *argv[])
{
/* ===================================================================== ;*/
/*
This Main Program's Logic is next Chart
1)Parameters Decleare(4 Blocks)
2)Input Some Parameters(5 Blocks)
3)Graphics Parameters Tuning(AUTO)
4)All Celestron Random Settings
5)Main Core Calculation
	5-1)Select one Point Celestron and to go Next Celestron
	5-2)Calculation of Gravity-Influence from All element Celestron for
	the (5-1's)one Celestron .
	5-3)Real Moving Operation of 5-1's Celestron.
	5-4)Goto 5-1's step and until Timelimit of inputted.
6)Graphic Subroutine (From Old Basic-Language's fossel)
7)END decleare;*/
/* ===================================================================== ;*/

/* ===================================================================== ;*/
/* 1)Parameters Decleare(4 Blocks) ;*/
/* ===================================================================== ;*/

/* == 1/4 Astro-Constant parameters ==================================== ;*/
double earthspeed;
	earthspeed = 2*3.141592/365.251;
	/*printf("%lf",earthspeed   ,"\n");*/
	/*Unit   is 1 AU;*/

double gconst,sungravity;
	 gconst=5.2/(pow(10,8));
	 sungravity=3.3*pow(10,5)*gconst*earthspeed;

/* 5.2:the local number in this machines.envilonment.;*/
/* 3.3*10^5 times of Earth, Sun Mass.;*/

/* == 2/4 Main-Using parameters ======================================= ;*/
/*V:Vector of Speed;*/
/*A:Accelation;*/
double x[1001],vx[1001],ax,sumax;
double y[1001],vy[1001],ay,sumay;
/*double z[1001],vz[1001],az,sumaz  (2-Demention Version,No-Use.) ;*/
double dx,dy;
/*temprary number.*/;


/*Time and Gravity Parameters;*/
double dtd     ;/*dtd:delta time days::calcuration of micro time unit;*/

double jsumdtd;
double pereachgrvmassave, eachgrvmassave,speedpara  ;
/*-ave:Average;*/

/* ==3/4 Help Prameters =============================================== ;*/
double howlong,timecounter;
double  n;
int random,i,j;
double frame;
int flag98;
/* n:preparation of n x-y points (Max is 50 points).                     ;*/
/* random :for random parameters of location and speed making,sub numbers.  ;*/
/* i :for culculation help parameter,to select one point and total accel-;*/
/* influation calculate,                                                 ;*/
/* j :for culculation help parameter,for one point,search other gravity- ;*/
/* points next to next.                                                  ;*/
/* flag98 is 1 for 640*400 machine                                       ;*/

double d ;/*(Roof length of Dynamic triangle X-Y's Diagram)'s ^2;*/

double accel,sumaccel,accelcheck1 ,accelcheck2;
/* accelcheck1-2:1 is one by one Over accelation Check.(No-Using);
		:2 is Total SUM accelation Checking. ;*/

/* == 4/4 Borland/Inprize Graphic Driver Locatiuon ===================== ;*/
int gd=DETECT,gm;
initgraph(&gd,&gm,"")     ;

/* ===================================================================== ;*/
/*  2)Input Some Parameters(5 Blocks)  ;*/
/* ===================================================================== ;*/

if (argc==1){
printf("Multi -Celestrons Moving Simulation 999 Version 1.05 \n");
printf("CopyRight 1999 Kazuki Kimijima Fujiwara/Kaztomi Miyama\n");
printf("VYE00525@Nifty.ne.jp/81966682@People.or.jp\n");
printf("URL:http://home.interlink.or.jp/~kimijima/\n");
printf("\nUsage:mc n eachgrvave speed howlong accelcheck frame flag98\n");
goto processend;}

if (argc <5){printf("Input Parameter Error!!\n");goto processend;}
if (argc <6){accelcheck2=0.001;}
if (argc <7){frame=5;}
if (argc <8){flag98=0;}

if (8<=argc){flag98     =atoi(argv[7]);}

if (7<=argc){frame      =atof(argv[6]);}
if (6<=argc){accelcheck2=atof(argv[5]);}

howlong                 =atof(argv[4]);
speedpara               =atof(argv[3]);
pereachgrvmassave       =atof(argv[2]);
n                       =atof(argv[1]);

/* 1/5 Prepare Celestron Points Number.;*/
printf ("Points Number.\n");
printf ("%lf\n",n);
/* Already 50 Points Preparated ;*/

/* 2/5 Prepare Celestron Average Gravity ;*/
printf ("How Mass of each points \n? of Sun.\n");
printf ("%lf\n" ,pereachgrvmassave);
eachgrvmassave = sungravity * pereachgrvmassave;

/* 3/5 Prepare Average First Speed for each Celestrons for Rumdom ;*/
printf ("Average of First Speed \nVector bigth \n:Unit Earthspeed around SUN.\n");
printf ("%lf\n",speedpara);
speedpara = speedpara * 2;

/* 4/5 Input Howlong times on Calculation (Unit : Years) ;*/
printf("HOWLONG by years.      \n");
printf ("%lf\n",howlong);
howlong = howlong * 365.251;

/* 5/5 Over Accelation Check for Uncorrectness (Only Total parameter ;*/
/* Accelcheck2                                                   ;*/
printf ("Correctionness \nOver Accel Check\n:Unit Eeath Lotation speed.\n");
printf ("Around SUN :Usually 0.001 \nEeath Speed. \n");
/*;*/
printf ("%lf\n",accelcheck2);
/* ==================================================================== ;*/
/* 3)Guide Circle Drawing     ;*/
/* ==================================================================== ;*/
if(flag98==98){circle(320,200,int(320*3/frame+.5));
		circle(320,200,int(320*5/frame+.5));}
else{circle (320, 240,int(320*3/frame+.5));
circle (320, 240,int(320*5/frame+.5));}

printf ("3 AU radius \n");
printf ("5 AU radius \n");

printf("Hit Any KEY:Then,start.\n");
getchar();
/* ==================================================================== ;*/
/* 4)All Celestron Ramdom Setting ;*/
/* ==================================================================== ;*/
for(random=1 ;random <=n;random=random+1)
{

	/*radius : 3 AU;*/

	x[random] = 3 * sqrt(actualrandom01()) * cosrandom();
	y[random] = 3 * sqrt(actualrandom01()) * sinrandom();

	vx[random] = earthspeed * speedpara * actualrandom01()* cosrandom();
	vy[random] = earthspeed * speedpara * actualrandom01()* sinrandom();

/*printf ("%lf %lf %lf %lf",x[random], y[random], vx[random], vy[random],"\n");
*/
pixelcircle(x[random],y[random],2,frame,flag98);
/* ==old NEXT random of BASIC === ;*/
}

/* ==================================================================== ;*/
/*5)Main Core Calculation
	5-1)Select one Point Celestron and to go Next Celestron.
	5-2)Calculation of Gravity-Influence from All element Celestron for
	the (5-1's)one Celestron .
	5-3)Real Moving Operation of 5-1's Celestron.
	5-4)Goto 5-1's step and until Timelimit of inputted.*/
/* ==================================================================== ;*/

jsumdtd=0;
dtd=.25;
for(timecounter=0;timecounter<=howlong;timecounter=timecounter+jsumdtd)
{
printf("Now Time %lf Years.\r",timecounter/365.251);
/*5-1*/
	for (i = 1;i<=n;i=i+1)
	{
		/*dtd=deltatimedays(first,1/4days.);*/
		jsumdtd=0;
		dtd=.25;

		/* 5-2*/
		secondroop:;

		sumax = 0;
		sumay = 0;

		/*printf ("\f","%lf %lf",timecounter,dtd,"\n");*/

		for (j = 1; j<=n;j=j+1)
		{
			/*if i==j, it's me to me calculation,
			then no calculation.;*/

			if (i != j)
			{

			dx=(x[i]-x[j]);
			dy=(y[i]-y[j]);

			d = pow(dx,2) + pow(dy,2);

			ax =  -dx/ (pow(d,1.5));
			ay =  -dy/ (pow(d,1.5));

			ax = eachgrvmassave * earthspeed * ax ;
			ay = eachgrvmassave * earthspeed * ay ;

			/*accel = sqrt(ax*ax + ay *ay);*/

			/* accelcheck1 = .1
			if (accelcheck1 < accel )
			{ dtd = dtd / 2goto secondroop;};*/

			sumax = sumax + ax;
			sumay = sumay + ay;

			sumaccel = sqrt(sumax*sumax + sumay *sumay);

			/* accelcheck2 = .1 tune this parameter--
			by manual-input;*/

			if (accelcheck2 < sumaccel*dtd ){
			dtd = dtd / 2;
			goto secondroop;}

			jsumdtd=jsumdtd+dtd;
			/*printf ("%lf",timecounter,
			"days", "%d", i ,"%d",j,"\n");*/
			}
		/* == old NEXT j of BASIC === ;*/
		}

		/*5-3*/
		vx[i] = vx[i] + sumax * dtd;
		vy[i] = vy[i] + sumay * dtd;
		x[i] = x[i] + vx[i] * dtd;
		y[i] = y[i] + vy[i] * dtd;

		/* No-check Version Vectorspeed.;*/

		pixelcircle(x[i],y[i],1,frame,flag98);

	}/* == old NXET i of BASIC === ; */


/*5-4*/



}

/* ============================================================ ;*/
/*7)processEND decleare;*/
/* ============================================================ ;*/
processend:;
printf("\nDemonstration Finished.:Hit Any KEY. \n");
getchar();
closegraph();
}