Re: [deal.II] Using FEInterfaceValues with FESystems?
Timo, Don't worry about that. I can't complain if someone else does the work no? :) On Sunday, March 14, 2021 at 1:07:49 p.m. UTC-4 Timo Heister wrote: > Bruno, I already have a working implementation of the extractors on my > machine. I just didn't get the time to clean up, write tests, and > submit the PR. Sorry. I will try to pick this up again in the next few > days. > > On Fri, Mar 12, 2021 at 12:26 PM Wolfgang Bangerth > wrote: > > > > On 3/11/21 11:10 PM, blais...@gmail.com wrote: > > > Let me read the documentation a bit and I could try to come up with a > PR. > > > Maybe this is something we could discuss? I'd be glad to help in > improving DG > > > methods :) > > > > Happy to discuss in whatever venue you propose! > > > > Start by looking at the FEValuesViews::Scalar class. We'd need a > corresponding > > FEInterfaceValuesViews::Scalar class. Once you know how the former > operates, > > you will probably be able to see how the latter needs to be implemented. > I > > foresee a lot of copy-paste of small functions where the only change is > in the > > 2-3 line bodies. > > > > Best > > Wolfgang > > > > -- > > > > 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+un...@googlegroups.com. > > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/5882752f-4530-16b6-41d4-b8b6c9d746f9%40colostate.edu > . > > > > -- > Timo Heister > http://www.math.clemson.edu/~heister/ > -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/633d1b6b-868d-4292-9b3a-0ef5269ba4d1n%40googlegroups.com.
Re: [deal.II] Using FEInterfaceValues with FESystems?
Bruno, I already have a working implementation of the extractors on my machine. I just didn't get the time to clean up, write tests, and submit the PR. Sorry. I will try to pick this up again in the next few days. On Fri, Mar 12, 2021 at 12:26 PM Wolfgang Bangerth wrote: > > On 3/11/21 11:10 PM, blais...@gmail.com wrote: > > Let me read the documentation a bit and I could try to come up with a PR. > > Maybe this is something we could discuss? I'd be glad to help in improving > > DG > > methods :) > > Happy to discuss in whatever venue you propose! > > Start by looking at the FEValuesViews::Scalar class. We'd need a corresponding > FEInterfaceValuesViews::Scalar class. Once you know how the former operates, > you will probably be able to see how the latter needs to be implemented. I > foresee a lot of copy-paste of small functions where the only change is in the > 2-3 line bodies. > > Best > Wolfgang > > -- > > 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. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/5882752f-4530-16b6-41d4-b8b6c9d746f9%40colostate.edu. -- Timo Heister http://www.math.clemson.edu/~heister/ -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAMRj59FfkCKePzNiUX5gqr2LeKO8gPC-h9MPCS05SEndJC%3DBdg%40mail.gmail.com.
Re: [deal.II] Using FEInterfaceValues with FESystems?
On 3/11/21 11:10 PM, blais...@gmail.com wrote: Let me read the documentation a bit and I could try to come up with a PR. Maybe this is something we could discuss? I'd be glad to help in improving DG methods :) Happy to discuss in whatever venue you propose! Start by looking at the FEValuesViews::Scalar class. We'd need a corresponding FEInterfaceValuesViews::Scalar class. Once you know how the former operates, you will probably be able to see how the latter needs to be implemented. I foresee a lot of copy-paste of small functions where the only change is in the 2-3 line bodies. Best Wolfgang -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/5882752f-4530-16b6-41d4-b8b6c9d746f9%40colostate.edu.
Re: [deal.II] Using FEInterfaceValues with FESystems?
Good to know. Let me read the documentation a bit and I could try to come up with a PR. Maybe this is something we could discuss? I'd be glad to help in improving DG methods :) On Thursday, March 11, 2021 at 4:29:04 p.m. UTC-5 Wolfgang Bangerth wrote: > On 3/11/21 10:14 PM, blais...@gmail.com wrote: > > I am currently trying to solve an advection-diffusion equation at very, > very > > high Péclet number. I first implemented a GLS-stabilized continuous > Galerkin > > solution, it works well, but clearly I think a DG approach would be > immensely > > better here. > > > > The issue is that the velocity field for my advection-diffusion stems > from a > > Continuous Galerkin solution with another dof handler, formulated using > a > > FESystem. Everything work well until them. The issue I have is that the > > FEInterfaceValues does not seem to want to work with FEValuesExtractors? > It > > seems like it cannot use the [] operator to extract my velocity values > out of > > my 4 components (velocity + pressure). > > > > I am sure i'm forgetting something stupid here, but how should I proceed > for > > vector valued problem using DG? > > No, you're not forgetting anything. FEInterfaceValues is still a pretty > new > class and nobody has ever written the equivalent of the FEValuesViews > classes > that you get from FEValues when you apply an extractor. > > I suspect that the effort would not be huge, and that one could just > copy-paste a lot of code from the existing views classes. If you were > interested in doing that, I'd be happy to discuss what I think is going to > be > necessary. But the short answer is really that FEInterfaceValues isn't > ready > for vector-valued problems at the moment. > > Best > Wolfgang > > -- > > 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/1189c3ca-63e1-4ebf-aace-395747dd4425n%40googlegroups.com.
Re: [deal.II] Using FEInterfaceValues with FESystems?
On 3/11/21 10:14 PM, blais...@gmail.com wrote: I am currently trying to solve an advection-diffusion equation at very, very high Péclet number. I first implemented a GLS-stabilized continuous Galerkin solution, it works well, but clearly I think a DG approach would be immensely better here. The issue is that the velocity field for my advection-diffusion stems from a Continuous Galerkin solution with another dof handler, formulated using a FESystem. Everything work well until them. The issue I have is that the FEInterfaceValues does not seem to want to work with FEValuesExtractors? It seems like it cannot use the [] operator to extract my velocity values out of my 4 components (velocity + pressure). I am sure i'm forgetting something stupid here, but how should I proceed for vector valued problem using DG? No, you're not forgetting anything. FEInterfaceValues is still a pretty new class and nobody has ever written the equivalent of the FEValuesViews classes that you get from FEValues when you apply an extractor. I suspect that the effort would not be huge, and that one could just copy-paste a lot of code from the existing views classes. If you were interested in doing that, I'd be happy to discuss what I think is going to be necessary. But the short answer is really that FEInterfaceValues isn't ready for vector-valued problems at the moment. Best Wolfgang -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/95dbdc73-3b8b-6402-98b9-4894446aa55b%40colostate.edu.
[deal.II] Using FEInterfaceValues with FESystems?
Dear all, I hope you are well. I am currently trying to solve an advection-diffusion equation at very, very high Péclet number. I first implemented a GLS-stabilized continuous Galerkin solution, it works well, but clearly I think a DG approach would be immensely better here. The issue is that the velocity field for my advection-diffusion stems from a Continuous Galerkin solution with another dof handler, formulated using a FESystem. Everything work well until them. The issue I have is that the FEInterfaceValues does not seem to want to work with FEValuesExtractors? It seems like it cannot use the [] operator to extract my velocity values out of my 4 components (velocity + pressure). I am sure i'm forgetting something stupid here, but how should I proceed for vector valued problem using DG? Best Bruno -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/26bd805e-4467-43f1-9132-1b61e1880964n%40googlegroups.com.