From af7b163040e76c5bbd64742738653531e24e77e6 Mon Sep 17 00:00:00 2001 From: sotech117 Date: Fri, 8 Dec 2023 12:33:50 -0500 Subject: add 4d rotation ops and view matrix calculation --- src/4dvecops/vectoroperations.cpp | 33 +++++++++++++++++++++++++++++++++ 1 file changed, 33 insertions(+) (limited to 'src/4dvecops/vectoroperations.cpp') diff --git a/src/4dvecops/vectoroperations.cpp b/src/4dvecops/vectoroperations.cpp index d41dad3..a54b758 100644 --- a/src/4dvecops/vectoroperations.cpp +++ b/src/4dvecops/vectoroperations.cpp @@ -27,4 +27,37 @@ glm::vec4 dot4( glm::vec4 u, glm::vec4 v) { return {u[0] * v[0], u[1] * v[1], u[2] * v[2], u[3] * v[3]}; +} + +glm::mat4 getViewMatrix4( + glm::vec4 fromPoint, + glm::vec4 toPoint, + glm::vec4 upVector, + glm::vec4 overVector) { + + // calculate e3 basis vector, the transformation col of view matrix + if (glm::distance(fromPoint, toPoint) < 0.0001f) { + throw std::runtime_error("fromPoint and toPoint are the same"); + } + glm::vec4 e3 = glm::normalize(fromPoint - toPoint); + + // calculate e2 basis vector, from the combinatory cross of up and over with e3 + glm::vec4 e2 = cross4(upVector, overVector, e3); + e2 = glm::normalize(e2); + if (glm::distance(e2, glm::vec4{0, 0, 0, 1}) < 0.0001f) { + throw std::runtime_error("invalid up vector"); + } + + // calculate e1 basis vector, from the cross of only the over vector + glm::vec4 e1 = cross4(overVector, e3, e2); + e1 = glm::normalize(e1); + if (glm::distance(e1, glm::vec4{0, 0, 0, 1}) < 0.0001f) { + throw std::runtime_error("invalid over vector"); + } + + // calculate e0 basis vector, the 4d orthogonal vector to the other 3 bases + glm::vec4 e0 = cross4(e3, e2, e1); + e0 = glm::normalize(e0); + + return {e2, e1, e0, e3}; } \ No newline at end of file -- cgit v1.2.3-70-g09d2 From f52fe6b21ae625fa325af4d78704abc8b2cfa689 Mon Sep 17 00:00:00 2001 From: sotech117 Date: Fri, 8 Dec 2023 12:36:52 -0500 Subject: small rename --- CMakeLists.txt | 2 +- src/4dvecops/vec4operations.cpp | 63 +++++++++++++++++++++++++++++++++++++++ src/4dvecops/vectoroperations.cpp | 63 --------------------------------------- 3 files changed, 64 insertions(+), 64 deletions(-) create mode 100644 src/4dvecops/vec4operations.cpp delete mode 100644 src/4dvecops/vectoroperations.cpp (limited to 'src/4dvecops/vectoroperations.cpp') diff --git a/CMakeLists.txt b/CMakeLists.txt index a55919f..c0bf031 100644 --- a/CMakeLists.txt +++ b/CMakeLists.txt @@ -63,7 +63,7 @@ add_executable(${PROJECT_NAME} src/accelerate/kdtree.h src/accelerate/bvh.cpp src/accelerate/bvh.h - src/4dvecops/vectoroperations.cpp + src/4dvecops/vec4operations.cpp src/4dvecops/rotations4d.cpp ) diff --git a/src/4dvecops/vec4operations.cpp b/src/4dvecops/vec4operations.cpp new file mode 100644 index 0000000..a54b758 --- /dev/null +++ b/src/4dvecops/vec4operations.cpp @@ -0,0 +1,63 @@ +#include "raytracer/raytracer.h" + +// vector operations on 4d vectors, +// reference: https://hollasch.github.io/ray4/Four-Space_Visualization_of_4D_Objects.html#chapter5 + +glm::vec4 cross4( + glm::vec4 u, + glm::vec4 v, + glm::vec4 w) { + float a = (v[0] * w[1]) - (v[1] * w[0]); + float b = (v[0] * w[2]) - (v[2] * w[0]); + float c = (v[0] * w[3]) - (v[3] * w[0]); + float d = (v[1] * w[2]) - (v[2] * w[1]); + float e = (v[1] * w[3]) - (v[3] * w[1]); + float f = (v[2] * w[3]) - (v[3] * w[2]); + + glm::vec4 result; + result[0] = (u[1] * f) - (u[2] * e) + (u[3] * d); + result[1] = -(u[0] * f) + (u[2] * c) - (u[3] * b); + result[2] = (u[0] * e) - (u[1] * c) + (u[3] * a); + result[3] = -(u[0] * d) + (u[1] * b) - (u[2] * a); + + return result; +} + +glm::vec4 dot4( + glm::vec4 u, + glm::vec4 v) { + return {u[0] * v[0], u[1] * v[1], u[2] * v[2], u[3] * v[3]}; +} + +glm::mat4 getViewMatrix4( + glm::vec4 fromPoint, + glm::vec4 toPoint, + glm::vec4 upVector, + glm::vec4 overVector) { + + // calculate e3 basis vector, the transformation col of view matrix + if (glm::distance(fromPoint, toPoint) < 0.0001f) { + throw std::runtime_error("fromPoint and toPoint are the same"); + } + glm::vec4 e3 = glm::normalize(fromPoint - toPoint); + + // calculate e2 basis vector, from the combinatory cross of up and over with e3 + glm::vec4 e2 = cross4(upVector, overVector, e3); + e2 = glm::normalize(e2); + if (glm::distance(e2, glm::vec4{0, 0, 0, 1}) < 0.0001f) { + throw std::runtime_error("invalid up vector"); + } + + // calculate e1 basis vector, from the cross of only the over vector + glm::vec4 e1 = cross4(overVector, e3, e2); + e1 = glm::normalize(e1); + if (glm::distance(e1, glm::vec4{0, 0, 0, 1}) < 0.0001f) { + throw std::runtime_error("invalid over vector"); + } + + // calculate e0 basis vector, the 4d orthogonal vector to the other 3 bases + glm::vec4 e0 = cross4(e3, e2, e1); + e0 = glm::normalize(e0); + + return {e2, e1, e0, e3}; +} \ No newline at end of file diff --git a/src/4dvecops/vectoroperations.cpp b/src/4dvecops/vectoroperations.cpp deleted file mode 100644 index a54b758..0000000 --- a/src/4dvecops/vectoroperations.cpp +++ /dev/null @@ -1,63 +0,0 @@ -#include "raytracer/raytracer.h" - -// vector operations on 4d vectors, -// reference: https://hollasch.github.io/ray4/Four-Space_Visualization_of_4D_Objects.html#chapter5 - -glm::vec4 cross4( - glm::vec4 u, - glm::vec4 v, - glm::vec4 w) { - float a = (v[0] * w[1]) - (v[1] * w[0]); - float b = (v[0] * w[2]) - (v[2] * w[0]); - float c = (v[0] * w[3]) - (v[3] * w[0]); - float d = (v[1] * w[2]) - (v[2] * w[1]); - float e = (v[1] * w[3]) - (v[3] * w[1]); - float f = (v[2] * w[3]) - (v[3] * w[2]); - - glm::vec4 result; - result[0] = (u[1] * f) - (u[2] * e) + (u[3] * d); - result[1] = -(u[0] * f) + (u[2] * c) - (u[3] * b); - result[2] = (u[0] * e) - (u[1] * c) + (u[3] * a); - result[3] = -(u[0] * d) + (u[1] * b) - (u[2] * a); - - return result; -} - -glm::vec4 dot4( - glm::vec4 u, - glm::vec4 v) { - return {u[0] * v[0], u[1] * v[1], u[2] * v[2], u[3] * v[3]}; -} - -glm::mat4 getViewMatrix4( - glm::vec4 fromPoint, - glm::vec4 toPoint, - glm::vec4 upVector, - glm::vec4 overVector) { - - // calculate e3 basis vector, the transformation col of view matrix - if (glm::distance(fromPoint, toPoint) < 0.0001f) { - throw std::runtime_error("fromPoint and toPoint are the same"); - } - glm::vec4 e3 = glm::normalize(fromPoint - toPoint); - - // calculate e2 basis vector, from the combinatory cross of up and over with e3 - glm::vec4 e2 = cross4(upVector, overVector, e3); - e2 = glm::normalize(e2); - if (glm::distance(e2, glm::vec4{0, 0, 0, 1}) < 0.0001f) { - throw std::runtime_error("invalid up vector"); - } - - // calculate e1 basis vector, from the cross of only the over vector - glm::vec4 e1 = cross4(overVector, e3, e2); - e1 = glm::normalize(e1); - if (glm::distance(e1, glm::vec4{0, 0, 0, 1}) < 0.0001f) { - throw std::runtime_error("invalid over vector"); - } - - // calculate e0 basis vector, the 4d orthogonal vector to the other 3 bases - glm::vec4 e0 = cross4(e3, e2, e1); - e0 = glm::normalize(e0); - - return {e2, e1, e0, e3}; -} \ No newline at end of file -- cgit v1.2.3-70-g09d2