[deal.II] Adaptive TimestepControl

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

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

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

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

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

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

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
> > 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

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