Move from using GaussQuadrature to FastGaussQuadrature.jl

Author: sriharshakandala

.buildkite/Manifest.toml

@@ -2,7 +2,7 @@
 
 julia_version = "1.10.8"
 manifest_format = "2.0"
-project_hash = "6ab89829ea190189b0319a6f8e22b3515e5283c2"
+project_hash = "08a56a97851a386e63bd497ee2d35065349e6442"
 
 [[deps.ADTypes]]
 git-tree-sha1 = "e2478490447631aedba0823d4d7a80b2cc8cdb32"
@@ -327,10 +327,10 @@ weakdeps = ["CUDA", "MPI"]
     ClimaCommsMPIExt = "MPI"
 
 [[deps.ClimaCore]]
-deps = ["Adapt", "BandedMatrices", "BlockArrays", "ClimaComms", "CubedSphere", "DataStructures", "ForwardDiff", "GaussQuadrature", "GilbertCurves", "HDF5", "InteractiveUtils", "IntervalSets", "KrylovKit", "LinearAlgebra", "MultiBroadcastFusion", "NVTX", "PkgVersion", "RecursiveArrayTools", "RootSolvers", "SparseArrays", "StaticArrays", "Statistics", "UnrolledUtilities"]
+deps = ["Adapt", "BandedMatrices", "BlockArrays", "ClimaComms", "CubedSphere", "DataStructures", "FastGaussQuadrature", "ForwardDiff", "GilbertCurves", "HDF5", "InteractiveUtils", "IntervalSets", "KrylovKit", "LinearAlgebra", "MultiBroadcastFusion", "NVTX", "PkgVersion", "RecursiveArrayTools", "RootSolvers", "SparseArrays", "StaticArrays", "Statistics", "UnrolledUtilities"]
 path = ".."
 uuid = "d414da3d-4745-48bb-8d80-42e94e092884"
-version = "0.14.27"
+version = "0.14.29"
 weakdeps = ["CUDA", "Krylov"]
 
     [deps.ClimaCore.extensions]
@@ -744,6 +744,12 @@ git-tree-sha1 = "acebe244d53ee1b461970f8910c235b259e772ef"
 uuid = "9aa1b823-49e4-5ca5-8b0f-3971ec8bab6a"
 version = "0.3.2"
 
+[[deps.FastGaussQuadrature]]
+deps = ["LinearAlgebra", "SpecialFunctions", "StaticArrays"]
+git-tree-sha1 = "fd923962364b645f3719855c88f7074413a6ad92"
+uuid = "442a2c76-b920-505d-bb47-c5924d526838"
+version = "1.0.2"
+
 [[deps.FastPower]]
 git-tree-sha1 = "58c3431137131577a7c379d00fea00be524338fb"
 uuid = "a4df4552-cc26-4903-aec0-212e50a0e84b"
@@ -896,12 +902,6 @@ git-tree-sha1 = "9cb62849057df859575fc1dda1e91b82f8609709"
 uuid = "d2c73de3-f751-5644-a686-071e5b155ba9"
 version = "0.73.13+0"
 
-[[deps.GaussQuadrature]]
-deps = ["SpecialFunctions"]
-git-tree-sha1 = "eb6f1f48aa994f3018cbd029a17863c6535a266d"
-uuid = "d54b0c1a-921d-58e0-8e36-89d8069c0969"
-version = "0.5.8"
-
 [[deps.Gettext_jll]]
 deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl", "Libiconv_jll", "Pkg", "XML2_jll"]
 git-tree-sha1 = "9b02998aba7bf074d14de89f9d37ca24a1a0b046"

.buildkite/Project.toml

@@ -21,7 +21,6 @@ DiffEqCallbacks = "459566f4-90b8-5000-8ac3-15dfb0a30def"
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
 FastBroadcast = "7034ab61-46d4-4ed7-9d0f-46aef9175898"
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
-GaussQuadrature = "d54b0c1a-921d-58e0-8e36-89d8069c0969"
 GilbertCurves = "88fa7841-ef32-4516-bb70-c6ec135699d9"
 IntervalSets = "8197267c-284f-5f27-9208-e0e47529a953"
 JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"

