build based on 790a208

dev/.documenter-siteinfo.json

@@ -1 +1 @@
-{"documenter":{"julia_version":"1.9.4","generation_timestamp":"2023-11-29T01:53:10","documenter_version":"1.1.1"}}
+{"documenter":{"julia_version":"1.9.4","generation_timestamp":"2023-12-06T02:45:52","documenter_version":"1.2.1"}}

dev/compare/index.html

@@ -567,4 +567,4 @@
   2.187066 seconds (8.27 k allocations: 76.855 MiB, 0.36% compilation time)
 
 julia> @time itp(xq, yq; method = Hiyoshi(2));
- 13.762652 seconds (9.26 k allocations: 76.920 MiB, 0.06% compilation time)</code></pre><h1 id="Conclusion"><a class="docs-heading-anchor" href="#Conclusion">Conclusion</a><a id="Conclusion-1"></a><a class="docs-heading-anchor-permalink" href="#Conclusion" title="Permalink"></a></h1><p>Overall, the smooth interpolants have the best performance, with <code>Farin(1)</code> and <code>Hiyoshi(2)</code> typically beating most interpolants. <code>Hiyoshi(2)</code> is much slower than the other interpolants, though, and <code>Farin(1)</code> may be a preferable interpolant if <span>$C^1$</span> continuity at the data sites is sufficient. For generating derivatives, the <code>Direct()</code> seems to beat the results with the <code>Iterative()</code> method in most situations.</p></article><nav class="docs-footer"><a class="docs-footer-prevpage" href="../swiss/">« Switzerland Elevation Data</a><a class="docs-footer-nextpage" href="../interpolation_math/">Interpolation Details »</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="auto">Automatic (OS)</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.1.1 on <span class="colophon-date" title="Wednesday 29 November 2023 01:53">Wednesday 29 November 2023</span>. Using Julia version 1.9.4.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>
+ 13.762652 seconds (9.26 k allocations: 76.920 MiB, 0.06% compilation time)</code></pre><h1 id="Conclusion"><a class="docs-heading-anchor" href="#Conclusion">Conclusion</a><a id="Conclusion-1"></a><a class="docs-heading-anchor-permalink" href="#Conclusion" title="Permalink"></a></h1><p>Overall, the smooth interpolants have the best performance, with <code>Farin(1)</code> and <code>Hiyoshi(2)</code> typically beating most interpolants. <code>Hiyoshi(2)</code> is much slower than the other interpolants, though, and <code>Farin(1)</code> may be a preferable interpolant if <span>$C^1$</span> continuity at the data sites is sufficient. For generating derivatives, the <code>Direct()</code> seems to beat the results with the <code>Iterative()</code> method in most situations.</p></article><nav class="docs-footer"><a class="docs-footer-prevpage" href="../swiss/">« Switzerland Elevation Data</a><a class="docs-footer-nextpage" href="../interpolation_math/">Interpolation Details »</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="auto">Automatic (OS)</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.2.1 on <span class="colophon-date" title="Wednesday 6 December 2023 02:45">Wednesday 6 December 2023</span>. Using Julia version 1.9.4.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>

dev/differentiation/index.html

@@ -48,7 +48,7 @@
 tri = triangulate([x'; y'])
 vorn = voronoi(tri)
 
-fig = Figure(fontsize=36, resolution=(1800, 600))
+fig = Figure(fontsize=36, size=(1800, 600))
 ax = Axis(fig[1, 1], xlabel="x", ylabel="y", width=600, height=600, title="(a): Data and triangulation", titlealign=:left)
 scatter!(ax, x, y, color=:black, markersize=9)
 triplot!(ax, tri, color=:black, linewidth=2, show_convex_hull=false)
@@ -83,7 +83,7 @@
 function plot_gradients(∇g, tri, f′, x, y)
     ∇g1 = first.(∇g)
     ∇g2 = last.(∇g)
-    fig = Figure(fontsize=36, resolution=(2400, 600))
+    fig = Figure(fontsize=36, size=(2400, 600))
     plot_f2(fig, x, y, ∇g1, "(a): ∂f̂/∂x", 1, tri, -3.5:0.5:3.0, true, "∂f̂/∂x")
     plot_f2(fig, x, y, ∇g2, "(b): ∂f̂/∂y", 3, tri, -3.5:0.5:3.0, true, "∂f̂/∂y")
     plot_f2(fig, x, y, getindex.(f′.(x, y), 1), "(c): ∂f/∂x", 2, tri, -3.5:0.5:3.0, true, "∂f/∂x")
@@ -111,7 +111,7 @@
     H₁₂ = getindex.(H, 3)
     H₂₂ = getindex.(H, 2)
 
