GMG for matrix-free VANS #1458
GMG for matrix-free VANS #1458
Conversation
| // AssertThrow( | ||
| // this->stabilization == | ||
| // Parameters::Stabilization::NavierStokesStabilization::pspg_supg, | ||
| // ExcMessage( | ||
| // "PSPG-SUPG stabilization is the only stabilization method currently |
There was a problem hiding this comment.
The variable this->stabilization is not initialized at this stage.
| const std::shared_ptr<Function<dim>> forcing_function, | ||
| const std::shared_ptr<SimulationControl> &simulation_control, | ||
| const std::shared_ptr<FESystem<dim>> fe); | ||
| MFNavierStokesPreconditionGMGBase( |
There was a problem hiding this comment.
The class hierarchy of the multigrid operators: MFNavierStokesPreconditionGMGBase -> MFNavierStokesPreconditionGMG -> MFNavierStokesVANSPreconditionGMG. The derivatives need to implement create_level_operator() and initialize(...).
| template <int dim> | ||
| void | ||
| FluidDynamicsMatrixFree<dim>::create_GMG() | ||
| { | ||
| gmg_preconditioner = std::make_shared<MFNavierStokesPreconditionGMG<dim>>( | ||
| this->simulation_parameters, | ||
| this->dof_handler, | ||
| this->dof_handler_fe_q_iso_q1); | ||
|
|
||
| gmg_preconditioner->reinit(this->mapping, | ||
| this->cell_quadrature, | ||
| this->forcing_function, | ||
| this->simulation_control, | ||
| this->fe); | ||
| } | ||
|
|
||
| template <int dim> | ||
| void | ||
| FluidDynamicsMatrixFree<dim>::initialize_GMG() | ||
| { | ||
| dynamic_cast<MFNavierStokesPreconditionGMG<dim> *>(gmg_preconditioner.get()) | ||
| ->initialize(this->simulation_control, | ||
| this->flow_control, | ||
| this->present_solution, | ||
| this->time_derivative_previous_solutions); | ||
| } |
There was a problem hiding this comment.
These are (new) function that can be overridden in the derivated classes (e.g., to create a different type of multgrid type and to pass additional information to the initialization methods).
|
At the moment, the void fraction is evaluated at every Newton step which is more that needed. However, we have the same issue with |
No description provided.