Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

On 07/26/2017 01:18 PM, sml.imf...@gmail.com wrote: I must admit that I don't really remember ... I think it was GridOut, since I specifically wanted to export only the grid and not some solution data. What dou you mean with "GridOut should work"? I was wrong thinking that the GridOut::write_v

Re: [deal.II] Component Mask for setting initial conditions

On 08/04/2017 09:49 AM, Maxi Miller wrote: Why is there only one value not equal to zero, and why is it moving? Shouldn't all the values be set already? Maxi -- this is not a question to which you can expect an answer on the mailing list. We are happy to help with conceptual questions, but y

Re: [deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

On 08/04/2017 07:07 AM, sml.imf...@gmail.com wrote: Yes, that makes sense to me. But now if we have neighboring faces with differrent refinement levels (e.g. the face element on one side is only half as long as the one on the other side) the situation gets a little more complicated, right? Sinc

Re: [deal.II] Component Mask for setting initial conditions

I am filling it up like that, and then (when filling system_matrix) I would like to take a look at the values: for (unsigned int q_point=0; q_point class InitialValues : public Function { private: equation_class eq_class; physics_equations equations; public:

Re: [deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

Yes, that makes sense to me. But now if we have neighboring faces with differrent refinement levels (e.g. the face element on one side is only half as long as the one on the other side) the situation gets a little more complicated, right? Since then you cannot just copy the coefficients 1 to 1,

Re: [deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

As indicated in the documentation , the finite elements FE_DGP and FE_FaceP are not based on nodal interpolation (as eg. FE_DGQ and FE_FaceQ), but on projection. Apparently, this prevents these elements from being used in combinat

[deal.II] FE_DGP/FE_FaceP and hp::DoFHandler

Dear all, As indicated in the documentation , the finite elements FE_DGP and FE_FaceP are not based on nodal interpolation (as eg. FE_DGQ and FE_FaceQ), but on projection. Apparently, this prevents these elements from being used in

Re: [deal.II] Component Mask for setting initial conditions

On 08/04/2017 05:41 AM, Maxi Miller wrote: Is an alternative to check the components, and set the conditions accordingly? I wrote it like | template doubleInitialValues::value(constPoint&p,constunsignedintcomponent)const { (void)p; Assert(component n_components,ExcIndexRange(component,0,this->n

[deal.II] Re: hp HDG in deal.II

Daniel, Now I have another question. I am trying to compute the L2 error over all >> the faces. For simplicity, I tested this starting from the tutorial >> step-51. I added another VectorTools::integrate_difference statement to the >> postprocess routine, with arguments dof_handler and solution

Re: [deal.II] Component Mask for setting initial conditions

Is an alternative to check the components, and set the conditions accordingly? I wrote it like template double InitialValues::value(const Point &p, const unsigned int component) const { (void) p; Assert(component < this->n_components, ExcIndexRange(component, 0, thi

Re: [deal.II] Component Mask for setting initial conditions

On 08/04/2017 04:23 AM, 'Maxi Miller' via deal.II User Group wrote: When creating boundary conditions, I can use the parameter component_mask in order to set boundary conditions only for a part of my system. Now I want to set initial conditions for the same system. One approach is using | Vecto

Re: [deal.II] Function for VectorTools::interpolate_boundary_values with different boundary values in multi-variable system

On 08/04/2017 04:19 AM, 'Maxi Miller' via deal.II User Group wrote: I changed the BoundaryCondition-class to | template classBoundaryValues:publicFunction { private: equation_class eq_class =carrier_density; physics_equations equations; public: BoundaryValues(equation_class equ

[deal.II] Component Mask for setting initial conditions

When creating boundary conditions, I can use the parameter component_mask in order to set boundary conditions only for a part of my system. Now I want to set initial conditions for the same system. One approach is using VectorTools::project but there I can not see an option for component_mask in

Re: [deal.II] Function for VectorTools::interpolate_boundary_values with different boundary values in multi-variable system

I changed the BoundaryCondition-class to template class BoundaryValues : public Function { private: equation_class eq_class = carrier_density; physics_equations equations; public: BoundaryValues(equation_class equ_class) : Function(3*dim), eq_class(equ

Re: [deal.II] Function for VectorTools::interpolate_boundary_values with different boundary values in multi-variable system

Dear Maxi, It would have been helpful if you'd posted the full error. However, I'm guessing that this is the problem: The function that you're trying to call is void VectorTools::interpol

Re: [deal.II] Function for VectorTools::interpolate_boundary_values with different boundary values in multi-variable system

On 08/04/2017 03:22 AM, Maxi Miller wrote: Then I get the following error: | error:nomatching functionforcall to ‘interpolate_boundary_values(dealii::DoFHandler<2>&,int,Step52::BoundaryValues<6>,dealii::ConstraintMatrix&)’ VectorTools::interpolate_boundary_values(dof_handler,1,BoundaryValues<3*

Re: [deal.II] Function for VectorTools::interpolate_boundary_values with different boundary values in multi-variable system

Then I get the following error: error: no matching function for call to ‘interpolate_boundary_values(dealii ::DoFHandler<2>&, int, Step52::BoundaryValues<6>, dealii::ConstraintMatrix &)’ VectorTools::interpolate_boundary_values(dof_handler,1, BoundaryValues<3*dim>(carrier_density),constrai

Re: [deal.II] Function for VectorTools::interpolate_boundary_values with different boundary values in multi-variable system

On 08/04/2017 01:38 AM, 'Maxi Miller' via deal.II User Group wrote: but how do I tell that in | VectorTools::interpolate_boundary_values(dof_handler,1,BoundaryValues<3*dim>(),constraint_matrix); | What happens when you do it? The function takes an additional ComponentMask argument that allow

Re: [deal.II] Need a deal.II user to solve diffusion reaction problem.

Yusheng, I am in the process to evaluate open source FEM software with a simple diffusion-reaction problem on the 3D cylinder. I am looking for a person who can solve this problem and be familiar Deal.II enough to compare difference modules. You may get more answers if you were more specifi

[deal.II] Function for VectorTools::interpolate_boundary_values with different boundary values in multi-variable system

I have a system with three 3d-variables, and would like to apply different boundary conditions on the constraint matrix, depending on the variable. According to the examples I have to create my own function: template class BoundaryValues : public Function { private: equatio