This is the mail archive of the
ecos-bugs@sourceware.org
mailing list for the eCos project.
[Bug 1001539] Single precision floating point math library
- From: bugzilla-daemon at bugs dot ecos dot sourceware dot org
- To: unassigned at bugs dot ecos dot sourceware dot org
- Date: Thu, 22 Mar 2012 13:42:23 +0000
- Subject: [Bug 1001539] Single precision floating point math library
- Auto-submitted: auto-generated
- References: <bug-1001539-777@http.bugs.ecos.sourceware.org/>
Please do not reply to this email. Use the web interface provided at:
http://bugs.ecos.sourceware.org/show_bug.cgi?id=1001539
--- Comment #19 from Sergei Gavrikov <sergei.gavrikov@gmail.com> 2012-03-22 13:42:17 GMT ---
(In reply to comment #18)
> (In reply to comment #12)
> > (In comment #9 I wrote)
> >
> > > All f-tests fail on eCos synthetic target
> >
> > My fault (that was wrong grep GDB output)! I'm sorry. Only below fails:
> >
> > FAIL:<sqrtf() failed tests> Line: 92, File: tests/vectors/sqrtf.c
> > FAIL:<log10f() failed tests> Line: 92, File: tests/vectors/log10f.c
> > FAIL:<frexpf() failed tests> Line: 92, File: tests/vectors/frexpf.c
> > FAIL:<acosf() failed tests> Line: 92, File: tests/vectors/acosf.c
> > FAIL:<asinf() failed tests> Line: 92, File: tests/vectors/asinf.c
> > FAIL:<logf() failed tests> Line: 92, File: tests/vectors/logf.c
> >
> > So, I got the same results for eCos sythetic target as for ARM7 target.
> >
> > Sergei
>
> Tests expect NaN for out-of-domain arguments, but functions return something
> else. For instance asinf(10) returns 0. Here is a little example:
IMO, it can depend on libm compatibility, look at libm compat. wrappers
(w_*.c), e.g. at w_asin{,f}.c. Though, I think in our case it should be
a result of call __ieee754_asin{,f}(). And if libm compat mode is not
IEEE? I'm looking at misc/standard.c (__kernel_standard)...
> int main(void)
> {
> volatile float fla, flb, flc;
> volatile double dla, dlb, dlc;
>
> printf("Enter floats: ");
> scanf("%f %f", &fla, &flb);
> printf("Enter doubles: ");
> scanf ("%lf %lf", &dla, &dlb);
>
> flc = fla + flb;
> dlc = dla + dlb;
> printf("float: %f + %f = %f\n", fla, flb, flc);
> printf("double: %f + %f = %f\n", dla, dlb, dlc);
>
> flc = fla / flb;
> dlc = dla / dlb;
> printf("float: %f / %f = %f\n", fla, flb, flc);
> printf("double: %f / %f = %f\n", dla, dlb, dlc);
>
> flc = asinf(fla);
> dlc = asin(dla);
> printf("float: asinf(%f) = %f\n", fla, flc);
> printf("double: asin(%f) = %f\n", dla, dlc);
>
> return 0;
> }
>
> And some printouts:
>
> Enter floats: 0 0
> Enter doubles: 0 0
> float: 0.000000 + 0.000000 = 0.000000
> double: 0.000000 + 0.000000 = 0.000000
> float: 0.000000 / 0.000000 = nan
> double: 0.000000 / 0.000000 = nan
> float: asinf(0.000000) = 0.000000
> double: asin(0.000000) = 0.000000
>
> Enter floats: 10 10
> Enter doubles: 10 10
> float: 10.000000 + 10.000000 = 20.000000
> double: 10.000000 + 10.000000 = 20.000000
> float: 10.000000 / 10.000000 = 1.000000
> double: 10.000000 / 10.000000 = 1.000000
> float: asinf(10.000000) = 0.000000
> double: asin(10.000000) = 0.000000
>
> As we can the same behavior is for double. I haven't check but I could
> imagine that double's tests either don't contain these cases or expect
> 0.
Thank you for the example.
> I wander if this is some compatibility issue?! FYI Linux man pages say
> they should return NaN.
Further I would refer to http://sourceware.org/newlib/libm.html and thus
to http://sourceware.org/newlib/libm.html#asin However, they also talk
about NaN, errno (EDOM), etc.
> Then next question: what to fix, functions or tests?
Do you trust your checkErrorAcceptableFloat() implementation? As I could
see asin*(10.0) should return N-a-N if libm compat is IEEE. IMHO, we have
to compare return values for d-implementations with f-ones and then test
(compare) behaviors of its checkers. I will try to investigate tonight.
Sergei
--
Configure bugmail: http://bugs.ecos.sourceware.org/userprefs.cgi?tab=email
------- You are receiving this mail because: -------
You are the assignee for the bug.