Writing floating point number to disk

[mathieu]

hi there,

I do not understand the syntax for ios_base::preci sion. Let say I
have a floating point:

const float f = 0.313244462;

How do I write is to a stream ?

const float f = 0.313244462;
std::cout << f << std::endl;

This truncate my number to an arbitrary number of digit after the
decimal point. all I am trying to do is write all meaningfull digit so
that when rereading this number from disk, I can find *exactly* the
same floating point number (I am using 'float' type and not 'double'
type).

Thanks
-Mathieu

[moongeegee]

Re: Writing floating point number to disk


const float f = 0.313244462;
std::cout << f << std::endl;
The output is 0.313244

const float f = 0.313244462;
cout.precision( 10);
std::cout << f << std::endl;
The out put is 0.313244462


On Sep 17, 8:27 am, mathieu <mathieu.malate ...@gmail.comwr ote:

hi there,
>
I do not understand the syntax for ios_base::preci sion. Let say I
have a floating point:
>
const float f = 0.313244462;
>
How do I write is to a stream ?
>
const float f = 0.313244462;
std::cout << f << std::endl;
>
This truncate my number to an arbitrary number of digit after the
decimal point. all I am trying to do is write all meaningfull digit so
that when rereading this number from disk, I can find *exactly* the
same floating point number (I am using 'float' type and not 'double'
type).
>
Thanks
-Mathieu

[mathieu]

Re: Writing floating point number to disk

[following initiated top post]

I do not want to specify an hardcoded '10' value in my program, this
should be deduced from the size of float and the number of significant
digit.

On Sep 17, 3:17 pm, moongeegee <moongee...@gma il.comwrote:

const float f = 0.313244462;
std::cout << f << std::endl;
The output is 0.313244
>
const float f = 0.313244462;
cout.precision( 10);
std::cout << f << std::endl;
The out put is 0.313244462
>
On Sep 17, 8:27 am, mathieu <mathieu.malate ...@gmail.comwr ote:
>

hi there,

>

I do not understand the syntax for ios_base::preci sion. Let say I
have a floating point:

>

const float f = 0.313244462;

>

How do I write is to a stream ?

>

const float f = 0.313244462;
std::cout << f << std::endl;

>

This truncate my number to an arbitrary number of digit after the
decimal point. all I am trying to do is write all meaningfull digit so
that when rereading this number from disk, I can find *exactly* the
same floating point number (I am using 'float' type and not 'double'
type).

>

Thanks
-Mathieu

[Juha Nieminen]

Re: Writing floating point number to disk

mathieu wrote:

I do not want to specify an hardcoded '10' value in my program

Then use '100' or whatever. AFAIK ostream only outputs at most that
many decimals. It won't output trailing zeros.

[James Kanze]

Re: Writing floating point number to disk

On Sep 17, 3:46 pm, Michael DOUBEZ <michael.dou... @free.frwrote:

mathieu a écrit :

[following initiated top post]

I do not want to specify an hardcoded '10' value in my
program, this should be deduced from the size of float and
the number of significant digit.

On 32 bit computer, float precision is +/- 5e-8 (i.e. 7 digits
precision) and double precision (64 bits) is 53 digits.

Really? On some 32 bit computers, maybe, but certainly not on
all. In C++0x, he can use std::numeric_li mits< float

>::max_digits 10 for the necessary precision, but this value

wasn't available in earlier versions of the standard (and thus
likely isn't available with his current compiler). There is an
std::numeric_li mits< float >::digits10, but this gives the
opposite value: the maximum number of digits you can read and
then output with no change. (For example, for an IEEE
float---the most common format on small and medium sized
machines, max_digits10 is 9, but digits10 is only 6.)

cout.setprecisi on(7);
should be enough for float but in doubt, you can set it to a
larger value if you want.

It's not enough on an Intel or AMD based PC, nor on a Sparc (the
machines I regularly use). It's probably not enough on a number
of other machines.

--
James Kanze (GABI Software) email:james.kan ze@gmail.com
Conseils en informatique orientée objet/
Beratung in objektorientier ter Datenverarbeitu ng
9 place Sémard, 78210 St.-Cyr-l'École, France, +33 (0)1 30 23 00 34

[James Kanze]

Re: Writing floating point number to disk

On Sep 17, 3:38 pm, Juha Nieminen <nos...@thanks. invalidwrote:

mathieu wrote:

I do not want to specify an hardcoded '10' value in my program

Then use '100' or whatever. AFAIK ostream only outputs at most
that many decimals. It won't output trailing zeros.

But he doesn't care about trailing zeros. His problem, as
originally stated, was to ensure round trip accuracy starting
from the internal representation: from an internal float to text
and back, ending up with exactly the same value. (This is, of
course, only possible if both the machine writing and the
machine reading have the same floating point format.) Ten is
sufficient for IEEE float; in fact, so is 9. For other formats,
you might need more or less.

And I would certainly not use 100. All you'll get is a lot of
extra useless and insignificant digits for most values. And
will output 100 digits if your using either fixed or scientific
format.

