PeriodicBC crashing on coordinates I don't have

Please:

*** include a *minimal* running example to reproduce the error.

Johan

On 06/29/2012 11:35 PM, Stav Gold wrote:
> New question #201835 on DOLFIN:
> https://answers.launchpad.net/dolfin/+question/201835
>
> At coordinate: x = 0 10416.7
> Traceback (most recent call last):
> File "tanhsbd.py", line 82, in <module>
> pbc1 = PeriodicBC(W.sub(0), hb)
> File "/usr/lib/python2.7/dist-packages/dolfin/cpp.py", line 18486, in __init__
> _cpp.PeriodicBC_swiginit(self,_cpp.new_PeriodicBC(*args))
> RuntimeError:
>
> *** -------------------------------------------------------------------------
> *** DOLFIN encountered an error. If you are not able to resolve this issue
> *** using the information listed below, you can ask for help at
> ***
> *** https://answers.launchpad.net/dolfin
> ***
> *** Remember to include the error message listed below and, if possible,
> *** include a *minimal* running example to reproduce the error.
> ***
> *** -------------------------------------------------------------------------
> *** Error: Unable to apply periodic boundary condition.
> *** Reason: Could not find a pair of matching degrees of freedom.
> *** Where: This error was encountered inside PeriodicBC.cpp.
> *** -------------------------------------------------------------------------
>
>
> Fairly straightforward--but this time, my mesh doesn't actually have a point at (0,10416.7)! It does have one at (0, 104166.666667), though, which you'll notice is bigger by a factor of 10.
>
> Any guesses as to what might be going on?
>

Never mind, it turns out I just hadn't symmetrized the mesh carefully enough.