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. -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
