I have the following function: function wind_factor(height) log1 = log(height / 0.0002) / log(6.0 / 0.0002) exp1 = exp(0.08163 * log(height/6.0)) return log1 + (log1 - exp1) end The parameter height is in the range 6…1000 . It is a bit slow for my purpose (about 23ns on my PC). Is th…

wind_factor(x)=evalpoly(1/x, (1.5949457323488772, 1518.919208108258, 104930.09322893005, -211416.6534637618, 1.3390365545830491e6, -5.705301564190215e6, 1.0384434802323304e7))/evalpoly(1/x, (1.0, 1092.3211847159544, 91063.22360243625)) here’s the absolute error: (and it runs in 3ns) [image] (foun…

How to approximate a function

General Usage

Oscar_Smith May 6, 2022, 1:17am 21

Since your range is (6,1000), the transform 1/x significantly reduces the domain you are approximating over.

Topic		Replies	Views
Is there a smarter way to evaluate an float power? Performance	27	1851	February 15, 2023
BetterExp 0.1.0 Package Announcements math	30	2894	September 9, 2020
A faster pow(x,1/12) function available in Fast12thRoot.jl Numerics	12	2649	June 3, 2019
Slow exp()? Numerics fast-math	12	2248	August 17, 2018
Speeding up evaluation of small exponential-like series Performance	19	591	May 17, 2024

How to approximate a function

Related topics