Move from using GaussQuadrature to FastGaussQuadrature.jl

Author: sriharshakandala

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},