<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: Times New Roman; font-size: 12pt; color: #000000'><style>p { margin: 0; }</style><div style="font-family: Times New Roman; font-size: 12pt; color: #000000">I'm not sure what you mean by 'coding by hand each algorithm for each thread'. <div><br></div><div>Do you want an entire Vec/Mat on each thread?</div><div>This is not how our model works currently and doing so will require considerable time and effort and may not be the worth it.</div><div><br></div><div>The way the pthread model in PETSc works is that every pthread vector/matrix uses threads from a thread pool for its operations. You can control how many and which threads you want from the thread pool.<br><div><br></div><div>Shri<br><br><hr id="zwchr"><blockquote style="border-left:2px solid rgb(16, 16, 255);margin-left:5px;padding-left:5px;">Jed, thanks for your response.<br><br><div class="gmail_quote">2012/3/19 Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im"><div class="gmail_quote">On Mon, Mar 19, 2012 at 15:47, Pedro Torres <span dir="ltr"><<a href="mailto:torres.pedrozpk@gmail.com" target="_blank">torres.pedrozpk@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>I has some experience with PETSc + MPI, but neither with PETSc+Pthreads nor Pthreads alone. Taking in to account that pthread was recenltly added to petsc, is it possible to code a program that solve the same linear system with two different solvers, one on each thread?. I guess that can be done calling different ksp object, one on each thread, but I don't have clear, first if this is possible and second how to do that. Sorry if my question is meaningless, I just started to study pthreads.</div>
<div><div><br></div><div>I will really appreciated any advice. Thanks in advance.</div></div></blockquote></div><br></div><div>No, the threading model does not work like this (and the performance of a system that did work that way would be a delicate beast).</div>
</blockquote><div> </div></div>Ok, so the only way to implement that is coding by hand each algorithm for each thread, rigth?. <div><br></div><div>Anyway, at least can I use petsc object like Vec/Mat on each thread?.</div>
<div><br></div><div>Thanks for your help!<br><br clear="all"><div><br></div>-- <br>Pedro Torres<br>
</div>
</blockquote><br></div></div></div></div></body></html>