$\newcommand{\Q}{{\mathbb Q}} $Let $f\in \Q[x]$ be a polynomial, and let $L/\Q$ be the finite Galois extension obtaining by adjoining to $\Q$ all roots of $f$. Magma knows how to compute $\Gamma:={\rm Gal}(L/\Q)$.
I am interested in algorithms and programs that can do the following:
(1) Find a finite set $S\ni\infty$ of places of $\Q$ such that the extension $L/\Q$ is unramified for all primes $p\in\overline S:=V(\Q)\smallsetminus S$ (where $V(\Q)$ denotes the set of all places of $\Q$).
(2) For each cyclic subgroup $\Delta\subseteq G$, somehow determine all primes $p\in\overline S$ with decomposition group $D_p\subseteq \Gamma$ conjugate to $\Delta$.
Example. Let $f(x)=x^2+1$. Then $L=\Q(\sqrt{-1})$ is unramified outside $\infty$ and $2$, and by Fermat's theorem, for $\Delta=\{1\}$ we obtain $p\equiv 1\pmod 4$, while for $\Delta=\Gamma={\rm Gal}(L/\Q)$ we have $p\equiv 3\pmod 4$.
EDIT:
Motivation. I want to describe using computer the Galois cohomology set $H^1(\Q,G)$ where $G$ is a connected reductive group over $\Q$. For this it more or less suffices to compute $H^1(\Q_v,G)$ for all places $v$ of $\Q$, and the Tate-Shafarevich kernel $Ш^1(\Q,G)$. There is a program computing $H^1({\mathbb R},G)$ in my paper with Willem A. de Graaf. The other guys can be computed from the decompositions groups $D_p$ for all primes p. This motivates my question.