From 0d4f3c3e2f016ff3f8c0055873fe2daa76210f69 Mon Sep 17 00:00:00 2001
From: Michael Reed <18372368+chakravala@users.noreply.github.com>
Date: Sat, 13 Apr 2024 18:56:58 -0400
Subject: [PATCH] fixed bug in expm1 evaluation

---
 Project.toml     |  2 +-
 src/composite.jl | 14 +++-----------
 2 files changed, 4 insertions(+), 12 deletions(-)

diff --git a/Project.toml b/Project.toml
index bcbca18..317b2a1 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "Grassmann"
 uuid = "4df31cd9-4c27-5bea-88d0-e6a7146666d8"
 authors = ["Michael Reed"]
-version = "0.8.9"
+version = "0.8.10"
 
 [deps]
 AbstractTensors = "a8e43f4a-99b7-5565-8bf1-0165161caaea"
diff --git a/src/composite.jl b/src/composite.jl
index 0492dd1..c3cf3a2 100644
--- a/src/composite.jl
+++ b/src/composite.jl
@@ -62,7 +62,6 @@ end
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             term .= out
-            out .= 0
             # term *= b/k
             $(loop[2])
             @inbounds norms .= (norms[2],norm(out),ns)
@@ -90,7 +89,6 @@ end
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             term .= out
-            out .= 0
             # term *= b/k
             $(loop[2])
             @inbounds norms .= (norms[2],norm(out),ns)
@@ -159,7 +157,7 @@ function Base.exp(t::Couple{V,B}) where {V,B}
     end
 end
 
-function Base.expm1(t::PseudoCouple{V}) where V
+function Base.expm1(t::PseudoCouple{V,B}) where {V,B}
     if isscalar(B)
         exp(t)-One(V)
     else
@@ -261,7 +259,6 @@ qlog_fast(b::AntiSpinor,x::Int=10000) = qlog_fast(Multivector(b),x)
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             prod .= out
-            out .= 0
             # prod *= w2
             $(loop[2])
             term .= out/k
@@ -293,7 +290,6 @@ end
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             prod .= out
-            out .= 0
             # prod *= w2
             $(loop[2])
             term .= out/k
@@ -437,7 +433,6 @@ end
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             term .= out
-            out .= 0
             # term *= τ/(k*(k-1))
             $(loop[2])
             @inbounds norms .= (norms[2],norm(out),ns)
@@ -467,7 +462,6 @@ end
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             term .= out
-            out .= 0
             # term *= τ/(k*(k-1))
             $(loop[2])
             @inbounds norms .= (norms[2],norm(out),ns)
@@ -523,7 +517,6 @@ end
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             term .= out
-            out .= 0
             # term *= τ/(k*(k-1))
             $(loop[2])
             @inbounds norms .= (norms[2],norm(out),ns)
@@ -552,7 +545,6 @@ end
             ns = norm(S)
             @inbounds ns ≈ norms[3] && break
             term .= out
-            out .= 0
             # term *= τ/(k*(k-1))
             $(loop[2])
             @inbounds norms .= (norms[2],norm(out),ns)
@@ -1084,9 +1076,9 @@ Base.rand(::AbstractRNG,::SamplerType{Couple}) = rand(Couple{rand(Manifold)})
 Base.rand(::AbstractRNG,::SamplerType{Couple{V}}) where V = rand(Couple{V,Submanifold{V}(UInt(rand(1:1<<mdims(V)-1)))})
 Base.rand(::AbstractRNG,::SamplerType{Couple{V,B}}) where {V,B} = Couple{V,B}(rand(Complex{Float64}))
 Base.rand(::AbstractRNG,::SamplerType{Couple{V,B,T}}) where {V,B,T} = Couple{V,B}(rand(Complex{T}))
-Base.rand(::AbstractRNG,::SamplerType{Couple{V,B,T} where B}) where {V,T} = rand(Couple{V,rand(Single{V}),T})
+Base.rand(::AbstractRNG,::SamplerType{Couple{V,B,T} where B}) where {V,T} = rand(Couple{V,Submanifold{V}(UInt(rand(1:1<<mdims(V)-1))),T})
 Base.rand(::AbstractRNG,::SamplerType{PseudoCouple}) = rand(PseudoCouple{rand(Manifold)})
 Base.rand(::AbstractRNG,::SamplerType{PseudoCouple{V}}) where V = rand(PseudoCouple{V,Submanifold{V}(UInt(rand(0:(1<<mdims(V)-1)-1)))})
 Base.rand(::AbstractRNG,::SamplerType{PseudoCouple{V,B}}) where {V,B} = PseudoCouple{V,B}(rand(Complex{Float64}))
 Base.rand(::AbstractRNG,::SamplerType{PseudoCouple{V,B,T}}) where {V,B,T} = PseudoCouple{V,B}(rand(Complex{T}))
-Base.rand(::AbstractRNG,::SamplerType{PseudoCouple{V,B,T} where B}) where {V,T} = rand(PseudoCouple{V,rand(Single{V}),T})
+Base.rand(::AbstractRNG,::SamplerType{PseudoCouple{V,B,T} where B}) where {V,T} = rand(PseudoCouple{V,Submanifold{V}(UInt(rand(0:(1<<mdims(V)-1)-1))),T})