[petsc-users] valgrind with petscmpiexec

2020-12-14 Thread Fande Kong
Hi All, I tried to use valgrind to check if the simulation is valgrind clean because I saw some random communication fails during the simulation. I tried this command-line petscmpiexec -valgrind -n 576 ../../../moose-app-oprof -i input.i -log_view -snes_view But I got the following error

Re: [petsc-users] Transform of algorithm containing zgemv/FFT/Slicing into PETSc-functions

2020-12-14 Thread Barry Smith via petsc-users
I think you can again use MatDenseGetArray() and do the copies directly respecting the shift that you desire. Each process will just do the local rows so you need not worry about parallelism. I think it may be as simple as get the array pointer for A, shift it by number of local rows *

[petsc-users] Transform of algorithm containing zgemv/FFT/Slicing into PETSc-functions

2020-12-14 Thread Roland Richter
Dear all, I am currently working on the transformation of an algorithm implemented using armadillo into PETSc. It is a forward/backward transformation, and boils down to the following steps (for the forward transformation): Assumed I have matrices A and B, defined as A =  |aa ab ac ad|