Hi, I recently learned about the class Algorithms::TimestepControl and how it is being used in the Algorithms::ThetaTimestepping class. It has a fairly limited set of strategies (uniform step size and doubling step size). I wonder if we can customize it to allow adaptive time stepping. Let me

### Re: [deal.II] Difference between BlockDynamicSparsityPattern and TrilinosWrappers::BlockSparsityPattern

Dear Wolfgang, Thank you, everything is clear now and i managed to accomplish what I wanted. Thanks! Bruno On Tuesday, 3 September 2019 19:09:00 UTC-4, Wolfgang Bangerth wrote: > > > Bruno, > > > Is there a different between how DynamicSparsityPatterns and > > BlockDynamicSparsityPatterns

### Re: [deal.II] Re: Little question about constructor

Dear Konrad, Appreciating for your detailed explanation, it really helps!! Thank you! I am learning c++ by the means of studying deal.ii, so actually I am a layman both in deal.ii and c++. I will read the reference page and books you recommended. Thanks again!! Best regards, Yuesu On Wed, Sep

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

Konrad, > 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 equation:

### [deal.II] Member description of MeshWorker::DoFInfo.matrix

Hello, I am reading step-12 of the tutorial. The following lines are from local integrator for interior face (dinfo is an alias to MeshWorker::DoFInfo). FullMatrix _v1_matrix = dinfo1.matrix(0, false).matrix; FullMatrix

### [deal.II] Re: Little question about constructor

Hi Yuesu, I think this is rather a general C++ question. Once you create an object all members are initialized by invoking a constructor. If you do not declare any constructor the compiler will silently generate a set of constructors. The default constructor that takes no arguments is one of

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

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)

> > boost::mpi is a nice idea in principle, but it's not particularly well > designed and few people use it. Some reasons are listed here: > > > https://scicomp.stackexchange.com/questions/3019/boostmpi-or-c-mpi-for-high-performance-scientific-applications > > > Best > W. > I think it

### Re: [deal.II] How to set material id with MPI

Hi Pham! >From your description I do not really get why you are specifically doing this, so maybe consider the following: I assume, you are flagging cells material ids on one locally owned part due to some custom condition - lets say some stress or function, you cannot formulate in the global