In this work we show the implementation of a linear scaling algorithm for the calculation of the Poisson integral. We use domain decomposition and non-uniform auxiliary grids (NGs) to calculate the electrostatic interaction. We demonstrate the approach within the PARSEC density functional theory code and perform calculations of long 1D carbon chains and other long molecules. Finally, we discuss possible applications to additional problems and geometries.