<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi!<br>
<br>
I am very interested by the ampl toolbox. However, there is a missing
feature I need for my work with truncated Newton algorithms: computing
the product of the hessian with a vector without evaluating the whole
hessian. This is available according to "hooking your solver..."<br>
<blockquote>
  <dl>
    <dd>
      <pre><tt>void hvcomp(real *HV, real *P, int nobj, real *OW, real *Y)
</tt></pre>
    </dd>
  </dl>
stores in <tt>HV</tt>
(a full vector of length <tt>n_var</tt>)
the Hessian of the Lagrangian times vector <tt>P</tt>.
In other words, <tt>hvcomp</tt>
computes <br>
  <img src="cid:part1.01070403.03090001@USherbrooke.ca"><br>
where <i>W</i>
is the Lagrangian Hessian, <a name="Wformula"></a><br>
  <img src="cid:part2.02070301.02090202@USherbrooke.ca"><br>
  <br>
</blockquote>
This is not in the toolbox. Would it be possible to add this function
to the toolbox?<br>
<br>
Thank you,<br>
<br>
JPD<br>
<br>
</body>
</html>