Update progressbar managers with existing fit results from ZarrTrace

Author: lucianopaz

pymc/sampling/mcmc.py

@@ -1157,13 +1157,25 @@ def _sample_many(
 
     with progress_manager:
         for i in range(chains):
+            trace = traces[i]
+            progress_manager.set_initial_state(
+                *trace.completed_draws_and_divergences(chain_specific=True)
+            )
+            progress_manager._progress.update(
+                progress_manager.tasks[i],
+                draws=progress_manager.completed_draws
+                if progress_manager.combined_progress
+                else progress_manager.draws,
+                divergences=progress_manager.divergences,
+                refresh=True,
+            )
             step.sampling_state = initial_step_state
             _sample(
                 draws=draws,
                 chain=i,
                 start=start[i],
                 step=step,
-                trace=traces[i],
+                trace=trace,
                 rng=rngs[i],
                 callback=callback,
                 progress_manager=progress_manager,

pymc/sampling/parallel.py

@@ -512,6 +512,9 @@ def __init__(
             progressbar=progressbar,
             progressbar_theme=progressbar_theme,
         )
+        if traces is not None:
+            for trace in traces:
+                self._progress.set_initial_state(*trace.completed_draws_and_divergences())
 
     def _make_active(self):
         while self._inactive and len(self._active) < self._max_active:

pymc/sampling/population.py

@@ -27,7 +27,6 @@
 from rich.progress import BarColumn, TextColumn, TimeElapsedColumn, TimeRemainingColumn
 
 from pymc.backends.base import BaseTrace
-from pymc.backends.zarr import ZarrChain
 from pymc.initial_point import PointType
 from pymc.model import Model, modelcontext
 from pymc.stats.convergence import log_warning_stats
@@ -110,6 +109,10 @@ def _sample_population(
 
     with CustomProgress(disable=not progressbar) as progress:
         task = progress.add_task("[red]Sampling...", total=draws)
+        for trace in traces:
+            progress.update(
+                task, completed=trace.completed_draws_and_divergences(chain_specific=True)[0]
+            )
         for _ in sampling:
             progress.update(task)
 
@@ -197,6 +200,7 @@ def __init__(
                         #     enumerate(progress_bar(steppers)) if progressbar else enumerate(steppers)
                         # ):
                         task = self._progress.add_task(description=f"Chain {c}")
+                        self._progress.update(task, completed=first_draw_idx)
                         secondary_end, primary_end = multiprocessing.Pipe()
                         stepper_dumps = cloudpickle.dumps(stepper, protocol=4)
                         process = multiprocessing.Process(

pymc/util.py

@@ -812,6 +812,7 @@ def __init__(
 
         self._show_progress = show_progress
         self.divergences = 0
+        self.draws = 0
         self.completed_draws = 0
         self.total_draws = draws + tune
         self.desc = "Sampling chain"
@@ -827,27 +828,35 @@ def __enter__(self):
     def __exit__(self, exc_type, exc_val, exc_tb):
         return self._progress.__exit__(exc_type, exc_val, exc_tb)
 
+    def set_initial_state(self, draws: int = 0, divergences: int = 0):
+        self.draws = draws
+        self.completed_draws += draws
+        self.divergences += divergences
+
     def _initialize_tasks(self):
         if self.combined_progress:
             self.tasks = [
                 self._progress.add_task(
                     self.desc.format(self),
-                    completed=0,
-                    draws=0,
+                    completed=self.completed_draws,
+                    draws=self.completed_draws,
                     total=self.total_draws * self.chains - 1,
                     chain_idx=0,
                     sampling_speed=0,
                     speed_unit="draws/s",
-                    **{stat: value[0] for stat, value in self.progress_stats.items()},
+                    **{
+                        stat: value[0] if stat != "diverging" else self.divergences
+                        for stat, value in self.progress_stats.items()
+                    },
                 )
             ]
 
         else:
             self.tasks = [
                 self._progress.add_task(
                     self.desc.format(self),
-                    completed=0,
-                    draws=0,
+                    completed=self.completed_draws,
+                    draws=self.draws,
                     total=self.total_draws - 1,
                     chain_idx=chain_idx,
                     sampling_speed=0,