Re: [deal.II] boundary matrix

2021-03-10 Thread Nikki Holtzer
I am using dealii version 9.1.1. The alternative you suggested has resolved my issue. Thank you, Nikki On Tuesday, March 9, 2021 at 4:16:36 AM UTC-5 Wolfgang Bangerth wrote: > On 3/7/21 10:16 PM, Nikki Holtzer wrote: > > Thank you for your input! That was my first attempt at tackling this. >

Re: [deal.II] boundary matrix

2021-03-09 Thread Wolfgang Bangerth
On 3/7/21 10:16 PM, Nikki Holtzer wrote: Thank you for your input! That was my first attempt at tackling this. However, I run into problems with the line for(const auto & face : cell-> face_iterators()). No matter what member function I choose off of the DoFCellAccessor class template referen

Re: [deal.II] boundary matrix

2021-03-07 Thread Nikki Holtzer
Thank you for your input! That was my first attempt at tackling this. However, I run into problems with the line for(const auto & face : cell-> face_iterators()). No matter what member function I choose off of the DoFCellAccessor class template reference, I get some version of the following:

Re: [deal.II] boundary matrix

2021-03-07 Thread Wolfgang Bangerth
On 3/7/21 9:00 PM, Nikki Holtzer wrote: Yes, I have tried it. The problem I am having is with the 'face->boundary_id() == 1' piece. I get an error which states that : 'face' was not declared in this scope. The only other place that 'face' appears in my code is where I set the boundary ids like

Re: [deal.II] boundary matrix

2021-03-07 Thread peter rum
I am not completely sure what you want to accomplish, but I think the following code should work for your purpose: ``` for(const auto & cell : triangulation.active_cell_iterators()) for(const auto & face : cell-> face_iterators()) if(face->at_boundary() && face->boundary_id() == 1) //

Re: [deal.II] boundary matrix

2021-03-07 Thread Nikki Holtzer
Yes, I have tried it. The problem I am having is with the 'face->boundary_id() == 1' piece. I get an error which states that : 'face' was not declared in this scope. The only other place that 'face' appears in my code is where I set the boundary ids like the following: triangulation.begin_acti

Re: [deal.II] boundary matrix

2021-03-07 Thread Wolfgang Bangerth
On 3/7/21 5:30 AM, Nikki Holtzer wrote: For a 1d problem, I would like to implement a "boundary matrix" for an absorbing boundary condition such as the one presented in Step 24 for only the left boundary. Is it necessary to have the 'if' statement be a conditional over faces like in Step 24 e

[deal.II] boundary matrix

2021-03-06 Thread Nikki Holtzer
Hello everyone! For a 1d problem, I would like to implement a "boundary matrix" for an absorbing boundary condition such as the one presented in Step 24 for only the left boundary. Is it necessary to have the 'if' statement be a conditional over faces like in Step 24 even though it is a 1d prob