How do I calculate the (hdiv-semi) norm only on certain components of my (vector-valued) solution?

57 views
Skip to first unread message

georgios....@gmail.com

unread,
Aug 11, 2018, 11:38:00 AM8/11/18
to deal.II User Group

Dear all, 

I would be most grateful for your help on the following matter.  I am solving a combined stokes-laplace problem (have a vector valued solution).  I am trying to calculate the hdiv norm on the first two components of my solution (i.e. the two components of the velocity), cellwise.   The implementation of my combined problem is from tutorial 46.  The integrate_difference is from tutorial 20 since I needed the mask for the velocity components.  However, I get an error saying 'no matching function for the call to integrate difference'.  Admittedly, the version of the function i am using (i.e. with the component mask as the last included argument) is not in vector_tools.h  which is where all the other versions of the integrate_differnce are I think.  I include the file for your consideration.  I would be most grateful for your input.  For succesful implementations of this version of integrate_difference please refer to tutorials 20, 55 and 56.

Kind regards,
Georgios
stokes_laplace_again.cpp

Wolfgang Bangerth

unread,
Aug 12, 2018, 5:05:25 PM8/12/18
to dea...@googlegroups.com
On 08/11/2018 09:37 AM, georgios....@gmail.com wrote:
>
> I would be most grateful for your help on the following matter.  I am solving
> a combined stokes-laplace problem (have a vector valued solution).  I am
> trying to calculate the hdiv norm on the first two components of my solution
> (i.e. the two components of the velocity), cellwise.   The implementation of
> my combined problem is from tutorial 46.  The integrate_difference is from
> tutorial 20 since I needed the mask for the velocity components.  However, I
> get an error saying 'no matching function for the call to integrate difference'.

The problem is that the functions that take an hp::DoFHandler object as first
argument want a hp::QCollection as fifth argument, but you just provide a
Quadrature object. If you pass a hp::QCollection instead, it compiles just fine.

Best
W.

--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

Georgios Sialounas

unread,
Aug 21, 2018, 12:38:34 PM8/21/18
to dea...@googlegroups.com
Dear Prof. Bangerth,

Thank you very much!  

Kind regards,
Georgios



--
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 a topic in the Google Groups "deal.II User Group" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/dealii/eaDIPSvvP48/unsubscribe.
To unsubscribe from this group and all its topics, send an email to dealii+un...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply all
Reply to author
Forward
0 new messages