Add Riemann solution calculator

Project.toml

@@ -1,7 +1,7 @@
 name = "KitBase"
 uuid = "86eed249-3a28-466f-8d3a-596821e1af9a"
 authors = ["Tianbai Xiao <[email protected]>"]
-version = "0.9.21"
+version = "0.9.22"
 
 [deps]
 BSON = "fbb218c0-5317-5bc6-957e-2ee96dd4b1f0"
@@ -25,6 +25,7 @@ Parameters = "d96e819e-fc66-5662-9728-84c9c7592b0a"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 RecipesBase = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
+Roots = "f2b01f46-fcfa-551c-844a-d8ac1e96c665"
 SciMLNLSolve = "e9a6253c-8580-4d32-9898-8661bb511710"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
@@ -38,18 +39,19 @@ CSV = "0.8, 0.9, 0.10"
 CUDA = "^2, ^3, ^4, ^5"
 Distributions = "0.21, 0.22, 0.23, 0.24, 0.25"
 FFTW = "^1"
-FastGaussQuadrature = "0.2, 0.3, 0.4, 0.5, 1"
+FastGaussQuadrature = "0.2, 0.3, 0.4, 0.5, ^1"
 FileIO = "^1"
 FiniteMesh = "0.3"
 ForwardDiff = "0.10"
 JLD2 = "0.4"
 MultivariatePolynomials = "0.3, 0.4, 0.5"
 NonlinearSolve = "^1, ^2, ^3"
-OffsetArrays = "1.8"
+OffsetArrays = "^1"
 Optim = "^1"
 Parameters = "0.12"
 RecipesBase = "^1"
 Reexport = "^1"
+Roots = "^2"
 SciMLNLSolve = "0.1"
 SpecialFunctions = "0.7, 0.8, 0.9, 0.10, ^1, ^2"
 StaticArrays = "^1"

src/KitBase.jl

@@ -16,6 +16,7 @@ import Base: *
 import BSON
 import JLD2
 import NonlinearSolve
+import Roots: Order1, find_zero
 import SciMLNLSolve: NLSolveJL
 using Base.Threads: @threads
 using CSV

src/Theory/theory.jl

@@ -38,3 +38,4 @@ include("theory_moments.jl")
 include("theory_particle.jl")
 include("theory_closure.jl")
 include("theory_quantum.jl")
+include("theory_solution.jl")