<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Faraz:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">I reviewed the code for how the matrix is being solved in Intel Paridso. It seems the matrix is getting reordered and sorted before calling Pardiso. I am not sure the reason for this. I know I need to ask the original developer of the application, but they are no longer responding.<br>
<br>
Is there something I am misunderstanding about how matrix's need to be re-ordered? I am assuming the metis or scotch algorithm does the reordering for me when using Mumps? Or am I supposed to reorder/sort it myself before calling Mumps?<br></blockquote><div> </div><div>Mumps uses parmetis or scotch for parallel symbolic factorization. For sequential symbolic factorization, it has several matrix orderings, which you can experiment with the option '-mat_mumps_icntl_7 <>'.</div><div>I doubt any of these ordering would match the performance of Pardiso.</div><div><br></div><div>Again, how large is your matrix? How do you run Pardiso in parallel? Can you use Pardiso on large matrices as mumps?</div><div><br></div><div>Hong</div></div><br></div></div>