Hi, Tim<br><br>Thank you for your help. I am really glad to get your help.<br>According to what you said, if the matrix A has been divided into several nodes in the cluster, you may use your parallel code to inverse A? <br>
My problem is that what is the distribution of the results?<br> thanks a lot.<br><br>Regards,<br>Yujie<br><br><div><span class="gmail_quote">On 2/5/08, <b class="gmail_sendername">Timothy Stitt</b> <<a href="mailto:tstitt@cscs.ch">tstitt@cscs.ch</a>> wrote:</span><blockquote class="gmail_quote" style="margin-top: 0; margin-right: 0; margin-bottom: 0; margin-left: 0; margin-left: 0.80ex; border-left-color: #cccccc; border-left-width: 1px; border-left-style: solid; padding-left: 1ex">
Yes Yujie, I was able to put together a parallel code to invert a large<br>sparse matrix with the help of the PETSc developers. If you need any<br>help or maybe a Fortran code template just let me know.<br><br>Best,<br><br>
Tim.<br><br>Waad Subber wrote:<br>> Hi<br>> There was a discussion between Tim Stitt and petsc developers about<br>> matrix inversion, and it was really helpful. That was in last Nov. You<br>> can check the emails archive<br>
><br>> <a href="http://www-unix.mcs.anl.gov/web-mail-archive/lists/petsc-users/2007/11/threads.html">http://www-unix.mcs.anl.gov/web-mail-archive/lists/petsc-users/2007/11/threads.html</a><br>><br>> Waad<br>><br>
> */Yujie <<a href="mailto:recrusader@gmail.com">recrusader@gmail.com</a>>/* wrote:<br>><br>> what is the difference between sequantial and parallel AIJ matrix?<br>> Assuming there is a matrix A, if<br>
> I partitaion this matrix into A1, A2, Ai... An.<br>> A is a parallel AIJ matrix at the whole view, Ai<br>> 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<br>
> MatMatSolve() if the matrix is not square? Thanks a lot.<br>><br>> Regards,<br>> Yujie<br>><br>><br>> On 2/4/08, *Barry Smith* <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a><br>
> <mailto:<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>>> wrote:<br>><br>><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<br>> 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<br>> manual,<br>> > however, I can't find some information. could you give me<br>> some advice?<br>
> ><br>> > thanks a lot.<br>> ><br>> > Regards,<br>> > Yujie<br>> ><br>><br>><br>><br>> ------------------------------------------------------------------------<br>
> Looking for last minute shopping deals? Find them fast with Yahoo!<br>> Search.<br>> <<a href="http://us.rd.yahoo.com/evt=51734/*http://tools.search.yahoo.com/newsearch/category.php?category=shopping">http://us.rd.yahoo.com/evt=51734/*http://tools.search.yahoo.com/newsearch/category.php?category=shopping</a>><br>
<br><br><br>--<br>Timothy Stitt<br>HPC Applications Analyst<br><br>Swiss National Supercomputing Centre (CSCS)<br>Galleria 2 - Via Cantonale<br>CH-6928 Manno, Switzerland<br><br>+41 (0) 91 610 8233<br><a href="mailto:stitt@cscs.ch">stitt@cscs.ch</a><br>
<br></blockquote></div><br>