On Tue, Oct 13, 2015 at 9:16 AM, Roy Stogner <royst...@ices.utexas.edu>
wrote:

>
> On Mon, 12 Oct 2015, Paul T. Bauman wrote:
>>
>>
>> OK, so we think this is the problem with the
>> over/underrefined_boundary_limit. In that case, we want "0" to
>> really mean no mismatch
>>
>
> Indeed we do.  Again from mesh_refinement.h:
>
>    * If \p overrefined_boundary_limit is negative, then level
>    * differences will be unlimited.
>    *
>    * \p overrefined_boundary_limit is 0 by default.  This implies that
>    * adaptive coarsening can only be done on an interior element if
>    * any boundary elements on its sides are simultaneously coarsened.
>    */
>
> So here "0" is documented to mean "0", a convention whose
> intuitiveness I would congratulate myself about if it wasn't
> immediately preceded by "negative numbers"-mean-"infinity".


So, the derpiness continues. Sorry.


> and, IMHO, we should support that case for
>> over/underrefined_boundary_limit. We can/will cook up a PR, but any
>> comments before we do so?
>>
>
> "IMHO, we should support that case" makes it sound like "Roy decided
> not to support that case but we should change that decision"; the
> reality is "that was the primary case Roy intended to support but we
> need a real programmer to revisit his cut-and-pasted hack job".
>

You're right. Sorry for the erroneous implication. That stemmed from our
apparent inability to adequately read documentation.


> It looks like the fix might just be two lines: change e.g.
>           if (_overrefined_boundary_limit)
> to
>           if (_overrefined_boundary_limit >= 0)
> in mesh_refinement.C, likewise for _underrefined_boundary_limit.
>
> But the hard work is that someone should do what I failed to do in
> the first place: add some more mixed dimension mesh refinement unit
> tests to actually *exercise* this option.
>

Yeah, once we've got a fix, we'll add tests on the libMesh and GRINS side.

Thanks,
------------------------------------------------------------------------------
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel

Reply via email to