  Precise and detailed knowledge of the internal structure of hadrons is one of
the most actual problems in elementary particle physics. In view of the planned
high energy physics facilities, in particular, the Electron-Ion Collider
constructed in Brookhaven National Laboratory, the Chinese Electron-Ion
Collider of China, or upgraded versions of CERN's LHC experiments, it is
important to prepare adequate theoretical tools to compare and correctly
interpret experimental results. One of the model frameworks allowing to
estimate hadron structure functions is the combination of the
McLerran-Venugopalan initial condition model together with the JIMWLK equation
which describes the evolution in rapidity of the initial distribution. In this
package, we present a parallel C++ implementation of both these ingredients. In
order to allow a thorough assessment of systematic effects, several
discretizations of the JIMWLK kernel are implemented both in position and
momentum spaces. The effects of the running coupling in three different
definitions are provided. The main code is supplemented with test and check
programs for all main functionalities. The clear structure of the code allows
easy implementation of further improvements such as the collinear constraint.
