domains?). For
traction BCs it is even more complicated since that means you need to
prescribe stress tensors. This unfortunately makes it difficult to
prescribe such BCs.
The strategy that I chose works for the lowest order Raviart-Thomas element
but what about higher order (because the dofs are e
be used here? Any ideas or
suggestions where I should look?
Best regards,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google
grid classes for
that? As far as I understand there are dofhandlers on each level of the
mesh. Can I get dofhandlers for all active subcells of say coarse cell
number n in level 5?
Best regards,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
On Wednesday, August 29, 2018 at 11:50:57 AM UTC+2, Jean-Paul Pelteret
wrote:
>
> Hi Konrad,
>
> Do you think it would be possible to (ab)use the multigrid classes for
> that? As far as I understand there are dofhandlers on each level of the
> mesh. Can I get dofhandlers for a
1 and the active children
> cells on mesh 2.
>
> Best,
> J-P
This is exactly what I am doing. Works. :-)
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received
).
On Thursday, August 30, 2018 at 6:41:51 PM UTC+2, Wolfgang Bangerth wrote:
>
>
> Jane: Konrad already gave the correct answer. There is also the function
> VectorTools::project_boundary_values_div_conforming:
>
> https://www.dealii.org/9.0.0/doxygen/deal.II/group
solutions for writing output
efficiently but
is there also a good solution solution for shared memory parallelism? I
would like to write vtk-files.
Maybe I missed it while browsing the documentation. Can anyone please point
me to the right spot?
Best,
Konrad
--
The deal.II project is located
Hi Jane Lee,
I recently came across a similar problem.
On Thursday, August 30, 2018 at 12:35:47 PM UTC+2, Jane Lee wrote:
>
> I believe the Neumann conditions are strongly imposed.
>
> And yes - I realised that inhomogeneous Neumann bc is ambiguous phrasing.
>
> I mean that I have a conditions
n easy to bring on
different nodes on a cluster). Thanks again for your comments.
Cheers,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscrib
not work so well
if you do not have scale separation (at least not in their standard
version).
Let me know about questions.
Cheers,
Konrad
On Wednesday, November 14, 2018 at 3:50:27 AM UTC+1, Ronghong Fan wrote:
>
> Dear Konrad,
> Do you have some progress about Multiscale Finite Eleme
let's get in
touch by email. This is my homepage with contact info
<https://www.clisap.de/research/a:-climate-dynamics-and-variability/crg-numerical-methods-in-geosciences/team-members/konrad-simon>
.
Best,
Konrad
On Monday, November 19, 2018 at 2:33:39 PM UTC+1, Ronghong Fan wrote:
&g
straints on degrees of freedom such as hanging nodes or periodic
boundary conditions. In such cases, building the sparsity pattern will
succeed, but you will get errors such as the current one at one point or
other when trying to write into the entries of the matrix.
It seems like I do not hand
Issue resolved. I used the wrong sparsity pattern. That was just a typo in
the code. Sorry for having bothered you with this but anyway:
Many thanks,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum
ribute_local_to_global(cell_matrix,
> local_dof_indices, advection_matrix);
> } // for number of quadrature points
>
> } // for each cell
>
> } // Simulator::create_advection_matrix()
>
> Sorry if that's too many details, I hope its clear. Thanks again for your
> he
> Hi Konrad,
>
> I figured out my error. I was imposing constraints differently on the
> advection matrix than from the other matrices. I fixed it by not using the
> constraints.local_to_global() function and just computing the full matrix
> and using constraints.conden
to -\Delta u + \epsilon u and setting f to zero (compatibility
is not necessary then) I see something that is close to what I would expect
in the pure Laplace case.)
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https
sides the rotation of the shape functions that I just
mentioned.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Gro
Hi Wolfgang,
Thank you, I guess I will give it a try. Do you have a recommendation for
literature about FEs on quad meshes?
Btw, the first de Rham complex I wrote down (see above) does not seem to be
correct on quadrilaterals.
Best,
Konrad
--
The deal.II project is located at http
without having to rescale.
Thanks again :-)
Konrad
On Monday, March 18, 2019 at 11:17:40 AM UTC+1, Daniel Arndt wrote:
>
> Konrad,
>
> I implemented a class derived from the Function class that evaluates
> a scalar or vector shape function at a given (set of) point(s) in a
part this will be CFL (note this is not a
sufficient condition) which scales like v\delta t / \delta x <= const. For
the diffusive part the time step constraint is much more severe since it
scales like D\delta t / (\delta x)^2 <= const. I implemented a similar
(also periodic scheme) and I did not run
starting point?
>
Thank you, I will see if I can make it.
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
&
here?
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from
for such
problems. I am not an expert in linear solvers but their efficiency is (as
so often) crucial here.
I would be grateful for any idea...
Cheers,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
>
> Yes, there are such interfaces. Take a look at SparseDirectUMPFACK (no
> PETSc required) and
>
>
> https://dealii.org/developer/doxygen/deal.II/classPETScWrappers_1_1SparseDirectMUMPS.html
>
>
> Best
> W.
>
> Thanks! The second one I was looking
corresponding
basis object are on the same cluster node?
Is there something like a distributed STL-like container whose elements I
can initialize on different cluster nodes (for example by querying
*cell->is_locally_owned()*) when looping over the triangulation?
Best,
Konrad
--
The deal
lii.org/current/doxygen/deal.II/classBlockVector.html>
velocity(solution.block(0));
Every second entry (in 2D) should then be your x-component of the velocity. ...
at least as far as I recall (so better double check that)
Best, Konrad
Hey everyone,
>
> I have the following problem
>
://www.boost.org/doc/libs/1_48_0/doc/html/boost/mpi/communicator.html
Never tried it but according to the documentation it seems possible.
Best,
Konrad
On Friday, August 23, 2019 at 6:35:18 PM UTC+2, Wolfgang Bangerth wrote:
>
> On 8/23/19 7:33 AM, Konrad wrote:
> >
> > Thanks
entify cells globally in a
distributed triangulation. I don't see though how the objects in the
std::map are distributed among cluster nodes as the cells are. Am I missing
something here?
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://grou
if you do not
provide them explicitly (and when you must do so). I really recommend the
book by Scott Meyers "Effective C++: 55 ways ..."
Hope that helps.
Best,
Konrad
On Wednesday, September 4, 2019 at 4:44:24 AM UTC+2, yuesu jin wrote:
>
> Hi all,
>
> I have a question
exactly when an
object contains pointers to whatever distributed data somewhere in memory
etc.
At any rate, Thanks for the hints.
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You
Dear deal.ii community,
I was wondering if there is a simple way to parse input (e.g., from a
parameter file) and use it similarly to the FunctionParser class (which
works for vector valued functions) just for a TensorFunction.
Anyone did that or something similar?
Best,
Konrad
Hi Charlie,
I can compile and run it. Could there be a problem with your deal.ii setup?
I am using v9.1.1.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received
::numbers::invalid_unsigned_int);
#endif
I guess that was the bug.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google
e relevant
tutorials? Parallelization is an issue since I would like to solve very
large problems in 3D.
Thanks in advance and best regards,
Konrad :-)
Here the Schur complement code:
template
class SchurComplementMPI : public Subscriptor
{
private:
using BlockType = typename BlockMatrixType:
he destructor, i.e.,
when the object goes out of scope. Only then
2. When you implement MPI code it is important to get into the right
mindset: You do not program code for one compute node only. You are
programming code for all nodes (at the same time). In your example the for
loop gets executed on
Hi Konrad,
I may add as a comment (without going into details): Your problem is - as
you correctly describe above - a mixed Poisson problem. You essentially
have two options to solve it and the way you choose influences the way you
need to treat the boundary conditions: in primal form
essage is actually quite clear. What I
found out in addition is that it is advantageous and helpful to read the
documentation ... ;-)
Best ,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
not tell what I am doing wrong). I do not invoke any
threading deliberately.
Sorry for posting this if inappropriate - you are of course not responsible
for debugging my code - but any hint would be appreciated.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mail
Thank you, David, so I will see if I can maybe still use it in the Schur
complement somehow.
Best,
Konrad
On Friday, September 6, 2019 at 7:46:23 PM UTC+2, David Wells wrote:
>
> Hi Konrad,
>
> I don't think that it is possible to use MUMPS with a block matrix for
> exactly
Dear deal.ii community,
is it possible to use MUMPS with a PETScWrappers::MPI::BlockSparseMatrix?
Don't find anything but I see that PETScWrappers::MPI::BlockSparseMatrix
does not inherit from PETScWrappers::MatrixBase.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org
articular the code in the "Possibilities for
> extensions" section) expands on this: If you have a preconditioner for
> the Schur complement, then you can just use that as part of the
> preconditioner for the overall system. It's the best approach we have
> for the Stokes eq
discourage
in-source-builds (for the tutorials that might be ok but not for more
complex code).
https://github.com/konsim83/Deal.ii-9.1.1_Source_code_organization
(just a suggestion that I use for my stuff, no setup of unit tests
implemented so far but that is easy to add)
Best,
Konrad
Hi dear deal.ii community!
I am trying to solve 2D transient problem in semiconductor devices and I'm
stuck with application of Dirichlet boundary values via ConstrainMatrix or
AffineMatrix (in the last version of this library)
The program (which I am rewriting) solves consecutively two
it is proper way to deal with this problem.
W dniu środa, 4 września 2019 18:17:23 UTC+2 użytkownik Wolfgang Bangerth
> napisał:
>>
>>
>> Konrad,
>>
>> > I am trying to solve 2D transient problem in semiconductor devices and
>> > I'm stuck with
Dear deal.ii community,
I have a quick question: Suppose I know a vector that is close to the
solution of my linear system. Is there a method to initialize an iterative
solver with this vector?
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum
solver_cg.h. I really hoped that would speed up things a bit when solving
time dependent problems.. Hm.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this mess
equation so it is not excluded that I
will ask for help in the near future on this forum (although I hope it
won't be necessary :).
best regards,
Konrad
W dniu czwartek, 5 września 2019 18:40:20 UTC+2 użytkownik Wolfgang
Bangerth napisał:
>
>
> Konrad,
>
> > I tried to gi
Hi Konrad,
Thank you very much for your comment! I will try to go into the detail of
it in next days. I must admit that to this point I've just used
::make_zero_boundary_constraints with regard to the conditions on E field.
I didn't have a better idea, and your solutions seem much more better
of
deal.ii? When looking at it I see that some version dependent boost patches
are applied but I am not familiar with the details.
Anyone knows this? Is there a workaround?
Thanks and best regards,
Konrad
-- Include
/tmp/u290231/spack-stage/dealii-9.1.1-gfk33pgt5rojhujhmuruaahuuzyzq2zm/spack-src/cmake
Hi!
I'm looking for quick advise. In my calculations I solve two equations (in
2D) consecutively: the Poisson equation where I obtain electric field *E*,
and then I solve the continuity equation where I obtain e.g. density of
electrons *n*. I want to get current that is *j* = C*nE* where
I need to do any special configuration steps for PETSC? I followed the
instructions that are documented on the deal.ii pages on "how to configure
Petsc". https://www.dealii.org/current/external-libs/petsc.html
Best,
Konrad
This is the error:
Running using PETSc.
Number of active cel
rong
> MPIEXEC pickedup by deal.II config.
> I wanted to fix that in Spack https://github.com/spack/spack/pull/11142
> but apparently this solution may not be fully functional for Slurm.
>
I guess I have a clue why I get the error. My backend nodes run on a
different architectu
d a solution the question is
jsut how much time it will eat...
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
Hi Denis,
I don't have the build folders any more so I can not post the error log.
But the error (using spack) occurred with both boost versions. I will post
something once I will find a solution.
At any rate, thank you for your help.
Best,
Konrad
--
The deal.II project is located at http
it on. But once I use
slurm to distribute the job across nodes I get "illegal instruction" erros.
Frustrating.
Nevertheless, thanks for your hint.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/de
so using deal.II? If you have
> any suggestions, please feel free to include them.
>
Look at tutorial step-44. I think you will find what you are looking for.
Hope that helps.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
htt
Hi Toni,
Seems like I missed that little note in the documentation. Thank you :-)
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you
little but annoying thing. I
get
make[4]: warning: jobserver unavailable: using -j1. Add '+' to parent make
rule.
when typing make -j8.
Do you have an idea what I could do?
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https
take the curl but divergence will be zero
by construction. Same if u is in H(div) with Raviart-Thomas approximation
(then the curl is zero by construction). How is the computation of the
gradient done internally?
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing
lls) as a key then
and connect to it an object that, for example, contains specific
information about your cell(s).
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You recei
all shape functions but is slow and cumbersome.
I actually only need it for Q1 elements and lowest order Nedelec and
Raviart-Thomas elements.
Anyone has an idea?
Cheers,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.c
Many thanks, Wolfgang.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from
Hi Zhidong,
On Monday, October 21, 2019 at 1:42:30 AM UTC+2, Zhidong Brian Zhang wrote:
>
> Thank you very much for your prompt reply, Konrad!
>
> My confusion is the output of cell->id(), for example,
> 0_3:000
> 0_3:200
> 0_3:003
> 0_3:006
> 0_3:406
> 0_3:606
&g
somehow. The funny thing is that when I step through the
code in debug mode I see that exactly the cast above fails. Funnily, the
cast dynamic_cast *>(
&(dof.get_triangulation())) works.
Now I am asking myself why? Am I missing something here?
Best regards,
Konrad
This
###
So essentially I want to build a library and then link two executables to
it. It compiles correctly but very slow since I can not invoke *make -j8*
Google unfortunately does not hel
Hi all!
I've encountered a problem with very simple two-cells merging in 2D
(deal.ii v.9.0.0).
Lets say that I have two rectangular cells (for two different material) and
I want to merge it together and then refine those cells properly.
I don't have a problem with this when I want to merge
:06:48 UTC+1 użytkownik Bruno Turcksin
napisał:
>
> Konrad,
>
> What does it mean you cannot merge the cells? Do you get the wrong
> results? Does the code segfault? Do you get an error message?
>
> Best,
>
> Bruno
>
--
The deal.II project is located at http://ww
Oh! Sorry for bothered you, problem lies somewhere else - this topic can be
deleted!
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the
with this?
A minimal example is attached. (I used deal.ii 9.2.0 but it should work
with 9.1.1, too.)
Best regards,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because
Add-on: The problem is in the function
RefineInterpolate::refine_and_interpolate_on_distributed_mesh()
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are
Add-on: The problem is in the function
RefineInterpolate::refine_and_interpolate_on_distributed_mesh()
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are
Dear deal.ii community,
I solved it, sorry for bugging you with it but simple mistakes can bug you
for long...
I simply forgot to re-distribute the dofs for the finite element after
refining the mesh. Ooofff :-/
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org
Hi Wolfgang,
On Tuesday, December 24, 2019 at 5:59:31 PM UTC+1, Wolfgang Bangerth wrote:
>
>
> Konrad,
> your email has no question :-) Is your problem that you can't call 'make
> -j8'
> and your question how to make that possible? If so, what happens if you
> try?
Many thanks, Matthias!
Works!
Best,
Konrad
On Wednesday, December 25, 2019 at 2:07:41 PM UTC+1, Matthias Maier wrote:
>
>
> On Fri, Dec 20, 2019, at 13:07 CST, Konrad Simon
and will I
have to manually MPI communicate such points to any processor and test if
they own the relevant cell?
Any help would be much appreciated.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum
Thank you, Bruno. :-)
On Tuesday, May 12, 2020 at 2:50:35 PM UTC+2, Bruno Turcksin wrote:
>
> Konrad,
>
> There is nothing out of the box. However, deal.II uses p4est which can
> use more that one layer of ghost cells. So you should take a look
> there to see how hard it is
Jean-Paul, thanks a lot. Works! :-)
Cheers,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User
v9.1.1
really depend on the development version of adol-c? Other versions <= 2.6.3
are on github.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because
Oh, seems like I missed that. Thanks you!
Konrad
On Sunday, October 11, 2020 at 11:35:34 PM UTC+2 Wolfgang Bangerth wrote:
> On 10/11/20 3:26 PM, Daniel Arndt wrote:
> >
> > have a look at
> >
> https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how
Thank you, Wolfgang and Daniel. Seems like I will go with the command line
then. I was just wondering if people here use Eclipse's PTP which sounded
like a good graphical tool. In my case it frequently crashes or simply gets
stuck.
Best,
Konrad
On Sunday, October 11, 2020 at 11:35:34 PM UTC+2
and/or is super slow.
What tools are you using? Any recommendations or hints?
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed
(but not always) p will be in a ghost layer. )
Any way to get the actual MPI rank of the owner?
Any help would be much appreciated.
Best,
Konrad
On Thursday, May 14, 2020 at 4:10:33 PM UTC+2 Konrad Simon wrote:
> Thank you, Bruno. :-)
>
>
> On Tuesday, May 12, 2020 at 2:50:35 PM UTC+2, Bruno Tu
.ii is trying to setup
constraints on dofs that are not owned by the MPI process.
I read a little bit in the discussions referring to Nedelec and refinement
- seems a difficult issue. Did anyone run into this, too?
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mail
>
> I suspect you are passing a fully distributed vector to that function, but
> it
> needs to read ghost elements of the vector. Have you tried copying the
> vector
> into a locally_relevant vector, and passing that to the function in
> question?
>
Thank you, Wolfgang, that was the issue.
> That looks like an unrelated error. Can you create a small testcase for
> this
> issue here?
>
I will try to come up with an example.
Thank you again and best regards,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options
rhs on
the orthogonal complement of the kernel of the kernel)
5. use the wrapped matrix in the solver.
Hope that helps.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received t
scaling+90 degree rotation and two
independent shifts) and dimension 6 in 3D. Once you project the rhs a
Krylov solver can deal with your singular problem.
Cheers,
Konrad
On Saturday, December 26, 2020 at 11:06:36 AM UTC+1 Konrad Simon wrote:
> Hi,
>
> On Saturday, December 26, 2020 at 6
Hi David,
Many thanks for the hint. After some research I believe I stumbled over this
issue#7970 <https://github.com/dealii/dealii/issues/7970>? Let's see what I
can do.
Best,
Konrad
On Tuesday, December 8, 2020 at 10:56:44 PM UTC+1 Wells, David wrote:
> Hi Konrad,
>
> Th
ould work for any H(div)-conformal shape function. Am I correct?
I am willing to fix that (in deal.ii itself) if necessary.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You
hasing the bug and hopefully come up with a merge request
soon.
Best,
Konrad
On Wednesday, December 9, 2020 at 3:44:58 PM UTC+1 Konrad Simon wrote:
> Hi David,
>
> Many thanks for the hint. After some research I believe I stumbled over this
> issue#7970 <https://github.com/deal
Hi Jean-Paul,
On Thursday, December 10, 2020 at 11:39:08 PM UTC+1 Jean-Paul Pelteret
wrote:
> HI Konrad,
>
> I have no familiarity with the H-div elements, so I could be wrong with
> this suggestion...
>
> The Fe_Nedelec element suffered from a similar issue, where adj
lement it in the current structures without interfering
with them (too much).
Best,
Konrad
On Wednesday, December 16, 2020 at 9:33:57 PM UTC+1 Konrad Simon wrote:
> Dear Jean-Paul, dear Deal.II community,
>
> I partially solved the problem of sign flipping and permuting the degrees
> of freedom and
ent
such a table then. This is a lot of work and more complicated for Nedelec
elements.
Best,
Konrad
On Saturday, December 12, 2020 at 8:07:19 PM UTC+1 Konrad Simon wrote:
> Hi Jean-Paul,
>
> On Thursday, December 10, 2020 at 11:39:08 PM UTC+1 Jean-Paul Pelteret
> wrote:
>
>
). This is the code snippet that I suspect to be
buggy. I hope the names are self explanatory (essentially it is a
modification of step-32). Note that dim=3 here.
Any help would be appreciated.
Best,
Konrad
/
// System and dof setup
what that means. Is there anything I can do? A workaround?
Thanks in advance and best regards,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to t
have
> to be the integral mean of the pressure that is zero. It could just be the
> arithmetic mean, which you can compute without access to ghost elements.
>
Many thanks, Wolfgang, that works! :-)
Konrad
--
The deal.II project is located at http://www.dealii.org/
For ma
dofs or vertex dofs). Now, lines can have
different orientations but lines also permute within a cell upon face
orientation/rotation changes.
Can anyone point me to the place in Deal.II where this is being taken care
of for FE_Q?
Best,
Konrad
--
The deal.II project is located at http
Dear Jean-Paul,
Many thanks for your reply.
On Tuesday, December 29, 2020 at 9:31:49 PM UTC+1 Jean-Paul Pelteret wrote:
> Hi Konrad,
>
> I'm sorry for taking some time to reply. To be honest, the inner working
> of the FE classes is not something that I've ever had the time or
&
ove to have fixed this long ago, and
> would have loved to help you more with it as well, but time availability
> is
> not always on my side :-((
>
I actually need all these vector elements so we both have an interest to
fix it.
Looking at the current implementation it seems like
ul).
Cheers,
Konrad
On Saturday, June 19, 2021 at 10:23:37 PM UTC+2 aael...@ncsu.edu wrote:
> Hi
> I am Abdo, Ph.D. student at NC state university and was looking for a FE
> package library that is flexible and can be used in modeling and analyzing
> 2D and 3D elastic wave equ
1 - 100 of 118 matches
Mail list logo