That's very bad news for me. This workaround will work for Laplace problem
but not for my case. I'm dealing with Stokes problem and my method requires
that second block of right hand side is zero, otherwise MinRes crashes and
GMRES have to be used (that is far less effective).
Can You give me
I think that the problem can be avoided by splitting operator A into two
parts: dirichlet nodes and other. The structure will be:
A1 0
0I
and the operator A1 applies matrix-vector product using fixed values on
Dirichlet nodes.
If I got it right, the problem is with this part of code:
Line:
Dear Michal,
You mean because once you apply an inhomogeneous Dirichlet condition on
the velocity you also get a contribution because the divergence is not
zero? You could work around that by applying a function in the whole
domain that satisfies the Dirichlet constraints and is divergence fre
Hello Prof. Bangerth,
I have attached a test case.
Best,
Yiyang
On Friday, October 27, 2017 at 5:12:49 PM UTC-5, Wolfgang Bangerth wrote:
>
> On 10/27/2017 02:13 PM, Yiyang Zhang wrote:
> >
> > Yes I think I am setting them in a way that is exactly same throughout
> > all processes.
> > Sinc
Michael,
are actually compiling with a MPI compiler?
Does replacing NULL by nullptr help?
Best,
Daniel
Am Samstag, 28. Oktober 2017 01:52:24 UTC+2 schrieb Michael:
>
> Thanks for your reply. I did install mpich_3.2-7_amd64.deb before. When I
> tried your simple mpi example, I can not even compi
W dniu sobota, 28 października 2017 18:52:57 UTC+2 użytkownik Martin
Kronbichler napisał:
>
> Dear Michal,
>
> You mean because once you apply an inhomogeneous Dirichlet condition on
> the velocity you also get a contribution because the divergence is not
> zero?
Exactly.
> You could work