 			p[v1 * V + v2] = v1;
 	}
 
+	int v;
 	FOR (v, 0, V - 1)
 		G[v * V + v] = 0;
 

	// cudaDeviceSynchronize waits for the kernel to finish, and returns
        cudaStatus = cudaDeviceSynchronize();
        HANDLE_ERROR(cudaStatus);
	
	FOR(block, 0, numOfBlock) 
	{
		fw_kernel_phase_one<VIRTUAL_BLOCK_SIZE><<<1, dimBlockP1>>>(block, n, pitch_int, dev_d, dev_p);

*/
template <int BLOCK_SIZE> __global__ void fw_kernel_phase_one(const unsigned int block, const unsigned int n, const size_t pitch, int * const d, int * const p)
{
	int i;
	int newPath;
	int newPred;

