[deal.II] Re: Periodic boundary conditions seem to be not applied

Bastian, can you provide us with a minimal example that shows your problem? What you are trying to do should be possible even if it might not be the right thing to do in your situation. Best, Daniel Am Mittwoch, 1. Juni 2016 02:10:36 UTC+2 schrieb Bastien Lauras: > > Hi, > > I am working on a 2

[deal.II] Re: Doubt in FEvalues

Hello Jean, Thanks for your help. It would have been difficult for me to find out without your support. Also, just to clarify, if the bug is indeed with Dof_handler then, why making a new fe_values object (after mesh movement) in my code gave correct gradients ? On Wednesday, June 1, 2016 a

[deal.II] avoid file system lock when calling P::D::SolutionTransfer::prepare_serialization() ?

Hi all, I am running on cluster with LUSTRE file system but it hasn't the 'flock' option turned on. So when I try to save a checkpoint for a future restart I get ``` File locking failed in ADIOI_Set_lock(fd 13,cmd F_SETLKW/7,type F_WRLCK/1,whence 0) with return value and errno 26. - I

[deal.II] Re: Doubt in FEvalues

You're welcome. As to why creating a new FEValues object solves the problem - perhaps it has something to do with the way that FEValues caches data. When you have just a single FEValues, it references "state" information, but when you create a new one the equivalent data is completely reevalua

[deal.II] Re: Doubt in FEvalues

Hi Jean, I am just thinking about it because there may be possibility of bug there as well if reinit(cell) and new fe_values object have something in common like catching data or anything. Even I am not sure about the internal working of the library, so I am just making a guess. Also, in the b

[deal.II] Re: Doubt in FEvalues

I've directed the developers towards this thread, but I'll make a comment about it. On Wednesday, June 1, 2016 at 4:17:22 PM UTC+2, RAJAT ARORA wrote: > > Hi Jean, > > I am just thinking about it because there may be possibility of bug there > as well if reinit(cell) and new fe_values object hav

[deal.II] The order in FEValues.get_quadrature_points and DoFTools::extract_boundary_dofs, are they the same with the dof indices?

Hi, all: I want to get coordinate information of particular boundary dofs. As stated in the title, I will use "FEValues.get_quadrature_points" and "DoFTools::extract_boundary_dofs" so that I get two vectors of points and dof respectively. But I don't know how the components of the vectors is or

[deal.II] Re: Periodic boundary conditions seem to be not applied

Hi Jean-Paul, Thank you for your answer. I am using deal.II 8.4.0. I've seen this add_periodicity() function, but it appears to me, in step-45

[deal.II] Re: Periodic boundary conditions seem to be not applied

Hi Daniel, So, I'm modelling a unit square, with an incompressible neo-hookean material. The formulation is (for the density of energy) : c_1 * (I_1 - 3) - p (I_3 - 1) where c_1 = mu / 2, p is a lagrangian parameter, and I_1 and I_3 are the first and third invariant of C = F^T * F. For me to se

[deal.II] Re: Periodic boundary conditions seem to be not applied

Ah, sorry, I misread what you said. So no, you don't need to use the add_periodicity() function that I suggested. But I would highly recommend that you take a look at Miehe's paper (or some equivalent, maybe from Prof Wrigger's, Prof Steinmann's or Prof Schroder's groups for example) before sp

[deal.II] Re: Periodic boundary conditions seem to be not applied

It is quite difficult to say what is going wrong without seeing the code. You probably want to check first if you really want to impose that kind of boundary conditions and follow Jean-Paul's advice. If the problem persists, it would be really helpful if you can give us a minimal compiling code

Re: [deal.II] avoid file system lock when calling P::D::SolutionTransfer::prepare_serialization() ?

Alberto, I am running on cluster with LUSTRE file system but it hasn't the 'flock' option turned on. So when I try to save a checkpoint for a future restart I get ``` File locking failed in ADIOI_Set_lock(fd 13,cmd F_SETLKW/7,type F_WRLCK/1,whence 0) with return value and errno 26. -

[deal.II] Read and write PETScWrappers::MPI::Vector to binary file

Hi, Wondering if there is any existing implementation for reading and writing PETScWrappers::MPI::Vector to files, essentially to get a restart capability in time stepping computations. PETScWrappers::VectorBase::write_ascii() provides the equivalent of VecView. But looking for something whic

[deal.II] Re: Read and write PETScWrappers::MPI::Vector to binary file

Hi Shiva, I do believe that the serialize() function is the one that you're looking for. Its a function thats inherited from the Subscriptor base class. J-P On Thursday, June 2, 2016 at 12: