On Tue, May 31, 2016 at 10:18 AM, Lorenzo Zanon <za...@aices.rwth-aachen.de> wrote:
> > > On 5/31/16 5:41 PM, John Peterson wrote: > > > > On Tue, May 31, 2016 at 3:23 AM, Lorenzo Zanon <za...@aices.rwth-aachen.de > > wrote: > >> Hello, >> >> I solved the problem. During the "online" run, I was writing down the >> coefficients of the reduced solution u_N in data files, which were then >> read in a separate run to compute the error bound. I converted this >> writing/reading from data files to binary; the precision is now preserved. >> > > The ASCII writing in legacy_write_offline_data_to_files() is done through > the Xdr class, but that should be writing 16 digits: > > xdr_cxx.C, line 746: > > *out << std::scientific >> << std::setprecision(16); > > > It looks like these lines were added in aff6b5ea, back on 2015-10-26, > though, so if you have al libmesh from before that, the precision might > have been lower. > > My version includes these lines, but still writing down the data file > would deliver 8-digit precision. The .dat file looks like: > > 0 > 0.202137 > -6.50782e-14 > 1.64219e-14 > 5.42741e-15 > Hmm... do you know if those are solution values, node (x,y,z) locations, or something else? I don't see a place in RBEIMEvaluation or its parent, RBEvaluation, where floating point values are written without going through an Xdr object. What compiler are you using currently? It's possible (but highly unlikely) that your compiler does not properly support std::setprecision()... -- John ------------------------------------------------------------------------------ What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic patterns at an interface-level. Reveals which users, apps, and protocols are consuming the most bandwidth. Provides multi-vendor support for NetFlow, J-Flow, sFlow and other flows. Make informed decisions using capacity planning reports. https://ad.doubleclick.net/ddm/clk/305295220;132659582;e _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
