Hi Michał,
It would be helpful if you could sent us a small, minimal example that
shows the problem (in particular, what type is "LevelVectorType" and
"SCPreconditoinMGs", "solver").
The example should be as short as possible and doesn't need to be
functional / or do anything at all - we essential
On 02/16/2018 07:34 PM, Moraad Biagooi wrote:
As you suggested, I tried a 2D mesh. The plots shows that the sign difference
happens at the edge of two neighbor symmetric cells answers. Two gradient
functions may calculate on different cells, so it looks like the signs are
different.
So it loo
Dear Prof. Wick,
I have used isotrope formulation for miehe shear loading (without local
refinement). I ca not see two crack branching which is described in the
article (A review on phase-field models of brittle fractureand a new fast
hybrid formulation)
I attach the test results.
Thanks for
Hi All,
I'm having problem when applying periodicity on meshes generated with Gmsh
or Abaqus. The problem arises when I run in parallel (in serial it works)
getting the following error:
void
dealii::GridTools::match_periodic_face_pairs(std::set >&, std::set::type, unsigned int> >&, int,
std:
Thank you for you reply,
The issue with the boundaries was indeed solved when I assigned the
boundary ids on the ghost cells as well.
As for the other problem with the refinement around the processor borders,
I was using a relatively complex way to move the triangulation vertices
each iteratio
On 02/15/2018 07:29 AM, Moraad Biagooi wrote:
I have some problems with the signs of VectorTools::point_gradient.
Here's an example:
I modified Step-4 of tutorial so that the mesh is a
GridGenerator::hyperball with half of its boundary_id equal to one.
After solving the problem, a function
Dear Wolfgang,
Thank you for your reply.
This problem is solved by adding
constraints.distribute (solution);
in
void MixedLaplaceProblem::solve ()
in my code.
Kind regards,
Jie
On Thursday, February 15, 2018 at 5:24:44 PM UTC+1, jie liu wrote:
>
>
> Dear all,
>
> I want to solve an 1D Poisso
> Dear David,
>
> I studied the step- 44 Tutorial Program and read already a lot about how
> to apply a Neumann contribution. But I have still two questions:
> In my case, I like to apply a force (not a force per unit reference area)
> as a Neumann contribution on the RHS. But for the Integral,
Dear David,
> I studied the step- 44 Tutorial Program and read already a lot about how to
> apply a Neumann contribution. But I have still two questions:
> In my case, I like to apply a force (not a force per unit reference area) as
> a Neumann contribution on the RHS. But for the Integral, Ni *
Dear David,
> So my first question is: Can I do it like this or is there a simpler way
of doing it?
I can't really comment on this. Hopefully someone with more experience than
me will chip in.
> And the second question is: How can I access this sort of Information?
The area might for exampl
10 matches
Mail list logo
