  We propose a geometric numerical analysis of SDEs admitting Lie symmetries
which allows us to individuate a symmetry adapted coordinates system where the
given SDE has notable invariant properties. An approximation scheme preserving
the symmetry properties of the equation is introduced. Our algorithmic
procedure is applied to the family of general linear SDEs for which two
theoretical estimates of the numerical forward error are established.
