Chuyi,
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?
>
Yes, you're right ConstraintMatrix::ge
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:
Chuyi,
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
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
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
Chuyi,
> 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_