Hi again,<div><br></div><div>Here's another question.</div><div>Let our function be:</div><div><br></div><div>y = myfunction(double x, double aux){</div><div> trace_on(0);</div><div> adouble xa, auxa;</div><div>
auxa <<= aux;</div><div> xa <<= x;</div><div> auxa * xa >>= y;</div><div> trace_off();</div><div>}</div><div><br></div><div>where "aux" is an extra parameter (and therefore we'll only want dy/dx).</div>
<div>As far as I know, dauxa/dx is assumed to be zero. Is there any way of giving dauxa/dx an initial value (different from zero)? Ideally, it would be something like:</div><div><br></div><div>auxa <<= aux;
</div><div>auxa.derivative() = whatever</div>
<div><br></div><div>Thanks a lot.</div><div>Alfonso </div><div><br><div class="gmail_quote">El 6 de marzo de 2012 03:02, Andrea Walther <span dir="ltr"><<a href="mailto:awalther@math.uni-paderborn.de" target="_blank">awalther@math.uni-paderborn.de</a>></span> escribió:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Alfonso,<div><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
The problem with the approach you suggest is that the evaluation of<br>
jacobian() and the other drivers would return, in addition to the "true"<br>
derivatives, the derivatives of the dependent variables with respect to<br>
those "fake" independent variables, right?<br>
</blockquote>
<br></div>
yes that is true. However, usually this would not harm the computing time<br>
too much and one can just ignore the additional derivatives.<br>
<br>
Currently, this is the easiest way to take additional parameters into account.<br>
<br>
Best regards,<div><div><br>
<br>
Andrea<br>
<br>
-- <br>
Prof. Dr. Andrea Walther<br>
Lehrstuhl fuer Mathematik und ihre Anwendungen<br>
Institut fuer Mathematik<br>
Universitaet Paderborn<br>
Warburger Str. 100<br>
33098 Paderborn<br>
<br>
Email: <a href="mailto:andrea.walther@uni-paderborn.de" target="_blank">andrea.walther@uni-paderborn.<u></u>de</a><br>
Phone: <a href="tel:%2B%2B49%205251%20602721" value="+495251602721" target="_blank">++49 5251 602721</a><br>
<a href="tel:%2B%2B49%205251%20602724" value="+495251602724" target="_blank">++49 5251 602724</a> (sekr.)<br>
Fax: <a href="tel:%2B%2B49%205251%20603728" value="+495251603728" target="_blank">++49 5251 603728</a><br>
<br>
**********<br>
<br>
<br>
<br>
</div></div></blockquote></div><br></div>