[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] Re: Evaluating FE-solution on distributed mesh, semi-Lagrangian method

2020-05-12 Thread Konrad Simon
for information from previous time steps. Is there a way in Deal.ii to control the size of the halo region. In my case I have pretty uniform grids if it helps. Best, Konrad On Monday, May 11, 2020 at 2:49:21 PM UTC+2, Bruno Turcksin wrote: > > Konrad, > > Unfortunately, you wil

[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

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

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

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

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

Re: [deal.II] projecting boundary values for H(div) [and maybe H(curl)]

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

[deal.II] projecting boundary values for H(div) [and maybe H(curl)]

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

Re: [deal.II] Nedelec elements of second kind and Brezzi-Douglas-Marini spaces

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

Re: [deal.II] Nedelec elements of second kind and Brezzi-Douglas-Marini spaces

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

[deal.II] Nedelec elements of second kind and Brezzi-Douglas-Marini spaces

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

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-11-20 Thread Konrad
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

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-11-16 Thread Konrad
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

[deal.II] Threading vtk output?

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

Re: [deal.II] Inhomogeneous neumann bc for step-20

2018-08-31 Thread Konrad
). 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

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-30 Thread Konrad
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

Re: [deal.II] Inhomogeneous neumann bc for step-20

2018-08-30 Thread Konrad
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

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Konrad
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

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Konrad
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

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Konrad
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

[deal.II] Design question about multiscale FEM with deal.ii

2018-08-28 Thread Konrad
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

  1   2   >