We develop a solution theory for singular elliptic stochastic PDEs with fractional
Laplacian, additive white noise and cubic nonlinearity. The method covers the whole
subcritical regime. It is based on the Wilsonian renormalization group theory and the
Polchinski flow equation.