#include <math.h>
#include <string.h>
+#include <assert.h>
#include "sysstuff.h"
#include "typedefs.h"
#include "macros.h"
case eelP3M_AD:
if (cr->duty & DUTY_PME)
{
+ assert(fr->n_tpi >= 0);
if (fr->n_tpi == 0 || (flags & GMX_FORCE_STATECHANGED))
{
pme_flags = GMX_PME_SPREAD_Q | GMX_PME_SOLVE;