fix

putianyi889 · Aug 11, 2024 · 5e8680c · 5e8680c
1 parent d22cb7e
commit 5e8680c
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 4 deletions.
diff --git a/src/BandedSylvesters.jl b/src/BandedSylvesters.jl
@@ -28,7 +28,7 @@ end
 
 buffer(R::BandedSylvesterRecurrence) = R.buffer
 function rdiv!(P::BandedSylvesterRecurrence, x)
-    rdiv!(P.buffer)
+    rdiv!(P.buffer, x)
     P
 end
 

diff --git a/src/PseudostableRecurrences.jl b/src/PseudostableRecurrences.jl
@@ -20,9 +20,10 @@ export stable_recurrence, precision_shift
 Estimates `log2` of the amplification of `P` by performing a full recurrence based on random initial conditions.
 """
 function precision_shift(P::AbstractLinearRecurrencePlan)
-    test, _ = init(P; init=:rand)
-    shift = 1
-    testmax = 1.0
+    test, start = init(P; init=:rand)
+    shift = 4 # safe choice from experiments
+    testmax = norm(start, Inf)
+    rdiv!(test, testmax)
     while true
         v = step!(test)
         if isnothing(v)