| 109 | } |
| 110 | |
| 111 | friend T angle(const vec4<T>& v1, const vec4<T>& v2) |
| 112 | { |
| 113 | T numerator = dot(v1, v2); |
| 114 | T denominator = v1.length() * v2.length(); |
| 115 | // assert(!isZero<T>(denominator)); |
| 116 | return std::acos(numerator/denominator); |
| 117 | } |
| 118 | |
| 119 | friend T distance(const vec4<T>& pt1, const vec4<T>& pt2) |
| 120 | { |