<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7654.12">
<TITLE>very small number math</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->
<P DIR=LTR><SPAN LANG="en-ca"><FONT FACE="Calibri">Hello,</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT FACE="Calibri">I think I’ve hit a precision limit</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri">, looking for a solution along the same sort of lines as posfun</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri"> or mfexp</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri">. My program stops (no error message) in a loop</FONT></SPAN><SPAN LANG="en-ca"> <FONT FACE="Calibri">when a particular equation attempts to divide 19616.8 by 7.57411e-154.</FONT></SPAN><SPAN LANG="en-ca"> <FONT FACE="Calibri">The preceding iteration worked with 21666/5.73085e-132. These are s</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri">illy</FONT></SPAN><SPAN LANG="en-ca"> <FONT FACE="Calibri">starting</FONT></SPAN><SPAN LANG="en-ca"> <FONT FACE="Calibri">estimates of a value which will get</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri"> reasonable with a little fitting, I just</FONT></SPAN><SPAN LANG="en-ca"> <FONT FACE="Calibri">want</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri"> to accept</FONT></SPAN><SPAN LANG="en-ca"> <FONT FACE="Calibri">a maximal result instead of bombing.</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri"></FONT></SPAN><SPAN LANG="en-ca"> <FONT FACE="Calibri">In this case I can just use a min function on the denominator</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri">, but wonder if there isn’t a conventional approach for this situation</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri">, say defining a minimum number threshold</FONT></SPAN><SPAN LANG="en-ca"><FONT FACE="Calibri">.</FONT></SPAN><SPAN LANG="en-ca"></SPAN></P>
<UL DIR=LTR><UL DIR=LTR>
<P DIR=LTR><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Mark Fowler</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Population Ecology Division</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Bedford Inst of Oceanography</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Dept Fisheries & Oceans</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Dartmouth NS Canada</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">B2Y 4A2</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Tel. (902) 426-3529</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Fax (902) 426-9710</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Email Mark.Fowler@dfo-mpo.gc.ca</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Home Tel. (902) 461-0708</FONT></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"><FONT SIZE=2 FACE="Arial">Home Email mark.fowler@ns.sympatico.ca</FONT></SPAN></P>
</UL></UL>
<P DIR=LTR><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN><SPAN LANG="en-ca"></SPAN></P>
<P DIR=LTR><SPAN LANG="en-ca"></SPAN></P>
</BODY>
</HTML>