format

README.md

@@ -31,14 +31,14 @@ To perform one-dimensional quadrature, we can simply construct an `IntegralProbl
 into the problem as the fourth argument of `IntegralProblem`.
 
 ```julia
-using Integrals 
-f(x, p) = sin(x*p)
+using Integrals
+f(x, p) = sin(x * p)
 p = 1.7
 prob = IntegralProblem(f, -2, 5, p)
 sol = solve(prob, QuadGKJL())
 ```
 
-For basic multidimensional quadrature we can construct and solve a `IntegralProblem`. Since we are using no arguments `p` in this example, we omit the fourth argument of `IntegralProblem` 
+For basic multidimensional quadrature we can construct and solve a `IntegralProblem`. Since we are using no arguments `p` in this example, we omit the fourth argument of `IntegralProblem`
 from above. The lower and upper bounds are now passed as vectors, with the `i`th elements of
 the bounds giving the interval of integration for `x[i]`.
 

ext/IntegralsFastGaussQuadratureExt.jl

@@ -48,4 +48,4 @@ function Integrals.__solvebp_call(prob::IntegralProblem, alg::Integrals.GaussLeg
     err = nothing
     SciMLBase.build_solution(prob, alg, val, err, retcode = ReturnCode.Success)
 end
-end
+end

src/algorithms.jl

@@ -121,4 +121,4 @@ function GaussLegendre(; n = 250, subintervals = 1, nodes = nothing, weights = n
         nodes, weights = gausslegendre(n)
     end
     return GaussLegendre(nodes, weights, subintervals)
-end
+end

test/gaussian_quadrature_tests.jl

@@ -59,36 +59,36 @@ alg = GaussLegendre(subintervals = 7)
 sol = solve(prob, alg)
 @test sol.u ≈ -exp(3) * 3.3 + 3.3 / exp(5) - 40 + cos(5) - cos(3)
 
-f = (x, p) -> exp(-x^2) 
+f = (x, p) -> exp(-x^2)
 prob = IntegralProblem(f, 0.0, Inf)
 alg = GaussLegendre()
 sol = solve(prob, alg)
-@test sol.u ≈ sqrt(π)/2
-alg = GaussLegendre(subintervals=1)
-@test sol.u ≈ sqrt(π)/2
-alg = GaussLegendre(subintervals=17)
-@test sol.u ≈ sqrt(π)/2
+@test sol.u ≈ sqrt(π) / 2
+alg = GaussLegendre(subintervals = 1)
+@test sol.u ≈ sqrt(π) / 2
+alg = GaussLegendre(subintervals = 17)
+@test sol.u ≈ sqrt(π) / 2
 
 prob = IntegralProblem(f, -Inf, Inf)
 alg = GaussLegendre()
 sol = solve(prob, alg)
 @test sol.u ≈ sqrt(π)
-alg = GaussLegendre(subintervals=1)
+alg = GaussLegendre(subintervals = 1)
 @test sol.u ≈ sqrt(π)
-alg = GaussLegendre(subintervals=17)
+alg = GaussLegendre(subintervals = 17)
 @test sol.u ≈ sqrt(π)
 
 prob = IntegralProblem(f, -Inf, 0.0)
 alg = GaussLegendre()
 sol = solve(prob, alg)
-@test sol.u ≈ sqrt(π)/2
-alg = GaussLegendre(subintervals=1)
-@test sol.u ≈ sqrt(π)/2
-alg = GaussLegendre(subintervals=17)
-@test sol.u ≈ sqrt(π)/2
+@test sol.u ≈ sqrt(π) / 2
+alg = GaussLegendre(subintervals = 1)
+@test sol.u ≈ sqrt(π) / 2
+alg = GaussLegendre(subintervals = 17)
+@test sol.u ≈ sqrt(π) / 2
 
 # Make sure broadcasting correctly handles the argument p 
-f = (x, p) -> 1 + x + x^p[1] - cos(x*p[2]) + exp(x)*p[3] 
+f = (x, p) -> 1 + x + x^p[1] - cos(x * p[2]) + exp(x) * p[3]
 p = [0.3, 1.3, -0.5]
 prob = IntegralProblem(f, 2, 6.3, p)
 alg = GaussLegendre()