Gromacs
2018.8
|
#include <cstdio>
#include "gromacs/topology/ifunc.h"
Declares functions for handling orientation restraints.
Functions | |
void | init_orires (FILE *fplog, const gmx_mtop_t *mtop, const t_inputrec *ir, const t_commrec *cr, t_state *globalState, t_oriresdata *od) |
Decides whether orientation restraints can work, and initializes all the orientation restraint stuff in *od (and assumes *od is already allocated. If orientation restraint are used, globalState is read and modified on the master rank (which is the only rank, since orientation restraints can not run in parallel). | |
real | calc_orires_dev (const gmx_multisim_t *ms, int nfa, const t_iatom fa[], const t_iparams ip[], const t_mdatoms *md, const rvec x[], const t_pbc *pbc, t_fcdata *fcd, history_t *hist) |
Calculates the time averaged D matrices, the S matrix for each experiment. More... | |
void | diagonalize_orires_tensors (t_oriresdata *od) |
Diagonalizes the order tensor(s) of the orienation restraints. More... | |
void | print_orires_log (FILE *log, t_oriresdata *od) |
Prints order parameter, eigenvalues and eigenvectors to the log file. | |
void | update_orires_history (t_fcdata *fcd, history_t *hist) |
Copies the new time averages that have been calculated in calc_orires_dev. | |
Variables | |
t_ifunc | orires |
Calculates the orientation restraint forces. | |
real calc_orires_dev | ( | const gmx_multisim_t * | ms, |
int | nfa, | ||
const t_iatom | fa[], | ||
const t_iparams | ip[], | ||
const t_mdatoms * | md, | ||
const rvec | x[], | ||
const t_pbc * | pbc, | ||
t_fcdata * | fcd, | ||
history_t * | hist | ||
) |
Calculates the time averaged D matrices, the S matrix for each experiment.
Returns the weighted RMS deviation of the orientation restraints.
void diagonalize_orires_tensors | ( | t_oriresdata * | od | ) |
Diagonalizes the order tensor(s) of the orienation restraints.
For each experiment eig containts first 3 eigenvalues and then the 3 eigenvectors. The eigenvalues are ordered on magnitude.