submodule(specialmatrices_circulant) circulant_inverse implicit none(type, external) contains module procedure inv_rdp B = circulant(real(ifft(1.0_dp/A%c_hat)) / A%n) end procedure end submodule
