I am looking into leveraging GPU's for two body orbit propagation. I am doing basic Kepler propagation, and only dealing with circular and elliptical orbits. I took a first shot at using single precision (8x speed for GPU's) and noticed some significant error in the orbit.
I wanted to know if it is known to be a problem dealing with single precision, or should it be ok? I am only concerned with error being in the region of 1km or so.
The precision error I am talking about is not to do with integration, like was mentioned I am using the closed form analytical method (very similar to the linked method) and using Newton-Rhapson integration technique for elliptical orbits.
I ended up finding out that the precision issue i was having was when calculating GMST and using that for the conversion from ECI to ECEF. I even tried using MJD, instead of JD, but it turns out if I just pre-compute the GMST for each time step using double precision, that solves my issue.