Carlo,

what is the error that you get? The program crashes is too vague.

Best,

Bruno

On Wednesday, October 18, 2017 at 6:46:11 AM UTC-4, Carlo Marcati wrote:
>
>
> Hi, 
> I am trying to use a SolutionTrasfer object with 
> PETScWrappers::MPI::Vector. When trying to do 
> solution_trasfer.interpolate(in, out), the program crashes at the last line 
> of the 
> void SolutionTransfer<dim, VectorType, DoFHandlerType>::interpolate
> (const VectorType &in,
>  VectorType       &out) const
> function, i.e., at
>     out=all_out[0];
> This happens because all_out[0] is in the wrong state. 
>
> Writing
> all_out[0].compress (VectorOperation::insert);
> before that line would probably solve the issue, though it is obviously 
> specific to PETSc vectors. 
>
> The version of the interpolate function that takes vectors of vectors as 
> input, i.e.,
> void SolutionTransfer<dim, VectorType, DoFHandlerType>::
> interpolate (const std::vector<VectorType> &all_in,
>              std::vector<VectorType>       &all_out) const
>  works fine.
>
> Best regards,
> Carlo
>
>
>

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
For more options, visit https://groups.google.com/d/optout.

Reply via email to