[petsc-users] R: local/global DMPlex Vec output

2024-10-24 Thread Semplice Matteo via petsc-users
Hi,   I tried again today and have (re?)discovered this example https://urldefense.us/v3/__https://petsc.org/release/src/dm/impls/plex/tutorials/ex14.c.html__;!!G_uCfscf7eWS!fJpyioZNC3oN15Qm8-PTF_nxD_NYc-Ou_-moVFROjwZyZREL7e0PwP-rDwjCcynulDAeLGL1WqkN7tyuZepLZ0L7BfI88y7favtI0g$ , but I cannot

Re: [petsc-users] Block sparse matrix

2024-10-24 Thread Barry Smith
Perhaps this is what you are looking for https://urldefense.us/v3/__https://petsc.org/release/manualpages/Mat/MATNEST/*matnest__;Iw!!G_uCfscf7eWS!e9L1JRLW7O-w_LDeaQomJxREvfqcm0FrZtlpjeE46ux3HjMJELyJmp9jx4uTJkK7WyQc57DtJFcC3nFe6jEJpOU$ ? > On Oct 24, 2024, at 8:48 PM, Diego Magela Lemos via p

Re: [petsc-users] Block sparse matrix

2024-10-24 Thread Barry Smith
> On Oct 24, 2024, at 9:02 PM, Matthew Knepley wrote: > > On Thu, Oct 24, 2024 at 8:49 PM Diego Magela Lemos via petsc-users > mailto:petsc-users@mcs.anl.gov>> wrote: >> Is there something similar to SciPy sparse block_array >> (https://urldefense.us/v3/__https://docs.scipy.org/doc/scipy/refe

Re: [petsc-users] Block sparse matrix

2024-10-24 Thread Matthew Knepley
On Thu, Oct 24, 2024 at 8:49 PM Diego Magela Lemos via petsc-users < petsc-users@mcs.anl.gov> wrote: > Is there something similar to SciPy sparse block_array ( > https://urldefense.us/v3/__https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.block_array.html*scipy.sparse.block_array__

[petsc-users] Block sparse matrix

2024-10-24 Thread Diego Magela Lemos via petsc-users
Is there something similar to SciPy sparse block_array ( https://urldefense.us/v3/__https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.block_array.html*scipy.sparse.block_array__;Iw!!G_uCfscf7eWS!YANe4IC4s7Ki7CYmRbt0miQwrrFq7e-6_btfizr9d4vir0dzWIncCQBkH3Cbr94Y_jxGNSAhJh59aSOHAHRESglg

Re: [petsc-users] local/global DMPlex Vec output

2024-10-24 Thread Matthew Knepley
On Thu, Oct 24, 2024 at 6:04 PM Matteo Semplice < matteo.sempl...@uninsubria.it> wrote: > Hi. The HDF5 solution looks good to me, but I now get this error > > Okay, I can make a workaround for this. Here is what is happening. When you output solutions, you really want the essential boundary condi

Re: [petsc-users] local/global DMPlex Vec output

2024-10-24 Thread Matteo Semplice via petsc-users
Hi. The HDF5 solution looks good to me, but I now get this error $ ../src/saveDemo Creating mesh with (10,10) faces [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Object is in wrong state [0]PETSC ERROR: Need to

Re: [petsc-users] local/global DMPlex Vec output

2024-10-24 Thread Matthew Knepley
I just looked at the code. The VTK code is very old, and does not check for cell overlap. We have been recommending that people use either HDF5 or CGNS, both of which work in this case I believe. I can fix VTK if that is what you want, but it might take me a little while as it is very busy at work

Re: [petsc-users] Problem running geoclaw example with petsc 3.22

2024-10-24 Thread Barry Smith
Good, some helpful information with the runs you made. In the crash below it made more progress it was able to allocate multiple regions of shared memory and access them. I don't know why it would crash later. Can you tell me all the steps with miniforge (which seems to be related to