Changes in libcfa/src/vec/vec.hfa [fd54fef:250dbae]

File:

• libcfa/src/vec/vec.hfa (modified) (7 diffs)

libcfa/src/vec/vec.hfa

@@ -1 @@
-//
-// Cforall Version 1.0.0 Copyright (C) 2021 University of Waterloo
-//
-// The contents of this file are covered under the licence agreement in the
-// file "LICENCE" distributed with Cforall.
-//
-// io/types.hfa --
-//
-// Author           : Dimitry Kobets
-// Created On       :
-// Last Modified By :
-// Last Modified On :
-// Update Count     :
-//
-
 #pragma once
 
 #include <math.hfa>
 
-trait fromint(T) {
+trait fromint(otype T) {
     void ?{}(T&, int);
 };
-trait zeroinit(T) {
+trait zeroinit(otype T) {
     void ?{}(T&, zero_t);
 };
-trait zero_assign(T) {
+trait zero_assign(otype T) {
     T ?=?(T&, zero_t);
 };
-trait subtract(T) {
+trait subtract(otype T) {
     T ?-?(T, T);
 };
-trait negate(T) {
+trait negate(otype T) {
     T -?(T);
 };
-trait add(T) {
+trait add(otype T) {
     T ?+?(T, T);
 };
-trait multiply(T) {
+trait multiply(otype T) {
     T ?*?(T, T);
 };
-trait divide(T) {
+trait divide(otype T) {
     T ?/?(T, T);
 };
-trait lessthan(T) {
+trait lessthan(otype T) {
     int ?<?(T, T);
 };
-trait equality(T) {
+trait equality(otype T) {
     int ?==?(T, T);
 };
-trait sqrt(T) {
+trait sqrt(otype T) {
     T sqrt(T);
 };
@@ -68 +53 @@
 }
 
-trait dottable(V, T) {
+trait dottable(otype V, otype T) {
     T dot(V, V);
 };
@@ -74 +59 @@
 static inline {
 
-forall(T | sqrt(T), V | dottable(V, T))
+forall(otype T | sqrt(T), otype V | dottable(V, T))
 T length(V v) {
    return sqrt(dot(v, v));
 }
 
-forall(T, V | dottable(V, T))
+forall(otype T, otype V | dottable(V, T))
 T length_squared(V v) {
    return dot(v, v);
 }
 
-forall(T, V | { T length(V); } | subtract(V))
+forall(otype T, otype V | { T length(V); } | subtract(V))
 T distance(V v1, V v2) {
     return length(v1 - v2);
 }
 
-forall(T, V | { T length(V); V ?/?(V, T); })
+forall(otype T, otype V | { T length(V); V ?/?(V, T); })
 V normalize(V v) {
     return v / length(v);
@@ -95 +80 @@
 
 // Project vector u onto vector v
-forall(T, V | dottable(V, T) | { V normalize(V); V ?*?(V, T); })
+forall(otype T, otype V | dottable(V, T) | { V normalize(V); V ?*?(V, T); })
 V project(V u, V v) {
     V v_norm = normalize(v);
@@ -102 +87 @@
 
 // Reflect incident vector v with respect to surface with normal n
-forall(T | fromint(T), V | { V project(V, V); V ?*?(T, V); V ?-?(V,V); })
+forall(otype T | fromint(T), otype V | { V project(V, V); V ?*?(T, V); V ?-?(V,V); })
 V reflect(V v, V n) {
     return v - (T){2} * project(v, n);
@@ -111 +96 @@
 // entering material (i.e., from air to water, eta = 1/1.33)
 // v and n must already be normalized
-forall(T | fromint(T) | subtract(T) | multiply(T) | add(T) | lessthan(T) | sqrt(T),
-       V | dottable(V, T) | { V ?*?(T, V); V ?-?(V,V); void ?{}(V&, zero_t); })
+forall(otype T | fromint(T) | subtract(T) | multiply(T) | add(T) | lessthan(T) | sqrt(T),
+       otype V | dottable(V, T) | { V ?*?(T, V); V ?-?(V,V); void ?{}(V&, zero_t); })
 V refract(V v, V n, T eta) {
     T dotValue = dot(n, v);
@@ -128 +113 @@
 // i is the incident vector
 // ng is the geometric normal of the surface
-forall(T | lessthan(T) | zeroinit(T), V | dottable(V, T) | negate(V))
+forall(otype T | lessthan(T) | zeroinit(T), otype V | dottable(V, T) | negate(V))
 V faceforward(V n, V i, V ng) {
     return dot(ng, i) < (T){0} ? n : -n;