dynamically using LDL, LU or Cholesky solver in dense matrix solve to speed up calculations (relates #194)