imt-ring 1.4.1__py3-none-any.whl → 1.5.1__py3-none-any.whl

ring/algorithms/generator/base.py

@@ -1,10 +1,8 @@
-from functools import partial
-from typing import Callable, Optional, Sequence
+from typing import Callable, Optional
 import warnings
 
 import jax
 import jax.numpy as jnp
-import tqdm
 import tree_utils
 
 from ring import base
@@ -12,9 +10,9 @@ from ring import utils
 from ring.algorithms import jcalc
 from ring.algorithms import kinematics
 from ring.algorithms.generator import batch
+from ring.algorithms.generator import finalize_fns
 from ring.algorithms.generator import motion_artifacts
-from ring.algorithms.generator import randomize
-from ring.algorithms.generator import transforms
+from ring.algorithms.generator import setup_fns
 from ring.algorithms.generator import types
 
 
@@ -26,92 +24,137 @@ class RCMG:
         setup_fn: Optional[types.SETUP_FN] = None,
         finalize_fn: Optional[types.FINALIZE_FN] = None,
         add_X_imus: bool = False,
-        add_X_imus_kwargs: Optional[dict] = None,
+        add_X_imus_kwargs: dict = dict(),
         add_X_jointaxes: bool = False,
-        add_X_jointaxes_kwargs: Optional[dict] = None,
+        add_X_jointaxes_kwargs: dict = dict(),
         add_y_relpose: bool = False,
         add_y_rootincl: bool = False,
         sys_ml: Optional[base.System] = None,
         randomize_positions: bool = False,
         randomize_motion_artifacts: bool = False,
         randomize_joint_params: bool = False,
-        randomize_anchors: bool = False,
-        randomize_anchors_kwargs: Optional[dict] = None,
-        randomize_hz: bool = False,
-        randomize_hz_kwargs: Optional[dict] = None,
         imu_motion_artifacts: bool = False,
-        imu_motion_artifacts_kwargs: Optional[dict] = None,
+        imu_motion_artifacts_kwargs: dict = dict(),
         dynamic_simulation: bool = False,
-        dynamic_simulation_kwargs: Optional[dict] = None,
+        dynamic_simulation_kwargs: dict = dict(),
         output_transform: Optional[Callable] = None,
         keep_output_extras: bool = False,
         use_link_number_in_Xy: bool = False,
+        cor: bool = False,
+        disable_tqdm: bool = False,
     ) -> None:
 
-        randomize_anchors_kwargs = _copy_kwargs(randomize_anchors_kwargs)
-        randomize_hz_kwargs = _copy_kwargs(randomize_hz_kwargs)
-
-        if randomize_hz:
-            finalize_fn = randomize.randomize_hz_finalize_fn_factory(finalize_fn)
-
-        partial_build_gen = partial(
-            _build_generator_lazy,
-            setup_fn=setup_fn,
-            finalize_fn=finalize_fn,
-            add_X_imus=add_X_imus,
-            add_X_imus_kwargs=add_X_imus_kwargs,
-            add_X_jointaxes=add_X_jointaxes,
-            add_X_jointaxes_kwargs=add_X_jointaxes_kwargs,
-            add_y_relpose=add_y_relpose,
-            add_y_rootincl=add_y_rootincl,
-            randomize_positions=randomize_positions,
-            randomize_motion_artifacts=randomize_motion_artifacts,
-            randomize_joint_params=randomize_joint_params,
-            imu_motion_artifacts=imu_motion_artifacts,
-            imu_motion_artifacts_kwargs=imu_motion_artifacts_kwargs,
-            dynamic_simulation=dynamic_simulation,
-            dynamic_simulation_kwargs=dynamic_simulation_kwargs,
-            output_transform=output_transform,
-            keep_output_extras=keep_output_extras,
-            use_link_number_in_Xy=use_link_number_in_Xy,
-        )
-
         sys, config = utils.to_list(sys), utils.to_list(config)
+        sys_ml = sys[0] if sys_ml is None else sys_ml
 
-        if randomize_anchors:
-            assert (
-                len(sys) == 1
-            ), "If `randomize_anchors`, then only one system is expected"
-            sys = randomize.randomize_anchors(sys[0], **randomize_anchors_kwargs)
+        for c in config:
+            assert c.is_feasible()
 
