Test Polak-Ribire form conjugate gradient method (of_method = cg1) in OFDFT, symmetry=on
