[deal.II] Re: Compute jacobian at point (no quadrature point)

Felix, Would you take a look at if [1] is what you want to have? [1]: https://www.dealii.org/8.4.0/doxygen/deal.II/classFEValuesBase.html#a665ab681011424cb1b0ba11c4155a121 Best, Weixiong 在 2017年6月13日星期二 UTC-7上午9:35:09，Felix Lorenz写道： > > > Hi everyone, > > I want to compute the transformation

[deal.II] Re: get_function_value error with hp class

Thanks, Daniel, I used FESystem to wrap the scalar FiniteElement and the FE_Nothing together. Here are the hp-relevant variables: FESystem fe_continuous; FESystem fe_discontinuous; hp::FECollection fe_c; hp::QCollection q_c; hp::QCollection q_face_c; This is how I do initialization

[deal.II] Re: Any easy how-to: petsc and solver preconditioner selection from input?

Daniel, Thanks for the reply. I guess I'd prefer to dig a bit into SolverSelector so everything can be wrapped into the input file. This could be a fun fix to our code. Best, Weixiong 在 2017年6月13日星期二 UTC-7上午3:11:02，Daniel Arndt写道： > > Weixiong, > > Am Samstag, 10. Juni 2017 12:54:50 UTC+2

[deal.II] Compute jacobian at point (no quadrature point)

Hi everyone, I want to compute the transformation of a gradient at a point (no quadrature point) on the reference cell to the real space. Mathematically formulated I want to compute:

[deal.II] Compute jacobian at a point (no quadrature point)

Hi everyone, I want to compute the transformation of the gradient of the shape function at a point (no quadrature point) on the reference cell to the real space. Mathematically written I want to compute: : point in the reference cell : jacobian : gradient of the i-th shape

[deal.II] Re: Any easy how-to: petsc and solver preconditioner selection from input?

Weixiong, Am Samstag, 10. Juni 2017 12:54:50 UTC+2 schrieb Weixiong Zheng: > > Hi all, > > I am developing something based on petsc solvers and preconditioners. We > actually expect to invoke user-requested solvers and preconditioners from > the input file. > > I saw selector classes for

[deal.II] Re: get_function_value error with hp class

Weixiong, [...] > Then I got confused: > 1. I used hp class a year ago or so, fv.get_function_values works > directly. Now I have to use the above change. Did I do something wrong? > 2. Then why fv.shape_value still works as it did? Isn't fv still a > FESystem object? > How did you initialize

[deal.II] Re: About applying periodic boundary condition on the same mesh by different methods

Hello Daniel, I think i've found where the problem was. It's because I used set_boundary_id with active_cell_iterator after global_refinement. But I found in source codes of collect_periodic_faces function it's just cell_iterator that has been used. I've instead called set_boundary_id at