Skip to content

Seeding not working for single chains #5378

Closed
@ricardoV94

Description

@ricardoV94

Seeding seems to not be respected when using a single chain:

import numpy as np
import pymc as pm

with pm.Model(rng_seeder=3):
    x = pm.Normal("x")

    tr1 = pm.sample(
        chains=1,
        random_seed=1,
        tune=0,
        draws=10,
        return_inferencedata=False,
        compute_convergence_checks=False,
    )
    tr2 = pm.sample(
        chains=1,
        random_seed=1,
        tune=0,
        draws=10,
        return_inferencedata=False,
        compute_convergence_checks=False,
    )

    assert np.allclose(tr1["x"], tr2["x"])  # Fails

There is also a similar failure when passing multiple seeds, but working with a single core (i.e., sequential sampling)

#5377 Includes a test which covers the failing cases

Metadata

Metadata

Assignees

No one assigned

    Type

    No type

    Projects

    No projects

    Milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions