Hi Daniel, > On 3 Jan 2017, at 15:09, Daniel Arndt <[email protected]> wrote: > > Do you have a minimal working example that shows this problem?
I don’t have a MWE for this, not yet at least. The tricky part is that this happens for some usage cases, but the same code runs fine in other, more-or-less as complicated, usage cases. > > Some ideas: > Can you confirm that the output vector you are using is "clean" in the sense > that there are no spurious ghost values before handing it to interpolate? I call interpolate() right after calling setup_system() in my application, which in turn does NOT assign any values to those vectors. Within the setup_system() the initialization of output vectors is done via MatrixFree<dim,double>::initialize_dof_vector(solution_vectors[i]); I added solution_vectors[i] = 0.; right after initialization but this does not change anything, neither does solution_vectors[i].zero_out_ghosts() > In other words: Does output.zero_out_ghosts() or output=0. help? Do you get > the same error if you are using a Trilinos or PETSc vector as output? Luckily I have a version of the main class which uses Trilinos (same p::d::Tria, 4 MPI cores and 2 threads each), when using it I do not have this error. BUT it does not mean this is not happening, maybe Trilinos parallel vectors don’t have this assert at all? > > Given that we test SolutionTransfer in the testsuite, I would suspect that > this is related to your setup of the output vector. So far I don’t see any indication that this is the case. I also tried running with a single thread, but this does not influence the results. p.s. the nature of application is such that I can not simply switch from deal.ii to Trilinos vectors. So the only way to have some MWE is to serialize p::d::Tria and write out vectors in ASCII or something. Regards, Denis. -- 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.
