On 8/3/23 11:25, Simon Wiesheier wrote:
The reason why I wanted to use an
try/catch block is to surrogate the use of an if-statement.
Consider these two variants:
Variant 1:
disable_abort_on_exception();
try
{
feValues.reinit(cell,...);
}
catch(dealii::Exception Base exception & )
{
Dear Simon,
Thank you very much for your reply. The wrapper code looks
embarrassingly simple!
> If you are working with a perfectly Cartesian background I would expect that
> it would actually be faster to generate the quadrature in real space and then
> transform the quadrature to reference
Hi.
I have done this in the past. I dug up the wrapper function I had for it, in
case we decide that we should add it to the library:
https://github.com/dealii/dealii/pull/15838
The problem with this approach is that evaluating the reference space level set
function gets very expensive.