added array dot prod, L2 norm, sum and difference function, also addeded the option of passing a force array to the integrator as an optional argument