[deal.II] Re: AMR , how to pass solution vector to refined mesh

2017-03-16 Thread Jaekwang Kim
Thank you, Dr. Bangerth! Step-15 tutorial was really helpful for me and I tried to used solution transfer class. I refined grid as follow. Because I am using Picard iteration, I use it as a previous_solution vector (in first refined mesh) when it is interpolated into refined meshes. temp

[deal.II] Re: Installation using Spack fails during 'ncurses' installation

2017-03-16 Thread Jean-Paul Pelteret
Dear Stephen, Since this is an issue related not directly to deal.II but rather Spack itself (ncurses is only a dependency of a dealii related package), I would suggest that you post your question on their forum . You'll generally get a quick resp

[deal.II] Re: Installation using Spack fails during 'ncurses' installation

2017-03-16 Thread Bruno Turcksin
Steve, I have tried to use spack several times on two clusters and it never worked for me (but it works fine on my own machine). I usually have to patch a bunch of things and at the end, I still have problems when I load the modules. I find it a lot easier to install everything manually. You ca

[deal.II] Installation using Spack fails during 'ncurses' installation

2017-03-16 Thread Stephen DeWitt
Hello, I'm trying to install dealii on a shared AFS file system. Originally, I tried to install everything manually, which I had done successfully on our HPC cluster, but ran into a PETSc link error and decided to try Spack. I followed the instructions on the deal.II wiki, but installation stopp

Re: [deal.II] Re: Internal instability of the GMRES Solver / Trilinos

2017-03-16 Thread Martin Kronbichler
Dear Pascal, No, you do not need to try the other solution. I'm glad I could help. (This asserts the approach that we need to be careful with the vector pool between different calls.) Best, Martin On 16.03.2017 15:21, Pascal Kraft wrote: > Hi Martin, > > I have tried a version > with > Growin

Re: [deal.II] Re: Internal instability of the GMRES Solver / Trilinos

2017-03-16 Thread Pascal Kraft
Hi Martin, I have tried a version with GrowingVectorMemory::release_unused_memory() at the end of each step and removed my change to trilinos_vector.cc l.247 (back to the version from dealii source) and it seems to work fine. I have not tried the other solution you proposed, should I? Would

Re: [deal.II] Re: Internal instability of the GMRES Solver / Trilinos

2017-03-16 Thread Pascal Kraft
Dear Martin, my local machine is dying to a Valgrind run at the moment, but as soon as that is done with one step I will put these changes in right away and post the results here (<6 hrs). >From what I make of the call stacks on process somehow gets out of the SameAs() call without being MPI-bl

Re: [deal.II] Re: Internal instability of the GMRES Solver / Trilinos

2017-03-16 Thread Martin Kronbichler
Dear Pascal, You are right, in your case one needs to call GrowingVectorMemory::release_unused_memory() rather than for the vector. Can you try that as well? The problem appears to be that the call to SameAs returns different results for different processors, which it should not be, which is why