From 5cce67439b488cc62c97ee5ec39251ad0aeff6d2 Mon Sep 17 00:00:00 2001 From: Kyle Bayes Date: Wed, 3 Jul 2024 12:37:48 -0700 Subject: [PATCH] Clean up user_util. PiperOrigin-RevId: 649166858 Change-Id: I79423dda628fb4bf3ceae438dd8dc5b96d69d2a9 --- src/user/user_util.cc | 56 +++++++++++++++++++------------------------ src/user/user_util.h | 32 ++++++++++++------------- 2 files changed, 41 insertions(+), 47 deletions(-) diff --git a/src/user/user_util.cc b/src/user/user_util.cc index 4d640091..ed3a7158 100644 --- a/src/user/user_util.cc +++ b/src/user/user_util.cc @@ -17,28 +17,21 @@ #include #include #include -#include #include #include #include -#include #include #include "engine/engine_crossplatform.h" -using std::isnan; -using std::string; -using std::numeric_limits; - - // check if numeric variable is defined -bool mjuu_defined(const double num) { - return !isnan(num); +bool mjuu_defined(double num) { + return !std::isnan(num); } // compute address of M[g1][g2] where M is triangular n-by-n -int mjuu_matadr(int g1, int g2, const int n) { +int mjuu_matadr(int g1, int g2, int n) { if (g1<0 || g2<0 || g1>=n || g2>=n) { return -1; } @@ -54,13 +47,13 @@ int mjuu_matadr(int g1, int g2, const int n) { // set 4D vector -void mjuu_setvec(double* dest, const double x, const double y, const double z, const double w) { +void mjuu_setvec(double* dest, double x, double y, double z, double w) { dest[0] = x; dest[1] = y; dest[2] = z; dest[3] = w; } -void mjuu_setvec(float* dest, const double x, const double y, const double z, const double w) { +void mjuu_setvec(float* dest, double x, double y, double z, double w) { dest[0] = (float)x; dest[1] = (float)y; dest[2] = (float)z; @@ -69,12 +62,12 @@ void mjuu_setvec(float* dest, const double x, const double y, const double z, co // set 3D vector -void mjuu_setvec(double* dest, const double x, const double y, const double z) { +void mjuu_setvec(double* dest, double x, double y, double z) { dest[0] = x; dest[1] = y; dest[2] = z; } -void mjuu_setvec(float* dest, const double x, const double y, const double z) { +void mjuu_setvec(float* dest, double x, double y, double z) { dest[0] = (float)x; dest[1] = (float)y; dest[2] = (float)z; @@ -82,28 +75,28 @@ void mjuu_setvec(float* dest, const double x, const double y, const double z) { // set 2D vector -void mjuu_setvec(double* dest, const double x, const double y) { +void mjuu_setvec(double* dest, double x, double y) { dest[0] = x; dest[1] = y; } // add to double array -void mjuu_addtovec(double* dest, const double* src, const int n) { +void mjuu_addtovec(double* dest, const double* src, int n) { for (int i=0; i -const double mjEPS = 1E-14; // minimum value in various calculations -const double mjMINMASS = 1E-6; // minimum mass allowed +const double mjEPS = 1E-14; // minimum value in various calculations +const double mjMINMASS = 1E-6; // minimum mass allowed // check if numeric variable is defined: !_isnan(num) -bool mjuu_defined(const double num); +bool mjuu_defined(double num); // compute linear address of M[g1][g2] where M is triangular n-by-n // return -1 if inputs are invalid -int mjuu_matadr(int g1, int g2, const int n); +int mjuu_matadr(int g1, int g2, int n); // set 4D vector -void mjuu_setvec(double* dest, const double x, const double y, const double z, const double w); -void mjuu_setvec(float* dest, const double x, const double y, const double z, const double w); +void mjuu_setvec(double* dest, double x, double y, double z, double w); +void mjuu_setvec(float* dest, double x, double y, double z, double w); // set 3D vector -void mjuu_setvec(double* dest, const double x, const double y, const double z); -void mjuu_setvec(float* dest, const double x, const double y, const double z); +void mjuu_setvec(double* dest, double x, double y, double z); +void mjuu_setvec(float* dest, double x, double y, double z); // set 2D vector -void mjuu_setvec(double* dest, const double x, const double y); +void mjuu_setvec(double* dest, double x, double y); // copy real-valued vector template @@ -50,18 +50,18 @@ void mjuu_copyvec(T1* dest, const T2* src, int n) { } // add to double array -void mjuu_addtovec(double* dest, const double* src, const int n); +void mjuu_addtovec(double* dest, const double* src, int n); // zero array -void mjuu_zerovec(double* dest, const int n); +void mjuu_zerovec(double* dest, int n); // zero float array -void mjuu_zerovec(float* dest, const int n); +void mjuu_zerovec(float* dest, int n); // dot-product in 3D double mjuu_dot3(const double* a, const double* b); -// distance beween 3D points +// distance between 3D points double mjuu_dist3(const double* a, const double* b); // L1 norm between vectors @@ -69,8 +69,8 @@ double mjuu_L1(const double* a, const double* b, int n); // normalize vector to unit length, return previous length // if norm(vec)