 	hostvfield = (float4***) malloc(N*sizeof(float4**));
 	hostvfield[0] = (float4**) malloc(N*N*sizeof(float4*));
 	hostvfield[0][0] = (float4*) malloc(N*N*N*sizeof(float4));
-	std::cout << 2 << std::endl;
+//	std::cout << 2 << std::endl;
 	for (int i=1; i < N; i++) {
 		hostvfield[i] = (hostvfield[0] + i*N);
 	}

}

__global__ void RK4line(float4* lineoutput, double dt, unsigned int steps, float4 loc) {
	lineoutput[0] = tex3D(dataTex, make_float3(127.5,127.5,0));
	for (unsigned int i=1; i < steps; i++) {
		loc = loc + RK4step(loc,dt);
		lineoutput[i] = loc;


int main(void) {
	cudaArray* dataArray;
//	std::cout << 1 << std::endl;
	cudaChannelFormatDesc channelDesc = cudaCreateChannelDesc(32,32,32,32,cudaChannelFormatKindFloat);
	cudaExtent extent = make_cudaExtent(N , N, N);
	checkCudaErrors(cudaMalloc3DArray(&dataArray, &channelDesc,extent));
