Adaptive Plot

f18cb960-2268-4668-b695-e831c8372bcd 83 4 5 false false 0 decimal

&[DATE] &[TIME] - &[FILENAME]

Adaptive Plot

z w v CAngle a w z - : b v w - : 12 mat a Re i a Im * - (b * Re a abs b abs * / acos Re eval 0 try 21 line :

1 f a b APlot ξ f a b APlotAngle APlotDist MaxIters APlotNG / APlotNG APlot :

1 f a b α p Μ Ν APlot δ b a - Μ / : k 1 Μ 1 + (range T k el a k 1 - (δ * + : for Z T f vectorize eval : d p Z max Z min - abs * eval : ng 1 Ν range N T length : k 1 N 2 - (range A k el Z k el Z k 1 + el Z k 2 + el CAngle : for k 1 N 1 - (range D k el Z k 1 + el Z k el - abs : for A max α ≤ (D max d ≤ (& eval break N T length : k 0 : t 0 : κ 0 : τ 0 : 15 mat n 0 N 1 - 2 / 1 - (range t k k 1 + : el T 2 n * 1 + el : θ Z 2 n * 1 + el Z 2 n * 2 + el Z 2 n * 3 + el CAngle : θ α > (Z 2 n * 2 + el Z 2 n * 1 + el - abs d > (| t k k 1 + : el 0.5 T 2 n * 1 + el T 2 n * 2 + el + (* : if t k k 1 + : el T 2 n * 2 + el : θ α > (Z 2 n * 2 + el Z 2 n * 3 + el - abs d > (| t k k 1 + : el 0.5 T 2 n * 2 + el T 2 n * 3 + el + (* : if 51 line for N 2 mod 0 ≡ t k k 1 + : el T N 1 - el : θ α > (Z N el Z N 1 - el - abs d > (| t k k 1 + : el 0.5 T N 1 - el T N el + (* : if 21 line if T t T N el stack eval : Z T f vectorize eval : 51 line if 41 line for Z 61 line :

CSspline

Index in X for the interval containing x, assuming that  
X is in ascending order. x\032E\ X\032E\ CS.L n\0324\ X\032E\ length eval : x\032E\ X\032E\ 1 el X\032E\ n\0324\ el lele x\032E\ 2 X\032E\ 1 el * X\032E\ 2 el - X\032E\ 2 el ltlt 1 x\032E\ X\032E\ n\0324\ 1 - el 2 X\032E\ n\0324\ el * X\032E\ n\0324\ 1 - el - lelt n\0324\ 1 - n\0324\ x\032E\ X\032E\ ≤ vectorize 11 findrows rows - if if CS.L: Value outside of range. error if 21 line :

t\032E\ h\032E\ y\032E\.0 y\032E\.1 y\032E\''.0 y\032E\''.1 CS.1 t\032E\ t\032E\ 1 - 2 / y\032E\''.0 t\032E\ 1 + 3 / y\032E\''.1 y\032E\''.0 - (* + (* h\032E\ * h\032E\ * y\032E\.1 + y\032E\.0 - (* y\032E\.0 + eval :

Cubic Spline interpolation function. M\032E\ x\032E\ CS X\032E\ Y\032E\ V\032E\ T\032E\ 14 mat M\032E\ eval : k\0324\ x\032E\ X\032E\ CS.L : h\0324\ X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - : x\032E\ X\032E\ k\0324\ el - h\0324\ / h\0324\ Y\032E\ k\0324\ el Y\032E\ k\0324\ 1 + el V\032E\ k\0324\ el V\032E\ k\0324\ 1 + el CS.1 41 line :

Cubic Spline interpolation vector of coefficients. X\032E\ Y\032E\ CS.V n\0324\ X\032E\ length : c\0324\ 1 el X\032E\ 3 el X\032E\ 1 el - X\032E\ 3 el X\032E\ 2 el - / - : d\0324\ 1 el 0 : b\0324\ X\032E\ 3 el X\032E\ 1 el - 3 / c\0324\ 1 el X\032E\ 2 el X\032E\ 1 el - 6 / * - : c\0324\ 2 el X\032E\ 2 el X\032E\ 1 el - 6 b\0324\ * / X\032E\ 3 el X\032E\ 2 el - X\032E\ 2 el X\032E\ 1 el - / X\032E\ 2 el X\032E\ 1 el - X\032E\ 3 el X\032E\ 2 el - / - (* : d\0324\ 2 el 1 b\0324\ / Y\032E\ 3 el Y\032E\ 2 el - X\032E\ 3 el X\032E\ 2 el - / Y\032E\ 2 el Y\032E\ 1 el - X\032E\ 2 el X\032E\ 1 el - / - (* : n\0324\ 3 > k\0324\ 3 n\0324\ 1 - (range b\0324\ X\032E\ k\0324\ 1 + el X\032E\ k\0324\ 1 - el - 3 / c\0324\ k\0324\ 1 - el X\032E\ k\0324\ el X\032E\ k\0324\ 1 - el - 6 / * - : c\0324\ k\0324\ el X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - 6 b\0324\ * / : d\0324\ k\0324\ el 1 b\0324\ / Y\032E\ k\0324\ 1 + el Y\032E\ k\0324\ el - X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - / Y\032E\ k\0324\ el Y\032E\ k\0324\ 1 - el - X\032E\ k\0324\ el X\032E\ k\0324\ 1 - el - / - d\0324\ k\0324\ 1 - el X\032E\ k\0324\ el X\032E\ k\0324\ 1 - el - 6 / * - (* : 31 line for if a\0324\ X\032E\ n\0324\ el X\032E\ n\0324\ 2 - el - (- X\032E\ n\0324\ el X\032E\ n\0324\ 1 - el - (c\0324\ n\0324\ 2 - el * - : b\0324\ X\032E\ n\0324\ 1 - el X\032E\ n\0324\ 2 - el - c\0324\ n\0324\ 1 - el a\0324\ * - : d\0324\ n\0324\ el X\032E\ n\0324\ el X\032E\ n\0324\ 1 - el - (d\0324\ n\0324\ 2 - el * - d\0324\ n\0324\ 1 - el a\0324\ * - b\0324\ / : k\0324\ n\0324\ 1 - (1 range d\0324\ k\0324\ el d\0324\ k\0324\ el c\0324\ k\0324\ el d\0324\ k\0324\ 1 + el * - : for d\0324\ 1 el d\0324\ 1 el X\032E\ 2 el X\032E\ 1 el - X\032E\ 3 el X\032E\ 2 el - / d\0324\ 3 el * - : X\032E\ Y\032E\ d\0324\ 032 stack 14 mat eval 131 line :

Adaptive Cubic Spline interpolation. 1 f a b CS.A x g x i x f * + : Z x g a b APlot : z Z 1 el : k 2 Z length range Z k 1 - el Re Z k el Re ≡ k Z k el 12 mat trace z z Z k el stack eval : if for z Re z Im CS.V eval 51 line :

Setup

MaxIters 100 :

APlotAngle 15 ° * :

APlotDist 10 % * :

APlotNG 4 :

For the cloud:

a b c b a ≤ (a c ≤ (& 11 line : lele

a b c ltlt b a < (a c < (& 11 line :

a b c b a ≤ (a c < (& 11 line : lelt

CSpline Test

CS Interpolation

X 125 range : x f x 5^4 + ln sin : 21 sys

M X X f vectorize CS.V : x p M x CS : 21 sys

x f x p 2 1 sys x f M 1 el M 2 el . 12 green augment 2 1 sys

Adaptive CS interpolation

a 1 : b 25 : 12 mat x f x 5^4 + ln sin : 21 sys

M x f a b CS.A : x p M x CS : 21 sys

x f x p 2 1 sys x f M 1 el M 2 el . 12 green augment 2 1 sys

X 1015 -^m * (12 m * (range : x f x x 2 m * / cos * : 21 sys

M X X f vectorize CS.V : x p M x CS : 21 sys

CS interpolation with units: fails if starts at 0*'m

x m * f x m * p 2 1 sys x m * f M 1 el M 2 el . 12 green augment 2 1 sys

Adaptive CS interp. Can start at 0*'m

a 0 m * : b 12 m * : 12 mat x f x x 2 m * / cos * : 21 sys

M x f a b CS.A : x p M x CS : 21 sys

x m * f x m * p 2 1 sys x m * f M 1 el M 2 el . 12 green augment 2 1 sys Adaptive Plot Test

Adaptive Plot

Cartesians:

f\0028\x\0029\\0020\\002D\ x i x f * + >

Polars:

ρ\0028\θ\0029\\0020\\002D\ θ ρ e i θ *^* >

Parametrics:

z\0028\t\0029\\0020\\002D\ t x i t y * + >

t.o 1 time :

t z 25 t * cos * i 7 t * sin * + :

Z t z 02 π * APlot :

Z length 247

Z Re Z Im . 8 green augment Z Re Z Im augment 2 1 sys

x f x i x * 18 x * cos * + :

Z x f 05 APlot :

Z length 276

Z Re Z Im . 8 green augment Z Re Z Im augment 2 1 sys

θ ρ 27 θ * cos * e i θ *^* :

Z θ ρ 02 π * APlot :

Z length 257

1 time t.o - 32.514 s *

Z Re Z Im . 8 green augment Z Re Z Im augment 2 1 sys

Alvaro