char PullName[4][256];
double UmbPos[4][3];
double UmbCons[4][3];
- bool Flipped[4];
+ gmx_bool Flipped[4];
} t_UmbrellaHeader;
typedef struct
double *pos;
double *z;
double * N, *Ntot;
- bool * Flipped;
+ gmx_bool * Flipped;
double dt;
- bool **bContrib;
+ gmx_bool **bContrib;
} t_UmbrellaWindow;
typedef struct
{
const char *fnTpr,*fnPullf,*fnPdo,*fnPullx;
- bool bTpr,bPullf,bPdo,bPullx;
+ gmx_bool bTpr,bPullf,bPdo,bPullx;
int bins,cycl;
- bool verbose,bShift,bAuto,bBoundsOnly;
- bool bFlipProf;
+ gmx_bool verbose,bShift,bAuto,bBoundsOnly;
+ gmx_bool bFlipProf;
real tmin, tmax, dt;
real Temperature,Tolerance;
int nBootStrap,histBootStrapBlockLength;
real dtBootStrap,zProfZero,alpha;
int bsSeed,stepchange;
- bool bHistBootStrap,bWeightedCycl,bHistOutOnly;
- bool bAutobounds,bNoprof;
+ gmx_bool bHistBootStrap,bWeightedCycl,bHistOutOnly;
+ gmx_bool bAutobounds,bNoprof;
real min,max,dz;
- bool bLog;
+ gmx_bool bLog;
int unit;
real zProf0;
- bool bProf0Set,bs_verbose;
- bool bHistEq, bTab;
+ gmx_bool bProf0Set,bs_verbose;
+ gmx_bool bHistEq, bTab;
double *tabX,*tabY,tabMin,tabMax,tabDz;
int tabNbins;
} t_UmbrellaOptions;
void read_pdo_data(FILE * file, t_UmbrellaHeader * header,
int fileno, t_UmbrellaWindow * win,
t_UmbrellaOptions *opt,
- bool bGetMinMax,real *mintmp,real *maxtmp)
+ gmx_bool bGetMinMax,real *mintmp,real *maxtmp)
{
int i,inttemp,bins,count;
real min,max,minfound,maxfound;
double temp,time,time0=0,dt;
char *ptr;
t_UmbrellaWindow * window=0;
- bool timeok,dt_ok=1;
+ gmx_bool timeok,dt_ok=1;
char *tmpbuf,fmt[256],fmtign[256];
int len=STRLEN,dstep=1;
{
int i,j,k;
double U,min=opt->min,dz=opt->dz,temp,ztot_half,distance,ztot,contrib;
- bool bAnyContrib;
+ gmx_bool bAnyContrib;
ztot=opt->max-opt->min;
void calc_profile(double *profile,t_UmbrellaWindow * window, int nWindows, t_UmbrellaOptions *opt,
- bool bExact)
+ gmx_bool bExact)
{
int i,k,j;
double num,ztot_half,ztot,distance,min=opt->min,dz=opt->dz;
double calc_z(double * profile,t_UmbrellaWindow * window, int nWindows, t_UmbrellaOptions *opt,
- bool bExact)
+ gmx_bool bExact)
{
int i,j,k;
double U=0,min=opt->min,dz=opt->dz,temp,ztot_half,distance,ztot;
void cyclicProfByWeightedCorr(double *profile,t_UmbrellaWindow *window,
int nWindows, t_UmbrellaOptions * opt,
- bool bAppendCorr2File, const char *fn,
+ gmx_bool bAppendCorr2File, const char *fn,
const output_env_t oenv)
{
int i,j,k,bins=opt->bins;
{
int nsynth,N,i,nbins,r_index;
double r;
- static bool bWarnout=0;
+ static gmx_bool bWarnout=0;
N=thisWindow->N[pullid];
int i,j,*randomArray=0,winid,pullid,ib;
int iAllPull,nAllPull,*allPull_winId,*allPull_pullId;
FILE *fp;
- bool bExact=FALSE;
+ gmx_bool bExact=FALSE;
/* init random */
void read_pull_xf(const char *fn, const char *fntpr,
t_UmbrellaHeader * header, t_UmbrellaWindow * window,
- t_UmbrellaOptions *opt, bool bGetMinMax,real *mintmp,
+ t_UmbrellaOptions *opt, gmx_bool bGetMinMax,real *mintmp,
real *maxtmp)
{
double **y,pos=0.,t,force,time0=0.,dt;
int ny,nt,bins,ibin,i,g,dstep=1,nColPerGrp,nColRef,nColExpect;
real min,max,minfound,maxfound;
- bool dt_ok,timeok,bHaveForce;
+ gmx_bool dt_ok,timeok,bHaveForce;
const char *quantity;
minfound=1e20;
};
static t_UmbrellaOptions opt;
- static bool bHistOnly=FALSE;
+ static gmx_bool bHistOnly=FALSE;
const char *en_unit[]={NULL,"kJ","kCal","kT",NULL};
const char *en_unit_label[]={"","E (kJ mol\\S-1\\N)","E (kcal mol\\S-1\\N)","E (kT)",};
t_UmbrellaHeader header;
t_UmbrellaWindow * window=NULL;
double *profile,maxchange=1e20;
- bool bMinSet,bMaxSet,bAutoSet,bExact=FALSE;
+ gmx_bool bMinSet,bMaxSet,bAutoSet,bExact=FALSE;
char **fninTpr,**fninPull,**fninPdo;
const char *fnPull;
FILE *histout,*profout;