-parallel option causes differing answers with different cores. Problem still exist but now mitigated
fix numerical instability; pcxa was not declared private and was being modified by the threads