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
MGTwoLevelTransfer: enable FE_Nothing #12891
Conversation
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.
Some random comments, looks like a useful addition! 👍
for (unsigned int i = 0; i < fe_collection.size(); ++i) | ||
result &= fu(fe_collection[i]); |
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.
Wouldn't it be easier or at least more readable with
for (unsigned int i = 0; i < fe_collection.size(); ++i) | |
result &= fu(fe_collection[i]); | |
for (unsigned int i = 0; i < fe_collection.size(); ++i) | |
if (!fu(fe_collection[i])) | |
return false; |
and return true
otherwise?
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.
Done.
Assert(reference_cell == | ||
dof_handler_coarse.get_fe(fe_index_pair.first.second) | ||
dof_handler_coarse.get_fe(fe_index_pair.first.first) | ||
.reference_cell(), | ||
ExcNotImplemented()); |
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.
So this was a bug before, but we never actually triggered the alternative element?
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. It was "just" an assert ;)
ddc346e
to
237310b
Compare
@kronbichler I have updated the PR. Ready from my side. |
/rebuild |
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.
Looks good to me, apart from some comments regarding the test.
|
||
|
||
/** | ||
* TODO |
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.
Can you describe what the test is doing?
template <int dim> | ||
class AnalyticalSolution : public Function<dim> | ||
{ | ||
public: | ||
AnalyticalSolution() | ||
{} | ||
|
||
virtual double | ||
value(const Point<dim> &p, const unsigned int component = 0) const override | ||
{ | ||
(void)component; | ||
|
||
return p[0]; | ||
} | ||
|
||
private: | ||
}; |
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.
How about using
dealii/include/deal.II/base/function_lib.h
Lines 1329 to 1330 in 21e7779
Monomial(const Tensor<1, dim, Number> &exponents, | |
const unsigned int n_components = 1); |
exponents = {1, 0, 0}
?
237310b
to
392a312
Compare
No description provided.