Re: [petsc-users] PETSC Matrix debugging

2024-04-01 Thread Shatanawi, Sawsan Muhammad via petsc-users
Hello Barry and Matthew, Yes MatView is what I was looking for. Thank you for your help. Bests, Sawsan From: Barry Smith Sent: Monday, April 1, 2024 1:14 PM To: Matthew Knepley Cc: Shatanawi, Sawsan Muhammad ; petsc-users@mcs.anl.gov Subject: Re:

Re: [petsc-users] ex19: Segmentation Violation when run with MUMPS on MacOS (arm64)

2024-04-01 Thread Satish Balay via petsc-users
On Mon, 1 Apr 2024, Zongze Yang wrote: > Thank you for your update. > > I found some links that suggest this issue is related to the Apple linker, > which is causing problems with Fortran linking. > > 1. >

Re: [petsc-users] PETSC Matrix debugging

2024-04-01 Thread Barry Smith
Note, you can also run with the option -mat_view and it will print each matrix that gets assembled. Also in the debugger you can do call MatView(mat,0) > On Apr 1, 2024, at 2:18 PM, Matthew Knepley wrote: > > This Message Is From an External Sender > This message came from outside your

Re: [petsc-users] PETSC Matrix debugging

2024-04-01 Thread Matthew Knepley
On Mon, Apr 1, 2024 at 1:57 PM Shatanawi, Sawsan Muhammad via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hello everyone, I hope this email finds you well. Is there a way we can > check how the matrix looks like after setting it. I have tried debugging it > with gdb- break points- and print

[petsc-users] PETSC Matrix debugging

2024-04-01 Thread Shatanawi, Sawsan Muhammad via petsc-users
Hello everyone, I hope this email finds you well. Is there a way we can check how the matrix looks like after setting it. I have tried debugging it with gdb- break points- and print statements, but it only gave me one value instead of a matrix. Thank you in advance for your time and

Re: [petsc-users] ex19: Segmentation Violation when run with MUMPS on MacOS (arm64)

2024-04-01 Thread Zongze Yang
Thank you for your update. I found some links that suggest this issue is related to the Apple linker, which is causing problems with Fortran linking. 1.

Re: [petsc-users] ex19: Segmentation Violation when run with MUMPS on MacOS (arm64)

2024-04-01 Thread Satish Balay via petsc-users
On Mon, 1 Apr 2024, Zongze Yang wrote: > > I noticed this in the config.log of OpenMPI: > ``` > configure:30230: checking to see if mpifort compiler needs additional linker > flags > configure:30247: gfortran -o conftest -fPIC -ffree-line-length-none > -ffree-line-length-0

Re: [petsc-users] ex19: Segmentation Violation when run with MUMPS on MacOS (arm64)

2024-04-01 Thread Zongze Yang
> On 1 Apr 2024, at 23: 38, Satish Balay wrote: > > On Sun, 31 Mar 2024, Zongze Yang wrote: >>> --- >>> petsc@ npro petsc % ./configure --download-bison --download-chaco --download-ctetgen ZjQcmQRYFpfptBannerStart This Message Is From an External Sender

Re: [petsc-users] Error loading data coming from a .hdf5 file into a DMSwarm

2024-04-01 Thread MIGUEL MOLINOS PEREZ
I see. Thank you for the feedback Matthew! Thanks, Miguel On 1 Apr 2024, at 18:32, Matthew Knepley wrote: On Mon, Apr 1, 2024 at 11:13 AM MIGUEL MOLINOS PEREZ mailto:mmoli...@us.es>> wrote: Dear Matthew, Thank you for your suggestion. I tried to update the vector with the information coming

Re: [petsc-users] Error loading data coming from a .hdf5 file into a DMSwarm

2024-04-01 Thread Matthew Knepley
On Mon, Apr 1, 2024 at 11:13 AM MIGUEL MOLINOS PEREZ wrote: > Dear Matthew, > > Thank you for your suggestion. I tried to update the vector with the > information coming from the hdf5 file inside the main function. Then I > print the vector two times (see the lines below), the first time it has

Re: [petsc-users] Error loading data coming from a .hdf5 file into a DMSwarm

2024-04-01 Thread MIGUEL MOLINOS PEREZ
Dear Matthew, I came up with a workaround for the problem. I duplicate the vector and use the duplicated copy to read the information from the hdf5 file. Then I swap both vectors and delete the copy. If I invoke VecView outside of the function, the value has been modified properly. However,

Re: [petsc-users] ex19: Segmentation Violation when run with MUMPS on MacOS (arm64)

2024-04-01 Thread Satish Balay via petsc-users
On Sun, 31 Mar 2024, Zongze Yang wrote: > > --- > > petsc@npro petsc % ./configure --download-bison --download-chaco > > --download-ctetgen --download-eigen --download-fftw --download-hdf5 > > --download-hpddm --download-hwloc --download-hypre --download-libpng > > --download-metis

Re: [petsc-users] Error loading data coming from a .hdf5 file into a DMSwarm

2024-04-01 Thread MIGUEL MOLINOS PEREZ
Dear Matthew, Thank you for your suggestion. I tried to update the vector with the information coming from the hdf5 file inside the main function. Then I print the vector two times (see the lines below), the first time it has the correct data. However, the second time, it has the same values

Re: [petsc-users] Error loading data coming from a .hdf5 file into a DMSwarm

2024-04-01 Thread Matthew Knepley
On Sun, Mar 31, 2024 at 4:08 PM MIGUEL MOLINOS PEREZ wrote: > Dear all, I am writing a function which store datasets (Vectors) coming > from a DMSwarm structure into a hdf5 file. This step is done nicely > write_function(){ PetscViewerHDF5Open(…) PetscViewerHDF5PushTimestepping(…) >