#ifndef GMX_MATH_VECTYPES_H
#define GMX_MATH_VECTYPES_H
+#include <cassert>
#include <cmath>
#include <algorithm>
#include <type_traits>
-#include "gromacs/utility/gmxassert.h"
#include "gromacs/utility/real.h"
#define XX 0 /* Defines for indexing in */
//! Allow vector scalar division
BasicVector<ValueType> operator/(const ValueType& right) const
{
- GMX_ASSERT(right != 0, "Cannot divide by zero");
+ assert((right != 0 && "Cannot divide by zero"));
return *this * (1 / right);
}
//! Divide vector by a scalar
BasicVector<ValueType>& operator/=(const ValueType& right)
{
- GMX_ASSERT(right != 0, "Cannot divide by zero");
+ assert((right != 0 && "Cannot divide by zero"));
return *this *= 1 / right;
}
BasicVector<ValueType> unitVector() const
{
const ValueType vectorNorm = norm();
- GMX_ASSERT(vectorNorm != 0, "unitVector() should not be called with a zero vector");
+ assert((vectorNorm != 0 && "unitVector() should not be called with a zero vector"));
return *this / vectorNorm;
}
int nsig,
real* sig,
int* totalNumberOfBondedInteractions,
- gmx_bool bSumEkinhOld,
+ bool bSumEkinhOld,
int flags)
/* instead of current system, gmx_booleans for summing virial, kinetic energy, and other terms */
{