On 4/25/19 10:13 AM, 'Maxi Miller' via deal.II User Group wrote: > After running some tests, I ended up reducing the problem to the transfer to > and from the Epetra-vectors. I have to write an interface to the model > (according to the instructions), and as shown in the code above. There I have > Epetra-Vectors created by my interface, with a size of > locally_relevant_dofs.size(), and TrilinosWrappers::MPI::Vectors. Transfer > from the Epetra-Vectors to the MPI::Vectors works fine, but the way back does > not work (The MPI::Vectors are larger than the Epetra_Vectors). Are there any > hints in how I still could fit the data from the MPI::Vector into the > Epetra_Vector? Or should I rather ask this on the Trilinos mailing list?
Good question. I think it would probably be very useful to have small testcase others could look at. The program you have attached is still 1,300 lines, which is far too much for anyone to understand. But since you have an idea of what the problem is, do you think you could come up with a small testcase that illustrates exactly the issue you have? It doesn't have to do anything useful at all -- in your case, I think all that's necessary is to create a vector, convert it as you describe above, and then output the sizes to show that they are not the same. Run this on two processors, and you should have something that could be done in 50 or 100 lines, and I think that might be small enough for someone who doesn't know your code to understand what is necessary. I've built these sorts of testcases either from scratch in the past, or by just keep removing things from a program that (i) are run after the time the problem happens, and (ii) then removing everything that is not necessary. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.colostate.edu/~bangerth/ -- 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.
