[deal.II] Re: Data exchanges in IO

Bruce, The p_local_solution vector is filled with values as the result of solving > a linear equation, p_local_rho is filled with values by looping over > individual elements and assigning values calculated from other quantities. > Part of the reason the p_local_rho vector doesn't have ghost

[deal.II] error when trying to initialize PETSc preconditioners

Hi, everybody. I am just trying to initialize either ICC, Jacobi or ParaSalis preconditioner from PETSc, using a PETSc matrix in parallel. So far I am running just one mpi process. After the call to ''.initialize'' I get: [0]PETSC ERROR: No support for this operation for this object type

Re: [deal.II] error when trying to initialize PETSc preconditioners

Hey Ignacio, good to hear from you! Sadly this is not documented well, but the ICC preconditioner is only implemented for BAIJ matrices and not MPIAIJ (which we use even if you run on one processor), see https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCICC.html You should be

[deal.II] Re: Vector-valued gradient of solution vector

Dear Daniel, Thank you for a fast answer. I think I got the matrix and right-hand side correct now, please see test.cc for the example code. What I want to from the start was to get the vector-valued gradient from a scalar field, but I could not get this to work. See section "Do not work" in

Re: [deal.II] error when trying to initialize PETSc preconditioners

Many thanks Timo. I tried PreconditionBlockJacobi and it worked. However, I am running a validation loop of progressively smaller meshes, and stops again with a error when I make the call to initialize for the second time (I have to do it, since the matrix is the same ... but has changed, I

[deal.II] Daniel Arndt and Jean-Paul Pelteret appointed deal.II Developers

All, I want to share the good news that Daniel Arndt and Jean-Paul Pelteret have accepted our invitation to become deal.II Developers. They are now listed at https://dealii.org/authors.html Both of them have provided patches throughout the library for several years already, but you all

[deal.II] Re: Vector-valued gradient of solution vector

Joel, Yes the matrix you are assembling is a vector-valued mass matrix now. For me your code is failing with void dealii::FEValuesBase::get_function_gradients(const InputVector&, std::vector >&) const [with

[deal.II] Data exchanges in IO

Hi, I'm trying to write out some data from a calculation being run in parallel. I've got two vectors, p_local_solution and p_local_rho of type PetscWrappers::MPI::Vector. The p_local_solution vector is initialized using a statement