[petsc-users] Hypre AMS usage

[Pierre Jolivet]

> On 6 Sep 2024, at 6:56 PM, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Fri, Sep 6, 2024 at 11:37 AM Karthikeyan Chockalingam - STFC UKRI via petsc-users <petsc-users at mcs.anl.gov <mailto:petsc-users at mcs.anl.gov>> 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
>> 
>> karthikeyan.chockalingam at stfc.ac.uk <mailto:karthikeyan.chockalingam at stfc.ac.uk>
>>   
>> 
>>  <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$>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240906/8b755a54/attachment.html>