On Fri, Jul 24, 2020 at 9:58 PM Wolfgang Bangerth
wrote:
> On 7/23/20 12:07 PM, Xuefeng Li wrote:
> >
> > Well, the above function calculates the gradients of a finite element at
> the
> > quadrature points of a cell, not at the nodal points of a cell.
> > Such a need arises in the following situ
On 7/23/20 12:07 PM, Xuefeng Li wrote:
Well, the above function calculates the gradients of a finite element at the
quadrature points of a cell, not at the nodal points of a cell.
Such a need arises in the following situation.
for ( x in vector_of_nodal_points )
v(x) = g(x, u(x), grad u(x)
On Thu, Jul 23, 2020 at 5:14 PM Daniel Arndt wrote:
>
> You can do similarly,
>
> Quadrature q(fe.get_unit_support_points());
> FEValues fe_values (..., q, update_q_points);
> for (const auto& cell)
> ...
> points = fe_values.get_quadrature_points();
> fe_values.get_function_values(values);
>
> We need to update function v based on solution u in the following loop.
>
> for ( x in vector_of_nodal_points )
> v(x) = f(x, u(x))
>
> where v(x) and u(x) are the nodal point value for functions v and u,
> respectively, and f() is some function depending on function u as well as
> the locati
On Thu, Jul 23, 2020 at 12:43 PM Daniel Arndt
wrote:
> However, I have two more related questions. BTW, I am a newbie in C++
>> programming. So my questions may seem absurd.
>>
>>1. Now that we have a vector holding all nodal point coordinates, and
>>another vector holding all nodal point
>
> However, I have two more related questions. BTW, I am a newbie in C++
> programming. So my questions may seem absurd.
>
>1. Now that we have a vector holding all nodal point coordinates, and
>another vector holding all nodal point values of the solution. How do we
>access every noda
On Sun, Jul 19, 2020 at 7:36 PM Wolfgang Bangerth
wrote:
> On 7/19/20 6:28 PM, Daniel Arndt wrote:
> >
> > The tutorial examples show only how to access values of the solution
> at
> > the quadrature points within each cell.
> >
> >
> https://github.com/dealii/dealii/wiki/Frequently-Asked
On 7/19/20 6:28 PM, Daniel Arndt wrote:
The tutorial examples show only how to access values of the solution at
the quadrature points within each cell.
https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-to-get-the-mapped-position-of-support-points-of-my-element
Is it possible to access nodal point values of a solution in deal.ii? And
> is it possible to revise nodal values of a solution pointwise?
>
Yes, of course, the elements of your solution vector correspond to degrees
of freedoms which are nodal values in your case.
You simply access them (both for
Hi, there!
I have some general questions about deal.ii, and here is the background
info concerning my question.
There are two functions u and v, defined over domain \Omega in 1D/2D/3D. We
use deal.ii to solve for the numerical approximation of function u, using
1st degree polynomial for finite
10 matches
Mail list logo