The HOC-Krylov Solver is a Fortran-based computational tool designed to efficiently solve Poisson’s equation using high-order compact (HOC) finite difference schemes integrated with Krylov subspace iterative methods. This combination enhances both the accuracy and convergence speed of solutions in two-dimensional domains. The solver is particularly suited for applications in computational fluid dynamics and electrostatics, where precise and rapid solutions are essential. • Integrates high-order compact finite difference schemes with multiple Krylov solvers to deliver highly accurate and efficient solutions of large Poisson systems. • Supports flexible sparse matrix storage (CSR, MSR, ELLPACK-ITPACK, DIA) with optimized MVM routines, enabling systematic performance evaluation across solver–storage– preconditioner combinations. • Modular Fortran architecture enables easy extension, facilitating research in CFD, electrostatics, and vorticity–streamfunction formulations while providing benchmarking and post-processing capabilities.