<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body dir="auto">
<div>Make sure the arguments of the function are also df1b2variable<br>
<br>
Sent from my iPad</div>
<div><br>
On Nov 24, 2014, at 2:51 AM, Stirrup, Oliver <<a href="mailto:oliver.stirrup.13@ucl.ac.uk">oliver.stirrup.13@ucl.ac.uk</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div><style>
<!--
@font-face
{font-family:Calibri}
@font-face
{font-family:Verdana}
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif"}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline}
span.EmailStyle17
{font-family:"Calibri","sans-serif";
color:windowtext}
.MsoChpDefault
{font-family:"Calibri","sans-serif"}
@page WordSection1
{margin:72.0pt 72.0pt 72.0pt 72.0pt}
div.WordSection1
{}
-->
</style>
<div class="WordSection1">
<p class="MsoNormal">Hello all</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I have found that use of the ‘pow’ function causes executables created by ADMB to crash when the first argument is zero, i.e. when a calculation of 0^x is attempted.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">In my template files, data relating to a time variable are read into a dvar_vector and (within a loop) each element of this vector is then raised to the power of a model parameter. The time variable is positive in most cases, but occasionally
takes the value of zero, leading to the problem described.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Using the standard mode for ADMB, I have overcome this problem by defining the following function to replace ‘pow’:</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><<< >>></p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">GLOBALS_SECTION</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"> #include <admodel.h></p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"> dvariable pow2(const dvariable & time, const dvariable & T)</p>
<p class="MsoNormal"> {</p>
<p class="MsoNormal"> if(time==0)</p>
<p class="MsoNormal"> return 0;</p>
<p class="MsoNormal"> else return pow(time,T);</p>
<p class="MsoNormal"> }</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><<< >>></p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I would now like to implement a similar solution for the random effects mode of ADMB, but have so far been unsuccessful in my attempts, with a large number of error messages on compilation. I think that my problem lies in the fact that
I do not really understand how to define and manipulate df1b2variable objects.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">If anyone can offer a solution or any advice regarding this problem, I would be very grateful.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Best wishes</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Oliver</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><b><span style="font-size:8.0pt; font-family:"Verdana","sans-serif"">Oliver Stirrup</span></b></p>
<p class="MsoNormal"><span style="font-size:8.0pt; font-family:"Verdana","sans-serif"">PhD Student</span></p>
<p class="MsoNormal"><span style="font-size:8.0pt; font-family:"Verdana","sans-serif""> </span></p>
<p class="MsoNormal"><span style="font-size:8.0pt; font-family:"Verdana","sans-serif"">MRC Clinical Trials Unit at UCL</span></p>
<p class="MsoNormal"><span style="font-size:8.0pt; font-family:"Verdana","sans-serif"">Institute of Clinical Trials & Methodology</span></p>
<p class="MsoNormal"><span style="font-size:8.0pt; font-family:"Verdana","sans-serif"">London UK</span></p>
<p class="MsoNormal"><span style="font-size:8.0pt; font-family:"Verdana","sans-serif"">e-mail:
<a href="oliver.stirrup.13@ucl.ac.uk"><span style="color:blue">oliver.stirrup.13@ucl.ac.uk</span></a></span></p>
<p class="MsoNormal"> </p>
</div>
</div>
</blockquote>
<blockquote type="cite">
<div><span>_______________________________________________</span><br>
<span>Users mailing list</span><br>
<span><a href="mailto:Users@admb-project.org">Users@admb-project.org</a></span><br>
<span><a href="http://lists.admb-project.org/mailman/listinfo/users">http://lists.admb-project.org/mailman/listinfo/users</a></span><br>
</div>
</blockquote>
<br>
<hr>
<font face="Arial" color="Gray" size="1"><br>
This internet e-mail message, and any files transmitted with it, contains confidential, privileged information that is intended only for the addressee. If you have received this e-mail message in error, please call us at (206) 634-1838 collect if necessary)
and ask to speak to the message sender. Nothing in this e-mail or the act of transmitting it, is to be construed as a waiver of any rights or privileges enjoyed by the sender or the International Pacific Halibut Commission pursuant to the International Organizations
Immunities Act, 22 U.S.C. Sec. 288 et seq.<br>
</font>
</body>
</html>