On 05/14/2015 11:42 AM, Hazen Babcock wrote: > On 05/13/2015 05:50 PM, Alan W. Irwin wrote: >> >> Hi Hazen: >> >> I have reviewed plfill.c, and the returned intersection from >> notcrossed is only used when status is zero. So I have (commit ID = >> b916d4b) considerably simplified noncrossed to do the PLINT >> transformation of the intersection and return those values _only_ when >> status = 0. >> >> This change solved all fpe problem for examples 25 and 30, but there >> are remaining fpe issues for fortran examples 21 and 33. The logic in >> noncrossed is pretty bullet-proof now so my prediction is those >> remaining fpe issues are completely independent of notcrossed or else >> there is a non-initialized variable somewhere in fortran examples 21 >> and 33. But I look forward to your definitive conclusion concerning >> those remaining fpe issues. > > The issue in example 21 seems to be in the qhull library. I will look > into building a qhull with debug flag to try and track this one down. > > The remaining issue in example 33 is that we are multiplying -2.0e-200 > times 9.9e-201 when calculating iedge[i] in pldrawcn() at line 695 of > src/plcont.c. The result of this calculation, ~1.0e-400, is smaller than > what can be represented using floating point numbers. I changed the > login in pldrawcn() to do the sign determination without doing this > multiplication. I've just pushed this fix.
Hm, got the signs backwards, hopefully it is correct now. -Hazen ------------------------------------------------------------------------------ One dashboard for servers and applications across Physical-Virtual-Cloud Widest out-of-the-box monitoring support with 50+ applications Performance metrics, stats and reports that give you Actionable Insights Deep dive visibility with transaction tracing using APM Insight. http://ad.doubleclick.net/ddm/clk/290420510;117567292;y _______________________________________________ Plplot-devel mailing list Plplot-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/plplot-devel