--
James Kanze (GABI Software) email:james.kan ze@gmail.com
Conseils en informatique orientée objet/
Beratung in objektorientier ter Datenverarbeitu ng
9 place Sémard, 78210 St.-Cyr-l'École, France, +33 (0)1 30 23 00 34

[Kai-Uwe Bux]

Re: Writing floating point number to disk

mathieu wrote:

hi there,
>
I do not understand the syntax for ios_base::preci sion. Let say I
have a floating point:
>
const float f = 0.313244462;
>
How do I write is to a stream ?
>
const float f = 0.313244462;
std::cout << f << std::endl;
>
This truncate my number to an arbitrary number of digit after the
decimal point. all I am trying to do is write all meaningfull digit so
that when rereading this number from disk, I can find *exactly* the
same floating point number (I am using 'float' type and not 'double'
type).

Maybe something like the following will help:

template < typename Float >
std::string to_string ( Float data ) {
std::stringstre am in;
unsigned long const digits =
static_cast< unsigned long >(
- std::log( std::numeric_li mits<Float>::ep silon() )
/ std::log( 10.0 ) );
if ( in << std::dec << std::setprecisi on(2+digits) << data ) {
return ( in.str() );
} else {
throw some_thing;
}
}


Best

Kai-Uwe Bux

[mathieu]

Re: Writing floating point number to disk

On Sep 17, 5:37 pm, Kai-Uwe Bux <jkherci...@gmx .netwrote:

mathieu wrote:

hi there,

>

I do not understand the syntax for ios_base::preci sion. Let say I
have a floating point:

>

const float f = 0.313244462;

>

How do I write is to a stream ?

>

const float f = 0.313244462;
std::cout << f << std::endl;

>

This truncate my number to an arbitrary number of digit after the
decimal point. all I am trying to do is write all meaningfull digit so
that when rereading this number from disk, I can find *exactly* the
same floating point number (I am using 'float' type and not 'double'
type).

>
Maybe something like the following will help:
>
template < typename Float >
std::string to_string ( Float data ) {
std::stringstre am in;
unsigned long const digits =
static_cast< unsigned long >(
- std::log( std::numeric_li mits<Float>::ep silon() )
/ std::log( 10.0 ) );
if ( in << std::dec << std::setprecisi on(2+digits) << data ) {
return ( in.str() );
} else {
throw some_thing;
}
>
}
>
Best
>
Kai-Uwe Bux

Where is the "2+" coming from in setprecision ?

I found one of your earlier post:

http://www.archivum.info/comp.lang.c++/2005-10/msg03220.html

where this offset does not appear.

Thanks

[Kai-Uwe Bux]

Re: Writing floating point number to disk

mathieu wrote:

On Sep 17, 5:37 pm, Kai-Uwe Bux <jkherci...@gmx .netwrote:

>mathieu wrote:

hi there,

>>

I do not understand the syntax for ios_base::preci sion. Let say I
have a floating point:

>>

const float f = 0.313244462;

>>

How do I write is to a stream ?

>>

const float f = 0.313244462;
std::cout << f << std::endl;

>>

This truncate my number to an arbitrary number of digit after the
decimal point. all I am trying to do is write all meaningfull digit so
that when rereading this number from disk, I can find *exactly* the
same floating point number (I am using 'float' type and not 'double'
type).

>>
>Maybe something like the following will help:
>>
>template < typename Float >
>std::string to_string ( Float data ) {
> std::stringstre am in;
> unsigned long const digits =
> static_cast< unsigned long >(
> - std::log( std::numeric_li mits<Float>::ep silon() )
> / std::log( 10.0 ) );
> if ( in << std::dec << std::setprecisi on(2+digits) << data ) {
> return ( in.str() );
> } else {
> throw some_thing;
> }
>>
>}
>>
>Best
>>
>Kai-Uwe Bux

>
Where is the "2+" coming from in setprecision ?
>
I found one of your earlier post:

http://www.archivum.info/comp.lang.c++/2005-10/msg03220.html

>
where this offset does not appear.

It's just a safeguard against typical off-by-one errors (so 1 would probably
suffice). In short, I refuse to think when CPU-cycles are much cheaper that
brain-cycles.


Best

Kai-Uwe Bux

[Juha Nieminen]

Re: Writing floating point number to disk

Michael DOUBEZ wrote:

On 32 bit computer, float precision is +/- 5e-8 (i.e. 7 digits
precision) and double precision (64 bits) is 53 digits.

You mean that's not the case with 64-bit computers?

(Or, in other words: Why mention "32 bit computer" explicitly, given
that the exact same floating point precision is used in most 64-bit
computers as well. In fact, it was also used in 16-bit Intel-based
computers too.)

[James Kanze]

Re: Writing floating point number to disk

On Sep 18, 12:12 pm, Michael DOUBEZ <michael.dou... @free.frwrote:

James Kanze a écrit :

On Sep 17, 3:46 pm, Michael DOUBEZ <michael.dou... @free.frwrote:

mathieu a écrit :

