<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Dec 11, 2014 at 9:54 PM, Fande Kong <span dir="ltr"><<a href="mailto:fdkong.jd@gmail.com" target="_blank">fdkong.jd@gmail.com</a>></span> wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><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></blockquote><div><br></div><div>This is from dt.c:</div><div><br></div><div><div>#include <petscblaslapack.h></div></div><div><br></div><div>  PetscBLASInt LDZ, N;</div><div><br></div><div><div>  ierr = PetscBLASIntCast(npoints,&N);CHKERRQ(ierr);</div><div>  LDZ  = N;</div><div>  ierr = PetscFPTrapPush(PETSC_FP_TRAP_OFF);CHKERRQ(ierr);</div><div>  PetscStackCallBLAS("LAPACKsteqr",LAPACKsteqr_("I",&N,x,w,Z,&LDZ,work,&info));</div><div>  ierr = PetscFPTrapPop();CHKERRQ(ierr);</div><div>  if (info) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_PLIB,"xSTEQR error");</div></div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><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:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span><br>
> On Dec 11, 2014, at 9:33 PM, Fande Kong <<a href="mailto:fdkong.jd@gmail.com" target="_blank">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>
</blockquote></div><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>