| 218 | } |
| 219 | |
| 220 | friend T angle(const vec3<T>& v1, const vec3<T>& v2) |
| 221 | { |
| 222 | T numerator = dot(v1, v2); |
| 223 | T denominator = v1.length() * v2.length(); |
| 224 | // assert(!isZero<T>(denominator)); |
| 225 | return std::acos(numerator/denominator); |
| 226 | } |
| 227 | |
| 228 | friend T distance(const vec3<T>& pt1, const vec3<T>& pt2) |
| 229 | { |