Re: [deal.II] Averaging derived solution quantities (mises stress)

2023-03-24 Thread Konrad Schneider
>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(),

Re: [deal.II] Averaging derived solution quantities (mises stress)

2023-03-17 Thread Konrad Schneider
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

Re: [deal.II] Averaging derived solution quantities (mises stress)

2023-03-17 Thread Konrad Schneider
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

Re: [deal.II] output derived quanteties (gradients of vector valued solution like strain)

2023-03-10 Thread Konrad Schneider
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

[deal.II] output derived quanteties (gradients of vector valued solution like strain)

2023-03-10 Thread Konrad Schneider
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

Re: [deal.II] Required mapping member initialization of type MappingQ

2023-03-08 Thread Konrad Schneider
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

[deal.II] Required mapping member initialization of type MappingQ

2023-03-08 Thread Konrad Schneider
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

[deal.II] MPI: Evaluating several distributed functions on each node

2022-04-01 Thread Konrad Simon
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

[deal.II] Triangulation from lower-dimensional object

2022-02-17 Thread Konrad Simon
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

[deal.II] Elasticity Preconditioner

2021-09-24 Thread Konrad Simon
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

[deal.II] Re: Question about deal.II flexibility

2021-06-19 Thread Konrad Simon
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

Re: [deal.II] p4est user pointer and distributed::triangulation

2021-05-14 Thread Konrad Simon
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

[deal.II] p4est user pointer and distributed::triangulation

2021-05-13 Thread Konrad Simon
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

Re: [deal.II] interpolate FE_Nelelec

2021-04-23 Thread Konrad Simon
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:/

Re: [deal.II] interpolate FE_Nelelec

2021-04-20 Thread Konrad Simon
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

Re: [deal.II] Re: interpolate FE_Nelelec

2021-04-19 Thread Konrad Simon
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

Re: [deal.II] 2d meshes and orientation

2021-04-06 Thread Konrad Simon
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

[deal.II] 2d meshes and orientation

2021-04-05 Thread Konrad Simon
(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

[deal.II] Re: Need help installing deal.ii with p4est and Trilinos

2021-02-19 Thread Konrad Simon
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: >

Re: [deal.II] Face and line orientation dependent dof permutations in FE_Q

2021-01-30 Thread Konrad Simon
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

[deal.II] Face and line orientation dependent dof permutations in FE_Q

2021-01-13 Thread Konrad Simon
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

Re: [deal.II] Mean value of component of TrilinosWrappers::MPI::BlockVector

2021-01-13 Thread Konrad Simon
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

Re: [deal.II] Re: Mean value of component of TrilinosWrappers::MPI::BlockVector

2021-01-12 Thread Konrad Simon
> 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

Re: [deal.II] Mean value of component of TrilinosWrappers::MPI::BlockVector

2021-01-12 Thread Konrad Simon
> > 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.

[deal.II] Re: Mean value of component of TrilinosWrappers::MPI::BlockVector

2021-01-11 Thread Konrad Simon
.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

[deal.II] Mean value of component of TrilinosWrappers::MPI::BlockVector

2021-01-10 Thread Konrad Simon
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

Re: [deal.II] Div-conforming elements in 3D

2020-12-29 Thread Konrad Simon
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 &

[deal.II] Re: Mean Value Constraints

2020-12-26 Thread Konrad Simon
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

[deal.II] Re: Mean Value Constraints

2020-12-26 Thread Konrad Simon
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

Re: [deal.II] Div-conforming elements in 3D

2020-12-18 Thread Konrad Simon
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

Re: [deal.II] Div-conforming elements in 3D

2020-12-16 Thread Konrad Simon
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: > >

Re: [deal.II] Div-conforming elements in 3D

2020-12-12 Thread Konrad Simon
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

Re: [deal.II] Div-conforming elements in 3D

2020-12-10 Thread Konrad Simon
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

Re: [deal.II] Div-conforming elements in 3D

2020-12-09 Thread Konrad Simon
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

Re: [deal.II] Boundary conditions on components in with FESystem

2020-12-02 Thread Konrad Simon
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

[deal.II] Boundary conditions on components in with FESystem

2020-11-06 Thread Konrad Simon
). 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

Re: [deal.II] Tools for parallel debugging

2020-10-12 Thread Konrad Simon
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

Re: [deal.II] Tools for parallel debugging

2020-10-12 Thread Konrad Simon
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

[deal.II] Tools for parallel debugging

2020-10-11 Thread Konrad Simon
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

Re: [deal.II] Re: Evaluating FE-solution on distributed mesh, semi-Lagrangian method

