anchor.diff:    !! @param[in] InverseSquareRoot of the overlap matrix.
     !! @param[in] nel the number of electrons.
     !! @param[out] Density the density matrix computed by this routine.
  +  !! @param[out] energy_value out the energy of the system (optional). 
     !! @param[out] chemical_potential_out the chemical potential (optional).
     !! @param[in] solver_parameters_in parameters for the solver (optional).
     SUBROUTINE HPCPPlus(Hamiltonian, InverseSquareRoot, nel, Density, &
  
positive.code:   !! @param[in] InverseSquareRoot of the overlap matrix.
    !! @param[in] nel the number of electrons.
    !! @param[out] Density the density matrix computed by this routine.
    !! @param[out] chemical_potential_out the chemical potential (optional).
    !! @param[in] solver_parameters_in parameters for the solver (optional)
    SUBROUTINE TRS2(Hamiltonian, InverseSquareRoot, nel, Density, &
  
negative.code:   !! @param[in] InverseSquareRoot of the overlap matrix.
    !! @param[in] nel the number of electrons.
    !! @param[out] Density the density matrix computed by this routine.
    !! @param[out] energy_value out the energy of the system (optional).
    !! @param[out] chemical_potential_out the chemical potential (optional).
    !! @param[in] solver_parameters_in parameters for the solver (optional).
    SUBROUTINE HPCP(Hamiltonian, InverseSquareRoot, nel, Density, &
  