return (m[XX][XX]+m[YY][YY]+m[ZZ][ZZ]);
}
-static inline real _divide(real a,real b,const char *file,int line)
+static inline real _divide_err(real a,real b,const char *file,int line)
{
if (fabs(b) <= GMX_REAL_MIN)
gmx_fatal(FARGS,"Dividing by zero, file %s, line %d",file,line);
-box[ZZ][XX]*box[ZZ][XX]-box[ZZ][YY]*box[ZZ][YY]);
}
-#define divide(a,b) _divide((a),(b),__FILE__,__LINE__)
+#define divide_err(a,b) _divide_err((a),(b),__FILE__,__LINE__)
#define mod(a,b) _mod((a),(b),__FILE__,__LINE__)
#ifdef __cplusplus
{
for(m=0; (m<DIM); m++)
{
- b_inv[m] = divide(1.0,box_size[m]);
+ b_inv[m] = divide_err(1.0,box_size[m]);
}
bTriclinic = TRICLINIC(box);
}