On Wed, 2015-07-08 at 13:09 +0200, Jan Blechta wrote:
> The mentioned error seems to be fixed in development version of 
> DOLFIN.

I'm glad to hear I wasn't imagining it!


> Your code passes ok - it just refines elsewhere than necessary.
> Switching the inequality sign in
> 
>         cell_markers[c] = gamma[c.index()] > gamma_0
> 
> seems to behave better.


I'm not quite sure what you mean here.  Do you mean using
gamma[c.index()] < gamma_0 instead?  Is it better in the fixed
development version?  Both inequalities give more or less the same
error on my 1.5.0 installation 
(as does "cell_markers[c] =  not(gamma[c.index()] <= gamma_0)")

> FEniCS 1.6.0 including Debian/Ubuntu binary packages should be 
> released
> soon.

Thanks, I look forward to that.

Drew

_______________________________________________
fenics mailing list
fenics@fenicsproject.org
http://fenicsproject.org/mailman/listinfo/fenics

Reply via email to