Re: [deal.II] Why is this asserted to be impossible?
On 09/18/2016 03:45 AM, Alex Zimmerman wrote: I question all of the asserts in this file: https://github.com/dealii/dealii/blob/master/source/numerics/vector_tools_rhs.cc My particular case is the one on line 36, for dim = 1 and spacedim = 1. Thanks for looking into this. This looks like we disabled these functions back in the day when computing boundary integrals in 1d did not work. I think this ought to work now. What happens if you just remove these 4 functions from the file, recompile the deal.II library, and test whether the result works for your code? Best W. -- Wolfgang Bangerth email: bange...@colostate.edu 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
Re: [deal.II] Why is this asserted to be impossible?
I question all of the asserts in this file: https://github.com/dealii/dealii/blob/master/source/numerics/vector_tools_rhs.cc My particular case is the one on line 36, for dim = 1 and spacedim = 1. Thanks for looking into this. On Sunday, September 18, 2016 at 2:42:01 AM UTC+2, Wolfgang Bangerth wrote: > > On 09/17/2016 05:20 AM, Alex Zimmerman wrote: > > > > Not understanding why this is impossible, I began digging through the > code a > > bit and thought I would write my own 1D implementation. I noticed that > the > > code in the existing template should work correctly without any changes. > So > > instead I simply copied and renamed the template, without the > implementation > > asserting impossibility in 1D ( Assert(false, ExcImpossibleInDim(1)); ). > My > > first test ran exactly as expected. > > > > deal.II often asserts that something doesn't make sense or is impossible > in a > > certain dimension. Sometimes I agree and sometimes I don't. This is the > first > > time that it's just seemed completely wrong to me. What am I missing? > > Quite possibly nothing -- we may simply be wrong. Can you point out the > file > and line where the offending Assert is located? > > Best > W. > > -- > > Wolfgang Bangerth email: bang...@colostate.edu > > 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
Re: [deal.II] Why is this asserted to be impossible?
On 09/17/2016 05:20 AM, Alex Zimmerman wrote: Not understanding why this is impossible, I began digging through the code a bit and thought I would write my own 1D implementation. I noticed that the code in the existing template should work correctly without any changes. So instead I simply copied and renamed the template, without the implementation asserting impossibility in 1D ( Assert(false, ExcImpossibleInDim(1)); ). My first test ran exactly as expected. deal.II often asserts that something doesn't make sense or is impossible in a certain dimension. Sometimes I agree and sometimes I don't. This is the first time that it's just seemed completely wrong to me. What am I missing? Quite possibly nothing -- we may simply be wrong. Can you point out the file and line where the offending Assert is located? Best W. -- Wolfgang Bangerth email: bange...@colostate.edu 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
[deal.II] Why is this asserted to be impossible?
I have a code that is an extension of step-26 with non-homogeneous Neumann boundary conditions and zero forcing function, among other changes. Until recently I had only actually run the code in 2D, though I have been continuing the "dimension independent" programming. Yesterday I tried running a simple 1D case for the first time. I was already using VectorTools::create_boundary_rhs to apply the Neumann boundary conditions. Unfortunately I discovered that the 1D implementations of create_boundary_rhs all assert that this is an impossible operation in 1D. Not understanding why this is impossible, I began digging through the code a bit and thought I would write my own 1D implementation. I noticed that the code in the existing template should work correctly without any changes. So instead I simply copied and renamed the template, without the implementation asserting impossibility in 1D ( Assert (false, ExcImpossibleInDim(1)); ). My first test ran exactly as expected. deal.II often asserts that something doesn't make sense or is impossible in a certain dimension. Sometimes I agree and sometimes I don't. This is the first time that it's just seemed completely wrong to me. What am I missing? -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
