__global__ void addKernel(float *POSX, float *POSY, float *POSR, float *AX_M, float *AY_M)
{
	int first = threadIdx.x;
	for (int next = 0; next < N; next++) {
		if (first != next) {
			int den;
			den = 10 * sqrt((POSX[next] - POSX[first])*(POSX[next] - POSX[first]) + (POSY[next] - POSY[first])*(POSY[next] - POSY[first]));
			AX_M[first] = (POSX[next] - POSX[first]) / den;
			AY_M[first] = (POSY[next] - POSY[first]) / den;
		}
	}
	POSR[first] += AX_M[first];
	POSR[first + N] += AY_M[first];
}

void GPU_Math()
{
	float *CU_POSX, *CU_POSY, *CU_POSR, *CU_AX_M, *CU_AY_M;
	int size = N;

	cudaMalloc((void**)&CU_POSX, size);
	cudaMalloc((void**)&CU_POSY, size);
	cudaMalloc((void**)&CU_POSR, size*2);
	cudaMalloc((void**)&CU_AX_M, size);
	cudaMalloc((void**)&CU_AY_M, size);

	cudaMemcpy(CU_POSX, POSX, size, cudaMemcpyKind::cudaMemcpyHostToDevice);
	cudaMemcpy(CU_POSY, POSY, size, cudaMemcpyKind::cudaMemcpyHostToDevice);
	cudaMemcpy(CU_POSR, POSR, size*2, cudaMemcpyKind::cudaMemcpyHostToDevice);
	cudaMemcpy(CU_AX_M, AX_M, size, cudaMemcpyKind::cudaMemcpyHostToDevice);
	cudaMemcpy(CU_AY_M, AY_M, size, cudaMemcpyKind::cudaMemcpyHostToDevice);

	addKernel<<<1, size>>>(CU_POSX, CU_POSY, CU_POSR, CU_AX_M, CU_AY_M);

	cudaMemcpy(POSR, CU_POSR, size*2, cudaMemcpyKind::cudaMemcpyDeviceToHost);
	cudaMemcpy(AX_M, CU_AX_M, size, cudaMemcpyKind::cudaMemcpyDeviceToHost);
	cudaMemcpy(AY_M, CU_AY_M, size, cudaMemcpyKind::cudaMemcpyDeviceToHost);

	cudaFree(CU_POSX);
	cudaFree(CU_POSY);
	cudaFree(CU_POSR);
	cudaFree(CU_AX_M);
	cudaFree(CU_AY_M);
}