Gabriel Weymouth

Looks like the dependencies need to be updated so newer versions of StaticArrays are allowed.

TwoPhase is unstable

1

Migrated from issue #10. The attached LilyPad.pde has a circle moving up and down in a two phase flow. If QUICK advection is used, the flow is fine. If the...

Added functions to get distances to primatives and meshes

14

Created new file distances.jl which computes the absolute_distance and the signed_distance to triangle meshes, HyperSpheres and HyperRects. These functions have been exported. The sphere and rect signed functions are straightforward,...

Large negative real number arguments return NaN for expintx

1

`expintx(complex(-x))` returns `NaN` for x more than around 700. It's almost like this is being handled by `-exp(-x)*expinti(x)` and we're getting an overflow error.

Problem with sum over OffsetArray

1

It's me again. That guy who keeps finding bugs when using `sum`. In [Polyester](https://github.com/JuliaSIMD/Polyester.jl/issues/108), in [KernelAbstractions](https://discourse.julialang.org/t/gpu-sum-closure-throwing-an-error/96658), and now in OffsetArrays. Sigh... Here is a case which slows down a loop...

Undefined nested iterating variable

7

This simple example throws an error ``` using Polyester x = collect(1:12) y = zeros(6) @batch for i in eachindex(y) y[i] = sum(x[j] for j in 2i-oneunit(i):2i) end y ```...

Extensive use of abs for complex tolerances

2

I just profiled `expintx(1,z::Complex)` and found the majority of the time is spent checking the convergence using `abs(stuff) < ϵ`. Would there be a problem with simply squaring all of...

I'm trying to recreate a 3D Chebyshev surrogate model from an old paper [Newman 1987](https://watermark.silverchair.com/sname-jsr-1987-31-2-79.pdf?token=AQECAHi208BE49Ooan9kkhW_Ercy7Dm3ZL_9Cf3qfKAc485ysgAAA8IwggO-BgkqhkiG9w0BBwagggOvMIIDqwIBADCCA6QGCSqGSIb3DQEHATAeBglghkgBZQMEAS4wEQQMtepSki4L2l02qssrAgEQgIIDdWtmV6N4FvoJOCVWjtQ38RGCwH9PVY4Kp1lqgd-jshuxCRGMYFUw_BgfmqlYRPD4ulrfLlEvjYLmAMwp9LUBDS1mku7tQACOwcXtjZzsfkR1vpn9C4NOlqRA7gHg4dk5XNKHET9uVeNnEgq5Dbp8W81F6RZoGZCeHfz1JfwmfG_eRJycNw80NJ7M1xXVOoFclip21u9vD0jAskiTpj1XevbU_kDq0pByifAPlS0HNgZ9Y7Js0DFs8caQ_DK9kyprzDXHEKClk23rCBnzkSshyVompRtOfKEhKLeJeVBTDCAw86JZ1FPPB4oR4GjcpH7tO4tEG_U0YqWlSxOzJoQGYTDinwIDHHHfOOONJa7ku-bGbIFsssrWtubqfh6UnEmem4vQRiAWKirzaLWBFVG2rjFGXajsv2AMt7-2Okwg88lusacbvcMuZoXpYtAO0IQejsN7qDbm5d6vuJNEl31-ncAV-r3ROAhPNWcDbAQOvKnr_YcNd6yvkCQ5wQkhiSK_frDfT1cmnfWwbY4oRqsAh921ZKSIOQdfe9B4uYtyMQ_1crRJ8zgwmzZWeZzK5nzR7FxLJUNxd7PmAS_O1xSfxxGTeDuY1WtWTCE6-IBFkRPvkCRY0BW_swBIKD3_ZdVWx1cMTNF_jm1VXCrzsxSZCpDms_GIVbpbOvU93V3ZTh0SZB7SZTQdHz-yJWH_Hyp9vNdKUNAx8mWls96uCLmWpg-HddtGMq0YnUF43m0JsvBkWIH9hNe19u34Hr1UdTdouQH_MnmCjFA37CsP4mBhjFT0WGgZpi_th9NxnV5JPBjDTQSi6B68YpUUq924qvK3lw58EQKosZGKWYArzGohkJsSIHmRr0JsTLAVoyKxp330fHyjJgsyOkmV0-X9G2x9bUQWcrTClMEh9a6vWc4pM2X4WaSNggk15UAqOX60WROVsDJpL2seoaFzta22AoWcZxLXEsWOeKgbIKLO_h2q09ZChlYq-DNqdqOBzjus49nDZWUW4SgGcPcuTSG56P-QvRlqBuVNxcUkxpMmIbgiiVpBLaj-N5ckIKFPQkLa2Ckubxm-mMtdPtBDm7Zkk_zAQYjLD3_Hq6LfSUMB8HhW0Ve8xeHU4JsOTV20LX49s5wQHwkss_HVg3SnmSKQiKwsYVYh2fIBPtorwbb-yn0OlaKxRj5OVw). ```julia using SpecialFunctions,QuadGK using Base.MathConstants: γ function bruteN(x,y,z;withS=true) ζ(t) = (z*sqrt(1-t^2)+y*t+im*abs(x))*sqrt(1-t^2) Ni(t) = imag(expintx(ζ(t))+log(ζ(t))+γ) 2/π*quadgk(Ni,-1,0,1,atol=eps())[1]+ifelse(withS,S(x,y,z),0)...

derived types

1

There is nothing in here about derived types. This is a big feature of modern Fortran and I'm surprised it isn't included.

Throw exception or warning when user puts in a bad geometry: 1. Exception: Geometry is "inside out" (panels running clockwise) when computing `properties`. 2. Warning: Using very few panels (maybe