     float frustum_shift_rz_ry_rx[16];
     matrixMul4(frustum_shift_rz_ry_rx, frustum, shift_rz_ry_rx);
 
+#if 0
     glEnable(GL_CULL_FACE);
     glEnable(GL_DEPTH_TEST);
     glUseProgram(simplePrg);

  glfwMakeContextCurrent(win);
  gladLoadGL(glfwGetProcAddress);

  int deviceCount;
  cudaGetDeviceCount(&deviceCount);
  if (deviceCount == 0) {
    fprintf(stderr, "No CUDA-enabled devices available.");
    return EXIT_FAILURE;
  }
  for (int i = 0; i < deviceCount; i++) {
    cudaDeviceProp dev_prop;
    cudaGetDeviceProperties(&dev_prop, i);
    fprintf(stderr, "%c[%i] %s cap=%d.%d\n", i == deviceIndex ? '*' : ' ', i, dev_prop.name, dev_prop.major, dev_prop.minor);
  }
  if (deviceIndex < 0 || deviceCount <= deviceIndex) {
    fprintf(stderr, "Illegal CUDA device index %d\n", deviceIndex);
    return EXIT_FAILURE;
  }

    glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
    CHECK_GL;

#endif
    glUseProgram(simplePrg);
    glBindVertexArray(cudaVbo);
    glUniformMatrix4fv(0, 1, GL_FALSE, rz_ry_rx);
    glUniformMatrix4fv(1, 1, GL_FALSE, frustum_shift_rz_ry_rx);
    glUniform4f(2, 1.f, 1.f, 1.f, 1.f);
    glDrawArrays(GL_POINTS, 0, N);


    glfwSwapBuffers(win);
    glfwPollEvents();
