Re: [deal.II] solving stabilized Stokes

2020-02-18 Thread Wolfgang Bangerth
On 2/14/20 5:34 AM, Richard Schussnig wrote: Could the observed behaviour be caused by applying an iterative solver*? (using ReductionControl and a reduction factor of 1e-13 which is really low) * Block-triangular Schur-complement-based approach, similar as in step-55 suggested in the further

[deal.II] Re: solving stabilized Stokes

2020-02-18 Thread Bruno Blais
Dear Richard, We implement a quasi identical scheme in lethe (PSPG for continuity and SUPG for momentum). You can find the code on the github, it might help you : https://github.com/lethe-cfd/lethe/blob/master/source/solvers/gls_navier_stokes.cc I can confirm that you should be getting order p+1

[deal.II] Re: Interpolating to globally refined distributed mesh

2020-02-18 Thread Konrad Simon
Add-on: The problem is in the function RefineInterpolate::refine_and_interpolate_on_distributed_mesh() -- 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 su

[deal.II] Re: Interpolating to globally refined distributed mesh

2020-02-18 Thread Konrad Simon
Add-on: The problem is in the function RefineInterpolate::refine_and_interpolate_on_distributed_mesh() -- 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 su

[deal.II] Interpolating to globally refined distributed mesh

2020-02-18 Thread Konrad Simon
Dear deal.ii community, I am trying to interpolate a distributed solution vector onto a globally refined distributed mesh. Actually I need to globally refine more than once but I believe I can iterate this procedure (is there a more efficient way? The finite element can be vector or scalar valu

Re: [deal.II] Question on how MatrixTools::MatrixCreator sets the appropriate flags of FeValues?

2020-02-18 Thread Daniel Arndt
Krishna, Adding to David's answer, mass matrix and laplace matrix are built for multiple components separately, i.e. there is no coupling between them. Best, Daniel Am Di., 18. Feb. 2020 um 10:17 Uhr schrieb David Wells : > Hi Krishna, > > I'm not sure that I understand your question - those tw

Re: [deal.II] Question on how MatrixTools::MatrixCreator sets the appropriate flags of FeValues?

2020-02-18 Thread David Wells
Hi Krishna, I'm not sure that I understand your question - those two functions create a mass matrix and a laplace corresponding to the given parameters (dof handler, mapping, quadrature, etc) - in the first case only shape values are needed and in the second only shape gradients. Hence the functio

[deal.II] Question on how MatrixTools::MatrixCreator sets the appropriate flags of FeValues?

2020-02-18 Thread Krishnakumar Gopalakrishnan
In the tutorials, Step-23 suddenly gets rid of the assemble_system() function and introduces MatrixTools::MatrixCreator::create_mass_matrix and MatrixTools::MatrixCreator::create_laplace_matrix However, there is no discussion whatsoever of how the appropriate update flags (bitfields) of the Fev