The purpose of our work is to identify the key features of algorithms in porting codes for calculating of essentially nonlinear processes to a modern cluster of hybrid architecture. The KPI equation with a source term was numerically modeled to reveal the features of the occurrence of extreme waves. A methodology for the implementation of boundary conditions for the modified KPI equation was also proposed.