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
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 *
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|