Loss of precision of parameters in a model file

Double precision floating-point number parameters in svm-train, such as gamma, are written to a model file in a format of %g by svm_save_model in svm.cpp. This causes loss of precision of parameters, hence a slight difference between the kernel function in svm-train and that of svm-predict. The effect on precision of prediction may be small, especially when Qfloat is defined as float, but I think the parameters should be written in a format like %.16g.

[cjlin1]

FYI, in the libsvm 3.23 released 2 days ago this has been corrected.

Tavian Barnes writes:

%.16 should be enough due to the normalized representation
(i.e., 17 digits in total)

%.16g is only 16 digits total. Here’s an example that shows it’s not enough:

#include <stdio.h>

int main() { double a = 18014398509481982.0; double b = 18014398509481980.0; printf(“%.16g\n%.16g\n%s\n”, a, b, a == b ? “true” : “false”); return 0; }

Output:

1.801439850948198e+16 1.801439850948198e+16 false

From the famous “What Every Computer Scientist Should Know About Floating-Point Arithmetic” paper, Theorem 15:

... The same argument applied to double precision shows
that 17 decimal digits are required to recover a double
precision number.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub, or mute the thread.*

[cjlin1]

I think you are right… We will change the code later Tavian Barnes writes:

%.16 should be enough due to the normalized representation
(i.e., 17 digits in total)

%.16g is only 16 digits total. Here’s an example that shows it’s not enough:

#include <stdio.h>

int main() { double a = 18014398509481982.0; double b = 18014398509481980.0; printf(“%.16g\n%.16g\n%s\n”, a, b, a == b ? “true” : “false”); return 0; }

Output:

1.801439850948198e+16 1.801439850948198e+16 false

From the famous “What Every Computer Scientist Should Know About Floating-Point Arithmetic” paper, Theorem 15:

... The same argument applied to double precision shows that
17 decimal digits are required to recover a double precision
number.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub, or mute the thread.*