diff --git a/src/rotations.f90 b/src/rotations.f90
index cf6f66af8..3f6778e1a 100644
--- a/src/rotations.f90
+++ b/src/rotations.f90
@@ -552,14 +552,8 @@ function om2ax(om) result(ax)
    LWORK = 20   
  
  ! call the eigenvalue solver
-#if (FLOAT==8)
    call dgeev('N','V',3,o,3,Wr,Wi,devNull,3,VR,3,WORK,LWORK,INFO)
-#elif (FLOAT==4)
-   call sgeev('N','V',3,o,3,Wr,Wi,devNull,3,VR,3,WORK,LWORK,INFO)
-#else
-   NO SUITABLE PRECISION FOR REAL SELECTED, STOPPING COMPILATION
-#endif
-   if (INFO /= 0) call IO_error(0_pInt,ext_msg='Error in om2ax/(s/d)geev: (S/D)GEEV return not zero')
+   if (INFO /= 0) call IO_error(0_pInt,ext_msg='Error in om2ax DGEEV return not zero')
    i = maxloc(merge(1.0_pReal,0.0_pReal,cEq(cmplx(Wr,Wi,pReal),cmplx(1.0_pReal,0.0_pReal,pReal),tol=1.0e-14_pReal)),dim=1) ! poor substitute for findloc
    ax(1:3) = VR(1:3,i)
    where (                dNeq0([om(2,3)-om(3,2), om(3,1)-om(1,3), om(1,2)-om(2,1)])) &