[deal.II] Re: Possibility to output stuffs on the integration points?

2017-08-25 Thread Chuyi Duan
to write structure a VTK file and fill it with some data to > be visualised at the quadrature points. You can modify and extend the code > there to suit your needs. > > Regards, > Jean-Paul > > On Friday, July 14, 2017 at 11:10:26 AM UTC+2, Chuyi Duan wrote: >> >> D

[deal.II] Possibility to output stuffs on the integration points?

2017-07-14 Thread Chuyi Duan
Dear all, I'm using a struct to store informations on each integration point, by generating an output file I can only get a mean value of the points in a cell (.vtk file), it looks so "coarse" but it is actually not. So what I want to do is to output the stored stuffs on each points. I've only

[deal.II] Re: About applying periodic boundary condition on the same mesh by different methods

2017-07-05 Thread Chuyi Duan
Thanks Daniel, i finished it. But i think in the second step ConstraintMatrix::get_constraint_lines() should be (*constraint_entries)[0].first by calling constraint_entries = ConstraintMatrix::get_constraint_entries(* IndexSet::ElementIterator), right? Best, Chuyi 在 2017年6月27日星期二 UTC+2上午12:

[deal.II] Re: About applying periodic boundary condition on the same mesh by different methods

2017-06-26 Thread Chuyi Duan
Hi Daniel, I still have a question about how to apply periodic boundary conditons: Is there any possibility now to apply PBC with a value 1.0 in x_component and y_component while -1.0 in z_component? I've read an older version of step-45 before (like 8.3.0), which showed we could change that wh

[deal.II] Re: About applying periodic boundary condition on the same mesh by different methods

2017-06-13 Thread Chuyi Duan
Hello Daniel, I think i've found where the problem was. It's because I used set_boundary_id with active_cell_iterator after global_refinement. But I found in source codes of collect_periodic_faces function it's just cell_iterator that has been used. I've instead called set_boundary_id at first

[deal.II] About applying periodic boundary condition on the same mesh by different methods

2017-06-01 Thread Chuyi Duan
Dear all, i was trying to use Dealii to solve a 3-dimensional elastic-plastic question on a cubic by applying periodic boundary condition like in step-45 (v. 8.5.0). std::vector:: cell_iterator> > periodicity_vector; GridTools::collect_periodic_faces(dof_handler, 1, 2, 2, periodicity_vector);/