2 * This file is part of the GROMACS molecular simulation package.
4 * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
5 * Copyright (c) 2001-2004, The GROMACS development team.
6 * Copyright (c) 2013,2014,2015,2016,2017 by the GROMACS development team.
7 * Copyright (c) 2018,2019,2020, by the GROMACS development team, led by
8 * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
9 * and including many others, as listed in the AUTHORS file in the
10 * top-level source directory and at http://www.gromacs.org.
12 * GROMACS is free software; you can redistribute it and/or
13 * modify it under the terms of the GNU Lesser General Public License
14 * as published by the Free Software Foundation; either version 2.1
15 * of the License, or (at your option) any later version.
17 * GROMACS is distributed in the hope that it will be useful,
18 * but WITHOUT ANY WARRANTY; without even the implied warranty of
19 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
20 * Lesser General Public License for more details.
22 * You should have received a copy of the GNU Lesser General Public
23 * License along with GROMACS; if not, see
24 * http://www.gnu.org/licenses, or write to the Free Software Foundation,
25 * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
27 * If you want to redistribute modifications to GROMACS, please
28 * consider that scientific software is very special. Version
29 * control is crucial - bugs must be traceable. We will be happy to
30 * consider code for inclusion in the official distribution, but
31 * derived work must not be called official GROMACS. Details are found
32 * in the README & COPYING files - if they are missing, get the
33 * official version at http://www.gromacs.org.
35 * To help us fund GROMACS development, we humbly ask that you cite
36 * the research papers on the package. Check out http://www.gromacs.org.
38 /* This file is completely threadsafe - keep it that way! */
48 #include "gromacs/utility/strconvert.h"
50 void cmp_int(FILE* fp, const char* s, int index, int i1, int i2)
56 fprintf(fp, "%s[%d] (%d - %d)\n", s, index, i1, i2);
60 fprintf(fp, "%s (%d - %d)\n", s, i1, i2);
65 void cmp_int64(FILE* fp, const char* s, int64_t i1, int64_t i2)
69 fprintf(fp, "%s (", s);
70 fprintf(fp, "%" PRId64, i1);
72 fprintf(fp, "%" PRId64, i2);
77 void cmp_us(FILE* fp, const char* s, int index, unsigned short i1, unsigned short i2)
83 fprintf(fp, "%s[%d] (%hu - %hu)\n", s, index, i1, i2);
87 fprintf(fp, "%s (%hu - %hu)\n", s, i1, i2);
92 void cmp_uc(FILE* fp, const char* s, int index, unsigned char i1, unsigned char i2)
98 fprintf(fp, "%s[%d] (%d - %d)\n", s, index, int{ i1 }, int{ i2 });
102 fprintf(fp, "%s (%d - %d)\n", s, int{ i1 }, int{ i2 });
107 gmx_bool cmp_bool(FILE* fp, const char* s, int index, gmx_bool b1, gmx_bool b2)
113 fprintf(fp, "%s[%d] (%s - %s)\n", s, index, gmx::boolToString(b1), gmx::boolToString(b2));
117 fprintf(fp, "%s (%s - %s)\n", s, gmx::boolToString(b1), gmx::boolToString(b2));
123 void cmp_str(FILE* fp, const char* s, int index, const char* s1, const char* s2)
125 if (std::strcmp(s1, s2) != 0)
129 fprintf(fp, "%s[%d] (%s - %s)\n", s, index, s1, s2);
133 fprintf(fp, "%s (%s - %s)\n", s, s1, s2);
138 gmx_bool equal_real(real i1, real i2, real ftol, real abstol)
140 return ((2 * std::fabs(i1 - i2) <= (fabs(i1) + fabs(i2)) * ftol) || std::fabs(i1 - i2) <= abstol);
143 gmx_bool equal_float(float i1, float i2, float ftol, float abstol)
145 return ((2 * std::fabs(i1 - i2) <= (std::fabs(i1) + std::fabs(i2)) * ftol)
146 || std::fabs(i1 - i2) <= abstol);
149 gmx_bool equal_double(double i1, double i2, real ftol, real abstol)
151 return ((2 * fabs(i1 - i2) <= (fabs(i1) + fabs(i2)) * ftol) || fabs(i1 - i2) <= abstol);
154 void cmp_real(FILE* fp, const char* s, int index, real i1, real i2, real ftol, real abstol)
156 if (!equal_real(i1, i2, ftol, abstol))
160 fprintf(fp, "%s[%2d] (%e - %e)\n", s, index, i1, i2);
164 fprintf(fp, "%s (%e - %e)\n", s, i1, i2);
169 void cmp_float(FILE* fp, const char* s, int index, float i1, float i2, float ftol, float abstol)
171 if (!equal_float(i1, i2, ftol, abstol))
175 fprintf(fp, "%s[%2d] (%e - %e)\n", s, index, i1, i2);
179 fprintf(fp, "%s (%e - %e)\n", s, i1, i2);
184 void cmp_double(FILE* fp, const char* s, int index, double i1, double i2, double ftol, double abstol)
186 if (!equal_double(i1, i2, ftol, abstol))
190 fprintf(fp, "%s[%2d] (%16.9e - %16.9e)\n", s, index, i1, i2);
194 fprintf(fp, "%s (%16.9e - %16.9e)\n", s, i1, i2);