Hi,
I have written this code for testing purpose.
In case of a rotated cube, the shape functions seemed to be good as they
were parallel to the edge, although there was some small round-off errors
at the value 0.
However, when I tested with the grid which was loaded from GMSH, the shape
functions f
For a weak solution, you can split the boundary integral (coming from
integration by parts) into a normal and a tangential part. With
no-normal flux contraints you remove the normal part only, so the
tangential part needs to be written as a boundary integral on the RHS
of your PDE unless it is zero
Hi Wolfgang
Thanks for your reply.
That was actually what I had done previously. I've tried all sorts by
looking at the values on different boundaries, indicators, tried multiple
problems, etc before I had asked the question, which is why I came to
thinking it had something with the implemen
Hi Wolfgang,
I can get a table if it would be useful.
I see what you mean in terms of convergence. I guess I was looking for the
accuracy pointwise on a boundary where the Dirichlet condition for the
pressure is imposed weakly. In my case, the value of the output on the
boundary was important
Dear Wolfgang,
> It is hard to imagine situations in which the mass matrix would be
singular.
> It is a positive definite form that gives rise to the mass matrix and so
it
> really shouldn't be singular at all. Can you show the code again with
which
> you build it?
It seems that my mesh is neithe
> The mesh has over 10 cells in it, it is super refined. And oddly, when
> the refinement level is less, it doesn't blow up. It's only after a
> certain point.
> it is an even global refinement, starting from a hyper divided rectangle.
> no fancy refinement.
>
> I onl
On 4/14/19 11:59 PM, Robert Spartus wrote:
>
> Thanks for the insightful discussion on the integrating issue. Wolfgang, I
> guess your last argument is the same as you gave in one of your fantastic
> lectures?
Yes. (Also, thanks for the compliment :-) )
> Incidentally, do you have any ideas o
Jane,
> I continued to find out why I wasn't getting the correct applied Dirichlet
> values on the boundary for a code very similar to step-20, where the
> Dirichlet
> condition is applied weakly using
>
> for (unsigned int face_no=0;
> face_no::faces_per_cell;
> ++face_no)
> if (cell->at_bo
On 4/11/19 9:17 PM, Phạm Ngọc Kiên wrote:
> Testing for an edge whose global vertices located from (0,0,0) to (0,0,1) in
> real coordinates.
> With a cube I get the shape function vectors at the dof related to the edge,
> for examples, (0,0,0), (0,0,-0.25), (0,0,-0.5), (0,0,-1), which are paralle
On 4/14/19 4:29 AM, illi wrote:
> I have the following code snippet for computing Eigenvalues using Power
> Method:
> |
> Vector x;
> x = solution;
> double v = 0.0;
> PrimitiveVectorMemory> mem;
> const EigenPower>::AdditionalData data(0.);
> EigenPower<> ep(solver_control, mem, data);
> ep.solve
10 matches
Mail list logo