[deal.II] Received signal SIGSEGV, Segmentation fault

2023-12-29 Thread Smith Jack
Hello, I am trying to compile DEAL.II on a server. There is no error in the compilation, but when I try to run an example, segmentation fault occurs. I create a environment by conda, in this environment, I installed gcc, gxx, gfortan, make, cmake, lapack, blas, scalapack, petsc, symengine I also

Re: [deal.II] access all elements of a row of a trilinos sparse matrix

2023-12-29 Thread Wolfgang Bangerth
On 12/29/23 11:39, Simon Wiesheier wrote: For efficiency reasons, I decided to loop manually over the rows I am interested in. In a serial program, that worked smoothly. But if I get the gist of your message, vmult and extracting the relevant elements is the only way to do this in parallel, r

Re: [deal.II] access all elements of a row of a trilinos sparse matrix

2023-12-29 Thread Simon Wiesheier
Yes, basically I want to compute a matrix-vector product. But if my system has say 500k DoFs, I only need say 1k elements of the result of vmult, because the number of rows I am interested in is quite small -- they correspond to a small set of DoFs living on the boundary. For efficiency reasons, I

Re: [deal.II] access all elements of a row of a trilinos sparse matrix

2023-12-29 Thread Wolfgang Bangerth
On 12/28/23 12:22, Simon wrote: I want to build the scalar product between rows of a TrilinosWrappers::SparseMatrix and a TrilinosWrappers::MPI::Vector. This sounds like you're computing a matrix-vector product -- the dot products you compute are the elements of the resulting vector. As you