<div dir="ltr">Hi Barry,<div><br></div><div>Thanks. </div><div><br></div><div>I know how to solve a dense linear system in the petsc, but I was wondering how to call a Lapack routine in the petsc.</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 11, 2014 at 8:47 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
> On Dec 11, 2014, at 9:33 PM, Fande Kong <<a href="mailto:fdkong.jd@gmail.com">fdkong.jd@gmail.com</a>> wrote:<br>
><br>
> Hi all,<br>
><br>
> How to call a Lapack routine to solve a dense linear system? Any simple example?<br>
<br>
</span>   Create the Mat with MatCreateSeqDense() then create the usual KSP and use for the solver options -pc_type lu -ksp_type preonly<br>
<br>
    See src/ksp/ksp/examples/tutorials/ex30.c<br>
<br>
<br>
   Barry<br>
<br>
><br>
> Thanks,<br>
><br>
<br>
</blockquote></div></div>