>triangulation.active_cell_iterators()){
typename DoFHandler::active_cell_iterator cell (>triangulation,
tria_cell->level(), tria_cell->index(), _handler_recovery);
typename DoFHandler::active_cell_iterator stress_cell (>
triangulation, tria_cell->level(), tria_cell->index(),
there was a small bug in my code I posted above. the following code is
running fine.
Konrad Schneider schrieb am Freitag, 17. März 2023 um 17:03:18 UTC+1:
> Dear Wolfgang,
> yes the mesh is quite coarse, but I want to investigate the functionality
> of how deal.ii outputs stresse
rix, since for that we need the an local_dof_indices object
like std::vector that points only to every second
dof index or so.
I would be grateful for some tips or directions.
Best
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https
is pitfall
that we have to use the vtu-format for tensor data.
But thanks anyway for the quick response (again).
Best regards
Konrad
Wolfgang Bangerth schrieb am Freitag, 10. März 2023 um 20:35:09 UTC+1:
> On 3/10/23 07:41, Konrad Schneider wrote:
> >
> > I am using deal.ii ver
hether you may be able to implement the missing
functionality. If you do, please consider providing a patch to the
deal.II development sources (see the deal.II website on how to
contribute).
I am using deal.ii version 9.4 and am wondering what I did wrong here. Do
you have any suggestions?
Best rega
Dear Wolfgang,
thanks for clearing this up and responding so fast. Your explanation makes
sense to me.
Best
Konrad
Wolfgang Bangerth schrieb am Mittwoch, 8. März 2023 um 15:03:13 UTC+1:
> On 3/8/23 06:34, Konrad Schneider wrote:
> >
> > Only if I uncomment line 30 of my cod
Dear all,
I am fairly new to deal.ii and have a question concerning the
MappingQ-class of deal.ii
Why do I have to initialize the private mapping member of type
MappingQ to get a compilable code? To illustrated what I mean, I did
thin out the step-11 tutorial in the following way:
#include
style as implemented in Deal.II).
Thanks in advance :-)
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 Groups
somehow extract a Triangulation form a Triangulation by
providing, for example, the boundary id on the Triangulation level?
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
stem? Our boundary conditions fix all degrees of freedom in the kernel.
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 Go
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
Thank you, Wolfgang! This is exactly what I am looking for. I don‘t need to
access the pointer from outside, I am trying to extend the p4est-interface
itself.
Best,
Konrad
On Friday, May 14, 2021 at 12:42:48 AM UTC+2 Wolfgang Bangerth wrote:
> On 5/13/21 3:05 PM, Konrad Simon wr
for example
https://www.dealii.org/current/doxygen/deal.II/distributed_2tria_8cc_source.html#l2989
Is this pointer used in other interfaces somehow? p4est itself does not
touch it so I am wondering if I can reset it to anything I‘d like.
Any Ideas?
Best,
Konrad
--
The deal.II project is located
rientation issues will be
> irrelevant to the mesh.
>
Lowest order Nedelec and all other elements I mentioned are fine on the
meshes you need I guess.
Hope that helps.
Best,
Konrad
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https:/
der and 3D all orders), as Jean-Paul pointed out, are another
option. They are meant to by-pass certain mesh orientation issues on
complicated geometries (I can tell you a bit about that since currently I
am chasing some problems there).
Best,
Konrad
On Tuesday, April 20, 2021 at 9:23:23 PM UTC+2
context and will be happy to share my
experiences with you. Just email me: konrad.si...@uni-hamburg.de
Best,
Konrad
On Monday, April 19, 2021 at 5:51:17 AM UTC+2 Wolfgang Bangerth wrote:
> On 4/18/21 12:24 PM, John Smith wrote:
> >
> > Thank you for your reply. It is a bit di
Dear Wolfgang,
On Tuesday, April 6, 2021 at 5:53:15 AM UTC+2 Wolfgang Bangerth wrote:
> It's possible this is indeed a bug. In most cases, we run the 2d meshes
> through the mesh orienter, so we rarely see 2d meshes that are not
> correctly
> oriented and it wouldn't surprise me if we have
(3d RT I could fix).
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 Groups
"deal.II User Group&q
r your system
and can be customized nicely (you can also link against MPI libraries and
other stuff provided by the target system's admin). This is the option that
I use.
Hope that helps.
Best regards,
Konrad
On Friday, February 19, 2021 at 10:07:00 AM UTC+1 venepalli...@gmail.com
wrote:
>
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
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
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
> 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
>
> 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.
.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
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
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
&
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
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
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:
>
>
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
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 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
). 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
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
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
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
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
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
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
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
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
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
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
: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
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
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
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?
###
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
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
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
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
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
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
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
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
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
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
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
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
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 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
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
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
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
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
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
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
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
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:
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
::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
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
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
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
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
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
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
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
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
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
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
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
://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
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
>
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
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
>
> 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
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
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
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
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
> 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
1 - 100 of 118 matches
Mail list logo