Allow subspace parameters to be arrays

Author: emstoudenmire

src/solvers/subspace/densitymatrix.jl

@@ -6,9 +6,9 @@ function subspace_expand(
   problem,
   local_state::ITensor,
   region_iterator;
+  expansion_factor,
+  max_expand,
   north_pass=1,
-  expansion_factor=default_expansion_factor(),
-  max_expand=default_max_expand(),
   trunc,
   kws...,
 )

src/solvers/subspace/subspace.jl

@@ -5,10 +5,27 @@ default_expansion_factor() = 1.5
 default_max_expand() = typemax(Int)
 
 function subspace_expand(
-  problem, local_state, region_iterator; subspace_algorithm=nothing, sweep, trunc, kws...
+  problem,
+  local_state,
+  region_iterator;
+  expansion_factor=default_expansion_factor(),
+  max_expand=default_max_expand(),
+  subspace_algorithm=nothing,
+  sweep,
+  trunc,
+  kws...,
 )
+  expansion_factor = get_or_last(expansion_factor, sweep)
+  max_expand = get_or_last(max_expand, sweep)
   return subspace_expand(
-    Backend(subspace_algorithm), problem, local_state, region_iterator; trunc, kws...
+    Backend(subspace_algorithm),
+    problem,
+    local_state,
+    region_iterator;
+    expansion_factor,
+    max_expand,
+    trunc,
+    kws...,
   )
 end