Dear Bruno,
Thanks a lot for the positive feedback! I’m really glad to know that those
changes worked well for you. Its always great to hear that the improvements
that we make to deal.II (be them small or large) are impactful and help our
users, and it certainly encourages us to keep at it.
On 4/3/19 8:04 AM, Muhammad Mashhood wrote:
>
> I am trying to apply the Neumann (force)
> BC on the edge using the right_hand_side() function which I modified as
> following:
>
>
> void right_hand_side (const std::vector > , // the active
> cell is also
On 4/3/19 6:53 AM, solee wrote:
> I am a beginner of deal.ii and some problem happened. I want to solve the
> Maxwell equations in the 2D cylindrical coordinates with axis-symmetry.
> When I create a weak form of the equation (see attached formula), I want to
> remove 1 / r from the first term.
In addition to the excellent idea by Juan Carlos:
> When I initialize the PETSc sparsity matrix, I don't know where the nonzeros
> are.
> Then I loop over the locally owned cells on each processor, and guadually
> know
> for each row of the matrix the locations and values of the nonzeros. Do
Dear group members,
I am trying to apply the Neumann (force)
BC on the edge using the right_hand_side() function which I modified as
following:
void right_hand_side (const std::vector > , // the
active cell is also imported when right hand side is called
Hello there
I am a beginner of deal.ii and some problem happened. I want to solve the
Maxwell equations in the 2D cylindrical coordinates with axis-symmetry. When I
create a weak form of the equation (see attached formula), I want to remove 1 /
r from the first term.
When I refer to other
Phani,
*[...]*
>
> *Please verify that both the operating system and the processor support
> Intel(R) AVX512DQ, AVX512F, AVX512CD, AVX512BW and AVX512VL instructions.*
>
> *source/numerics/CMakeFiles/obj_numerics_inst.dir/build.make:91: recipe
> for target
Hey Pai Liu,
At some point in my research I had to do something similar to what you want
to do. In the application of the so-called Dirichlet-to-Neumann maps for
Helmholtz problems, one has to assemble a dense block corresponding to the
boundary degrees of freedom, so I had to include extra