<div dir="ltr">Hi Andre,<div><br></div><div style>TAO specializes in unconstrained and bound-constrained optimization, there is not a lot of support for linear constrained optimization.</div><div style>There is an interior point solver (ipm) that can accept general constraints, there are new functions for setting up these constraints, they aren't quite solid yet, and I've not very experienced in using the petsc4py package, but I can give you some general help.</div>
<div style><br></div><div style>There is a simple C example using nonlinear constraints in src/tao/examples/tutorials/toy.c, you should be able to easily  modify this for your example, where the equality constraint function will evaluate Ax-b and the equality jacobian will be the A matrix (you won't need to set the inequality constraint or jacobian). Because the ipm method builds a KKT system and solves it, it doesn't work well with iterative methods, a direct solver like superlu may be necessary.</div>
<div style><br></div><div style>I don't know enough about the actual python bindings to give you an example program in python, but it should follow pretty directly from the C example. Please let me know if you have any specific questions.</div>
<div style><br></div><div style>Jason Sarich</div><div style><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, May 14, 2014 at 9:29 AM, André Timótheo <span dir="ltr"><<a href="mailto:atmmachado@gmail.com" target="_blank">atmmachado@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div>
<div dir="ltr">p.s. I am not quite sure of the need of petsc-dev. I Just read about the inclusion of tao solvers in it and installed petsc-dev, cython and petsc4py on my ubuntu 12.04. 
<div><br>
</div>
<div>Unfortunately petsc4py (and the old tao4py) documentation does not have this type of python documented examples. After some correspondence in petsc4py mailing list they suggested that I might find some help on petsc-users mailing list.</div>

</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">2014-05-14 10:07 GMT-03:00 André Timótheo <span dir="ltr">
<<a href="mailto:atmmachado@gmail.com" target="_blank">atmmachado@gmail.com</a>></span>:<div><div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div style="font-family:arial,sans-serif;font-size:13px">
<div>
<div>
<div> I read about the merger of the TAO solvers on the PETSC-DEV. </div>
<div><br>
</div>
<div>    How can I use the TAO's constrainded optimization solver on </div>
<div>the PETSC-DEV (via petsc4py)? </div>
<div><br>
</div>
<div>    Can you show me some simple python script to deal with classical linear </div>
<div>constrainded optimization problems like:</div>
</div>
<br>
</div>
        minimize sum(x) subject to x >= 0 and Ax = b      </div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">    Thanks for your time.<br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Andre</div>
</div>
</blockquote>
</div></div></div>
<br>
</div>
</div>

</blockquote></div><br></div>