New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Improvement of the matrix sparsity pattern of DG on single core #1100
Conversation
if (parameters.use_discontinuous_temperature_discretization | ||
&& parameters.use_limiter_for_discontinuous_temperature_solution) | ||
apply_limiter_to_dg_solutions(AdvectionField::temperature()); | ||
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This change is independent of the other one. Can you separate this into a different pull request?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
OK. I will do.
/run-tests |
constraints, false, | ||
Utilities::MPI:: | ||
this_mpi_process(mpi_communicator)); | ||
DoFTools::make_flux_sparsity_pattern (dof_handler, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
you don't need to run both functions, just make_flux_sparsity_pattern should do both.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
oh, or do you need this because of the constraints?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Yes, for the constraints. Otherwise, it gives run-error messages.
This pull request contains two changes:
for single core, call the make_flux_sparsity_pattern for DG with face coupling, which has less nonzero numbers; for multiple cores, call the make_flux_sparsity_pattern for DG with face coupling at every where, but it allows to use multiple processor.
move the DG BP limiter to the solver.cc inside the function call of solve_advection.
@egpuckett