Hi<br>There was a discussion between Tim Stitt and petsc developers about matrix inversion, and it was really helpful. That was in last Nov. You can check the emails archive<br><br>http://www-unix.mcs.anl.gov/web-mail-archive/lists/petsc-users/2007/11/threads.html<br><br>Waad<br><br><b><i>Yujie <recrusader@gmail.com></i></b> wrote:<blockquote class="replbq" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"> what is the difference between sequantial and parallel AIJ matrix? Assuming there is a matrix A, if I partitaion this matrix into A1, A2, Ai... An.<br>A is a parallel AIJ matrix at the whole view, Ai is a sequential AIJ matrix? I want to operate Ai at each node. <br> In addition, whether is it possible to get general inverse using MatMatSolve() if the matrix is not square? Thanks a lot.<br> <br>Regards,<br>Yujie<br> <br><br><div><span class="gmail_quote">On
2/4/08, <b class="gmail_sendername">Barry Smith</b> <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"> <br> For sequential AIJ matrices you can fill the B matrix with the<br>identity and then use<br>MatMatSolve().<br><br> Note since the inverse of a sparse matrix is dense the B matrix is<br>a SeqDense matrix.<br><br> Barry<br><br>On Feb 4, 2008, at 12:37 AM, Yujie wrote:<br><br>> Hi,<br>> Now, I want to inverse a sparse matrix. I have browsed the manual,<br>> however, I can't find some information. could you give me some advice?<br> ><br>> thanks a lot.<br>><br>> Regards,<br>> Yujie<br>><br><br></blockquote></div><br> </blockquote><br><p> 
<hr size=1>Looking for last minute shopping deals? <a href="http://us.rd.yahoo.com/evt=51734/*http://tools.search.yahoo.com/newsearch/category.php?category=shopping">
Find them fast with Yahoo! Search.</a>