definite integration fixed; ISSAC'24 paper added

Author: shahriariravanian

Project.toml

@@ -1,7 +1,7 @@
 name = "SymbolicNumericIntegration"
 uuid = "78aadeae-fbc0-11eb-17b6-c7ec0477ba9e"
 authors = ["Shahriar Iravanian <[email protected]>"]
-version = "1.9.0"
+version = "1.9.1"
 
 [deps]
 DataDrivenDiffEq = "2445eb08-9709-466a-b3fc-47e12bd697a2"

README.md

@@ -38,10 +38,26 @@ integrate(sin(a * x) * cos(b * x), x; symbolic = true, detailed = false)
 
 # Citation
 
-If you use SymbolicNumericIntegration.jl in your research, please cite [this paper](https://arxiv.org/abs/2201.12468):
+If you use SymbolicNumericIntegration.jl in your research, please cite the [arxiv](https://arxiv.org/abs/2201.12468) and/or [ISSAC'24](https://dl.acm.org/doi/abs/10.1145/3666000.3669714) papers:
 
 ```
 
+@inproceedings{10.1145/3666000.3669714,
+author = {Iravanian, Shahriar and Gowda, Shashi and Rackauckas, Christopher},
+title = {Hybrid Symbolic-Numeric and Numerically-Assisted Symbolic Integration},
+year = {2024},
+isbn = {9798400706967},
+publisher = {Association for Computing Machinery},
+address = {New York, NY, USA},
+url = {https://doi.org/10.1145/3666000.3669714},
+doi = {10.1145/3666000.3669714},
+pages = {410–418},
+numpages = {9},
+keywords = {sparse regression, symbolic integration, symbolic-numeric computation},
+location = {Raleigh, NC, USA},
+series = {ISSAC '24}
+}
+
 @article{Iravanian2022,
 author = {Shahriar Iravanian and Carl Julius Martensen and Alessandro Cheli and Shashi Gowda and Anand Jain and Yingbo Ma and Chris Rackauckas},
 doi = {10.48550/arxiv.2201.12468},

src/integral.jl

@@ -49,7 +49,7 @@ julia> integrate(x * sin(a * x), (x, 0, 1); symbolic = true, detailed = false)
 
 Returns a tuple of (solved, unsolved, err) if `detailed == true`, where
 
-    solved: the solved integral 
+    solved: the solved integral
     unsolved: the residual unsolved portion of the input
     err: the numerical error in reaching the solution
 
@@ -122,7 +122,7 @@ function integrate(eq, xx::Tuple; kwargs...)
     sol = integrate(eq, x; kwargs...)
 
     if sol isa Tuple
-        if first(sol) != 0 && sol[2] == 0
+        if !isequal(sol, 0) && sol[2] == 0
             return substitute(first(sol), Dict(x => hi)) -
                    substitute(first(sol), Dict(x => lo))
         else
@@ -200,7 +200,7 @@ function integrate_sum(eq, x; bypass = false, kwargs...)
     return expand(solved), unsolved, ε₀
 end
 
-# integrate_term is the central part of the univariate integration code that 
+# integrate_term is the central part of the univariate integration code that
 # tries different methods to integrate `eq`.
 #
 # note: this function will be replaced with solver(prob::Problem) in symbolic.jl
@@ -295,8 +295,8 @@ function integrate_term(eq, x; kwargs...)
     end
 end
 
-# try_integrate is the main dispatch point to call different sparse solvers. 
-# It tries to find a linear combination of the basis, whose derivative is 
+# try_integrate is the main dispatch point to call different sparse solvers.
+# It tries to find a linear combination of the basis, whose derivative is
 # equal to eq
 function try_integrate(eq, x, basis; plan = default_plan())
     if isempty(basis)