Dear Matt:<br><br>I checked the codes about MatMatSolve(). However, currently, PETSc didn't realize its parallel version. Is it right? I want to inverse the matrix parallelly. could you give me some examples about it? thanks a lot.<br>
<br>Regards,<br>Yujie<br><br><div><span class="gmail_quote">On 2/27/08, <b class="gmail_sendername">Matthew Knepley</b> <<a href="mailto:knepley@gmail.com">knepley@gmail.com</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">
On Wed, Feb 27, 2008 at 11:05 AM, Yujie <<a href="mailto:recrusader@gmail.com">recrusader@gmail.com</a>> wrote:<br> > Dear Sanjay:<br> ><br> > Thank you for your reply. I don't understand what you said. Now, I want to<br>
> use spooles package to inverse a sparse SPD matrix. I have further checked<br> > the inferface about spooles in PETSc. I find although spooles can deal with<br> > AX=B (B may be a dense matrix) with parallel LU factorization.<br>
> However, PETSc only provide the following:<br> > 51: PetscErrorCode MatSolve_MPISpooles(Mat A,Vec b,Vec x)<br> > I don't set b to a matrix even if I use<br> > 178: PetscErrorCode MatFactorNumeric_MPISpooles(Mat A,MatFactorInfo<br>
> *info,Mat *F) for LU factorization.<br> ><br> > Could you have any suggestions about this? thanks a lot.<br><br><br>MatMatSolve()<br><br> <a href="http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatMatSolve.html">http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatMatSolve.html</a><br>
<br> Matt<br><br><br> > Regards,<br> > Yujie<br> ><br> > On 2/26/08, Sanjay Govindjee <<a href="mailto:sanjay@ce.berkeley.edu">sanjay@ce.berkeley.edu</a>> wrote:<br> > > from my make file<br> > ><br>
> > -@${MPIRUN} -s all -np $(NPROC) $(PROGNAME) -ksp_type preonly<br> > > -ksp_monitor -pc_type cholesky -mat_type mpisbaijspooles -log_summary<br> > > -on_error_attach_debugger -mat_spooles_symmetryflag 0 -options_left<br>
> ><br> > ><br> > > -sg<br> > ><br> > ><br> > > Yujie wrote:<br> > > > Hi, everyone<br> > > ><br> > > > I have compiled PETSc with spooles. However, I try to find how to use<br>
> > > this package in PETSc directory. I can't find any examples for it.<br> > > > Could you give me some advice? I want to use spooles to inverse a<br> > > > sparse matrix. thanks a lot.<br>
> > ><br> > > > Regards,<br> > > > Yujie<br> > ><br> > ><br> ><br> ><br><br><br><br><br>--<br> What most experimenters take for granted before they begin their<br> experiments is infinitely more interesting than any results to which<br>
their experiments lead.<br> -- Norbert Wiener<br><br></blockquote></div><br>