Why does GL divide `gl_Position` by W for you rather than letting you do it yourself?

by

I’d like to extent on BDL’s answer. It is not only about the perspective interpolation. It is also about the clipping. The space the value gl_Position is supposed to be provided in is called clip space, and this is before the division by w.

The (default) clip volume of OpenGL is defined in clip space as

-w <= x,y,z <= w   (with w varying per vertex)

After the division by w we get

-1 <= x,y,z <= 1   (in NDC coordinates).

However, if you try to do the clipping after the division by w, and would check against that cube in NDC, you get a problem, because all clip space points fullfilling this:

 w <= x,y,z <= -w (in clip space)

will also fullfill the NDC constraint.

The thing here is that points behind the camera will be transformed to somewhere in front of the camera, mirrored (since x/-1 is the same as -x/1). This also happens to the z coordinate. One might argue that this is irrelevant, because any point behind the camera is projected behind (in the sense of more far away than) the far plane, as per the construction of the typical projection matrix, so it will lie outside of the viewing volume in either case.

But if you have a primitive where at least one point is inside the view volume, and at least one point is behind the camera, you should have a primitive which intersects the near plane also. However, after the division by w, it will intersect the far plane now!. So clipping in NDC space, after the division, is much harder to get right. I tried to visualize this in this drawing:

top-down view of eye space and NDC with and without clipping
(the drawing is to-scale, the depth range of projection is much shorter than anyone would typically use, to better illustrate the issue).

The clipping is done as a fixed-function stage in hardware and it has to be done before the division, hence you should provide the correct clip-space coordinates to work on.

(Note: actual GPUs might not use an extra clipping stage at all, they actually might also use a clipless rasterizer, like it is speculated in Fabian Giesen’s blog article there. There are some algorithms like Olano and Greer (1997). However, this all works by doing the rasterization directly in homogenous coordinates, so we still need the w…)

More Related Contents:

  1. Should I ever use a `vec3` inside of a uniform buffer or shader storage buffer object?
  2. Frequency of shader invocations in rendering commands
  3. Rendering meshes with multiple indices
  4. How to render depth linearly in modern OpenGL with gl_FragCoord.z in fragment shader?
  5. Why transform normals with the transpose of the inverse of the modelview matrix?
  6. glVertexAttribPointer and glVertexAttribFormat: What’s the difference?
  7. How to implement 2D raycasting light effect in GLSL
  8. Do conditional statements slow down shaders?
  9. GLSL fixed function fragment program replacement
  10. Passing a list of values to fragment shader
  11. In a fragment shader, why can’t I use a flat input integer to index a uniform array of sampler2D?
  12. What is state-of-the-art for text rendering in OpenGL as of version 4.1? [closed]
  13. How to debug a GLSL shader?
  14. How to use Pivot Point in Transformations
  15. glLineStipple deprecated in OpenGL 3.1
  16. How to Using the #include in glsl support ARB_shading_language_include
  17. From RGB to HSV in OpenGL GLSL
  18. Current state and solutions for OpenGL over Windows Remote [closed]
  19. How do you get the modelview and projection matrices in OpenGL?
  20. Multiple textures in GLSL – only one works
  21. Is it possible to rotate an object around its own axis and not around the base coordinate’s axis?
  22. GLSL, semaphores?
  23. GLSL – Using custom output attribute instead of gl_Position
  24. Reflection and refraction impossible without recursive ray tracing?
  25. What are the Attribute locations for fixed function pipeline in OpenGL 4.0++ core profile?
  26. GCC, stringification, and inline GLSL?
  27. Opengl pixel perfect 2D drawing
  28. Drawing a border on a 2d polygon with a fragment shader
  29. OpenGL: scale then translate? and how?
  30. how to enable vertical sync in opengl?

Leave a Comment