[petsc-users] Question about the performance of KSP solver

[Jed Brown]

Hi, Sanjay. The method uses matrix-free p-multigrid with AMG coarse solves. Attached is one of the models we're using as a test problem (hyperelastic model on an extruded Schwarz-P surface). Docs for the code we're using.

https://ratel.micromorph.org/
https://gitlab.com/micromorph/ratel
https://doi.org/10.21105/joss.02945 (about libCEED's data structures)

We'd love to test these methods out on more problems. We're working on improving integration with PETSc -- the goal is that users of PETSc will be able to write libCEED qfunctions and run with these data structures (which are optimized for both CPU and GPU) with almost no other code modifications.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot from 2022-02-27 13-23-50.png
Type: image/png
Size: 2436118 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220227/3515a5a6/attachment-0001.png>
-------------- next part --------------

Sanjay Govindjee <s_g at berkeley.edu> writes:

> Hi Jed,
>  ? Do you have a reference for this?
> -sanjay
>
> On 2/27/22 7:16 AM, Jed Brown wrote:
>> If you're interested in using different data structures, our experience is that we can solve similar problem sizes using Q2 elements in a few seconds (2-10) on a single node.