-    fig = Figure(fontsize=36, resolution=(2400, 600))
+    fig = Figure(fontsize=36, size=(2400, 600))
     plot_f2(fig, x, y, H₁₁, "(a): ∂²f̂/∂x²", 1, tri, -35:5:30, true, "∂²f̂/∂x²")
     plot_f2(fig, x, y, H₂₂, "(c): ∂²f̂/∂y²", 3, tri, -35:5:30, true, "∂²f̂/∂y²")
     plot_f2(fig, x, y, H₁₂, "(e): ∂²f̂/∂x∂y", 5, tri, -35:5:30, true, "∂²f̂/∂x∂y")
@@ -164,7 +164,7 @@
     ∇g1 = first.(∇g)
     ∇g2 = last.(∇g)
     ∇f = [f′(x, y) for x in xg, y in yg]
-    fig = Figure(fontsize=36, resolution=(2400, 600))
+    fig = Figure(fontsize=36, size=(2400, 600))
     plot_f2(fig, xg, yg, ∇g1, "(a): ∂f̂/∂x", 1, -3.5:0.5:3.0, true, "∂f̂/∂x")
     plot_f2(fig, xg, yg, ∇g2, "(b): ∂f̂/∂y", 3, -3.5:0.5:3.0, true, "∂f̂/∂y")
     plot_f2(fig, xg, yg, first.(∇f), "(c): ∂f/∂x", 2, -3.5:0.5:3.0, true, "∂f/∂x")
@@ -191,7 +191,7 @@
     H₁₂ = getindex.(H, 3)
     H₂₂ = getindex.(H, 2)
     Hf = [f′′(x, y) for x in xg, y in yg]
-    fig = Figure(fontsize=36, resolution=(2400, 600))
+    fig = Figure(fontsize=36, size=(2400, 600))
     plot_f2(fig, xg, yg, H₁₁, "(a): ∂²f̂/∂x²", 1, -35:5:30, true, "∂²f̂/∂x²")
     plot_f2(fig, xg, yg, H₂₂, "(c): ∂²f̂/∂y²", 3, -35:5:30, true, "∂²f̂/∂y²")
     plot_f2(fig, xg, yg, H₁₂, "(e): ∂²f̂/∂x∂y", 5, -35:5:30, true, "∂²f̂/∂x∂y")
@@ -248,4 +248,4 @@
 7.479964687679311
 
 julia> εH_nohull = rrmserr(f′′.(_x, _y), to_mat.(Hg), ∂, _x, _y)
-38.884740966379056</code></pre><p>The errors are smaller, though not by much.</p></article><nav class="docs-footer"><a class="docs-footer-prevpage" href="../interpolation/">« Interpolation</a><a class="docs-footer-nextpage" href="../swiss/">Switzerland Elevation Data »</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="auto">Automatic (OS)</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.1.1 on <span class="colophon-date" title="Wednesday 29 November 2023 01:53">Wednesday 29 November 2023</span>. Using Julia version 1.9.4.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>
+38.884740966379056</code></pre><p>The errors are smaller, though not by much.</p></article><nav class="docs-footer"><a class="docs-footer-prevpage" href="../interpolation/">« Interpolation</a><a class="docs-footer-nextpage" href="../swiss/">Switzerland Elevation Data »</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="auto">Automatic (OS)</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.2.1 on <span class="colophon-date" title="Wednesday 6 December 2023 02:45">Wednesday 6 December 2023</span>. Using Julia version 1.9.4.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>

dev/differentiation_math/index.html

@@ -21,4 +21,4 @@
 \overline{\boldsymbol g}_1 &= \gamma_i^\prime g_{i1}, \\
 \overline{\boldsymbol g}_2 &= \gamma_i^\prime g_{i2}, \\
 \boldsymbol{\bm\theta} &= \begin{bmatrix} \frac{\partial f(\boldsymbol x_0)}{\partial x} & \frac{\partial f(\boldsymbol x_0)}{\partial y} & \frac{\partial^2 f(\boldsymbol x_0)}{\partial x^2} & \frac{\partial f(\boldsymbol x_0)}{\partial y^2} & \frac{\partial f(\boldsymbol x_0)}{\partial x\partial y} \end{bmatrix}^T.
