From 4040308baafbee50c4376960b99735e1901f8462 Mon Sep 17 00:00:00 2001
From: henry <henry>
Date: Tue, 27 Jan 2009 21:17:41 +0000
Subject: [PATCH] Simplified the pow-of-2 check.  See
 http://openfoam.cfd-online.com/cgi-bin/forum/show.cgi?126/10819

---
 src/randomProcesses/fft/fft.C | 5 ++---
 1 file changed, 2 insertions(+), 3 deletions(-)

diff --git a/src/randomProcesses/fft/fft.C b/src/randomProcesses/fft/fft.C
index 4531c52..9509e97 100644
--- a/src/randomProcesses/fft/fft.C
+++ b/src/randomProcesses/fft/fft.C
@@ -48,9 +48,8 @@ void fft::transform
 {
     forAll(nn, idim)
     {
-        // Check for power of two
-        unsigned int dimCount = nn[idim];
-        if (!dimCount || (dimCount & (dimCount - 1)))
+        scalar pow2 = log(scalar(nn[idim]))/log(scalar(2));
+        if ((pow2 - int(pow2 + 0.5)) > SMALL)
         {
             FatalErrorIn
             (
-- 
2.11.4.GIT