Precompute LU factors and add to intermediates.
Precompute LU factors and add to intermediates.

- This speeds up repeated calls of sht with the same intermediates,
  since the P array does not change across calls.

praveenv253/sht,praveenv253/sht