On 1/2/19 3:24 AM, 'Maxi Miller' via deal.II User Group wrote: > that would definitely be useful for me (else that function is not usable when > running multiple MPI threads). I can take a look at the function itself, but > am not sure if I can write a patch for it without additional help.
That's where we're happy to help! Take a look at the implementation of the function here: https://github.com/dealii/dealii/blob/master/include/deal.II/numerics/fe_field_function.templates.h#L263 The exception is thrown in line 297/298. In essence, what one would need to do is this: (i) instead of throwing an exception, just condition the code in lines 300-308 with an `if (cell->is_artificial() == false)` and just ignore all points for which the cell is artificial; (ii) add an `else` clause in case where the cell is artificial in which we fill the output `values` array with invalid values; (iii) *after* the loop, if one of the cells was artificial, throw the exception. Want to give this a try? A guide on how to contribute patches is here: https://github.com/dealii/dealii/wiki/Contributing Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
