Merge "Merge branch release-4-6 into master"
[alexxy/gromacs.git] / src / gromacs / mdlib / sim_util.c
index 46ada2a183b202e27ffa55720f2459c5e0550c53..98f8f4e0e66ada65c056956609138134ece76900 100644 (file)
 #include "nbnxn_cuda_data_mgmt.h"
 #include "nbnxn_cuda/nbnxn_cuda.h"
 
-#if 0
-typedef struct gmx_timeprint {
-
-} t_gmx_timeprint;
-#endif
-
-/* Portable version of ctime_r implemented in src/gmxlib/string2.c, but we do not want it declared in public installed headers */
-char *
-gmx_ctime_r(const time_t *clock, char *buf, int n);
-
-
 double
 gmx_gettime()
 {
@@ -426,8 +415,7 @@ static void posres_wrapper(FILE *fplog,
                   fr->rc_scaling, fr->ePBC, fr->posres_com, fr->posres_comB);
     if (bSepDVDL)
     {
-        fprintf(fplog, sepdvdlformat,
-                interaction_function[F_POSRES].longname, v, dvdl);
+        gmx_print_sepdvdl(fplog, interaction_function[F_POSRES].longname, v, dvdl);
     }
     enerd->term[F_POSRES] += v;
     /* If just the force constant changes, the FEP term is linear,
@@ -480,7 +468,7 @@ static void pull_potential_wrapper(FILE *fplog,
                        cr, t, lambda[efptRESTRAINT], x, f, vir_force, &dvdl);
     if (bSepDVDL)
     {
-        fprintf(fplog, sepdvdlformat, "Com pull", enerd->term[F_COM_PULL], dvdl);
+        gmx_print_sepdvdl(fplog, "Com pull", enerd->term[F_COM_PULL], dvdl);
     }
     enerd->dvdl_lin[efptRESTRAINT] += dvdl;
 }
@@ -507,7 +495,7 @@ static void pme_receive_force_ener(FILE           *fplog,
                       &cycles_seppme);
     if (bSepDVDL)
     {
-        fprintf(fplog, sepdvdlformat, "PME mesh", e, dvdl);
+        gmx_print_sepdvdl(fplog, "PME mesh", e, dvdl);
     }
     enerd->term[F_COUL_RECIP] += e;
     enerd->dvdl_lin[efptCOUL] += dvdl;
@@ -1698,12 +1686,6 @@ void do_force_cutsGROUP(FILE *fplog, t_commrec *cr,
            groups, &(inputrec->opts), top, mdatoms,
            cr, nrnb, lambda, dvdlambda, &enerd->grpp, bFillGrid,
            bDoLongRangeNS);
-        if (bSepDVDL)
-        {
-            fprintf(fplog, sepdvdlformat, "LR non-bonded", 0.0, dvdlambda);
-        }
-        enerd->dvdl_lin[efptVDW]  += dvdlambda[efptVDW];
-        enerd->dvdl_lin[efptCOUL] += dvdlambda[efptCOUL];
 
         wallcycle_stop(wcycle, ewcNS);
     }
@@ -2388,8 +2370,7 @@ void calc_dispcorr(FILE *fplog, t_inputrec *ir, t_forcerec *fr,
 
         if (fr->bSepDVDL && do_per_step(step, ir->nstlog))
         {
-            fprintf(fplog, sepdvdlformat, "Dispersion correction",
-                    *enercorr, dvdlambda);
+            gmx_print_sepdvdl(fplog, "Dispersion correction", *enercorr, dvdlambda);
         }
         if (fr->efep != efepNO)
         {