Changeset b545cad

Timestamp:

Nov 11, 2019, 2:41:41 PM (5 years ago)

Author:

Dmitry Kobets <dkobets@…>

Branches:

ADT, arm-eh, ast-experimental, enum, forall-pointer-decay, jacob/cs343-translation, master, new-ast-unique-expr, pthread-emulation, qualifiedEnum

Children:

a39fd1d

Parents:

44f41997

git-author:

Dmitry Kobets <dkobets@…> (10/04/19 18:51:57)

git-committer:

Dmitry Kobets <dkobets@…> (11/11/19 14:41:41)

Message:

Added refraction

File:

• libcfa/src/vector.hfa (modified) (4 diffs)

libcfa/src/vector.hfa

@@ -1 +1 @@
 #pragma once
-#include "math.hfa"
+#include <math.hfa>
 #include <iostream.hfa>
 
 //---------------------- Vector Types ----------------------
 // TODO: make generic, as per glm
-
 
 struct vec2 {
@@ -21 +20 @@
         return os;
     }
-    void ?|?( ostype & os, const vec2 v ) {
+    void ?|?( ostype & os, const vec2& v ) {
         (ostype &)(os | v); ends( os );
     }
+}
+
+void ?{}(vec2& vec, zero_t) with (vec) {
+    x = y = 0;
+}
+void ?{}(vec2& vec, vec2& other) with (vec) {
+    [x,y] = other.[x,y];
 }
 
@@ -31 +37 @@
 vec2 ?*?(const vec2& v, float scalar) with (v) {
     return [x * scalar, y * scalar];
+}
+vec2 ?*?(float scalar, const vec2& v) {
+    return v * scalar;
 }
 vec2 ?/?(const vec2& v, float scalar) with (v) {
@@ -58 +67 @@
     return v / sqrt(dot(v, v));
 }
+
+// project vector u onto vector v
+static inline vec2 project(const vec2& u, const vec2& v) {
+    vec2 v_norm = normalize(v);
+    return v_norm * dot(u, v_norm);
+}
+
+/* returns the reflection direction : v - 2.0 * project(v, n)
+ * for incident vector v and surface normal n
+ */
+static inline vec2 reflect(const vec2& v, const vec2& n) {
+    return v - 2 * project(v, n);
+}
+
+// incident vector v, surface normal n
+// eta = ratio of indices of refraction between starting material and
+// entering material (i.e., from air to water, eta = 1/1.33)
+static inline vec2 refract(const vec2& v, const vec2& n, float eta) {
+    float dotValue = dot(n, v);
+    float k = 1 - eta * eta * (1 - dotValue * dotValue);
+    if (k < 0) {
+        return 0;
+    }
+    return eta * v - (eta * dotValue + sqrt(k)) * n;
+}