On Wed, 2015-07-08 at 16:08 +0200, Jan Blechta wrote: > On Wed, 8 Jul 2015 11:43:50 +0000 > Drew Parsons <[email protected]> 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. >
Ok :) Thanks Jan. Drew _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
