RE: [deal.II] Computing the solution gradient at the quadrature point on a face

Hi Jean-Pau, You’re always helpful. The deformation dependent loads are discussed in Wriggers, Wood and some other books and papers. I’m studying that. It seems the formulation is complicated. I believe that the auto-differentiation tools is a good approach to this nonlinear problem and I’ll dig

[deal.II] Re: Boost version problem

Update: the problem is that, when I installed the newest version of dealii, the boost version that I had separately installed was 1.67.0. Instead of using the version bundled in dealii (version 1.70.0), dealii configured with the separately installed version. However, whenever it does the

[deal.II] Boost version problem

Hi folks, I'm getting the error "The version number of boost that you are compiling with does not match the version number of boost found during deal.II's configuration step." when I try to compile the first tutorial program. It is thrown at /usr/local/include/deal.II/base/config.h:508:17

[deal.II] finite element with shape functions having delta_ij property at quadrature points

Dear all, I am solving a problem in 2d using FE_Q(2) elements and a gauss quadrature rule with (fe.degree +1) quadrature points in each co-ordinate direction, that is, I have in total nine quadrature points. My question pertains to the following: At each cell, I need to approximate a field

Re: [deal.II] point_value, Real parts, step-29

Another thing: You probably also need to initialise with the right number of components. So something like Vector vecSol (fe.n_components()); > On 10. Aug 2021, at 19:40, Jean-Paul Pelteret wrote: > > Hi Hermes, > > You don’t say what errors you’re seeing, but my guess is that it now doesn’t

Re: [deal.II] point_value, Real parts, step-29

Hi Hermes, You don’t say what errors you’re seeing, but my guess is that it now doesn’t compile. This variant of the function (the one that Daniel linked to) returns void, so you should call it before outputting the result: Vector vecSol; VectorTools::point_value(dof_handler,

[deal.II] deal.II Newsletter #177

Hello everyone! This is deal.II newsletter #177. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #12639: Replace AssertDimension (proposed by gfcas)

Re: [deal.II] point_value, Real parts, step-29

Thank you for your answer. I am not sure if I fully understand your suggestion. Do you mean something like that: Vector<*double*> vecSol; std::cout << "Solution at (0.2,0.2): "<< VectorTools::point_value(dof_handler, solution,Point<2>(0.2, 0.2),vecSol)<< std::endl; I still get some errors.

Re: [deal.II] point_value, Real parts, step-29

Hermes, Use another overload. The one returning the solution as a parameter should work: https://www.dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#acd358e9b110ccbf4a7f76796d206b9c7 Best, Daniel Am Di., 10. Aug. 2021 um 09:41 Uhr schrieb Hermes Sampedro < hermesampe...@gmail.com>:

[deal.II] point_value, Real parts, step-29

Dear all, It is explained in Step-3 how to evaluate the solution in a point. I am trying to do the same for Step-29, to evaluate the real and imaginary parts separately in a single point: *std::cout << "Solution at (0.2,0.2): "<< VectorTools::point_value(dof_handler, solution,Point<2>(0.2,

Re: [deal.II] Inaccurate convergence rate using ParsedConvergenceTable class

If you want to use symbolic calculations, you could also leverage SymEngine, and use Functions::SymbolicFunction https://www.dealii.org/current/doxygen/deal.II/classFunctions_1_1SymbolicFunction.html In your code, you could create such object by making a std::unique_ptr> fun; parse its