code: def _merge_scipy_results(results, same_t=False):
      keys = list(results[0].keys())
      merged_results = {key: [] for key in keys}
  
      for result in results:
          for key, value in merged_results.items():
              value.append(result[key])
  
      if same_t:
          merged_results["t"] = gs.moveaxis(
              gs.stack(merged_results["t"]), 0, 1
          )
          merged_results["y"] = gs.moveaxis(
              gs.stack(merged_results["y"]), 0, 1
          )
  
      return merged_results
  
imports: from abc import ABC, abstractmethod
  import scipy
  import geomstats.backend
  import geomstats.integrator
  from geomstats.errors import check_parameter_accepted_values
  from geomstats.numerics._common import result_to_backend_type
label_desc: Use image node_id 1 for calling a global function with example usage: _merge_scipy_results(results, same_t) and returns: merged_results