import random


def random_param_shift(vals, sigmas):
    """Add a random (normal) shift to a parameter set, for testing"""

    assert len(vals) == len(sigmas)
    shifts = [random.gauss(0, sd) for sd in sigmas]
    newvals = [(x + y) for x, y in zip(vals, shifts)]

    return newvals