fprintf(stderr,"Loaded with Money\n\n");
}
- if (PAR(cr) && !((Flags & MD_PARTDEC) || EI_TPI(inputrec->eI)))
+ if (PAR(cr) && !((Flags & MD_PARTDEC) ||
+ EI_TPI(inputrec->eI) ||
+ inputrec->eI == eiNM))
{
cr->dd = init_domain_decomposition(fplog,cr,Flags,ddxyz,rdd,rconstr,
dddlb_opt,dlb_scale,
snew(*enerd,1);
init_enerdata(top_global->groups.grps[egcENER].nr,ir->n_flambda,*enerd);
- /* Init bin for energy stuff */
- *mdebin = init_mdebin((*outf)->fp_ene,top_global,ir,NULL);
+ if (mdebin != NULL)
+ {
+ /* Init bin for energy stuff */
+ *mdebin = init_mdebin((*outf)->fp_ene,top_global,ir,NULL);
+ }
clear_rvec(mu_tot);
calc_shifts(ems->s.box,fr->shift_vec);
unsigned long Flags,
gmx_runtime_t *runtime)
{
- t_mdebin *mdebin;
const char *NM = "Normal Mode Analysis";
gmx_mdoutf_t *outf;
int natoms,atom,d;
state_global,top_global,state_work,&top,
&f,&f_global,
nrnb,mu_tot,fr,&enerd,&graph,mdatoms,&gstat,vsite,constr,
- nfile,fnm,&outf,&mdebin);
+ nfile,fnm,&outf,NULL);
natoms = top_global->natoms;
snew(fneg,natoms);
if (MASTER(cr))
{
- print_ebin(NULL,FALSE,FALSE,FALSE,fplog,atom,t,eprAVER,
- FALSE,mdebin,fcd,&(top_global->groups),&(inputrec->opts));
-
fprintf(stderr,"\n\nWriting Hessian...\n");
gmx_mtxio_write(ftp2fn(efMTX,nfile,fnm),sz,sz,full_matrix,sparse_matrix);
}
sprintf(filemode, bAppendFiles ? "a+" : "w+");
- if (ir->eI != eiNM
+ if ((EI_DYNAMICS(ir->eI) || EI_ENERGY_MINIMIZATION(ir->eI))
#ifndef GMX_FAHCORE
&&
!(EI_DYNAMICS(ir->eI) &&
{
of->fp_trn = open_trn(ftp2fn(efTRN,nfile,fnm), filemode);
}
- if (!EI_ENERGY_MINIMIZATION(ir->eI) &&
+ if (EI_DYNAMICS(ir->eI) &&
ir->nstxtcout > 0)
{
of->fp_xtc = open_xtc(ftp2fn(efXTC,nfile,fnm), filemode);
of->xtc_prec = ir->xtcprec;
}
- of->fp_ene = open_enx(ftp2fn(efEDR,nfile,fnm), filemode);
+ if (EI_DYNAMICS(ir->eI) || EI_ENERGY_MINIMIZATION(ir->eI))
+ {
+ of->fp_ene = open_enx(ftp2fn(efEDR,nfile,fnm), filemode);
+ }
of->fn_cpt = opt2fn("-cpo",nfile,fnm);
if (ir->efep != efepNO && ir->nstdhdl > 0 &&
(ir->separate_dhdl_file == sepdhdlfileYES ) &&
- !EI_ENERGY_MINIMIZATION(ir->eI))
+ EI_DYNAMICS(ir->eI))
{
if (bAppendFiles)
{
for(i=0;i<cb->nsamples;i++)
{
samples_t *s=cb->s[i];
- sample_range_t *r=&(ca->r[i]);
+ sample_range_t *r=&(cb->r[i]);
if (r->use)
{
if (!s->hist)
atom_id **index;
rvec *x,hbox;
matrix box;
- real t,ccut,dist,ang;
+ real t,ccut,dist=0.0,ang=0.0;
double max_nhb,aver_nhb,aver_dist;
- int h,i,j,k=0,l,start,end,id,ja,ogrp,nsel;
+ int h=0,i,j,k=0,l,start,end,id,ja,ogrp,nsel;
int xi,yi,zi,ai;
int xj,yj,zj,aj,xjj,yjj,zjj;
int xk,yk,zk,ak,xkk,ykk,zkk;