djairoh
/
cuda-based-raytrace
mirror of https://github.com/MartinOpat/cuda-based-raytrace.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Renamed var. for consistency with sphere.h

This commit is contained in:
Martin Opat 2024-12-23 22:16:22 +01:00
parent 3dc6c8e05f
commit 79d0587d97
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/linalg/mat.h
View File

@ -1,6 +1,6 @@
#pragma once
__device__ Vec3 computeGradient(float* volData, const int d_volumeWidth, const int d_volumeHeight, const int d_volumeDepth, int x, int y, int z) {
__device__ Vec3 computeGradient(float* volumeData, const int d_volumeWidth, const int d_volumeHeight, const int d_volumeDepth, int x, int y, int z) {
// Finite difference for partial derivatives.
// For boundary voxels - clamp to the boundary.
// Normal should point from higher to lower intensities
@ -13,12 +13,12 @@ __device__ Vec3 computeGradient(float* volData, const int d_volumeWidth, const i
int zp = min(z + 1, d_volumeDepth - 1);
// Note: Assuming data is linearized (idx = z*w*h + y*w + x) TODO: Unlinearize if data not linear
float gx = volData[z * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + xp]
- volData[z * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + xm];
float gy = volData[z * d_volumeWidth * d_volumeHeight + yp * d_volumeWidth + x ]
- volData[z * d_volumeWidth * d_volumeHeight + ym * d_volumeWidth + x ];
float gz = volData[zp * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + x ]
- volData[zm * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + x ];
float gx = volumeData[z * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + xp]
- volumeData[z * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + xm];
float gy = volumeData[z * d_volumeWidth * d_volumeHeight + yp * d_volumeWidth + x ]
- volumeData[z * d_volumeWidth * d_volumeHeight + ym * d_volumeWidth + x ];
float gz = volumeData[zp * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + x ]
- volumeData[zm * d_volumeWidth * d_volumeHeight + y * d_volumeWidth + x ];
return Vec3(gx, gy, gz);
}