Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-18 Thread jane . lee
This did exactly what I needed. Thanks!! On Saturday, March 16, 2019 at 5:48:14 PM UTC, Jean-Paul Pelteret wrote: > > There’s a function to compute the bounding box in GridTools: > > https://www.dealii.org/current/doxygen/deal.II/namespaceGridTools.html#ae1ec55abefa31cf001fd29d8d4d993f1 > > > On 1

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-16 Thread Jean-Paul Pelteret
There’s a function to compute the bounding box in GridTools: https://www.dealii.org/current/doxygen/deal.II/namespaceGridTools.html#ae1ec55abefa31cf001fd29d8d4d993f1 > On 16 Mar 2019, at 1

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-16 Thread jane . lee
Thanks for your suggestions. BoundingBox should for now be sufficient, but I am having some trouble using it properly. I am unsure what to put in to use it? triangulation? cell? I think i need something like BoundingBox mybox; mybox.get_boundary_points(); but I'm a little unsure where I put

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-15 Thread Wolfgang Bangerth
> Yes, what you put into much better words than mine is exactly what I am > needing - For a given quadrature point at > (x,y), find how far the domain extends above (x,y) in y-direction? > > So I am looking to find the y-coordinate of the point which is directly > above the (x,y) in question,

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-15 Thread jane . lee
Hi Jean Paul, Thanks for this - for now this should be sufficient, though for the rest I will have to do a long blurb. Do you have an example of how to use BoundingBox? I'm a little confused by the document. I don't know what i have to set for the argument before i do something (i am guessing

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-15 Thread Jean-Paul Pelteret
Dear Jane, There is the BoundingBox class that could tell you the overall extents of the domain. So if your geometry is rectangular and not rotated then you could work out the difference between the y-coordinate of the top

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-15 Thread jane . lee
Hi Wolfgang, Yes, what you put into much better words than mine is exactly what I am needing - For a given quadrature point at (x,y), find how far the domain extends above (x,y) in y-direction? So I am looking to find the y-coordinate of the point which is directly above the (x,y) in question,

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-14 Thread Wolfgang Bangerth
On 3/14/19 5:50 PM, jane@jandj-ltd.com wrote: > > I've tried looking a cell->face(face_no)->at_boundary() type input within the > cell iterator and by using something like cell->vertex(v)(1) but I can't > quite > get my head around how to find the y coordinate of the top boundary for every

[deal.II] Finding the y-coordinate of the boundary during the assembly

2019-03-14 Thread jane . lee
Hi all, I am solving a problem where one of the functions in the assembly process requires the knowledge of the boundary coordinate of the mesh, which moves at each timestep. My solution should be only varying in the y-direction (it is indeed with the current output) and I have a known linear