Hi Giang, > I did the git bisect (14 iterations) and it reveals that the commit 0eeeec59a54d7385a8e862fa114201e91c0f3c53 gives the problem.
For convenience, I linked to the commit <https://github.com/dealii/dealii/commit/0eeeec59a54d7385a8e862fa114201e91c0f3c53> with that specific hash. This commit is not part of the deal.ii 9.2 release, so it can not be the cause of the diverging solution as you found it with version 9.2. Further, that commit only changes parts of the `MatrixFree` implementation, which the code-gallery program does not use. I'm confused: So your solution was fine right before that particular commit? > In addition somewhere between commit d5ab9e084ae9aae301b121c22a5548539df7032e and 51c7b294f5f81c5c0f1ac1aa09670ea1729f563e makes the code run significantly slower. Did you notice the slowdown only in `Debug` mode or also in `Release` mode? Best, Marc On Monday, June 14, 2021 at 11:06:02 AM UTC-6 hgbk...@gmail.com wrote: > Hi Wolfgang > > Thanks for the hint. Luckily the old code compiled and runs fine with > v9.0.0. > > I did the git bisect (14 iterations) and it reveals that the commit > 0eeeec59a54d7385a8e862fa114201e91c0f3c53 gives the problem. > > In addition somewhere between > commit d5ab9e084ae9aae301b121c22a5548539df7032e > and 51c7b294f5f81c5c0f1ac1aa09670ea1729f563e makes the code run > significantly slower. > > Unfortunately, I could not thoroughly go into the code and investigate the > problem. This is only to find out which commit has an issue. > > The test is performed with a slightly modified version of the original > example (see attached file) in which > > KellyErrorEstimator<dim>::estimate(dof_handler, > face_quad_formula, > typename FunctionMap<dim>::type(), > present_solution, > estimated_error_per_cell, > fe.component_mask(velocity)); > > is replaced by > > KellyErrorEstimator<dim>::estimate(dof_handler, > face_quad_formula, > {}, > present_solution, > estimated_error_per_cell, > fe.component_mask(velocity)); > > This is because the FunctionMap is deprecated. I hope that it does not > affect the results. > > Best > Giang > > On Wednesday, June 9, 2021 at 10:43:44 PM UTC+2 Wolfgang Bangerth wrote: > >> On 6/5/21 9:02 AM, hgbk...@gmail.com wrote: >> > >> > I have tried to run this example with deal.II 9.1.0 and this problem >> still >> > persists. With deal.II 9.0.0/9.0.1 the example fails to compile because >> it >> > requires affine_constraints.h, which is not yet introduced: >> > >> > >> /home/hbui/workspace2/deal.II/time_dependent_navier_stokes/time_dependent_navier_stokes.cc:9:10: >> >> >> > fatal error: deal.II/lac/affine_constraints.h: No such file or >> directory >> > #include <deal.II/lac/affine_constraints.h> >> >> This class used to be called ConstraintMatrix and was defined in >> deal.II/lac/constraint_matrix.h. >> >> You could go back in the history of the NS gallery program to match the >> 9.0 >> release. The github repository for the code gallery is here: >> https://github.com/dealii/code-gallery >> The history of the .cc file is here: >> >> >> https://github.com/dealii/code-gallery/commits/master/time_dependent_navier_stokes/time_dependent_navier_stokes.cc >> >> >> Best >> W. >> >> -- >> ------------------------------------------------------------------------ >> Wolfgang Bangerth email: bang...@colostate.edu >> 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/f35fb9c6-1c2c-4124-9425-ebb0e4ee30cbn%40googlegroups.com.
