Merge release-5-0 into release-5-1
[alexxy/gromacs.git] / src / gromacs / mdlib / perf_est.c
index 182ff35cd752b4a5a110bde75fbf30e9a1334de3..0951875dd3475965f04d473e4ae9362045d896f4 100644 (file)
@@ -3,7 +3,7 @@
  *
  * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
  * Copyright (c) 2001-2008, The GROMACS development team.
- * Copyright (c) 2012,2014, by the GROMACS development team, led by
+ * Copyright (c) 2012,2014,2015, by the GROMACS development team, led by
  * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
  * and including many others, as listed in the AUTHORS file in the
  * top-level source directory and at http://www.gromacs.org.
  * To help us fund GROMACS development, we humbly ask that you cite
  * the research papers on the package. Check out http://www.gromacs.org.
  */
-#include "config.h"
+#include "gmxpre.h"
+
+#include "gromacs/legacyheaders/perf_est.h"
 
 #include <math.h>
 
+#include "gromacs/legacyheaders/types/commrec.h"
 #include "gromacs/math/vec.h"
+#include "gromacs/mdlib/nbnxn_consts.h"
+#include "gromacs/mdlib/nbnxn_search.h"
 #include "gromacs/topology/topology.h"
 #include "gromacs/utility/fatalerror.h"
 
-#include "gromacs/legacyheaders/perf_est.h"
-#include "gromacs/legacyheaders/types/commrec.h"
-#include "nbnxn_search.h"
-#include "nbnxn_consts.h"
-
 /* Computational cost of bonded, non-bonded and PME calculations.
  * This will be machine dependent.
  * The numbers here are accurate for Intel Core2 and AMD Athlon 64
@@ -99,7 +99,7 @@ int n_bonded_dx(gmx_mtop_t *mtop, gmx_bool bExcl)
      */
     ndx      = 0;
     ndx_excl = 0;
-#if __ICC == 1400 || __ICL == 1400
+#if defined _ICC && __ICC == 1400 || defined __ICL && __ICL == 1400
 #pragma novector /* Work-around for incorrect vectorization */
 #endif
     for (mb = 0; mb < mtop->nmolblock; mb++)
@@ -178,6 +178,7 @@ static void pp_group_load(gmx_mtop_t *mtop, t_inputrec *ir, matrix box,
     nq                = 0;
     nlj               = 0;
     *bChargePerturbed = FALSE;
+    *bTypePerturbed   = FALSE;
     for (mb = 0; mb < mtop->nmolblock; mb++)
     {
         molt = &mtop->moltype[mtop->molblock[mb].type];
@@ -292,6 +293,7 @@ static void pp_verlet_load(gmx_mtop_t *mtop, t_inputrec *ir, matrix box,
     nqlj              = 0;
     nq                = 0;
     *bChargePerturbed = FALSE;
+    *bTypePerturbed   = FALSE;
     for (mb = 0; mb < mtop->nmolblock; mb++)
     {
         molt = &mtop->moltype[mtop->molblock[mb].type];