Sign in to follow this  
drowne

Different behaviour between opencl c99 and c++ with the same implementation of a realtime voxel raycaster

Recommended Posts

[color=#000000][font=Arial,]
I am working with opencl to develop a voxel raycasting engine. I am trying to do something similar to Gigavoxels by Crassin (http://maverick.inria.fr/Publications/2011/Cra11/). In this paper, they are using an octree to store the voxel data. For the moment I am just trying to descend inside the octree until I reach a leaf that contains rendering data.[/font][/color][color=#000000][font=Arial,]
I have made two implementations: one in OpenCl on the GPU and another C++ on the CPU. The problem I am experiencing is that on the GPU the algorithm is going through a wrong number of levels until reaching a leaf inside the octree. The CPU version gives correct results. The algorithm for both versions is the same and the code is almost similar.[/font][/color][color=#000000][font=Arial,]
Do you guys know what could be the problem? Could it be a hardware problem, an OpenCl problem or am I doing something wrong? I am experiencing the same result on three different nVidia GPUs.[/font][/color]

Share this post


Link to post
Share on other sites
[quote name='coderchris' timestamp='1339425432' post='4948171']
You said the code is similar, but I'm guessing not exactly the same. I'd bet the problem lies somewhere in the minor differences...
[/quote]


[color=#000000][font=Arial,]Here is the C++ Code:[/font][/color]
[color=#000000][font=Arial,][code]
// Calculate actual ray stepping position
glm::vec4 pos = eyeRay_o + eyeRay_d * t;

uint offset = 0;
//check if root is leaf
uint leafFlag = GetLeafBit(octreeNodes[0]);
//get children address of root
uint childrenAddress = GetChildAddress(octreeNodes[0]);

while (iterations < 30) {
iterations++;

// Calculate subdivision offset
offset = (uint)(pos.x * 2) + (uint)(pos.y * 2) * 2 + (uint)(pos.z * 2) * 4;

if (leafFlag == 1) {
//return some colour and exit the loop
break;
}
else
{
glm::uvec4 off = glm::uvec4(pos.x * 2, pos.y * 2, pos.z * 2, pos.w * 2);
pos.x = 2 * pos.x - off.x;
pos.y = 2 * pos.y - off.y;
pos.z = 2 * pos.z - off.z;
pos.w = 2 * pos.w - off.w;
}

// Extract node data from the children
finalAddress = childrenAddress + offset;
leafFlag = GetLeafBit(nodes[finalAddress]);
childrenAddress = GetChildAddress(nodes[finalAddress]);
}
[/code]

Here is the OpenCL Code:[/font][/color]

[code]
// Calculate actual ray stepping position
float4 position = rayOrigin + rayDirection * t;
uint offset = 0;
//check if root is leaf
uint leafFlag = extractOctreeNodeLeaf(octreeNodes[0]);
//get children address of root
uint childrenAddress = extractOctreeNodeAddress(octreeNodes[0]);

//position will be in the [0, 1] interval
//size of octree is 1
while (iterations < 30) {
iterations++;

//calculate the index of the next child based on the position in the current subdivision
offset = (uint)(position.x * 2) + (uint)(position.y * 2) * 2 + (uint)(position.z * 2) * 4;

if (leafFlag == 1) {
//return some colour and exit the loop
break;
}
else
{
//transform the position inside the parent
//to the position inside the child subdivision
//size of child will be considered to be 1
uint4 off;
off.x = floor(position.x * 2);
off.y = floor(position.y * 2);
off.z = floor(position.z * 2);
off.w = floor(position.w * 2);
position = 2 * position - off;
}

// Extract node data from the children
finalAddress = childrenAddress + offset;
leafFlag = extractOctreeNodeLeaf(octreeNodes[finalAddress]);
//each node has an index to an array of 8 children - the index points to the first child
childrenAddress = extractOctreeNodeAddress(octreeNodes[finalAddress]);
}
[/code] Edited by drowne

Share this post


Link to post
Share on other sites
Probably the error is in extractOctreeNodeLeaf or extractOctreeNodeAddress or similar. Try posting the entire code files, and perhaps someone will spot an error.

Also, floating point calculations are often not guaranteed to yield exactly the same results when run on different devices.

Share this post


Link to post
Share on other sites
[quote name='Erik Rufelt' timestamp='1339428887' post='4948188']
Probably the error is in extractOctreeNodeLeaf or extractOctreeNodeAddress or similar. Try posting the entire code files, and perhaps someone will spot an error.

Also, floating point calculations are often not guaranteed to yield exactly the same results when run on different devices.
[/quote]


Both functions just do some bit operations:

OpenCL version:

[code]
inline char extractOctreeNodeLeaf(uint value) {
value = value >> 1;
return value & 1;
}

inline uint extractOctreeNodeAddress(uint value) {
return value >> 2;
}
[/code]

C++ version:

[code]
inline byte GetLeafBit(uint value)
{
value = value >> 0x1;
return value & 0x1;
}

inline uint GetChildAddress(uint value)
{
return value >> 0x2;
}
[/code]

Share this post


Link to post
Share on other sites
You should try verifying that the "octreeNodes" array is the same (byte for byte) on the GPU version as it is on the CPU version.

If thats true, also make sure your test is using the same values for rayOrigin, rayDirection, and t.

And like Erik said, it could just be floating point playing tricks on you.

It probably wont make a difference, but you could try changing the opencl version from using floor() inside the else to casting (like you do in the CPU version).

Share this post


Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now

Sign in to follow this