DOUBLE PRECISION FUNCTION KLDSTL(X1, ETA1, X2, ETA2)
#endif
#if !defined(CLEO_DOUBLE)
      REAL FUNCTION KLDSTL(X1, ETA1, X2, ETA2)
#endif
 
C  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
C    Returns distance between two lines.
C
C  X1      real array
C          (x,y,z) of initial point defining line 1
C
C  ETA1    real array
C          direction cosines of line 1 (do not have to be normalized)
C
C  X2      real array
C          (x,y,z) of initial point defining line 2
C
C  ETA2    real array
C          direction cosines of line 2 (do not have to be normalized)
C  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>