Re: [deal.II] Re: Saving particles for checkpointing

On 12/9/20 5:32 PM, blais...@gmail.com wrote: I have one tiny question actually. In all tests, the std::ostringstream oss which contains some information about the number of particles and the data being serialized is kept. Generally, I guess this should be written to an auxiliary file, then

[deal.II] Re: Saving particles for checkpointing

I have one tiny question actually. In all tests, the std::ostringstream oss which contains some information about the number of particles and the data being serialized is kept. Generally, I guess this should be written to an auxiliary file, then re-read before deserialization can occur? On

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

Thank you, Mark. I just built dealii from the source (deal.II-9.3.0-pre). And my little test is passing now. Thanks for the help. -Kaushik On Wed, Dec 9, 2020 at 4:36 PM Kaushik Das wrote: > Thank you Mark. > I am using the dealii lib that I got from apt-get from > deal.ii-9.2.0-backports. > I

[deal.II] Re: Saving particles for checkpointing

I think they will! I feel kind-off dumb not having realized that :) Thanks! Bruno On Wednesday, December 9, 2020 at 6:06:17 p.m. UTC-5 peterrum wrote: > Do the following links help? > - > https://github.com/dealii/dealii/blob/master/tests/serialization/particle_handler_01.cc > - >

[deal.II] Re: Saving particles for checkpointing

Do the following links help? - https://github.com/dealii/dealii/blob/master/tests/serialization/particle_handler_01.cc - https://github.com/dealii/dealii/blob/master/tests/serialization/particle_handler_02.cc Peter On Wednesday, 9 December 2020 at 23:15:09 UTC+1 blais...@gmail.com wrote: >

[deal.II] Saving particles for checkpointing

Dear all, I hope you are well. I am trying to add checkpointing in our particle code. Particles do not have a save function, but they have a serialize function, which I think would be exactly what I need. However, the function is not documented. It takes as an argument something of class

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

Thank you Mark. I am using the dealii lib that I got from apt-get from deal.ii-9.2.0-backports. I used PETSc and the abort was on even 1 cpus. I tried 2, 3, 6 cpus and all aborted similarly. I will get the latest master branch and build that. Thanks, Kaushik On Wed, Dec 9, 2020 at 4:23 PM Marc

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

>From your stacktrace I can see you are using PETSc and deal.II 9.2.0 which already incorporates the specified patch. Would you try to build the actual master branch anyways? On Wednesday, December 9, 2020 at 2:11:59 PM UTC-7 Marc Fehling wrote: > Hi Kaushik, > > I am unable to reproduce your

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

Hi Kaushik, I am unable to reproduce your problem with the code you provided on the latest build of deal.II and Trilinos. - On how many processes did you run your program? - Did you use PETSc or Trilinos? - Could you try to build deal.II on the latest master branch? There is a

Re: [deal.II] initial values of different-type data

You have to be a bit careful - unfortunately, C++ has one of the most convoluted rules when objects are zero initialized or not. As a rule of thumb: * whenever you are working with a plain "double", "float", "int", etc. it is best to explicitly initialize the variable (because C++ in most

[deal.II] initial values of different-type data

Hi dealii developers, I would like to ask(confirm)whether zero is the initial value of different-type data, e.g. for double, int, Vector and Tensor. For example, if I just define double x;. Does x have an initial value of 0.0? Moreover, in my task, I set initial values of a variable, e.g. "a",

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

Hi David, Many thanks for the hint. After some research I believe I stumbled over this issue#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, > > This isn't a very

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

Hi Marc and Bruno, I was able to reproduce this abort on an even simpler test. Please see the attached file. Initial grid: /* * --- * | 0 | 0 | * --- * | 1 | 1 | 0 - FEQ, 1 - FE_Nothing * --- */ /* Set refine flags: * --- * | R | R |FEQ * --- *

[deal.II] dealii on Apple M1

Hi All, Has anyone tried to compile dealii for Apple M1 ? Can you share your experience? I've tried to do this for version 9.2 without success. -- 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 ---