#pragma OPENCL EXTENSION cl_khr_fp64 : enable




__kernel void n_body(__global  double *X, __global double *Y, __global double *Z,__global  double *UX, __global double *UY, __global double *UZ,__global  double *M,int N,__global  double *X_N, __global double *Y_N, __global double *Z_N) 



{
		  
int i =get_global_id(0);
#pragma OPENCL EXTENSION cl_khr_fp64 : enable

	 double ax,ay,az,f,vx,vy,vz,inv_dist,inv_dist_cube,t,dx,dy,dz,offset;
        ax=0;ay=0;az=0;f=0;vx=0;vy=0;vz=0;inv_dist=0;inv_dist_cube=0;t=3600;dx=0;dy=0;dz=0;offset=1;
		for (int j=0; j<N; j++)
			{

		dx=X[j]-X[i];
                dy=Y[j]-Y[i];
                dz=Z[j]-Z[i];


                inv_dist=1.0/(sqrt(dx*dx+dy*dy+dz*dz+offset));
                inv_dist_cube=inv_dist*inv_dist*inv_dist;

                f=M[j]*inv_dist_cube;

                ax+=f*dx;
                ay+=f*dy;
                az+=f*dz;
       


            }
       vx=UX[i] +ax*t;
        vy=UY[i] +ay*t;
        vz=UZ[i] +az*t;
//	vx=4;
//        vy=4;
//       vz=4;
//       barrier(CLK_LOCAL_MEM_FENCE | CLK_GLOBAL_MEM_FENCE);

     //   X_N[i]=((vx*vx -UX[i]*UX[i])/(2*ax))+X[i];
     //   Y_N[i]=((vy*vy -UY[i]*UY[i])/(2*ay))+Y[i];
     //   Z_N[i]=((vz*vz -UZ[i]*UZ[i])/(2*az))+Z[i];

	 X_N[i]=((vx*vx -UX[i]*UX[i])/(2*ax))+X[i]; 
	 Y_N[i]=((vy*vy -UY[i]*UY[i])/(2*ay))+Y[i];

//	X_N[i]=X[i];
//        Y_N[i]=3.14;

 //       Z_N[i]=((vz*vz -UZ[i]*UZ[i])/(2*az))+Z[i];
	Z_N[i]=Z[i]+2;

        UX[i]=vx;
        UY[i]=vy;
        UZ[i]=vz;
	

  ax=0;ay=0;az=0;f=0;vx=0;vy=0;vz=0;inv_dist=0;inv_dist_cube=0;t=1000;dx=0;dy=0;dz=0;offset=1;


}