Computing \(J\)-ideals of a matrix over a principal ideal domain (Q2359556)

If \(J\) is an ideal in a unique factorization domain \(D\) and \(B\) is an \(n\times n\) matrix over \(D\), then the \(J\)-ideal \(N_J(B)\) of \(B\) is defined by \[ N_J(B)=\{f\in D[X]:\;f(B)\equiv 0 \pmod{J}\}. \] This definition is due to the second author who showed earlier [Linear Algebra Appl. 494, 44--69 (2016; Zbl 1369.13030)] that to determine \(N_J(B)\) it suffices to find \(N_{p^t}(B)\) for \(t=1,2,\dots\) and a finite set of prime elements \(p\) of \(D\), and gave a formula for \(N_{p^t}(B)\) involving a certain finite set \(\Pi(p)\) of polynomials . In this paper the authors provide an algorithm determining \(\Pi(p)\).