Re: [deal.II] Using solution on from one FE problem as a boundary condition for another

2020-01-20 Thread Ernesto Ismail
Hi Daniel and J-P, Thanks for your advice. I will try what you suggest. I feel like a bit of an idiot here, but I can't seem to find the source code for the interpolate_boundary_values() function. I've rooted around various source folders on Git and most of them don't contain much e.g. https:/

Re: [deal.II] IterativeInverse function appears to have disappeared between v8.5.0 and v9.0.0?

2020-01-20 Thread Daniel Arndt
Krishna, [...] > Can someone please point me in the right direction regarding the > documentation of this class for the latest release of dealii? > Have a look at https://www.dealii.org/current/doxygen/deal.II/changes_between_8_5_0_and_9_0_0.html, item 70. Best, Daniel -- The deal.II project i

[deal.II] IterativeInverse function appears to have disappeared between v8.5.0 and v9.0.0?

2020-01-20 Thread Krishnakumar Gopalakrishnan
I was just finishing Prof Bangerth's video lecture 21 wherein the concept of Schur Complement is introduced. For computing the inverse of the mass matrix (M), it looks like Dealii uses the facilities of the IterativeInverse class. The latest notes in the webpage for the lecture, https://www.ma

Re: [deal.II] discontinous contour over elements

2020-01-20 Thread Bruno Blais
Do you use Tecplot to calculate vorticity from the velocity field or do you calculate the vorticity from your code, and then visualize it from tecplot? The way deal.II visualizes gradients (or vorticity in this case) is the correct way it should be done, because it is visualized on an "element b

Re: [deal.II] Re: Application of inverse mass matrix in matrix-free context using cell_loop instead of scale()

2020-01-20 Thread 'Maxi Miller' via deal.II User Group
Found a bug in the library, reported here: https://github.com/dealii/dealii/issues/9405, and therefore I got wrong results. Am Sonntag, 19. Januar 2020 13:50:58 UTC+1 schrieb Maxi Miller: > > Hei, > I attached a working MWE. Changing between the approach suggested above > and scale() is done by

[deal.II] Difference between laplace-application through a matrix and through get_function_gradients()

2020-01-20 Thread 'Maxi Miller' via deal.II User Group
When applying the laplace operator to my solution vector, what is the difference between forming an explicit matrix containing the laplace operator by using template void LaplaceProblem::assemble_mass_matrix (){ TimerOutput::Scope t(computing_timer, "Mass matrix assembly");

Re: [deal.II] Different shape representations with manifolds on the same triangulation

2020-01-20 Thread Juan Carlos Araujo Cabarcas
Dear Jean-Paul, thanks again for your support and kind suggestions. I have worked with MappingQEulerian some time before, and as I remember, I dropped it because it requires the use of a vector field defined in the whole domain in order to curve geometries. After this, I started using ChartManif

Re: [deal.II] Different shape representations with manifolds on the same triangulation

2020-01-20 Thread Juan Carlos Araujo Cabarcas
Dear Jean-Paul, thanks again for your support and kind suggestions. I have worked with MappingQEulerian some time before, and as I remember, I dropped it because it requires a vector field define in the whole domain in order to curve the geometry. After this, I started using ChartManifolds and T

Re: [deal.II] Different shape representations with manifolds on the same triangulation

2020-01-20 Thread Juan Carlos Araujo Cabarcas
Dear Jean-Paul, thanks again for your support and kind suggestions. I have worked with MappingQEulerian some time before, and as I remember, I dropped it because it requires a vector field define in the whole domain in order to curve the geometry. After this, I started using ChartManifolds and T