readme

Author: mscroggs

README.md

@@ -1,3 +1,4 @@
 # QuadratureRulesDotOrg
 
-[![Build Status](https://github.com/mscroggs/QuadratureRulesDotOrg.jl/actions/workflows/CI.yml/badge.svg?branch=main)](https://github.com/mscroggs/QuadratureRulesDotOrg.jl/actions/workflows/CI.yml?query=branch%3Amain)
+The code in this repo is generated by [the quadraturerules repo](https://github.com/quadraturerules/quadraturerules).
+If you want to open an issue, please open it [on the quadrature rules repo](https://github.com/quadraturerules/quadraturerules/issues).

src/QuadratureRulesDotOrg.jl

@@ -1,5 +1,65 @@
 module QuadratureRulesDotOrg
 
-# Write your package code here.
+include("domains.jl")
+include("gauss_legendre.jl")
+include("gauss_lobatto_legendre.jl")
+include("hammer_marlowe_stroud.jl")
+include("sauter_schwab.jl")
+include("xiao_gimbutas.jl")
+include("centroid_quadrature.jl")
+include("vertex_quadrature.jl")
+
+@enum QuadratureRule begin
+    QR_GaussLegendre = 1
+    QR_GaussLobattoLegendre = 3
+    QR_HammerMarloweStroud = 6
+    QR_SauterSchwab = 7
+    QR_XiaoGimbutas = 2
+    QR_CentroidQuadrature = 5
+    QR_VertexQuadrature = 4
+end
+
+"Get a quadrature rule for a single integral."
+function single_integral_quadrature(
+    rtype::QuadratureRule,
+    domain::Domain,
+    order::Integer,
+)
+    if rtype == QR_GaussLegendre
+        return gauss_legendre(domain, order)
+    end
+    if rtype == QR_GaussLobattoLegendre
+        return gauss_lobatto_legendre(domain, order)
+    end
+    if rtype == QR_HammerMarloweStroud
+        return hammer_marlowe_stroud(domain, order)
+    end
+    if rtype == QR_XiaoGimbutas
+        return xiao_gimbutas(domain, order)
+    end
+    if rtype == QR_CentroidQuadrature
+        return centroid_quadrature(domain, order)
+    end
+    if rtype == QR_VertexQuadrature
+        return vertex_quadrature(domain, order)
+    end
+    throw("Unsupported rule for single integral: $rtype")
+end
+
+"Get a quadrature rule for a double integral."
+function double_integral_quadrature(
+    rtype::QuadratureRule,
+    domain::Domain,
+    order::Integer,
+)
+    if rtype == QR_SauterSchwab
+        return sauter_schwab(domain, order)
+    end
+    throw("Unsupported rule for single integral: $rtype")
+end
+
+export QuadratureRule
+export Domain
+export single_integral_quadrature
 
 end

test/runtests.jl

@@ -2,5 +2,10 @@ using QuadratureRulesDotOrg
 using Test
 
 @testset "QuadratureRulesDotOrg.jl" begin
-    # Write your tests here.
+    points, weights = QuadratureRulesDotOrg.single_integral_quadrature(
+        QuadratureRulesDotOrg.QR_GaussLobattoLegendre,
+        QuadratureRulesDotOrg.Domain_Interval,
+        4
+    )
+    @test isapprox(sum(weights), 1.0)
 end