2020-09-27 Thread Konrad Simon
(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

Re: [deal.II] Re: Evaluating FE-solution on distributed mesh, semi-Lagrangian method

2020-05-14 Thread Konrad Simon
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

[deal.II] Evaluating FE-solution on distributed mesh, semi-Lagrangian method

2020-05-11 Thread 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

Re: [deal.II] Installation with spack - no access to adol-c

2020-03-18 Thread Konrad Simon
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

[deal.II] Installation with spack - no access to adol-c

2020-03-17 Thread Konrad Simon
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

[deal.II] Re: Interpolating to globally refined distributed mesh

2020-02-19 Thread Konrad Simon
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

[deal.II] Re: Interpolating to globally refined distributed mesh

2020-02-18 Thread Konrad Simon
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

[deal.II] Re: Interpolating to globally refined distributed mesh

2020-02-18 Thread Konrad Simon
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

[deal.II] Interpolating to globally refined distributed mesh

2020-02-18 Thread Konrad Simon
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

[deal.II] Re: Simple merge_triangulation problem

2020-01-23 Thread Konrad Wiśniewski
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

[deal.II] Re: Simple merge_triangulation problem

2020-01-23 Thread Konrad Wiśniewski
: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

[deal.II] Simple merge_triangulation problem

2020-01-23 Thread Konrad Wiśniewski
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

Re: [deal.II] Re: cmake with library and executables

2019-12-25 Thread Konrad Simon
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

Re: [deal.II] Re: cmake with library and executables

2019-12-24 Thread 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?

[deal.II] Re: cmake with library and executables

2019-12-20 Thread Konrad Simon
### 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

[deal.II] projecting function onto TrilinosWrappers::MPI::BlockVector

2019-11-24 Thread Konrad Simon
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

[deal.II] Question about postprocessing

2019-11-09 Thread Konrad Simon
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

[deal.II] cmake with library and executables

2019-11-08 Thread Konrad Simon
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

Re: [deal.II] Evaluate shape functions for an element on a given cell in a triangulation

2019-10-27 Thread Konrad Simon
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

[deal.II] Re: How to use CellId class

2019-10-22 Thread Konrad Simon
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

[deal.II] Evaluate shape functions for an element on a given cell in a triangulation

2019-10-20 Thread Konrad Simon
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

[deal.II] Re: How to use CellId class

2019-10-20 Thread Konrad Simon
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

[deal.II] Re: Issue with boost serialization and spack?

2019-10-02 Thread Konrad Simon
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

[deal.II] Re: Issue with boost serialization and spack?

2019-10-02 Thread Konrad Simon
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

[deal.II] Re: Issue with boost serialization and spack?

2019-10-02 Thread Konrad Simon
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

[deal.II] Data postprocessing with two dof handlers

2019-10-01 Thread Konrad Wiśniewski
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

[deal.II] Re: Issue with boost serialization and spack?

2019-09-30 Thread Konrad Simon
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

[deal.II] Issue with boost serialization and spack?

2019-09-28 Thread Konrad Simon
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

[deal.II] Re: Configuration of PETSc

2019-09-20 Thread Konrad Simon
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

[deal.II] Re: deal.II for Inverse problems in non-linear elasticity

2019-09-20 Thread Konrad Simon
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

[deal.II] Configuration of PETSc

2019-09-20 Thread Konrad Simon
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

Re: [deal.II] Initialization of iterative solver

2019-09-18 Thread Konrad Simon
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

[deal.II] Initialization of iterative solver

2019-09-18 Thread Konrad Simon
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

[deal.II] Re: On using project_boundary_values_div_conforming for (RT-DG) in parallel case

2019-09-11 Thread 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

Re: [deal.II] Question about (parallel) linear solver

2019-09-11 Thread Konrad
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

[deal.II] Re: On the usage of Utilities::MPI

2019-09-11 Thread Konrad
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

[deal.II] Question about (parallel) linear solver

2019-09-11 Thread Konrad
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:

Re: [deal.II] Re: Error with Petsc and and multithreading?

2019-09-09 Thread Konrad
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

[deal.II] Re: Error with Petsc and and multithreading?

2019-09-09 Thread Konrad
::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

Re: [deal.II] Poisson equation, implementation of boundary values, mixed finite element, semiconductor devices

2019-09-08 Thread Konrad Wiśniewski
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

Re: [deal.II] MUMPS with PETScWrappers::MPI::BlockSparseMatrix

2019-09-08 Thread Konrad
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

[deal.II] Error with Petsc and and multithreading?

2019-09-07 Thread Konrad
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

Re: [deal.II] Poisson equation, implementation of boundary values, mixed finite element, semiconductor devices

2019-09-07 Thread Konrad
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

Re: [deal.II] Poisson equation, implementation of boundary values, mixed finite element, semiconductor devices

2019-09-06 Thread Konrad Wiśniewski
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

[deal.II] MUMPS with PETScWrappers::MPI::BlockSparseMatrix

2019-09-06 Thread Konrad
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

[deal.II] Function parser for tensor function

2019-09-06 Thread Konrad
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

Re: [deal.II] Poisson equation, implementation of boundary values, mixed finite element, semiconductor devices

2019-09-05 Thread Konrad Wiśniewski
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

[deal.II] Re: Little question about constructor

2019-09-04 Thread Konrad
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

[deal.II] Poisson equation, implementation of boundary values, mixed finite element, semiconductor devices

2019-09-04 Thread Konrad Wiśniewski
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

Re: [deal.II] Distributing objects on cluster nodes according to distributed triangulation (MPI)

2019-09-04 Thread Konrad
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

Re: [deal.II] Distributing objects on cluster nodes according to distributed triangulation (MPI)

2019-09-03 Thread Konrad
://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

Re: [deal.II] Distributing objects on cluster nodes according to distributed triangulation (MPI)

2019-08-23 Thread Konrad
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

[deal.II] Re: getting components of Blockvector

2019-08-22 Thread Konrad
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 >

[deal.II] Distributing objects on cluster nodes according to distributed triangulation (MPI)

2019-08-22 Thread Konrad
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

Re: [deal.II] Some problems with Eclipse IDE

2019-07-22 Thread Konrad
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

Re: [deal.II] Parallel sparse direct solver?

2019-04-29 Thread 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

[deal.II] Parallel sparse direct solver?

2019-04-29 Thread Konrad
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

Re: [deal.II] Pure Neumann BCs in mixed form

2019-03-26 Thread Konrad
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

Re: [deal.II] Pure Neumann BCs in mixed form

2019-03-26 Thread Konrad
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

[deal.II] Pure Neumann BCs in mixed form

2019-03-23 Thread Konrad
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

Re: [deal.II] Re: advection diffusion with periodic boundaries

2019-03-23 Thread Konrad
> 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   2   >