ups

Author: FraCpl

src/JTools.jl

@@ -43,8 +43,8 @@ export plotBox, plotBox!, lines3!, scatter3!, multilines!, plotSphere!, plotSens
 export plotCorrelation, plotCorrelation!
 include("functions/makieTools.jl")
 
-export crossmat, logrange, unwrap!, modd, isMultiple, signum
-export polyfit, polyval, interp1, mul3x1!, mul3x3!
+export crossmat, unwrap!, modd, isMultiple, signum #, logrange
+export polyfit, polyval, interp1, mul3x1!, mul3x3!, evalpoly!
 include("functions/math.jl")
 
 export lsq, lsqWeight

src/functions/math.jl

@@ -1,5 +1,5 @@
 @inline crossmat(v) = [0.0 -v[3]  v[2]; v[3] 0.0 -v[1]; -v[2] v[1] 0.0];
-@inline logrange(a, b, length=100) = exp10.(range(a; stop=b, length=length))
+# @inline logrange(a, b, length=100) = exp10.(range(a; stop=b, length=length))
 
 @inline function unwrap!(x, period=2π)
 	y = convert(eltype(x), period)
@@ -46,6 +46,16 @@ end
     return y
 end
 
+function evalpoly!(y::Vector{T}, t::T, C) where T
+    y .= C[1]
+    tk = T(1)
+    @inbounds for i in 2:lastindex(C)
+        tk *= t
+        y .+= C[i].*tk
+    end
+    return
+end
+
 @inline function interp1(x, y, xi)
     if xi ≥ x[end]
         return y[end]

src/functions/units.jl

@@ -3,13 +3,15 @@ module Units
 m = 1.0
 cm = 1e-2
 mm = 1e-3
+μm = 1e-6
 km = 1e3
 ft = 0.3048
 inch = 0.0254
 
 # Time
 s = 1.0
 ms = 1e-3
+μs = 1e-6
 min = 60.0
 hr = 3600.0
 hrs = 3600.0
@@ -34,6 +36,7 @@ kg = 1.0
 N = 1.0
 kN = 1e3
 mN = 1e-3
+μN = 1e-6
 
 # Angular rate
 rpm = 2π/min
@@ -49,6 +52,7 @@ ng = g*1e-9
 # Frequency
 Hz = 1.0
 mHz = 1e-3
+μHz = 1e-6
 kHz = 1e3
 
 # Constants