Fixed small bug in pe's HCSITS solver and added overrelaxation parameter
There was a small bug in the orthogonal projection model of the HCSITS solver. I fixed it and made the overrelaxation parameter omega editable by the user. This is necessary because it must be adopted depending on the simulation scenario. I also added some comments and tried to improve readability in general.