• Advertisement
  • Popular Tags

  • Popular Now

  • Advertisement
  • Similar Content

    • By Fadey Duh
      Good evening everyone!

      I was wondering if there is something equivalent of  GL_NV_blend_equation_advanced for AMD?
      Basically I'm trying to find more compatible version of it.

      Thank you!
    • By Jens Eckervogt
      Hello guys, 
       
      Please tell me! 
      How do I know? Why does wavefront not show for me?
      I already checked I have non errors yet.
      using OpenTK; using System.Collections.Generic; using System.IO; using System.Text; namespace Tutorial_08.net.sourceskyboxer { public class WaveFrontLoader { private static List<Vector3> inPositions; private static List<Vector2> inTexcoords; private static List<Vector3> inNormals; private static List<float> positions; private static List<float> texcoords; private static List<int> indices; public static RawModel LoadObjModel(string filename, Loader loader) { inPositions = new List<Vector3>(); inTexcoords = new List<Vector2>(); inNormals = new List<Vector3>(); positions = new List<float>(); texcoords = new List<float>(); indices = new List<int>(); int nextIdx = 0; using (var reader = new StreamReader(File.Open("Contents/" + filename + ".obj", FileMode.Open), Encoding.UTF8)) { string line = reader.ReadLine(); int i = reader.Read(); while (true) { string[] currentLine = line.Split(); if (currentLine[0] == "v") { Vector3 pos = new Vector3(float.Parse(currentLine[1]), float.Parse(currentLine[2]), float.Parse(currentLine[3])); inPositions.Add(pos); if (currentLine[1] == "t") { Vector2 tex = new Vector2(float.Parse(currentLine[1]), float.Parse(currentLine[2])); inTexcoords.Add(tex); } if (currentLine[1] == "n") { Vector3 nom = new Vector3(float.Parse(currentLine[1]), float.Parse(currentLine[2]), float.Parse(currentLine[3])); inNormals.Add(nom); } } if (currentLine[0] == "f") { Vector3 pos = inPositions[0]; positions.Add(pos.X); positions.Add(pos.Y); positions.Add(pos.Z); Vector2 tc = inTexcoords[0]; texcoords.Add(tc.X); texcoords.Add(tc.Y); indices.Add(nextIdx); ++nextIdx; } reader.Close(); return loader.loadToVAO(positions.ToArray(), texcoords.ToArray(), indices.ToArray()); } } } } } And It have tried other method but it can't show for me.  I am mad now. Because any OpenTK developers won't help me.
      Please help me how do I fix.

      And my download (mega.nz) should it is original but I tried no success...
      - Add blend source and png file here I have tried tried,.....  
       
      PS: Why is our community not active? I wait very longer. Stop to lie me!
      Thanks !
    • By codelyoko373
      I wasn't sure if this would be the right place for a topic like this so sorry if it isn't.
      I'm currently working on a project for Uni using FreeGLUT to make a simple solar system simulation. I've got to the point where I've implemented all the planets and have used a Scene Graph to link them all together. The issue I'm having with now though is basically the planets and moons orbit correctly at their own orbit speeds.
      I'm not really experienced with using matrices for stuff like this so It's likely why I can't figure out how exactly to get it working. This is where I'm applying the transformation matrices, as well as pushing and popping them. This is within the Render function that every planet including the sun and moons will have and run.
      if (tag != "Sun") { glRotatef(orbitAngle, orbitRotation.X, orbitRotation.Y, orbitRotation.Z); } glPushMatrix(); glTranslatef(position.X, position.Y, position.Z); glRotatef(rotationAngle, rotation.X, rotation.Y, rotation.Z); glScalef(scale.X, scale.Y, scale.Z); glDrawElements(GL_TRIANGLES, mesh->indiceCount, GL_UNSIGNED_SHORT, mesh->indices); if (tag != "Sun") { glPopMatrix(); } The "If(tag != "Sun")" parts are my attempts are getting the planets to orbit correctly though it likely isn't the way I'm meant to be doing it. So I was wondering if someone would be able to help me? As I really don't have an idea on what I would do to get it working. Using the if statement is truthfully the closest I've got to it working but there are still weird effects like the planets orbiting faster then they should depending on the number of planets actually be updated/rendered.
    • By Jens Eckervogt
      Hello everyone, 
      I have problem with texture
      using System; using OpenTK; using OpenTK.Input; using OpenTK.Graphics; using OpenTK.Graphics.OpenGL4; using System.Drawing; using System.Reflection; namespace Tutorial_05 { class Game : GameWindow { private static int WIDTH = 1200; private static int HEIGHT = 720; private static KeyboardState keyState; private int vaoID; private int vboID; private int iboID; private Vector3[] vertices = { new Vector3(-0.5f, 0.5f, 0.0f), // V0 new Vector3(-0.5f, -0.5f, 0.0f), // V1 new Vector3(0.5f, -0.5f, 0.0f), // V2 new Vector3(0.5f, 0.5f, 0.0f) // V3 }; private Vector2[] texcoords = { new Vector2(0, 0), new Vector2(0, 1), new Vector2(1, 1), new Vector2(1, 0) }; private int[] indices = { 0, 1, 3, 3, 1, 2 }; private string vertsrc = @"#version 450 core in vec3 position; in vec2 textureCoords; out vec2 pass_textureCoords; void main(void) { gl_Position = vec4(position, 1.0); pass_textureCoords = textureCoords; }"; private string fragsrc = @"#version 450 core in vec2 pass_textureCoords; out vec4 out_color; uniform sampler2D textureSampler; void main(void) { out_color = texture(textureSampler, pass_textureCoords); }"; private int programID; private int vertexShaderID; private int fragmentShaderID; private int textureID; private Bitmap texsrc; public Game() : base(WIDTH, HEIGHT, GraphicsMode.Default, "Tutorial 05 - Texturing", GameWindowFlags.Default, DisplayDevice.Default, 4, 5, GraphicsContextFlags.Default) { } protected override void OnLoad(EventArgs e) { base.OnLoad(e); CursorVisible = true; GL.GenVertexArrays(1, out vaoID); GL.BindVertexArray(vaoID); GL.GenBuffers(1, out vboID); GL.BindBuffer(BufferTarget.ArrayBuffer, vboID); GL.BufferData(BufferTarget.ArrayBuffer, (IntPtr)(vertices.Length * Vector3.SizeInBytes), vertices, BufferUsageHint.StaticDraw); GL.GenBuffers(1, out iboID); GL.BindBuffer(BufferTarget.ElementArrayBuffer, iboID); GL.BufferData(BufferTarget.ElementArrayBuffer, (IntPtr)(indices.Length * sizeof(int)), indices, BufferUsageHint.StaticDraw); vertexShaderID = GL.CreateShader(ShaderType.VertexShader); GL.ShaderSource(vertexShaderID, vertsrc); GL.CompileShader(vertexShaderID); fragmentShaderID = GL.CreateShader(ShaderType.FragmentShader); GL.ShaderSource(fragmentShaderID, fragsrc); GL.CompileShader(fragmentShaderID); programID = GL.CreateProgram(); GL.AttachShader(programID, vertexShaderID); GL.AttachShader(programID, fragmentShaderID); GL.LinkProgram(programID); // Loading texture from embedded resource texsrc = new Bitmap(Assembly.GetEntryAssembly().GetManifestResourceStream("Tutorial_05.example.png")); textureID = GL.GenTexture(); GL.BindTexture(TextureTarget.Texture2D, textureID); GL.TexParameter(TextureTarget.Texture2D, TextureParameterName.TextureMagFilter, (int)All.Linear); GL.TexParameter(TextureTarget.Texture2D, TextureParameterName.TextureMinFilter, (int)All.Linear); GL.TexImage2D(TextureTarget.Texture2D, 0, PixelInternalFormat.Rgba, texsrc.Width, texsrc.Height, 0, PixelFormat.Bgra, PixelType.UnsignedByte, IntPtr.Zero); System.Drawing.Imaging.BitmapData bitmap_data = texsrc.LockBits(new Rectangle(0, 0, texsrc.Width, texsrc.Height), System.Drawing.Imaging.ImageLockMode.ReadOnly, System.Drawing.Imaging.PixelFormat.Format32bppRgb); GL.TexSubImage2D(TextureTarget.Texture2D, 0, 0, 0, texsrc.Width, texsrc.Height, PixelFormat.Bgra, PixelType.UnsignedByte, bitmap_data.Scan0); texsrc.UnlockBits(bitmap_data); GL.Enable(EnableCap.Texture2D); GL.BufferData(BufferTarget.TextureBuffer, (IntPtr)(texcoords.Length * Vector2.SizeInBytes), texcoords, BufferUsageHint.StaticDraw); GL.BindAttribLocation(programID, 0, "position"); GL.BindAttribLocation(programID, 1, "textureCoords"); } protected override void OnResize(EventArgs e) { base.OnResize(e); GL.Viewport(0, 0, ClientRectangle.Width, ClientRectangle.Height); } protected override void OnUpdateFrame(FrameEventArgs e) { base.OnUpdateFrame(e); keyState = Keyboard.GetState(); if (keyState.IsKeyDown(Key.Escape)) { Exit(); } } protected override void OnRenderFrame(FrameEventArgs e) { base.OnRenderFrame(e); // Prepare for background GL.Clear(ClearBufferMask.ColorBufferBit); GL.ClearColor(Color4.Red); // Draw traingles GL.EnableVertexAttribArray(0); GL.EnableVertexAttribArray(1); GL.BindVertexArray(vaoID); GL.UseProgram(programID); GL.BindBuffer(BufferTarget.ArrayBuffer, vboID); GL.VertexAttribPointer(0, 3, VertexAttribPointerType.Float, false, 0, IntPtr.Zero); GL.ActiveTexture(TextureUnit.Texture0); GL.BindTexture(TextureTarget.Texture3D, textureID); GL.BindBuffer(BufferTarget.ElementArrayBuffer, iboID); GL.DrawElements(BeginMode.Triangles, indices.Length, DrawElementsType.UnsignedInt, 0); GL.DisableVertexAttribArray(0); GL.DisableVertexAttribArray(1); SwapBuffers(); } protected override void OnClosed(EventArgs e) { base.OnClosed(e); GL.DeleteVertexArray(vaoID); GL.DeleteBuffer(vboID); } } } I can not remember where do I add GL.Uniform2();
    • By Jens Eckervogt
      Hello everyone
      For @80bserver8 nice job - I have found Google search. How did you port from Javascript WebGL to C# OpenTK.?
      I have been searched Google but it shows f***ing Unity 3D. I really want know how do I understand I want start with OpenTK But I want know where is porting of Javascript and C#?
       
      Thanks!
  • Advertisement
  • Advertisement
Sign in to follow this  

OpenGL problem with ray/box intersection in GLSL

This topic is 559 days old which is more than the 365 day threshold we allow for new replies. Please post a new topic.

If you intended to correct an error in the post then please contact us.

Recommended Posts

Hi all,
 
For part of an OpenGL ray-tracer I'm working on (as a hobby) I'm attempting to find the interval t in which a ray intersects an area-aligned-bounding-box (AABB).

I've created a GLSL function (left un-optimised for clarity) based on the slab method described here: https://tavianator.com/fast-branchless-raybounding-box-intersections/:

#version 330

const float INFINITY = 1.0 / 0.0;

struct AABB
{
    vec3 min;
    vec3 max;
};

struct Ray
{
    vec3 origin;
    vec3 direction;
};

float rayAABBIntersection (Ray ray, AABB aabb)
{
        vec3 inverseRayDirection = 1.0 / ray.direction;

	float t1;
	float t2;
	float tmin = -INFINITY;
	float tmax = INFINITY;

	if (ray.direction.x != 0.0) 
	{
	    t1 = (aabb.min.x - ray.origin.x) * inverseRayDirection.x;
	    t2 = (aabb.max.x - ray.origin.x) * inverseRayDirection.x;
	    tmin = max(tmin, min(t1, t2));
	    tmax = min(tmax, max(t1, t2));
	}
    
	if (ray.direction.y != 0.0) 
	{
	    t1 = (aabb.min.y - ray.origin.y) * inverseRayDirection.y;
	    t2 = (aabb.max.y - ray.origin.y) * inverseRayDirection.y;
	    tmin = max(tmin, min(t1, t2));
	    tmax = min(tmax, max(t1, t2));
	}
	
	if (ray.direction.z != 0.0) 
	{
	    t1 = (aabb.min.z - ray.origin.z) * inverseRayDirection.z;
	    t2 = (aabb.max.z - ray.origin.z) * inverseRayDirection.z;
	    tmin = max(tmin, min(t1, t2));
	    tmax = min(tmax, max(t1, t2));
	}
	
	if (tmin < 0)
	{
	    return tmax;
	}
	else
	{
	    return tmin;
	}
}

I've tested the function first in C by replacing `min` and `max` with `fmin` and `fmax` functions from <math.h>.

 

In C the code seems to work fine and correctly returns the t value for the ray intersection, however in the GLSL depending on the ray position and angle this gives me incorrect results.

 

If anyone has any insight about why the GLSL isn't functioning as expected, I'd be very grateful!

 

(Note: I'm debugging on OSX and have had problems building glslDebugger, so have had to resort to outputting colours to assist in debugging which proves very time consuming - if anyone knows of any Mac based tools which could help, I'd be grateful for those also).

 

 

Share this post


Link to post
Share on other sites
Advertisement

I could be mistaken but I guess that in 3.3 GLSL, division by zero would result in undefined behavior. From page 21 of the official spec (https://www.opengl.org/registry/doc/GLSLangSpec.3.30.6.pdf) :

 

"Similarly, treatment of conditions such as divide by 0 may lead to an unspecified result, but in no case should such a condition lead to the interruption or termination of processing. Generally, there are no signaling NaNs, and operating on NaNs (Not a Number) or infs (positive or negative infinities) gives undefined results."

Share this post


Link to post
Share on other sites

The division by zero is indeed a problem with this approach. E.g. the Ember raytracer just uses FLT_MIN if one of the direction elements is zero (at least in one of their older versions). Another option is to use the Separating Axis Test (SAT) for ray/box instead. This will not clip the ray though, but only provide a boolean answer whether the ray intersects the box. I use this in most of my tree traversals as I found it faster and more robust.

Share this post


Link to post
Share on other sites

Shameless self-promotion: http://www.gamedev.net/page/resources/_/technical/math-and-physics/intersection-math-algorithms-learn-to-derive-r3033 - this is probably more optimal algorithm for the GPU (I've been using (similar version to this) in OpenCL/CUDA for quite a bit of time, I'm sure GLSL is able to handle it too).

 

Don't try to build your code to avoid divisions by 0, especially in world of real time ray tracing, build it in such way so you can handle division by zeroes (use the properties of IEEE floating point math - like INF value or NAN value).

 

Let me try to dig something on my AABB intersection here on local machine.

 

EDIT:

 

To add here - this is my gpu-based Ray-AABB test done when working with KD-tree (and testing against AABB sorrounding whole tree):

float4 v1 = (aabbMin - o) * invRay;
float4 v2 = (aabbMax - o) * invRay;
float4 n = min(v1, v2);
float4 f = max(v1, v2);
float enter = max(n.x, max(n.y, n.z));
float exit = min(f.x, min(f.y, f.z));

if (exit > 0.0f && enter < exit)
{
    // We have a hit
}

Where:

  • aabbMin - is float4 which contains min point in X, Y and Z members and 1.0 in W member
  • aabbMax - is float4 which contains max point in X, Y and Z members and 1.0 in W member 
  • o - is ray origin point with 1.0 in W member
  • invRay - is float4(1.0 / rayDirection.x, 1.0 / rayDirection.y, 1.0 / rayDirection.z, 0.0)

Functions min (and respectively max) are working per component, like:

float4 min(float4 a, float4 b)
{
    return (float4)(min(a.x, b.x), min(a.y, b.y), min(a.z, b.z), 0.0);
}
Edited by Vilem Otte

Share this post


Link to post
Share on other sites

Thanks for the advice everyone!

 

@SapphireG and @Dirk the advice on the handling of infinite and NaN values was very helpful - it caused me to test fundamentally the output of min and max functions with NaN values...this has helped me confirm that my rendering error is not caused by this function.

 

@Vilem - thanks for sharing your code - your Ray AABB test looks very similar to an optimised one I had been testing:

float rayAABBIntersection (Ray ray, AABB aabb)
{
   float t[10];
   t[1] = (aabb.min.x - ray.origin.x) * ray.inverseDirection.x;
   t[2] = (aabb.max.x - ray.origin.x) * ray.inverseDirection.x;
   t[3] = (aabb.min.y - ray.origin.y) * ray.inverseDirection.y;
   t[4] = (aabb.max.y - ray.origin.y) * ray.inverseDirection.y;
   t[5] = (aabb.min.z - ray.origin.z) * ray.inverseDirection.z;
   t[6] = (aabb.max.z - ray.origin.z) * ray.inverseDirection.z;
   t[7] = max(max(min(t[1], t[2]), min(t[3], t[4])), min(t[5], t[6]));
   t[8] = min(min(max(t[1], t[2]), max(t[3], t[4])), max(t[5], t[6]));
   t[9] = (t[8] < 0 || t[7] > t[8]) ? INFINITY : t[7];
   return t[9];
}

but yours looks to be slightly more optimal (using vector operations to remove a few lines by the looks of it)!

I hope you don't mind if I use yours for inspiration! ;)

 

I've since moved to working in OpenGL 4.1 (GLSL #version 410) and by checking the results of `min` and `max` have confirmed they work the same way as my C code when presented with NaN values - (if one of the values is NaN, the other is always returned).

 

I think the error in my ray-tracing must be elsewhere....so I'd like to quickly break down the steps I'm taking to calculate a ray to see if my approach is valid:

  1. I upload a simple cube VBO to the GPU.
  2. In my vert shader I pass the vert-position (before applying a model-view-projection transform) to the fragment shader. This should give me the surface position for a fragment in model space. (Correct?)
  3. I multiply the eye position (in world space) by the inverse model-view matrix to get the eye position in model space and pass this to the fragment shader.
  4. To calculate a ray for each fragment, I use the eye-position as the origin, and the direction is calculated as the fragment surface position, minus the eye position, normalized.

Is this a valid approach for calculating a (model-space) ray per fragment?

 

Many thanks again for all your help, it's greatly appreciated!

Share this post


Link to post
Share on other sites

 

) I'm attempting to find the interval t in which a ray intersects an area-aligned-bounding-box (AABB).

 

I've created a GLSL function (left un-optimised for clarity) based on the slab method described here

 

 

you define z_far thats right the same z_far you define but you store it in shader

 

in vertex shader you pass to fragment vertex position as: aabb_world_matrix * v 

 

in fragment shader

float t = DistanceBetweenPoints(Camera_position, world_vertex_pos) / z_far;

 

where in vertex shader

world_vertex_pos.x = dp43(WM1, Vpos);
world_vertex_pos.y = dp43(WM2, Vpos);
world_vertex_pos.z = dp43(WM3, Vpos);
float dp43(vec4 matrow, vec3 p)
{
return ( (matrow.x*p.x) + (matrow.y*p.y) + (matrow.z*p.z) + matrow.w );
}

where WM1 is vec4 and stores world matrix row 1

where WM2 is vec4 and stores world matrix row 2

where WM3 is vec4 and stores world matrix row 3

where WM4 is vec4 and stores world matrix row 4

 

and Vpos is an attribute that i pass ass vertex buffer anyway its old glsl

 

hypnotized-smiley.gif

Share this post


Link to post
Share on other sites
Sign in to follow this  

  • Advertisement