Re: [deal.II] Re: Function transfer from coarse to fine grid

2019-04-07 Thread Ina
Thank you for your reply. Then I will just use VectorTools::interpolate_to_
different_mesh. 

Best,
Ina

-- 
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] project_boundary_values for components

2019-04-07 Thread Wolfgang Bangerth
On 4/6/19 12:29 PM, Anton wrote:
> 
> Basically if I interpolate a step function (which can be exactly presented in 
> for example DGQ1 space if the discontinuity is at the node) I get a 1-element 
> wide transition.  In ASCII pictures, I expect  ___  but I get ___/---  
> instead :)

Ah, yes -- ASCII art to the fore! Yes, the problem is that the DGQ(k) elements 
have support (interpolation) points at the vertices, and so you'll get the 
same value for both cells at these vertices. That's unavoidable. Projection is 
the solution.


> So if I understand you correctly, one would need to (1) obtain 
> ConstraintMatrix from project_boundary_values; (2) go over each row of the 
> constraint matrix and check which component it is; (3) if it is the "right" 
> component  then copy the row.  Could I ask how does one figure out the 
> component number from the row number in the ConstraintMatrix - I cannot 
> immediately figure this out from DoFHandler interface?

You can't do it this way, but there are functions in DoFTools that allow you 
to extract *all* DoF indices for a particular component from a DoFHandler, and 
then you can cross-check the two objects.

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] Poisson test code fails for a fe-degree >= 2

2019-04-07 Thread Wolfgang Bangerth
On 4/7/19 1:54 PM, 'Maxi Miller' via deal.II User Group wrote:
> Could do that, but not until Easter (Deadlines before that), if that is ok?

Any time is good for contributions!
Cheers
  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] Poisson test code fails for a fe-degree >= 2

2019-04-07 Thread 'Maxi Miller' via deal.II User Group
Could do that, but not until Easter (Deadlines before that), if that is ok?

Am Freitag, 5. April 2019 17:39:30 UTC+2 schrieb Wolfgang Bangerth:
>
> On 4/5/19 9:34 AM, 'Maxi Miller' via deal.II User Group wrote: 
> > Oh, thanks! That happens if I start from an example, and don't write 
> > everything from scratch... 
>
> Since I stepped in that trap myself a couple of weeks ago and had to 
> have Timo rescue me, I'm thinking that we ought to just change this in 
> all tutorial programs where the quadrature formula degree is not tied to 
> the FE degree. Would you want to write up a pull request to make such a 
> change? 
>
> 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.
For more options, visit https://groups.google.com/d/optout.