Index: src/TwoPunctures.c =================================================================== --- src/TwoPunctures.c (revision 102) +++ src/TwoPunctures.c (working copy) @@ -281,11 +281,11 @@ } char valbuf_p[100], valbuf_m[100]; - sprintf (valbuf_p,"%f", (float) *mp); + sprintf (valbuf_p,"%.17g", (double)*mp); CCTK_ParameterSet ("par_m_plus", "twopunctures", valbuf_p); - sprintf (valbuf_m,"%f", (float) *mm); + sprintf (valbuf_m,"%.17g", (double)*mm); CCTK_ParameterSet ("par_m_minus", "twopunctures", valbuf_m); - new_mass = (double)*mp + (double)*mm; + new_mass = *mp + *mm; } } admMass = (*mp + *mm @@ -297,8 +297,8 @@ F_of_v (cctkGH, nvar, n1, n2, n3, v, F, u); CCTK_VInfo (CCTK_THORNSTRING, - "The two puncture masses are %g and %g", - (double) par_m_minus, (double) par_m_plus); + "The two puncture masses are %.17g and %.17g", + (double) *mm, (double) *mp); /* print out ADM mass, eq.: \Delta M_ADM=2*r*u=4*b*V for A=1,B=0,phi=0 */ admMass = (*mp + *mm