>[following initiated top post]

>I do not want to specify an hardcoded '10' value in my
>program, this should be deduced from the size of float and
>the number of significant digit.

On 32 bit computer, float precision is +/- 5e-8 (i.e. 7 digits
precision) and double precision (64 bits) is 53 digits.

Really?

That's the values I take when computing the precision of a calculus.

On some 32 bit computers, maybe, but certainly not on all.

Provided they use the IEEE representation, they have 24 bits
to encode the fractional part and thus 7 digits precision. I
thought it was the most widespread norm in use.

The most widespread, but far from universal. Most mainframes
use something different.

In C++0x, he can use std::numeric_li mits<
float >::max_digits 10 for the necessary precision, but this
value wasn't available in earlier versions of the standard
(and thus likely isn't available with his current compiler).
There is an std::numeric_li mits< float >::digits10, but this
gives the opposite value: the maximum number of digits you
can read and then output with no change. (For example, for
an IEEE float---the most common format on small and medium
sized machines, max_digits10 is 9, but digits10 is only 6.)

cout.setprecisi on(7);
should be enough for float but in doubt, you can set it to a
larger value if you want.

It's not enough on an Intel or AMD based PC, nor on a Sparc (the
machines I regularly use). It's probably not enough on a number
of other machines.

I understand 8 digits since it could reduce an increased error
on the last digit ( ...88 written ...9 read ...91 ) but how
they came up with 9 or more ? (not that it is critical).

To tell the truth, I don't know. I just copied the value out of
the draft standard:-). (I'd have guessed 8 myself. But I know
that I'm no specialist in this domain. I only know enough to
know how much I don't know.)

--
James Kanze (GABI Software) email:james.kan ze@gmail.com
Conseils en informatique orientée objet/
Beratung in objektorientier ter Datenverarbeitu ng
9 place Sémard, 78210 St.-Cyr-l'École, France, +33 (0)1 30 23 00 34

[James Kanze]

Re: Writing floating point number to disk

On Sep 18, 3:40 pm, Juha Nieminen <nos...@thanks. invalidwrote:

Michael DOUBEZ wrote:

On 32 bit computer, float precision is +/- 5e-8 (i.e. 7 digits
precision) and double precision (64 bits) is 53 digits.

You mean that's not the case with 64-bit computers?

(Or, in other words: Why mention "32 bit computer" explicitly,
given that the exact same floating point precision is used in
most 64-bit computers as well. In fact, it was also used in
16-bit Intel-based computers too.)

Exactly. The problem isn't the number of bits, but the floating
point format (although I can't quite see an IEEE format of a
machine which isn't 8/16/32/64 bits). Thus, I'm aware of at
least four floating point formats in common use today, two of
which are used on 32 bit machines.

--
James Kanze (GABI Software) email:james.kan ze@gmail.com
Conseils en informatique orientée objet/
Beratung in objektorientier ter Datenverarbeitu ng
9 place Sémard, 78210 St.-Cyr-l'École, France, +33 (0)1 30 23 00 34

[Juha Nieminen]

Re: Writing floating point number to disk

James Kanze wrote:

Thus, I'm aware of at
least four floating point formats in common use today, two of
which are used on 32 bit machines.

Three, actually. Intel FPUs support 80-bit floating point numbers, and
usually you can create one with "long double".

[Rune Allnor]

Re: Writing floating point number to disk

On 18 Sep, 15:40, Juha Nieminen <nos...@thanks. invalidwrote:

Michael DOUBEZ wrote:

On 32 bit computer, float precision is +/- 5e-8 (i.e. 7 digits
precision) and double precision (64 bits) is 53 digits.

>
  You mean that's not the case with 64-bit computers?
>
  (Or, in other words: Why mention "32 bit computer" explicitly, given
that the exact same floating point precision is used in most 64-bit
computers as well. In fact, it was also used in 16-bit Intel-based
computers too.)

floats need not be 32 bit on 64-bit computers, just the same
way ints don't need to be 16 bit on a 32-bit computer.

Rune

[James Kanze]

Re: Writing floating point number to disk

On Sep 19, 12:28 pm, Juha Nieminen <nos...@thanks. invalidwrote:

James Kanze wrote:

Thus, I'm aware of at least four floating point formats in
common use today, two of which are used on 32 bit machines.

Three, actually. Intel FPUs support 80-bit floating point
numbers, and usually you can create one with "long double".

That's not what I meant. I was thinking of four different
formats for float, or for double. Off hand, for float, there's
IEEE (PC's and most Unix boxes), IBM floating point (also 32
bits, but base 16, for starters), and the floating point formats
for the Unisys mainframes (36 bits on one architecture, 48 on
the other, both base 8). There are certainly others as well,
but I don't know if you'll find them on any modern machines.

--
James Kanze (GABI Software) email:james.kan ze@gmail.com
Conseils en informatique orientée objet/
Beratung in objektorientier ter Datenverarbeitu ng
9 place Sémard, 78210 St.-Cyr-l'École, France, +33 (0)1 30 23 00 34