Re: [deal.II] Why is this asserted to be impossible?

2016-09-18 Thread Wolfgang Bangerth

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?

2016-09-18 Thread Alex Zimmerman
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?

2016-09-17 Thread Wolfgang Bangerth

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.