Code beautification with uncrustify
[alexxy/gromacs.git] / src / tools / polynomials.c
index 310d365365be31b75e80199a376eb5cd4f304f18..b311db86f6b0b48cabd9009add83358edcbdca44 100644 (file)
@@ -1,11 +1,11 @@
 /*
- * 
+ *
  *                This source code is part of
- * 
+ *
  *                 G   R   O   M   A   C   S
- * 
+ *
  *          GROningen MAchine for Chemical Simulations
- * 
+ *
  *                        VERSION 3.2.0
  * Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
  * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
  * modify it under the terms of the GNU General Public License
  * as published by the Free Software Foundation; either version 2
  * of the License, or (at your option) any later version.
- * 
+ *
  * If you want to redistribute modifications, please consider that
  * scientific software is very special. Version control is crucial -
  * bugs must be traceable. We will be happy to consider code for
  * inclusion in the official distribution, but derived work must not
  * be called official GROMACS. Details are found in the README & COPYING
  * files - if they are missing, get the official version at www.gromacs.org.
- * 
+ *
  * To help us fund GROMACS development, we humbly ask that you cite
  * the papers on the package - you can find them in the top README file.
- * 
+ *
  * For more info, check our website at http://www.gromacs.org
- * 
+ *
  * And Hey:
  * Green Red Orange Magenta Azure Cyan Skyblue
  */
 #include "gmx_fatal.h"
 #include "gstat.h"
 
-real LegendreP(real x,unsigned long m)
+real LegendreP(real x, unsigned long m)
 
 {
-  real polynomial=0,x2,x3;
+    real polynomial = 0, x2, x3;
 
-  switch (m) {
-  case eacP0:
-    polynomial=1.0;
-    break;
-  case eacP1:
-    polynomial=x;
-    break;
-  case eacP2:
-    x2=x*x;
-    polynomial=1.5*x2 - 0.5;
-    break;
-  case eacP3:
-    x2=x*x;
-    polynomial=(35*x2*x2 - 30*x2 + 3)/8;
-    break;
-  case eacP4:
-    x2=x*x;
-    x3=x2*x;
-    polynomial=(63*x3*x2 - 70*x3 + 15*x)/8;
-    break;
-  default:
-    gmx_fatal(FARGS,"Legendre polynomials of order %d are not supported, %s %d",
-               m,__FILE__,__LINE__);
-  }
-  return (polynomial);
+    switch (m)
+    {
+        case eacP0:
+            polynomial = 1.0;
+            break;
+        case eacP1:
+            polynomial = x;
+            break;
+        case eacP2:
+            x2         = x*x;
+            polynomial = 1.5*x2 - 0.5;
+            break;
+        case eacP3:
+            x2         = x*x;
+            polynomial = (35*x2*x2 - 30*x2 + 3)/8;
+            break;
+        case eacP4:
+            x2         = x*x;
+            x3         = x2*x;
+            polynomial = (63*x3*x2 - 70*x3 + 15*x)/8;
+            break;
+        default:
+            gmx_fatal(FARGS, "Legendre polynomials of order %d are not supported, %s %d",
+                      m, __FILE__, __LINE__);
+    }
+    return (polynomial);
 }