Thanks for the comments
On Tue, 27 Sep 2005, Wout Ruijter wrote:
> Three questions regarding periodic boundary conditions:
> 1 - Michael already showed how to implement periodic bc's of type u(l)=u(r),
> how do you add u(l)-u(r)=x? Is that a matter of modifying the righthandside?
> I mean, is the rhs taken into account for constrained dof?
You'll have to do so explicitly. I don't recall whether Michael
finally ended up using a penalty method or building libMesh
constraints to do his boundary conditions, but only the former method
will work for your case. I don't think the DofConstraintRow objects
can include right hand side terms.
> 2 - Is there any way in which of two elements with the same error estimate
> only one will be refined?
If you have a mesh on which 400 elements have error 2 and 600 have
error 1, then AFAIK telling libMesh to refine the worst 30% of
elements will cause a random 300 elements from the first group to be
flagged for refinement. Normally this is barely noticeable: any
elements "left behind" by one refinement step are likely to be the
first elements in line for refinement on the next step. If you need
two meshes to match exactly or one mesh to be perfectly symmetric,
though, this will be a problem.
> If this is not the case I would be able to ensure matching meshes by
> modifying the error estimates, or should I modify the refinement
> flags directly?
Either way is good; I'm not sure which way would be quicker to code.
To be safe make sure to modify the flags conservatively: only set new
refinement flags and remove coarsening flags, never set coarsening or
> 3- Am I right to presume that in the case that no interpolated constraints
> are used the _dof_coupling can be used for x==0?
> I guess this could speed up the whole matter in cases with many couplings.
I don't think I understand this question.