On Wed, 8 Jul 2015 11:43:50 +0000 Drew Parsons <d.pars...@murdoch.edu.au> wrote:
> 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)") I haven't examined what your code does at all. I just looked onto a plot of the solution and a refinement did not seem "correct". So I changed '>' into '<' and it got much more refined near non-zero second derivative. It's not probably the right thing to do. So never mind. Jan > > > 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
