Re: [deal.II] Re: fe_enriched and step-47

2017-04-09 Thread Simon Sticko
Hi. Regarding this: > On Thursday, April 6, 2017 at 1:51:41 PM UTC+2, Denis Davydov wrote: > > The issue with integrating and quadrature would still remain. I've been implementing an algorithm to generate bulk and surface quadrature given a level set function as boundary description. It's

Re: [deal.II] Re: Bug Report in GridTools::get_active_neighbors()

2019-11-29 Thread Simon Sticko
Hi, Fixed in #9108, already merged in. I should probably have written that here for completeness. Sorry. Best, Simon -- 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

Re: [deal.II] Re: FE_Q

2019-11-30 Thread Simon Sticko
Hi, Just to elaborate on what Daniel already said, Oliver, note that it is the constructor of the Step3-class that calls the constructor of FE_Q, that is Step3::Step3() : fe(1) , dof_handler(triangulation) {} If you miss this call in the Step3-constructor, you will get the error you

[deal.II] Re: Bug Report in GridTools::get_active_neighbors()

2019-11-28 Thread Simon Sticko
Hi, I think it's just the documentation that can be slightly improved here. The function get_active_child_cells has a similar signature but it has this note in its documentation: * @note Since in C++ the MeshType template argument can not be deduced from * a function call, you will have to

[deal.II] Re: Resolving with different RHS

2020-06-02 Thread Simon Sticko
Hi, a better option than computing the inverse is to factorize the matrix. This can be done using the SparseDirectUMFPACK solver: https://www.dealii.org/current/doxygen/deal.II/classSparseDirectUMFPACK.html You might want to take a look at step 22, which uses this solver:

[deal.II] Re: Error calculating shape gradient

2020-06-10 Thread Simon Sticko
Hi, from the error message you can see that the element you are using only has 1 component. You get an error because you are trying to access component 1, which doesn't exist. Since your element should have dim components, there is likely something wrong with how you create your element. It

[deal.II] Re: Error calculating shape gradient

2020-06-10 Thread Simon Sticko
{ > fe_collection.push_back(FE_Q(deg)); > quadrature_collection.push_back(QGauss(deg+1)); > face_quadrature_collection.push_back(QGauss(deg+1)); > } > > I tried to look into the fe_collection.n_components(), it yields 1. > is that what you meant? > > BR, > Ihsan > &

[deal.II] Re: How to get normal / tangential vectors at nodes, not at quadrature points?

2020-06-05 Thread Simon Sticko
Hi, I suspect that the condition you have would be easier to enforce weakly by adding terms to the weak formulation, as is done in e.g. discontinuous Galerkin. Nevertheless, if you still want the normals at the nodes you can create a "dummy" quadrature which has the nodes on the face as

Re: [deal.II] How to transform points on faces from unit cells to real cells?

2020-06-22 Thread Simon Sticko
gt;> Doug >> >> On Monday, June 22, 2020 at 5:05:57 PM UTC-4, Simon Sticko wrote: >>> >>> Hi, >>> >>> Lex, given your previous question: >>> >>> https://groups.google.com/forum/#!topic/dealii/xghVE7Km78o >>> >>> If

Re: [deal.II] How to transform points on faces from unit cells to real cells?

2020-06-22 Thread Simon Sticko
Hi, Lex, given your previous question: https://groups.google.com/forum/#!topic/dealii/xghVE7Km78o If you are already using an FEFaceValues object with a dummy-quadrature (created from FiniteElement::get_unit_face_support_points()) then you can use one of the following functions on FEFaceValues

[deal.II] Re: Non Homogenous Boundary Conditions on a Cylinder

2020-06-26 Thread Simon Sticko
Hi, I'm not sure I understand your question correctly, but if you want to impose an inhomogeneous Neumann boundary condition you would typically write a small class representing the boundary condition you want: class NeumannBoundaryCondition : public Function<3> { public: double value(const

[deal.II] Re: Non Homogenous Boundary Conditions on a Cylinder

2020-06-30 Thread Simon Sticko
int)); >} > } > > } > } > > > This seems correct to me, but the code is showing a cyldiner that looks > like this: > > > [image: Screen Shot 2020-06-29 at 5.50.19 PM.png] > which doesn't appear to be the right solution. Even when I change

[deal.II] Re: Non Homogenous Boundary Conditions on a Cylinder

2020-07-01 Thread Simon Sticko
0; j < dofs_per_cell; ++j) > system_matrix.add(local_dof_indices[i], > local_dof_indices[j], > cell_matrix(i, j)); > for (unsigned int i = 0; i < dofs_per_cell; ++i) > system_rhs(local_dof_indices[i]