initial import factorization tools

also worth to note is that this method creates a 3rd matrix, its not adding x to acc.  We could make our own addTo method which deals with the underlying arrays directly as this would be faster.  However "it aint broke" at least with 200 dims