Project.toml

@@ -10,8 +10,8 @@ BlockArrays = "8e7c35d0-a365-5155-bbbb-fb81a777f24e"
 ClimaComms = "3a4d1b5c-c61d-41fd-a00a-5873ba7a1b0d"
 CubedSphere = "7445602f-e544-4518-8976-18f8e8ae6cdb"
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
+FastGaussQuadrature = "442a2c76-b920-505d-bb47-c5924d526838"
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
-GaussQuadrature = "d54b0c1a-921d-58e0-8e36-89d8069c0969"
 GilbertCurves = "88fa7841-ef32-4516-bb70-c6ec135699d9"
 HDF5 = "f67ccb44-e63f-5c2f-98bd-6dc0ccc4ba2f"
 InteractiveUtils = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
@@ -51,8 +51,8 @@ CubedSphere = "0.2, 0.3"
 DataStructures = "0.18.13"
 Dates = "1"
 FastBroadcast = "0.3.1"
+FastGaussQuadrature = "1.0.2"
 ForwardDiff = "0.10.15"
-GaussQuadrature = "0.5.8"
 GilbertCurves = "0.1"
 HDF5 = "0.16.16, 0.17"
 InteractiveUtils = "1"

src/Quadratures/Quadratures.jl

@@ -1,7 +1,7 @@
 
 module Quadratures
 
-import GaussQuadrature
+import FastGaussQuadrature
 import StaticArrays: SVector, SMatrix, MMatrix
 import LinearAlgebra: Diagonal
 
@@ -53,6 +53,8 @@ function quadrature_points end
     GLL{Nq}()
 
 Gauss-Legendre-Lobatto quadrature using `Nq` quadrature points.
+
+https://mathworld.wolfram.com/LobattoQuadrature.html
 """
 struct GLL{Nq} <: QuadratureStyle{Nq} end
 
@@ -61,8 +63,8 @@ Base.show(io::IO, ::GLL{Nq}) where {Nq} =
 
 unique_degrees_of_freedom(::GLL{Nq}) where {Nq} = Nq - 1
 @generated function quadrature_points(::Type{FT}, ::GLL{Nq}) where {FT, Nq}
-    points, weights = GaussQuadrature.legendre(FT, Nq, GaussQuadrature.both)
-    :($(SVector{Nq}(points)), $(SVector{Nq}(weights)))
+    points, weights = FastGaussQuadrature.gausslobatto(Nq)
+    :($(SVector{Nq, FT}(points)), $(SVector{Nq, FT}(weights)))
 end
 
 """
@@ -77,8 +79,8 @@ Base.show(io::IO, ::GL{Nq}) where {Nq} =
 
 unique_degrees_of_freedom(::GL{Nq}) where {Nq} = Nq
 @generated function quadrature_points(::Type{FT}, ::GL{Nq}) where {FT, Nq}
-    points, weights = GaussQuadrature.legendre(FT, Nq, GaussQuadrature.neither)
-    :($(SVector{Nq}(points)), $(SVector{Nq}(weights)))
+    points, weights = FastGaussQuadrature.gausslegendre(Nq)
+    :($(SVector{Nq, FT}(points)), $(SVector{Nq, FT}(weights)))
 end
 
 """
@@ -220,7 +222,7 @@ end
         quadrature_points(FT, quadfrom())[1],
     )
 end
-
+#=
 """
     V = orthonormal_poly(points, quad)
 
@@ -247,7 +249,7 @@ function spectral_filter_matrix(
     V = orthonormal_poly(points, quad)
     return V * Diagonal(Σ) / V
 end
-
+=#
 function cutoff_filter_matrix(
     ::Type{FT},
     quad::GLL{Nq},