Doctest & ref fixes

Author: Ames

docs/src/shocktube.md

@@ -25,7 +25,7 @@ using PyThermo.ShockTube
 driven = Species("N2", P=100u"kPa", T=300u"K")
 
 # Calculate shock jump for Mach 2.0 shock
-shocked, velocity = shockjump(driven, 2.0)
+shocked, velocity = shockjump!(driven, 2.0)
 
 println("Post-shock pressure: ", pressure(shocked))
 println("Post-shock temperature: ", temperature(shocked))
@@ -51,7 +51,7 @@ driven = Mixture(["He" => 0.95, "acetone" => 0.05], T=18u"°C", P=85u"kPa")
 Ms = 2.2
 
 # Calculate required driver pressure
-driver_calc = driverpressure(driver, driven, Ms)
+driver_calc = driverpressure!(driver, driven, Ms)
 println("Required driver pressure: ", pressure(driver_calc))
 # Output: ~3.73 MPa
 ```
@@ -73,7 +73,7 @@ driver = Species("He")
 driven = Mixture(["He" => 0.95, "acetone" => 0.05], T=18u"°C", P=85u"kPa")
 Ms = 2.2
 
-result = shockcalc(driver, driven, Ms)
+result = shockcalc!(driver, driven, Ms)
 
 # Access all states
 println("Initial driven state: ", result.driven)
@@ -168,7 +168,7 @@ This is equivalent to:
 
 ```julia
 # Manual approach
-shocked, u_shocked = shockjump(driven, Ms)
+shocked, u_shocked = shockjump!(driven, Ms)
 sol = riemann_interface(shocked, test_gas, u_shocked, 0.0)
 ```
 
@@ -277,7 +277,7 @@ driven = Mixture(["He" => 0.95, "acetone" => 0.05], T=18u"°C", P=85u"kPa")
 Ms = 2.2
 
 # Complete calculation
-result = shockcalc(driver, driven, Ms)
+result = shockcalc!(driver, driven, Ms)
 
 println("=== Shock Tube Results ===")
 println("Incident shock Mach: ", result.Ms)
@@ -354,13 +354,12 @@ end
 ## API reference
 
 ```@docs
-shockjump
 shockjump!
-driverpressure
 driverpressure!
-shockcalc
 shockcalc!
+ShockCalc
 riemann_interface
+RiemannSolution
 ```
 
 ## Theory background

src/ShockTube.jl

@@ -53,8 +53,8 @@ julia> round(Int, ustrip(u"m/s", velocity))
 ```
 
 # See Also
-- [`shockcalc`](@ref): Comprehensive shock tube calculation
-- [`driverpressure`](@ref): Calculate required driver pressure
+- [`shockcalc!`](@ref): Comprehensive shock tube calculation
+- [`driverpressure!`](@ref): Calculate required driver pressure
 """
 function shockjump!(gas, Mach)
     α1 = (γ(gas) + 1) / (γ(gas) - 1)
@@ -107,8 +107,8 @@ julia> round(typeof(1.0u"MPa"), pressure(driver_calc), digits=2)
 ```
 
 # See Also
-- [`shockjump`](@ref): Calculate shock jump conditions
-- [`shockcalc`](@ref): Comprehensive shock tube calculation
+- [`shockjump!`](@ref): Calculate shock jump conditions
+- [`shockcalc!`](@ref): Comprehensive shock tube calculation
 """
 function driverpressure!(driver, driven, Ms)
     γ1, γ4 = γ(driven), γ(driver)
@@ -165,7 +165,7 @@ julia> round(typeof(1.0u"kPa"), pressure(result.shocked), digits=1)
 ```
 
 # See Also
-- [`shockcalc`](@ref): Function that creates `ShockCalc` objects
+- [`shockcalc!`](@ref): Function that creates `ShockCalc` objects
 """
 struct ShockCalc{DR<:PyThermo.Chemical,DV<:PyThermo.Chemical,U<:Unitful.Velocity}
     driver::DR
@@ -235,7 +235,7 @@ julia> driver = Species("He")
 Species(He, 298.1 K, 1.013e+05 Pa)
 
 julia> driven = Mixture(["He" => 0.95, "acetone" => 0.05], T = 18u"°C", P = 85u"kPa")
-Mixture(95.0% helium, 5.00% acetone, 291.1 K, 8.500e+04 Pa)
+Mixture(95% He, 5% acetone, 291.1 K, 8.500e+04 Pa)
 
 julia> result = shockcalc(driver, driven, 2.2);
 
@@ -251,8 +251,8 @@ The calculation follows standard shock tube theory:
 4. Returns complete state information for analysis
 
 # See Also
-- [`shockjump`](@ref): Calculate shock jump conditions only
-- [`driverpressure`](@ref): Calculate required driver pressure only
+- [`shockjump!`](@ref): Calculate shock jump conditions only
+- [`driverpressure!`](@ref): Calculate required driver pressure only
 - [`ShockCalc`](@ref): Result container struct
 """
 function shockcalc!(driver, driven, Ms)
@@ -302,7 +302,7 @@ julia> round(u"m/s", sol.u_star, digits=1)
 ```
 
 # See Also
-- [`riemann_interface`](@ref): Function that creates `RiemannSolution` objects
+- [`shockcalc!`](@ref): Function that creates `ShockCalc` objects
 """
 struct RiemannSolution{P<:Unitful.Pressure, V<:Unitful.Velocity, D<:Unitful.Density, T<:Unitful.Temperature}
     p_star::P
@@ -409,7 +409,7 @@ julia> round(ustrip(u"m/s", sol.S_contact), digits=1)
 
 # See Also
 - [`RiemannSolution`](@ref): Result container struct
-- [`shockjump`](@ref): Calculate shock jump conditions
+- [`shockjump!`](@ref): Calculate shock jump conditions
 """
 function riemann_interface(left::Chemical, right::Chemical, u_left=0.0, u_right=0.0)
     # Extract primitive variables