> On 6 Sep 2024, at 6:56 PM, Matthew Knepley <[email protected]> wrote: > > On Fri, Sep 6, 2024 at 11:37 AM Karthikeyan Chockalingam - STFC UKRI via > petsc-users <[email protected] <mailto:[email protected]>> wrote: >> Hello, >> >> >> >> I am trying to use the Hypre AMS preconditioner for the first time. >> >> >> >> I am following the example problem from Hypre >> >> https://urldefense.us/v3/__https://github.com/hypre-space/hypre/blob/3caa81955eb8d1b4e35d9b450e27cf6d07b50f6e/src/examples/ex15.c*L954__;Iw!!G_uCfscf7eWS!Z9zog8mGTxJtHDwPnjywf0ppfrJPq23lCLWeE3gO-FYd_QgXq8S0g41HEHnsEbQ0YOFVE8vmL-OlCwyXA0U3cw$ >> >> <https://urldefense.us/v3/__https://github.com/hypre-space/hypre/blob/3caa81955eb8d1b4e35d9b450e27cf6d07b50f6e/src/examples/ex15.c*L954__;Iw!!G_uCfscf7eWS!biuCOXkV-qD-e8iTXiwaD9XJ0EWuqr2PidZJdPqbnFzCshQFo-mclGSvGBOPuGSTkjPvMKnIfJzWos7QShcSitpwyrhL5ySgA0XM$> >> >> >> I have so far successfully set the discrete gradient operator and vertex >> co-ordinates, >> >> >> >> // Set discrete gradient >> >> petscErr = PCHYPRESetDiscreteGradient(pc, par_G); >> >> >> >> // Set vertex coordinates >> >> petscErr = PCHYPRESetEdgeConstantVectors(pc, par_xvec, par_yvec, par_zvec); >> >> >> >> Do I need to set the following smoothing options? >> >> https://urldefense.us/v3/__https://github.com/hypre-space/hypre/blob/3caa81955eb8d1b4e35d9b450e27cf6d07b50f6e/src/examples/ex15.c*L965__;Iw!!G_uCfscf7eWS!Z9zog8mGTxJtHDwPnjywf0ppfrJPq23lCLWeE3gO-FYd_QgXq8S0g41HEHnsEbQ0YOFVE8vmL-OlCwwsyo4fRw$ >> >> <https://urldefense.us/v3/__https://github.com/hypre-space/hypre/blob/3caa81955eb8d1b4e35d9b450e27cf6d07b50f6e/src/examples/ex15.c*L965__;Iw!!G_uCfscf7eWS!biuCOXkV-qD-e8iTXiwaD9XJ0EWuqr2PidZJdPqbnFzCshQFo-mclGSvGBOPuGSTkjPvMKnIfJzWos7QShcSitpwyrhL55yOMB-v$> > It cannot hurt. I would set them to begin with. > >> Also, do I need to convert from MATMPIAIJ to CSR? >> >> https://urldefense.us/v3/__https://github.com/hypre-space/hypre/blob/3caa81955eb8d1b4e35d9b450e27cf6d07b50f6e/src/examples/ex15.c*L984__;Iw!!G_uCfscf7eWS!Z9zog8mGTxJtHDwPnjywf0ppfrJPq23lCLWeE3gO-FYd_QgXq8S0g41HEHnsEbQ0YOFVE8vmL-OlCwzzrmDpiA$ >> >> <https://urldefense.us/v3/__https://github.com/hypre-space/hypre/blob/3caa81955eb8d1b4e35d9b450e27cf6d07b50f6e/src/examples/ex15.c*L984__;Iw!!G_uCfscf7eWS!biuCOXkV-qD-e8iTXiwaD9XJ0EWuqr2PidZJdPqbnFzCshQFo-mclGSvGBOPuGSTkjPvMKnIfJzWos7QShcSitpwyrhL53Trv1si$> > No. > >> What are the other PETSc calls to invoke AMS? Is there an example problem I >> can look at? >>
-pc_hypre_type ams or PetscCall(PCHYPRESetType(pc, "ams")) Thanks, Pierre > I do not know. I don't think we have an example. > > Thanks, > > Matt > >> Thank you. >> >> >> >> Karthik. >> >> >> >> -- >> >> Karthik Chockalingam, Ph.D. >> >> Senior Research Software Engineer >> >> High Performance Systems Engineering Group >> >> Hartree Centre | Science and Technology Facilities Council >> >> [email protected] >> <mailto:[email protected]> >> >> >> <image001.png> >> >> >> > > > -- > What most experimenters take for granted before they begin their experiments > is infinitely more interesting than any results to which their experiments > lead. > -- Norbert Wiener > > https://urldefense.us/v3/__https://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!Z9zog8mGTxJtHDwPnjywf0ppfrJPq23lCLWeE3gO-FYd_QgXq8S0g41HEHnsEbQ0YOFVE8vmL-OlCwxghWiarQ$ > > <https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!Z1YtcjIl3JraPfN9GbBHhno-LcLOXJurvRZ3MKzzbqSMDLtq5x7_DICKdBj5idJnKPOPq7kmcpCvA1RwP5cb$>
