more efficient PS1 calculation