The scheme you wrote looks fine. If a=1, then you should get q=0 or something
close to machine zero. What are you getting ? Here you are projecting the
gradient onto a discontinuous polynomial space. The blowup probably happens
when you solve your advection equation, and maybe you need to look t
My apologies for a second message. Right after posting, I realized my code
snippet is potentially wrong due to a copy and paste error. You have to
call MPI_Barrier before releasing the mutex or else rank 0 could release to
mutex early, depending on what the underlying call library call is doing.
I wanted to follow up on this since I was able to find a solution which
works for me. Hopefully it will help someone :)
Thank you very much for the help.
1) As suggested by Dr. Bangerth above, MPI_THREAD_MULTIPLE can work fine
with deal.ii if you are willing to handle the initialization of the
Lee:
I want to set 2 constraints (P = 0 and \delta P = 0) on the scalar pressure
variable on a single point (4, 0.1) at the outlet. Here, \delta P is the
Newton's iteration update value for P.
I used "constraints_hp_nonzero" to interpolate the real boundary values for my
physical problem a
Hello Deal.ii Friends,
I want to set 2 constraints (P = 0 and \delta P = 0) on the scalar pressure
variable on a single point (4, 0.1) at the outlet. Here, \delta P is the
Newton's iteration update value for P.
I used "constraints_hp_nonzero" to interpolate the real boundary values for
my phys
Hello everyone!
This is deal.II newsletter #292.
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:
#17334: Test examples/step-30: remove test (proposed by tamiko)
https:
