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, led by
7 * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
8 * and including many others, as listed in the AUTHORS file in the
9 * top-level source directory and at http://www.gromacs.org.
11 * GROMACS is free software; you can redistribute it and/or
12 * modify it under the terms of the GNU Lesser General Public License
13 * as published by the Free Software Foundation; either version 2.1
14 * of the License, or (at your option) any later version.
16 * GROMACS is distributed in the hope that it will be useful,
17 * but WITHOUT ANY WARRANTY; without even the implied warranty of
18 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
19 * Lesser General Public License for more details.
21 * You should have received a copy of the GNU Lesser General Public
22 * License along with GROMACS; if not, see
23 * http://www.gnu.org/licenses, or write to the Free Software Foundation,
24 * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
26 * If you want to redistribute modifications to GROMACS, please
27 * consider that scientific software is very special. Version
28 * control is crucial - bugs must be traceable. We will be happy to
29 * consider code for inclusion in the official distribution, but
30 * derived work must not be called official GROMACS. Details are found
31 * in the README & COPYING files - if they are missing, get the
32 * official version at http://www.gromacs.org.
34 * To help us fund GROMACS development, we humbly ask that you cite
35 * the research papers on the package. Check out http://www.gromacs.org.
41 #include "gromacs/utility/exceptions.h"
42 #include "gromacs/utility/programcontext.h"
43 #include "gromacs/utility/smalloc.h"
44 #include "gromacs/utility/stringutil.h"
46 struct gmx_output_env_t
48 explicit gmx_output_env_t(const gmx::IProgramContext
&context
)
49 : programContext(context
),
54 trajectory_io_verbosity(0) {}
57 const gmx::IProgramContext
&programContext
;
59 /* the time unit, enum defined in oenv.h */
60 time_unit_t time_unit
;
61 /* view of file requested */
63 /* xvg output format, enum defined in oenv.h */
64 xvg_format_t xvg_format
;
65 /* The level of verbosity for this program */
67 /* The level of verbosity during trajectory I/O. Default=1, quiet=0. */
68 int trajectory_io_verbosity
;
71 /* The source code in this file should be thread-safe.
72 Please keep it that way. */
74 /******************************************************************
76 * T R A J E C T O R Y S T U F F
78 ******************************************************************/
80 /* read only time names */
81 /* These must correspond to the time units type time_unit_t in oenv.h */
82 static const real timefactors
[] = { real(0), real(1e3
), real(1), real(1e-3), real(1e-6), real(1e-9), real(1e-12), real(0) };
83 static const real timeinvfactors
[] = { real(0), real(1e-3), real(1), real(1e3
), real(1e6
), real(1e9
), real(1e12
), real(0) };
84 static const char *time_units_str
[] = {
85 nullptr, "fs", "ps", "ns", "us",
88 static const char *time_units_xvgr
[] = {
89 nullptr, "fs", "ps", "ns",
94 /***** OUTPUT_ENV MEMBER FUNCTIONS ******/
96 void output_env_init(gmx_output_env_t
**oenvp
,
97 const gmx::IProgramContext
&context
,
98 time_unit_t tmu
, gmx_bool view
, xvg_format_t xvg_format
,
103 gmx_output_env_t
*oenv
= new gmx_output_env_t(context
);
105 oenv
->time_unit
= tmu
;
107 oenv
->xvg_format
= xvg_format
;
108 oenv
->verbosity
= verbosity
;
109 const char *env
= getenv("GMX_TRAJECTORY_IO_VERBOSITY");
110 oenv
->trajectory_io_verbosity
= (env
!= nullptr ? strtol(env
, NULL
, 10) : 1);
113 GMX_CATCH_ALL_AND_EXIT_WITH_FATAL_ERROR
;
116 void output_env_init_default(gmx_output_env_t
**oenvp
)
120 gmx_output_env_t
*oenv
= new gmx_output_env_t(gmx::getProgramContext());
123 GMX_CATCH_ALL_AND_EXIT_WITH_FATAL_ERROR
;
126 void output_env_done(gmx_output_env_t
*oenv
)
132 int output_env_get_verbosity(const gmx_output_env_t
*oenv
)
134 return oenv
->verbosity
;
137 int output_env_get_trajectory_io_verbosity(const gmx_output_env_t
*oenv
)
139 return oenv
->trajectory_io_verbosity
;
142 std::string
output_env_get_time_unit(const gmx_output_env_t
*oenv
)
144 return time_units_str
[oenv
->time_unit
];
147 std::string
output_env_get_time_label(const gmx_output_env_t
*oenv
)
149 return gmx::formatString("Time (%s)", time_units_str
[oenv
->time_unit
] ?
150 time_units_str
[oenv
->time_unit
] : "ps");
153 std::string
output_env_get_xvgr_tlabel(const gmx_output_env_t
*oenv
)
155 return gmx::formatString("Time (%s)", time_units_xvgr
[oenv
->time_unit
] ?
156 time_units_xvgr
[oenv
->time_unit
] : "ps");
159 real
output_env_get_time_factor(const gmx_output_env_t
*oenv
)
161 return timefactors
[oenv
->time_unit
];
164 real
output_env_get_time_invfactor(const gmx_output_env_t
*oenv
)
166 return timeinvfactors
[oenv
->time_unit
];
169 real
output_env_conv_time(const gmx_output_env_t
*oenv
, real time
)
171 return time
*timefactors
[oenv
->time_unit
];
174 void output_env_conv_times(const gmx_output_env_t
*oenv
, int n
, real
*time
)
177 double fact
= timefactors
[oenv
->time_unit
];
181 for (i
= 0; i
< n
; i
++)
188 gmx_bool
output_env_get_view(const gmx_output_env_t
*oenv
)
193 xvg_format_t
output_env_get_xvg_format(const gmx_output_env_t
*oenv
)
195 return oenv
->xvg_format
;
198 const char *output_env_get_program_display_name(const gmx_output_env_t
*oenv
)
200 const char *displayName
= nullptr;
204 displayName
= oenv
->programContext
.displayName();
206 GMX_CATCH_ALL_AND_EXIT_WITH_FATAL_ERROR
;
211 const gmx::IProgramContext
&
212 output_env_get_program_context(const gmx_output_env_t
*oenv
)
214 return oenv
->programContext
;