adding isotropic/axis_aligned keywords and fixing units of gradient computations
updating optimize to use vector instead of pars
optimize should subscript kernel not vector