Parallel computing for Navier-Stokes
Hi,
I have difficulty finding proper examples of parallel execution of coupled Navier-Stokes problems (like https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/tree/master/test/freeflow/navierstokes/channel/3d_nonuniform). When I apply the typical workflow for enabling parallel run (as described in the handbook, by including <dumux/linear/amgbackend.hh>
and replacing the sequential solver with using LinearSolver = AMGBiCGSTABBackend<LinearSolverTraits<GridGeometry>>; auto linearSolver = std::make_shared<LinearSolver>(leafGridView, gridGeometry->dofMapper());
), I face a couple of errors similar to this:
/dune-istl/dune/istl/novlpschwarz.hh:80:42: error: no type named ‘ConstColIterator’ in ‘class Dune::MultiTypeBlockMatrix<Dune::MultiTypeBlockVector<Dune::BCRSMatrix<Dune::FieldMatrix<double, 3, 3>, std::allocator<Dune::FieldMatrix<double, 3, 3> > >, Dune::BCRSMatrix<Dune::FieldMatrix<double, 3, 1>, std::allocator<Dune::FieldMatrix<double, 3, 1> > > >, Dune::MultiTypeBlockVector<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 3>, std::allocator<Dune::FieldMatrix<double, 1, 3> > >, Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > > > >’
I checked it with both MomentumGridGeometry
and MassGridGeometry
, both resulting in the same error. I see that there are quite a few examples and tests for AMGBiCGSTABBackend
for Stokes free flow and flow in porous media, but I couldn't find any relevant thing for NS where mass and momentum equations are coupled. Is parallel AMG backend not supported for such coupled problems?
I'm checking all these things with DuMux 3.6.