-\end{align*}\]</p><p>To solve this linear system, let</p><p class="math-container">\[\boldsymbol D = \begin{bmatrix} \overline{\boldsymbol A} \\ \overline{\boldsymbol B} \\ \overline{\boldsymbol C} \end{bmatrix}, \quad \boldsymbol c = \begin{bmatrix} \overline{\boldsymbol w} \\ \overline{\boldsymbol g}_1 \\ \overline{\boldsymbol g}_2 \end{bmatrix},\]</p><p>so that <span>$\boldsymbol D^T\boldsymbol D\boldsymbol\theta = \boldsymbol D^T\boldsymbol c$</span>. These are just the normal equations for <span>$\boldsymbol D\boldsymbol \theta = \boldsymbol c$</span>, thus we can estimate the gradients and Hessians by simply solving <span>$\boldsymbol D\boldsymbol \theta = \boldsymbol c$</span>.</p><h1 id="Generation-Away-from-the-Data-Sites"><a class="docs-heading-anchor" href="#Generation-Away-from-the-Data-Sites">Generation Away from the Data Sites</a><a id="Generation-Away-from-the-Data-Sites-1"></a><a class="docs-heading-anchor-permalink" href="#Generation-Away-from-the-Data-Sites" title="Permalink"></a></h1><p>It is possible to extend these ideas so that we can approximate the derivative at any point <span>$\boldsymbol x_0 \in \mathcal C(\boldsymbol X)$</span>. Using the associated interpolant, simply approximate <span>$z_0$</span> with the value of the interpolant at <span>$\boldsymbol x_0$</span>, and then replace <span>$W_i$</span> by <span>$\lambda_i/\|\boldsymbol x_i-\boldsymbol x_0\|$</span>, where <span>$\lambda_i$</span> is the Sibson coordinate at <span>$\boldsymbol x_i$</span> relative to <span>$\boldsymbol x_0$</span>. If using a direct approach to approximate gradients and Hessians, Sibson coordinates cannot be used (because you can&#39;t extend the weights out to <span>$N_0^2$</span>) and so <span>$W_i$</span> remains as is in that case. Note that the <span>$N_0$</span> neighbourhoods are now the sets of natural neighbours.</p></article><nav class="docs-footer"><a class="docs-footer-prevpage" href="../interpolation_math/">« Interpolation Details</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="auto">Automatic (OS)</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.1.1 on <span class="colophon-date" title="Wednesday 29 November 2023 01:53">Wednesday 29 November 2023</span>. Using Julia version 1.9.4.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>
+\end{align*}\]</p><p>To solve this linear system, let</p><p class="math-container">\[\boldsymbol D = \begin{bmatrix} \overline{\boldsymbol A} \\ \overline{\boldsymbol B} \\ \overline{\boldsymbol C} \end{bmatrix}, \quad \boldsymbol c = \begin{bmatrix} \overline{\boldsymbol w} \\ \overline{\boldsymbol g}_1 \\ \overline{\boldsymbol g}_2 \end{bmatrix},\]</p><p>so that <span>$\boldsymbol D^T\boldsymbol D\boldsymbol\theta = \boldsymbol D^T\boldsymbol c$</span>. These are just the normal equations for <span>$\boldsymbol D\boldsymbol \theta = \boldsymbol c$</span>, thus we can estimate the gradients and Hessians by simply solving <span>$\boldsymbol D\boldsymbol \theta = \boldsymbol c$</span>.</p><h1 id="Generation-Away-from-the-Data-Sites"><a class="docs-heading-anchor" href="#Generation-Away-from-the-Data-Sites">Generation Away from the Data Sites</a><a id="Generation-Away-from-the-Data-Sites-1"></a><a class="docs-heading-anchor-permalink" href="#Generation-Away-from-the-Data-Sites" title="Permalink"></a></h1><p>It is possible to extend these ideas so that we can approximate the derivative at any point <span>$\boldsymbol x_0 \in \mathcal C(\boldsymbol X)$</span>. Using the associated interpolant, simply approximate <span>$z_0$</span> with the value of the interpolant at <span>$\boldsymbol x_0$</span>, and then replace <span>$W_i$</span> by <span>$\lambda_i/\|\boldsymbol x_i-\boldsymbol x_0\|$</span>, where <span>$\lambda_i$</span> is the Sibson coordinate at <span>$\boldsymbol x_i$</span> relative to <span>$\boldsymbol x_0$</span>. If using a direct approach to approximate gradients and Hessians, Sibson coordinates cannot be used (because you can&#39;t extend the weights out to <span>$N_0^2$</span>) and so <span>$W_i$</span> remains as is in that case. Note that the <span>$N_0$</span> neighbourhoods are now the sets of natural neighbours.</p></article><nav class="docs-footer"><a class="docs-footer-prevpage" href="../interpolation_math/">« Interpolation Details</a><div class="flexbox-break"></div><p class="footer-message">Powered by <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> and the <a href="https://julialang.org/">Julia Programming Language</a>.</p></nav></div><div class="modal" id="documenter-settings"><div class="modal-background"></div><div class="modal-card"><header class="modal-card-head"><p class="modal-card-title">Settings</p><button class="delete"></button></header><section class="modal-card-body"><p><label class="label">Theme</label><div class="select"><select id="documenter-themepicker"><option value="documenter-light">documenter-light</option><option value="documenter-dark">documenter-dark</option><option value="auto">Automatic (OS)</option></select></div></p><hr/><p>This document was generated with <a href="https://github.com/JuliaDocs/Documenter.jl">Documenter.jl</a> version 1.2.1 on <span class="colophon-date" title="Wednesday 6 December 2023 02:45">Wednesday 6 December 2023</span>. Using Julia version 1.9.4.</p></section><footer class="modal-card-foot"></footer></div></div></div></body></html>