Re: [deal.II] Subtle issues in time-dependent Dirichlet boundary conditions

2019-05-16 Thread llf m
Dear Jie, In your conclusion here the point 2 may be ambiguous, > Apply nonzero_constraints at the first Newton iteration in every time step. Is that means: use nonzero_constraints to distribute newton update or solution or both? In my thought, in moving nonzero constraints case: In

Re: [deal.II] Separate components from the solution vector

2019-05-16 Thread llf m
Thanks Wolfgang, I get the point! Best M. Wolfgang Bangerth 于2019年5月17日周五 上午1:21写道: > On 5/13/19 8:45 PM, llf m wrote: > > I want to output the strain tensor and I use the > *DataPostprocessorTensor > > *which described in > > >

Re: [deal.II] Applying boundary masks when using project_boundary_values()

2019-05-16 Thread 'Maxi Miller' via deal.II User Group
If I understand that correctly, the correct order for doing that would be: -> Interpolate boundary condition without any component mask, giving me a vector -> Loop over all elements in the vector, and cross-check with the mask (using extract_dofs()) if they should be interpolated in that way ->

Re: [deal.II] Interpolate correct solution onto mesh using interpolate() with FE_Bernstein elements

2019-05-16 Thread 'Maxi Miller' via deal.II User Group
I will take a look, thanks for the help! Would it be useful to have it in a form which could be introduced in the library later (as drop-in replacement for interpolate(), for example)? Am Donnerstag, 16. Mai 2019 22:00:23 UTC+2 schrieb luca.heltai: > > > > On 16 May 2019, at 9:16, 'Maxi Miller'

Re: [deal.II] Interpolate correct solution onto mesh using interpolate() with FE_Bernstein elements

2019-05-16 Thread luca.heltai
> On 16 May 2019, at 9:16, 'Maxi Miller' via deal.II User Group > wrote: > > I would like to use the interpolate()-function to interpolate the correct > function onto my mesh for comparison during some simple tests. This works > fine when using FE_Q-elements, but fails when using

Re: [deal.II] Separate components from the solution vector

2019-05-16 Thread Wolfgang Bangerth
On 5/13/19 8:45 PM, llf m wrote: > I want to output the strain tensor and I use the *DataPostprocessorTensor > *which described in > https://dealii.org/developer/doxygen/deal.II/classDataPostprocessorTensor.html >  . > But as I am dealing with a vector-valued problem, > how can I pass the

Re: [deal.II] Interpolate correct solution onto mesh using interpolate() with FE_Bernstein elements

2019-05-16 Thread Wolfgang Bangerth
On 5/16/19 1:16 AM, 'Maxi Miller' via deal.II User Group wrote: > I would like to use the interpolate()-function to interpolate the correct > function onto my mesh for comparison during some simple tests. This works > fine > when using FE_Q-elements, but fails when using FE_Bernstein elements.

Re: [deal.II] Applying boundary masks when using project_boundary_values()

2019-05-16 Thread Wolfgang Bangerth
On 5/15/19 8:37 AM, 'Maxi Miller' via deal.II User Group wrote: > I have to replace VectorTools::interpolate_boundary_values with > VectorTools::project_boundary_values at several places in my project. I could > apply the boundary values only for certain parts in my FESystem using > component

[deal.II] Interpolate correct solution onto mesh using interpolate() with FE_Bernstein elements

2019-05-16 Thread 'Maxi Miller' via deal.II User Group
I would like to use the interpolate()-function to interpolate the correct function onto my mesh for comparison during some simple tests. This works fine when using FE_Q-elements, but fails when using FE_Bernstein elements. The error message I get is An error occurred in line <3098> of file