Re: [deal.II] Speeding up characteristic interpolation

2017-09-06 Thread Praveen C
Hello Martin > Regarding the performance of your code, there are several fundamental > problems: > - You run the constructor of Quadrature (allocates memory) and FEValues > (allocates lots of memory, evaluates a number of things you don't need) that > are both not made for use within the

Re: [deal.II] Output nodal scalar

2017-09-06 Thread Jie Cheng
Hi Wolfgang I tried DataPostprocessorTensor as soon as the commits merged in. It is easy to use, I really like it! Thank you! Unfortunately there is still a bug that I couldn't figure out, which is basically the same as my previous post: after successfully writing the output, the program

Re: [deal.II] Fail to set up a custom pretty-printing for TrilinosWrappers::MPI::BlockVector

2017-09-06 Thread Qing Yin
Dear Dr. Bangerth, Thanks for your reply. I also didn't find a good way to debug it. Best, Qing 在 2017年9月1日星期五 UTC-6下午5:52:48,Wolfgang Bangerth写道: > > > Qing, > > > I use the deal.II pretty-printers > > and it works well for > > template

Re: [deal.II] Wrong solution at the boundary in restarting the solution

2017-09-06 Thread Wolfgang Bangerth
* Check the following: Is sol_u a vector that has ghost elements before serialization? Is sol_u a vector that has ghost elements after serialization? Because if these vectors don't have ghost elements, then what you see visually may not be what you actually have. The vectors

Re: [deal.II] Wrong solution at the boundary in restarting the solution

2017-09-06 Thread Dante De Santis
Thanks for the answer. > I can't right away see what the issue may be, but to help debug this: > * Change your code in such a way that just before serialization, you set > all entries of the solution vector to one. This way, you know exactly > what you have. What you get back should match