             inv_quat_multiply(orientation, conjqm2, mbody);
             transpose_dot(ex_space, ey_space, ez_space, mbody, angmom2);
             
-            angmom2.x *= 0.5;
-            angmom2.y *= 0.5;
-            angmom2.z *= 0.5;
-            angmom2.w = 0.0;
+            angmom2.x *= 0.5f;
+            angmom2.y *= 0.5f;
+            angmom2.z *= 0.5f;
+            angmom2.w = 0.0f;
             
             // update angular velocity
             float4 angvel2;

            vel2.x = scale_t * vel.x + dtfm * force.x;
            vel2.y = scale_t * vel.y + dtfm * force.y;
            vel2.z = scale_t * vel.z + dtfm * force.z;
            vel2.w = 0.0;
            
            tmp = vel2.x * vel2.x + vel2.y * vel2.y + vel2.z * vel2.z;
            akin_t = body_mass * tmp;

    float x2, x4;
    x2 = x * x;
    x4 = x2 * x2;
    return (1.0f + (1.0f/6.0f) * x2 + (1.0f/120.0f) * x4 + (1.0f/5040.0f) * x2 * x4 + (1.0f/362880.0f) * x4 * x4);
    }

/*! Kernel to zero virial contribution from particles from rigid bodies
