The interpolation finite difference method (IFDM) allows free numerical analysis of elliptic partial differential equations over arbitrary domains. Conventionally, in the finite difference method (FDM), the calculation is performed using the second-order accuracy central difference. For engineering problems, second-order accuracy calculations are often sufficient. On the other hand, much research has been done on improving the accuracy of numerical calculations. Although there is much research in the field of FDM, the development of numerical calculations by the spectral method is decisive in improving calculation accuracy. Numerical calculations are usually performed by double precision calculations. If double precision calculations ensure 15 significant digits in floating point computing, such numerical calculations will be the ultimate goal to reach. A numerical calculation that does not seem to have an error even though it originally has an error is defined as a virtual error-zero (VE0) calculation. In this paper, we will examine an overall picture of high-accuracy numerical calculation by IFDM in the numerical calculation of the 1D Poisson equation. It becomes clear that a VE0 calculation is always possible in the numerical calculation method defined as the compact interpolation finite difference scheme (CIFD(m)).