2007/8/14, William Stein <[EMAIL PROTECTED]>: > Correct me if I'm wrong, but I don't think that makes sense. > If I do this: > > sage: a = RQDF(5) > sage: number_of_partitions(1000) > 24061467864032622473692149727991 > sage: del a > > then during the number_of_partitions call the CPU is set to > the wrong mode.
You are right, I did not think about this. 2007/8/13, Jonathan Bober <[EMAIL PROTECTED]>: > So -ffloat-store doesn't work. So on cpus without sse2, the quaddouble > wrapper needs to use the fpu fix, but probably should be rewritten so > that it doesn't affect other things. > One possibility is that one wrap _add_c_impl, _sub_c_impl, > etc., by setting the CPU mode, doing the op, then unsetting > the cpu mode. This is probably the robust solution (??), > but might have potential efficiency issues (which should be > investigated with careful benchmarking). It's looking like the *only* solution. Looks fail-safe, and of course, very ugly :). On a related note, it looks like GSL has a "-gsl" compilation flag that "Use[s] GSL to control floating point rounding and precision". didier --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~----------~----~----~----~------~----~------~--~---
