Viewing file:  pr58673-2.c (4.79 KB) -rw-r--r--
Select action/file-type:
 (+) |  (+) |  (+) | Code (+) | Session (+) |  (+) | SDB (+) |  (+) |  (+) |  (+) |  (+) |  (+) |
/* { dg-do compile { target { powerpc*-*-* && lp64 } } } */
/* { dg-skip-if "" { powerpc*-*-darwin* } } */
/* { dg-require-effective-target powerpc_p8vector_ok } */
/* { dg-options "-mdejagnu-cpu=power8 -O3 -funroll-loops" } */
#include <stddef.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
typedef long unsigned int size_t;
typedef struct _IO_FILE FILE;
typedef float real;
typedef real rvec[3];
typedef real matrix[3][3];
typedef real tensor[3][3];
enum
{
F_BONDS, F_G96BONDS, F_MORSE, F_CUBICBONDS, F_CONNBONDS, F_HARMONIC,
F_ANGLES, F_G96ANGLES, F_PDIHS, F_RBDIHS, F_IDIHS, F_LJ14, F_COUL14, F_LJ,
F_BHAM, F_LJLR, F_DISPCORR, F_SR, F_LR, F_WPOL, F_POSRES, F_DISRES,
F_DISRESVIOL, F_ORIRES, F_ORIRESDEV, F_ANGRES, F_ANGRESZ, F_SHAKE,
F_SHAKENC, F_SETTLE, F_DUMMY2, F_DUMMY3, F_DUMMY3FD, F_DUMMY3FAD,
F_DUMMY3OUT, F_DUMMY4FD, F_EQM, F_EPOT, F_EKIN, F_ETOT, F_TEMP, F_PRES,
F_DVDL, F_DVDLKIN, F_NRE
};
typedef union
{
struct
{
}
bham;
struct
{
real rA, krA, rB, krB;
}
harmonic;
}
t_iparams;
typedef struct
{
t_iparams *iparams;
}
t_idef;
typedef struct
{
}
t_inputrec;
typedef struct
{
}
t_commrec;
typedef struct
{
}
t_forcerec;
typedef struct
{
}
t_mdatoms;
typedef struct
{
}
t_filenm;
enum
{
eoPres, eoEpot, eoVir, eoDist, eoMu, eoForce, eoFx, eoFy, eoFz, eoPx, eoPy,
eoPz, eoPolarizability, eoDipole, eoObsNR, eoMemory =
eoObsNR, eoInter, eoUseVirial, eoNR
};
extern char *eoNames[eoNR];
typedef struct
{
int bPrint;
}
t_coupl_LJ;
typedef struct
{
int eObs;
t_iparams xi;
}
t_coupl_iparams;
typedef struct
{
real act_value[eoObsNR];
real av_value[eoObsNR];
real ref_value[eoObsNR];
int bObsUsed[eoObsNR];
int nLJ, nBU, nQ, nIP;
t_coupl_LJ *tcLJ;
}
t_coupl_rec;
void xvgr_legend ();
real calc_deviation ();
void pr_dev ();
static void
pr_ff (t_coupl_rec * tcr, real time, t_idef * idef, t_commrec * cr, int nfile,
t_filenm fnm[])
{
static FILE *prop;
static FILE **out = ((void *) 0);
static FILE **qq = ((void *) 0);
static FILE **ip = ((void *) 0);
char buf[256];
char *leg[] = {
"C12", "C6"
};
char **raleg;
int i, j, index;
if ((prop == ((void *) 0)) && (out == ((void *) 0)) && (qq == ((void *) 0))
&& (ip == ((void *) 0)))
{
for (i = j = 0; (i < eoObsNR); i++)
{
if (tcr->bObsUsed[i])
{
raleg[j++] =
(__extension__
(__builtin_constant_p (eoNames[i])
&& ((size_t) (const void *) ((eoNames[i]) + 1) -
(size_t) (const void *) (eoNames[i]) ==
1) ? (((const char *) (eoNames[i]))[0] ==
'\0' ? (char *) calloc ((size_t) 1,
(size_t) 1) : (
{
size_t
__len
=
strlen
(eoNames
[i])
+
1;
char
*__retval
=
(char
*)
malloc
(__len);
__retval;}
)): strdup (eoNames[i])));
raleg[j++] =
(__extension__
(__builtin_constant_p (buf)
&& ((size_t) (const void *) ((buf) + 1) -
(size_t) (const void *) (buf) ==
1) ? (((const char *) (buf))[0] ==
'\0' ? (char *) calloc ((size_t) 1,
(size_t) 1) : (
{
size_t
__len
=
strlen
(buf)
+
1;
char
*__retval
=
(char
*)
malloc
(__len);
__retval;}
)): strdup (buf)));
}
}
if (tcr->nLJ)
{
for (i = 0; (i < tcr->nLJ); i++)
{
if (tcr->tcLJ[i].bPrint)
{
xvgr_legend (out[i], (sizeof (leg) / sizeof ((leg)[0])),
leg);
}
}
}
}
}
void
do_coupling (FILE * log, int nfile, t_filenm fnm[], t_coupl_rec * tcr, real t,
int step, real ener[], t_forcerec * fr, t_inputrec * ir,
int bMaster, t_mdatoms * md, t_idef * idef, real mu_aver,
int nmols, t_commrec * cr, matrix box, tensor virial,
tensor pres, rvec mu_tot, rvec x[], rvec f[], int bDoIt)
{
int i, j, ati, atj, atnr2, type, ftype;
real deviation[eoObsNR], prdev[eoObsNR], epot0, dist, rmsf;
real ff6, ff12, ffa, ffb, ffc, ffq, factor, dt, mu_ind;
int bTest, bPrint;
t_coupl_iparams *tip;
if (bPrint)
{
pr_ff (tcr, t, idef, cr, nfile, fnm);
}
for (i = 0; (i < eoObsNR); i++)
{
deviation[i] =
calc_deviation (tcr->av_value[i], tcr->act_value[i],
tcr->ref_value[i]);
prdev[i] = tcr->ref_value[i] - tcr->act_value[i];
}
if (bPrint)
pr_dev (tcr, t, prdev, cr, nfile, fnm);
for (i = 0; (i < atnr2); i++)
{
factor = dt * deviation[tip->eObs];
switch (ftype)
{
case F_BONDS:
if (fabs (tip->xi.harmonic.krA) > 1.2e-38)
idef->iparams[type].harmonic.krA *=
(1 + factor / tip->xi.harmonic.krA);
}
}
}
|