-        if randomize_hz:
-            sys, config = randomize.randomize_hz(sys, config, **randomize_hz_kwargs)
-        else:
-            # create zip
-            N_sys = len(sys)
-            sys = sum([len(config) * [s] for s in sys], start=[])
-            config = N_sys * config
-            assert len(sys) == len(config)
-
-        if sys_ml is None:
-            # TODO
-            if False and len(sys) > 1:
-                warnings.warn(
-                    "Batched simulation with multiple systems but no explicit `sys_ml`"
+        self.gens = []
+        for _sys in sys:
+            self.gens.append(
+                _build_mconfig_batched_generator(
+                    sys=_sys,
+                    config=config,
+                    setup_fn=setup_fn,
+                    finalize_fn=finalize_fn,
+                    add_X_imus=add_X_imus,
+                    add_X_imus_kwargs=add_X_imus_kwargs,
+                    add_X_jointaxes=add_X_jointaxes,
+                    add_X_jointaxes_kwargs=add_X_jointaxes_kwargs,
+                    add_y_relpose=add_y_relpose,
+                    add_y_rootincl=add_y_rootincl,
+                    sys_ml=sys_ml,
+                    randomize_positions=randomize_positions,
+                    randomize_motion_artifacts=randomize_motion_artifacts,
+                    randomize_joint_params=randomize_joint_params,
+                    imu_motion_artifacts=imu_motion_artifacts,
+                    imu_motion_artifacts_kwargs=imu_motion_artifacts_kwargs,
+                    dynamic_simulation=dynamic_simulation,
+                    dynamic_simulation_kwargs=dynamic_simulation_kwargs,
+                    output_transform=output_transform,
+                    keep_output_extras=keep_output_extras,
+                    use_link_number_in_Xy=use_link_number_in_Xy,
+                    cor=cor,
                 )
-            sys_ml = sys[0]
+            )
 
-        self.gens = []
-        for _sys, _config in tqdm.tqdm(
-            zip(sys, config), desc="building generators", total=len(sys)
-        ):
-            self.gens.append(partial_build_gen(sys=_sys, config=_config, sys_ml=sys_ml))
+        self._n_mconfigs = len(config)
+        self._size_of_generators = [self._n_mconfigs] * len(self.gens)
+
+        self._disable_tqdm = disable_tqdm
+
+    def _compute_repeats(self, sizes: int | list[int]) -> list[int]:
+        "how many times the generators are repeated to create a batch of `sizes`"
+
+        S, L = sum(self._size_of_generators), len(self._size_of_generators)
+
+        def assert_size(size: int):
+            assert self._n_mconfigs in utils.primes(size), (
+                f"`size`={size} is not divisible by number of "
+                + f"`mconfigs`={self._n_mconfigs}"
+            )
+
+        if isinstance(sizes, int):
+            assert (sizes // S) > 0, f"Batchsize or size too small. {sizes} < {S}"
+            assert sizes % S == 0, f"`size`={sizes} not divisible by {S}"
+            repeats = L * [sizes // S]
+        else:
+            for size in sizes:
+                assert_size(size)
+
+            assert len(sizes) == len(
+                self.gens
+            ), f"len(`sizes`)={len(sizes)} != {len(self.gens)}"
 
-    def _to_data(self, sizes, seed):
-        return batch.batch_generators_eager_to_list(self.gens, sizes, seed=seed)
+            repeats = [
+                size // size_of_gen
+                for size, size_of_gen in zip(sizes, self._size_of_generators)
+            ]
+            assert 0 not in repeats
+
+        return repeats
+
+    def to_lazy_gen(
+        self, sizes: int | list[int] = 1, jit: bool = True
+    ) -> types.BatchedGenerator:
+        return batch.generators_lazy(self.gens, self._compute_repeats(sizes), jit)
+
+    @staticmethod
+    def _number_of_executions_required(size: int) -> int:
+        _, vmap = utils.distribute_batchsize(size)
+
+        eager_threshold = utils.batchsize_thresholds()[1]
+        primes = iter(utils.primes(vmap))
+        n_calls = 1
+        while vmap > eager_threshold:
+            prime = next(primes)
+            n_calls *= prime
+            vmap /= prime
+
+        return n_calls
 
     def to_list(self, sizes: int | list[int] = 1, seed: int = 1):
-        return self._to_data(sizes, seed)
+        "Returns list of unbatched sequences as numpy arrays."
+        repeats = self._compute_repeats(sizes)
+        sizes = list(jnp.array(repeats) * jnp.array(self._size_of_generators))
+
+        reduced_repeats = []
+        n_calls = []
+        for size, repeat in zip(sizes, repeats):
+            n_call = self._number_of_executions_required(size)
+            gcd = utils.gcd(n_call, repeat)
+            n_calls.append(gcd)
+            reduced_repeats.append(repeat // gcd)
+        jits = [N > 1 for N in n_calls]
+
+        gens = []
+        for i in range(len(repeats)):
+            gens.append(
+                batch.generators_lazy([self.gens[i]], [reduced_repeats[i]], jits[i])
+            )
+
+        return batch.generators_eager_to_list(gens, n_calls, seed, self._disable_tqdm)
 
     def to_pickle(
         self,
@@ -120,19 +163,9 @@ class RCMG:
         seed: int = 1,
         overwrite: bool = True,
     ) -> None:
-        data = tree_utils.tree_batch(self._to_data(sizes, seed))
+        data = tree_utils.tree_batch(self.to_list(sizes, seed))
         utils.pickle_save(data, path, overwrite=overwrite)
 
-    def to_hdf5(
-        self,
-        path: str,
-        sizes: int | list[int] = 1,
-        seed: int = 1,
-        overwrite: bool = True,
-    ) -> None:
-        data = tree_utils.tree_batch(self._to_data(sizes, seed))
-        utils.hdf5_save(path, data, overwrite=overwrite)
-
     def to_eager_gen(
         self,
         batchsize: int = 1,
@@ -140,14 +173,15 @@ class RCMG:
         seed: int = 1,
         shuffle: bool = True,
     ) -> types.BatchedGenerator:
-        return batch.batch_generators_eager(
-            self.gens, sizes, batchsize, seed=seed, shuffle=shuffle
-        )
+        data = self.to_list(sizes, seed)
+        assert len(data) >= batchsize
 
-    def to_lazy_gen(
-        self, sizes: int | list[int] = 1, jit: bool = True
-    ) -> types.BatchedGenerator:
-        return batch.batch_generators_lazy(self.gens, sizes, jit=jit)
+        def data_fn(indices: list[int]):
+            return tree_utils.tree_batch([data[i] for i in indices])
+
+        return batch.generator_from_data_fn(
+            data_fn, list(range(len(data))), shuffle, batchsize
+        )
 
     @staticmethod
     def eager_gen_from_paths(
@@ -159,7 +193,7 @@ class RCMG:
         tree_transform=None,
     ) -> tuple[types.BatchedGenerator, int]:
         paths = utils.to_list(paths)
-        return batch.batched_generator_from_paths(
+        return batch.generator_from_paths(
             paths,
             batchsize,
             include_samples,
@@ -169,19 +203,26 @@ class RCMG:
         )
 
 
-def _copy_kwargs(kwargs: dict | None) -> dict:
-    return dict() if kwargs is None else kwargs.copy()
+def _copy_dicts(f) -> dict:
+    def _f(*args, **kwargs):
+        _copy = lambda obj: obj.copy() if isinstance(obj, dict) else obj
+        args = tuple([_copy(ele) for ele in args])
+        kwargs = {k: _copy(v) for k, v in kwargs.items()}
+        return f(*args, **kwargs)
+
+    return _f
 
 
-def _build_generator_lazy(
+@_copy_dicts
+def _build_mconfig_batched_generator(
     sys: base.System,
-    config: jcalc.MotionConfig,
+    config: list[jcalc.MotionConfig],
     setup_fn: types.SETUP_FN | None,
     finalize_fn: types.FINALIZE_FN | None,
     add_X_imus: bool,
-    add_X_imus_kwargs: dict | None,
+    add_X_imus_kwargs: dict,
     add_X_jointaxes: bool,
-    add_X_jointaxes_kwargs: dict | None,
+    add_X_jointaxes_kwargs: dict,
     add_y_relpose: bool,
     add_y_rootincl: bool,
     sys_ml: base.System,
@@ -189,23 +230,14 @@ def _build_generator_lazy(
     randomize_motion_artifacts: bool,
     randomize_joint_params: bool,
     imu_motion_artifacts: bool,
-    imu_motion_artifacts_kwargs: dict | None,
+    imu_motion_artifacts_kwargs: dict,
     dynamic_simulation: bool,
-    dynamic_simulation_kwargs: dict | None,
+    dynamic_simulation_kwargs: dict,
     output_transform: Callable | None,
     keep_output_extras: bool,
     use_link_number_in_Xy: bool,
-) -> types.Generator | types.GeneratorWithOutputExtras:
-    assert config.is_feasible()
-
-    imu_motion_artifacts_kwargs = _copy_kwargs(imu_motion_artifacts_kwargs)
-    dynamic_simulation_kwargs = _copy_kwargs(dynamic_simulation_kwargs)
-    add_X_imus_kwargs = _copy_kwargs(add_X_imus_kwargs)
-    add_X_jointaxes_kwargs = _copy_kwargs(add_X_jointaxes_kwargs)
-
-    # default kwargs values
-    if "hide_injected_bodies" not in imu_motion_artifacts_kwargs:
-        imu_motion_artifacts_kwargs["hide_injected_bodies"] = True
+    cor: bool,
+) -> types.BatchedGenerator:
 
     if add_X_jointaxes or add_y_relpose or add_y_rootincl:
         if len(sys_ml.findall_imus()) > 0:
@@ -227,183 +259,120 @@ def _build_generator_lazy(
                 "`imu_motion_artifacts` is enabled but not `randomize_motion_artifacts`"
             )
 
-        if "hide_injected_bodies" in imu_motion_artifacts_kwargs:
-            if imu_motion_artifacts_kwargs["hide_injected_bodies"] and False:
-                warnings.warn(
-                    "The flag `hide_injected_bodies` in `imu_motion_artifacts_kwargs` "
-                    "is set. This will try to hide injected bodies. This feature is "
-                    "experimental."
-                )
-
         if "prob_rigid" in imu_motion_artifacts_kwargs:
             assert randomize_motion_artifacts, (
                 "`prob_rigid` works by overwriting damping and stiffness parameters "
                 "using the `randomize_motion_artifacts` flag, so it must be enabled."
             )
 
-    noop = lambda gen: gen
-    return GeneratorPipe(
-        transforms.GeneratorTrafoSetupFn(setup_fn) if setup_fn is not None else noop,
-        (
-            transforms.GeneratorTrafoSetupFn(jcalc._init_joint_params)
-            if randomize_joint_params
-            else noop
-        ),
-        transforms.GeneratorTrafoRandomizePositions() if randomize_positions else noop,
-        (
-            transforms.GeneratorTrafoSetupFn(
+    def _setup_fn(key: types.PRNGKey, sys: base.System) -> base.System:
+        pipe = []
+        if imu_motion_artifacts and randomize_motion_artifacts:
+            pipe.append(
                 motion_artifacts.setup_fn_randomize_damping_stiffness_factory(
-                    prob_rigid=imu_motion_artifacts_kwargs.get("prob_rigid", 0.0),
-                    all_imus_either_rigid_or_flex=imu_motion_artifacts_kwargs.get(
-                        "all_imus_either_rigid_or_flex", False
-                    ),
-                    imus_surely_rigid=imu_motion_artifacts_kwargs.get(
-                        "imus_surely_rigid", []
-                    ),
+                    **imu_motion_artifacts_kwargs
                 )
             )
-            if (imu_motion_artifacts and randomize_motion_artifacts)
-            else noop
-        ),
-        # all the generator trafors before this point execute in reverse order
-        # to see this, consider gen[0] and gen[1]
-        # the GeneratorPipe will unpack into the following:
-        # gen[1] will unfold into
-        # >>> sys = gen[1].setup_fn(sys)
-        # >>> return gen[0](sys)
-        # <-------------------- GENERATOR MIDDLE POINT ------------------------->
-        # all the generator trafos after this point execute in order
-        # >>> Xy, extras = gen[-2](*args)
-        # >>> return gen[-1].finalize_fn(extras)
-        (
-            transforms.GeneratorTrafoDynamicalSimulation(**dynamic_simulation_kwargs)
-            if dynamic_simulation
-            else noop
-        ),
-        (
-            motion_artifacts.GeneratorTrafoHideInjectedBodies()
-            if (
-                imu_motion_artifacts
-                and imu_motion_artifacts_kwargs["hide_injected_bodies"]
-            )
-            else noop
-        ),
-        (
-            transforms.GeneratorTrafoFinalizeFn(finalize_fn)
-            if finalize_fn is not None
-            else noop
-        ),
-        transforms.GeneratorTrafoIMU(**add_X_imus_kwargs) if add_X_imus else noop,
-        (
-            transforms.GeneratorTrafoJointAxisSensor(
-                sys_noimu, **add_X_jointaxes_kwargs
-            )
-            if add_X_jointaxes
-            else noop
-        ),
-        transforms.GeneratorTrafoRelPose(sys_noimu) if add_y_relpose else noop,
-        transforms.GeneratorTrafoRootIncl(sys_noimu) if add_y_rootincl else noop,
-        (
-            transforms.GeneratorTrafoNames2Indices(sys_noimu)
-            if use_link_number_in_Xy
-            else noop
-        ),
-        GeneratorTrafoRemoveInputExtras(sys),
-        noop if keep_output_extras else GeneratorTrafoRemoveOutputExtras(),
-        (
-            transforms.GeneratorTrafoLambda(output_transform, input=False)
-            if output_transform is not None
-            else noop
-        ),
-    )(config)
-
-
-def _generator_with_extras(
-    config: jcalc.MotionConfig,
-) -> types.GeneratorWithInputOutputExtras:
-    def generator(
-        key: types.PRNGKey, sys: base.System
-    ) -> tuple[types.Xy, types.OutputExtras]:
-        if config.cor:
+        if randomize_positions:
+            pipe.append(setup_fns._setup_fn_randomize_positions)
+        if randomize_joint_params:
+            pipe.append(jcalc._init_joint_params)
+        if setup_fn is not None:
+            pipe.append(setup_fn)
+
+        for f in pipe:
+            key, consume = jax.random.split(key)
+            sys = f(consume, sys)
+        if cor:
             sys = sys._replace_free_with_cor()
-
-        key_start = key
-        # build generalized coordintes vector `q`
-        q_list = []
-
-        def draw_q(key, __, link_type, link):
-            joint_params = link.joint_params
-            # limit scope
-            joint_params = (
-                joint_params[link_type]
-                if link_type in joint_params
-                else joint_params["default"]
+        return sys
+
+    def _finalize_fn(Xy: types.Xy, extras: types.OutputExtras):
+        pipe = []
+        if dynamic_simulation:
+            pipe.append(finalize_fns.DynamicalSimulation(**dynamic_simulation_kwargs))
+        if imu_motion_artifacts and imu_motion_artifacts_kwargs.get(
+            "hide_injected_bodies", True
+        ):
+            pipe.append(motion_artifacts.HideInjectedBodies())
+        if finalize_fn is not None:
+            pipe.append(finalize_fns.FinalizeFn(finalize_fn))
+        if add_X_imus:
+            pipe.append(finalize_fns.IMU(**add_X_imus_kwargs))
+        if add_X_jointaxes:
+            pipe.append(
+                finalize_fns.JointAxisSensor(sys_noimu, **add_X_jointaxes_kwargs)
             )
-            if key is None:
-                key = key_start
-            key, key_t, key_value = jax.random.split(key, 3)
-            draw_fn = jcalc.get_joint_model(link_type).rcmg_draw_fn
-            if draw_fn is None:
-                raise Exception(f"The joint type {link_type} has no draw fn specified.")
-            q_link = draw_fn(config, key_t, key_value, sys.dt, joint_params)
-            # even revolute and prismatic joints must be 2d arrays
-            q_link = q_link if q_link.ndim == 2 else q_link[:, None]
-            q_list.append(q_link)
-            return key
-
-        keys = sys.scan(draw_q, "ll", sys.link_types, sys.links)
-        # stack of keys; only the last key is unused
-        key = keys[-1]
-
-        q = jnp.concatenate(q_list, axis=1)
-
-        # do forward kinematics
-        x, _ = jax.vmap(kinematics.forward_kinematics_transforms, (None, 0))(sys, q)
-
-        Xy = ({}, {})
-        return Xy, (key, q, x, sys)
-
-    return generator
-
-
-class GeneratorPipe:
-    def __init__(self, *gen_trafos: Sequence[types.GeneratorTrafo]):
-        self._gen_trafos = gen_trafos
-
-    def __call__(
-        self, config: jcalc.MotionConfig
-    ) -> (
-        types.GeneratorWithInputOutputExtras
-        | types.GeneratorWithOutputExtras
-        | types.GeneratorWithInputExtras
-        | types.Generator
-    ):
-        gen = _generator_with_extras(config)
-        for trafo in self._gen_trafos:
-            gen = trafo(gen)
-        return gen
-
-
-class GeneratorTrafoRemoveInputExtras(types.GeneratorTrafo):
-    def __init__(self, sys: base.System):
-        self.sys = sys
-
-    def __call__(
-        self,
-        gen: types.GeneratorWithInputExtras | types.GeneratorWithInputOutputExtras,
-    ) -> types.Generator | types.GeneratorWithOutputExtras:
-        def _gen(key):
-            return gen(key, self.sys)
-
-        return _gen
-
-
-class GeneratorTrafoRemoveOutputExtras(types.GeneratorTrafo):
-    def __call__(
-        self,
-        gen: types.GeneratorWithOutputExtras | types.GeneratorWithInputOutputExtras,
-    ) -> types.Generator | types.GeneratorWithInputExtras:
-        def _gen(*args):
-            return gen(*args)[0]
-
-        return _gen
+        if add_y_relpose:
+            pipe.append(finalize_fns.RelPose(sys_noimu))
+        if add_y_rootincl:
+            pipe.append(finalize_fns.RootIncl(sys_noimu))
+        if use_link_number_in_Xy:
+            pipe.append(finalize_fns.Names2Indices(sys_noimu))
+
+        for f in pipe:
+            Xy, extras = f(Xy, extras)
+        return Xy, extras
+
+    def _gen(key: types.PRNGKey):
+        key, *consume = jax.random.split(key, len(config) + 1)
+        syss = jax.vmap(_setup_fn, (0, None))(jnp.array(consume), sys)
+
+        qs = []
+        for i, _config in enumerate(config):
+            key, _q = draw_random_q(key, syss[i], _config)
+            qs.append(_q)
+        qs = jnp.stack(qs)
+
+        @jax.vmap
+        def _vmapped_context(key, q, sys):
+            x, _ = jax.vmap(kinematics.forward_kinematics_transforms, (None, 0))(sys, q)
+            Xy, extras = ({}, {}), (key, q, x, sys)
+            return _finalize_fn(Xy, extras)
+
+        keys = jax.random.split(key, len(config))
+        Xy, extras = _vmapped_context(keys, qs, syss)
+        output = (Xy, extras) if keep_output_extras else Xy
+        output = output if output_transform is None else output_transform(output)
+        return output
+
+    return _gen
+
+
+def draw_random_q(
+    key: types.PRNGKey,
+    sys: base.System,
+    config: jcalc.MotionConfig,
+) -> tuple[types.Xy, types.OutputExtras]:
+
+    key_start = key
+    # build generalized coordintes vector `q`
+    q_list = []
+
+    def draw_q(key, __, link_type, link):
+        joint_params = link.joint_params
+        # limit scope
+        joint_params = (
+            joint_params[link_type]
+            if link_type in joint_params
+            else joint_params["default"]
+        )
+        if key is None:
+            key = key_start
+        key, key_t, key_value = jax.random.split(key, 3)
+        draw_fn = jcalc.get_joint_model(link_type).rcmg_draw_fn
+        if draw_fn is None:
+            raise Exception(f"The joint type {link_type} has no draw fn specified.")
+        q_link = draw_fn(config, key_t, key_value, sys.dt, joint_params)
+        # even revolute and prismatic joints must be 2d arrays
+        q_link = q_link if q_link.ndim == 2 else q_link[:, None]
+        q_list.append(q_link)
+        return key
+
+    keys = sys.scan(draw_q, "ll", sys.link_types, sys.links)
+    # stack of keys; only the last key is unused
+    key = keys[-1]
+
+    q = jnp.concatenate(q_list, axis=1)
+
+    return key, q