Let k be a field of characteristic 0. Given a polynomial mapping f = (f1, . . . , fp) from kn to kp, the local Bernstein–Sato ideal of f at a point a ∈ kn is defined as an ideal of the ring of polynomials in s = (s1, . . . , sp). We propose an algorithm for computing local Bernstein–Sato ideals by combining Gröbner bases in rings of differential operators with primary decomposition in a polynom...