Oh, sorry, I didn't express myself correctly. <br>I know how to install ColPack and having it recognised by ADOL-C (e.g. configuring ColPack to use the "/usr" prefix).<br>However there are some troubles when using ColPack for optimisation purposes (IPOPT), so I was wondering if enabling the sparse drivers without ColPack I could still benefit of some "extra" features, in my case e.g. if enabling sparsity without ColPack I would have been able to calculate and call the Jacobian only at a user-provided row/column points.<br>
In the povided code snipped the jacobian is called only over these points, but I don't know if it is stll calculated for all the points.<br><br>Thank you,<br>Antonello<br><br><br><div class="gmail_quote">2012/8/16 Kshitij Kulshreshtha <span dir="ltr"><<a href="mailto:kshitij@math.upb.de" target="_blank">kshitij@math.upb.de</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello, this means that you enabled sparse drivers in the configure<br>
script but it was not able to find the binary of libColPack or the<br>
ColPackHeaders.h file. Please make sure you used --with-colpack=DIR in<br>
the configure command and make sure ColPack is installed correctly in<br>
DIR. On 64-bit systems you may have to tell ColPack to install the<br>
binary in ${prefix}/lib64 by passing --libdir=${prefix}/lib64 to the<br>
configure command of ColPack.<br>
<div><div class="h5"><br>
As on 2012-08-16 10:34, Antonello Lobianco did write:<br>
> Hello, what does this exactly means?<br>
><br>
> Build sparse drivers: yes<br>
> Build with ColPack: no<br>
> Compressed sparse structures will not be available<br>
> Only sparsity patterns can be computed<br>
><br>
> I'm using ADOL-C as AD tool for the IPOPT solver and having problems<br>
> with ColPack when I need to run multiple optimisation problems (in<br>
> series) I implemented an algorithm for finding sparsity structure for my<br>
> specific model, at least for the Jacobian.<br>
><br>
> Then I replaced the eval_jac_g() function in the LuksanVlcek1 or<br>
> MittelmannDistCntrlNeumA examples<br>
> (ADOL-C/examples/additional_examples/ipopt) with:<br>
><br>
> eval_jac_g(){<br>
> if (values == NULL) {<br>
> // return the structure of the jacobian,<br>
> for(Index idx=0;idx<nzjelements.size();idx++){<br>
> iRow[idx] = nzjelements[idx][0];<br>
> jCol[idx] = nzjelements[idx][1];<br>
> }<br>
> }<br>
> else {<br>
> // return the values of the jacobian of the constraints<br>
> jacobian(tag_g,m,n,x,Jac);<br>
> for(Index idx=0;idx<nzjelements.size();idx++){<br>
> values[idx] = Jac[nzjelements[idx][0]][nzjelements[idx][1]];<br>
> }<br>
> }<br>
> return true;<br>
> }<br>
><br>
> where vector < vector <Index> > nzelements is the vector of non zero<br>
> elements in my model, with nzelements[i][0] the row position and<br>
> nzelements[i][1] the column one.<br>
><br>
> Is this the most efficient implementation or there could be a better way<br>
> to compute sparsity patterns?<br>
><br>
> Thank you,<br>
> Antonello<br>
><br>
><br>
><br>
><br>
> --<br>
> Antonello Lobianco<br>
> INRA, Laboratoire d'Economie Forestière<br>
> 14 Rue Girardet - 54000 Nancy, France<br>
> Tel: +33.652392310<br>
> Email: <a href="mailto:antonello.lobianco@nancy-engref.inra.fr">antonello.lobianco@nancy-engref.inra.fr</a><br>
</div></div>> <mailto:<a href="mailto:antonello.lobianco@nancy-engref.inra.fr">antonello.lobianco@nancy-engref.inra.fr</a>><br>
> <a href="http://antonello.lobianco.org" target="_blank">http://antonello.lobianco.org</a> <<a href="http://antonello.lobianco.org/" target="_blank">http://antonello.lobianco.org/</a>><br>
><br>
><br>
> _______________________________________________<br>
> ADOL-C mailing list<br>
> <a href="mailto:ADOL-C@list.coin-or.org">ADOL-C@list.coin-or.org</a><br>
> <a href="http://list.coin-or.org/mailman/listinfo/adol-c" target="_blank">http://list.coin-or.org/mailman/listinfo/adol-c</a><br>
><br>
<br>
--<br>
Dr. Kshitij Kulshreshtha<br>
<br>
Institut für Mathematik,<br>
Universität Paderborn,<br>
Warburger Straße 100,<br>
33098 Paderborn.<br>
<br>
Büro: A3.235<br>
<br>
Privatanschrift:<br>
Arnikaweg 62<br>
33100 Paderborn.<br>
_______________________________________________<br>
ADOL-C mailing list<br>
<a href="mailto:ADOL-C@list.coin-or.org">ADOL-C@list.coin-or.org</a><br>
<a href="http://list.coin-or.org/mailman/listinfo/adol-c" target="_blank">http://list.coin-or.org/mailman/listinfo/adol-c</a></blockquote></div><br><br clear="all"><br>-- <br>Antonello Lobianco<br>INRA, Laboratoire d'Economie Forestière<br>
14 Rue Girardet - 54000 Nancy, France<br>Tel: +33.652392310<br>Email: <a href="mailto:antonello.lobianco@nancy-engref.inra.fr" target="_blank">antonello.lobianco@nancy-engref.inra.fr</a><br><a href="http://antonello.lobianco.org/" target="_blank">http://antonello.lobianco.org</a><br>