QuantumKitHub · leburgel · Nov 20, 2025 · Oct 31, 2025 · Oct 31, 2025 · Oct 31, 2025
diff --git a/docs/make.jl b/docs/make.jl
@@ -31,7 +31,8 @@ links = InterLinks(
     "TensorOperations" => "https://quantumkithub.github.io/TensorOperations.jl/stable/",
     "KrylovKit" => "https://jutho.github.io/KrylovKit.jl/stable/",
     "BlockTensorKit" => "https://lkdvos.github.io/BlockTensorKit.jl/dev/",
-    "MatrixAlgebraKit" => "https://quantumkithub.github.io/MatrixAlgebraKit.jl/stable/"
+    "MatrixAlgebraKit" => "https://quantumkithub.github.io/MatrixAlgebraKit.jl/stable/",
+    "MPSKitModels" => "https://quantumkithub.github.io/MPSKitModels.jl/dev/"
 )
 
 # include MPSKit in all doctests

diff --git a/docs/src/examples/quantum1d/8.bose-hubbard/index.md b/docs/src/examples/quantum1d/8.bose-hubbard/index.md
diff --git a/docs/src/examples/quantum1d/8.bose-hubbard/main.ipynb b/docs/src/examples/quantum1d/8.bose-hubbard/main.ipynb
diff --git a/examples/Cache.toml b/examples/Cache.toml
@@ -5,6 +5,7 @@
 "2.haldane" = "804433b690faa1ce268a430edb4185af77a38de7a9f53e097795688a53c30c5e"
 "6.hubbard" = "af14e1df11b2392a69260ce061a716370a2ce50b5c907f0a7d2548e796d543fe"
 "7.xy-finiteT" = "7afb26fb9ff8ef84722fee3442eaed2ddffda4a5f70a7844b61ce5178093706c"
+"8.bose-hubbard" = "4ca414d4b76bbffbd88a7c153dce245c9fab5b7cee0fe0b706b9a10ed87e29e9"
 "3.ising-dqpt" = "2d314fd05a75c5c91ff81391c7bf166171b35a7b32933e9490756dc584fe8437"
 "5.haldane-spt" = "3de1b3baa1e4c5dc2252bbe8688d0c6ac8e5d5265deeb6ab66818bbfb02dd4aa"
 "4.xxz-heisenberg" = "1a2093a182f3ce070b53488484d1024e45496b291c1b74e3f370201d4c056be2"

diff --git a/examples/Project.toml b/examples/Project.toml
@@ -2,6 +2,7 @@
 BenchmarkFreeFermions = "5b68a00d-ca4d-4b58-ab0c-dc082ade624e"
 Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
 KrylovKit = "0b1a1467-8014-51b9-945f-bf0ae24f4b77"
+LaTeXStrings = "b964fa9f-0449-5b57-a5c2-d3ea65f4040f"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Literate = "98b081ad-f1c9-55d3-8b20-4c87d4299306"
 MPSKit = "bb1c41ca-d63c-52ed-829e-0820dda26502"

diff --git a/examples/quantum1d/8.bose-hubbard/main.jl b/examples/quantum1d/8.bose-hubbard/main.jl
diff --git a/src/algorithms/toolbox.jl b/src/algorithms/toolbox.jl
@@ -43,7 +43,7 @@ Calculate the Galerkin error, which is the error between the solution of the ori
 Concretely, this is the overlap of the current state with the single-site derivative, projected onto the nullspace of the current state:
 
 ```math
-\\epsilon = |VL * (VL' * \\frac{above}{\\partial AC_{pos}})|
+\\epsilon = \\left|VL ⋅ \\left(VL^{\\dagger} ⋅ \\frac{\\partial \\text{above}}{\\partial AC_{\\text{pos}}}\\right)\\right|
 ```
 """
 function calc_galerkin(