On Feb 8, 2011, at 5:35 AM, Fadoua El Moustaid wrote:
> > Hi all, > I have a problem with defining a mesh in the following way: > > normally: > nx = 100 > dx = 0.1 > mesh = Grid1D(nx = nx, dx = dx) > > it is a discretization of the the interval [0, 1]. In my case I need a > discretization for [0, 0.1] and for [0.9, 1] to be with space step dx = > 0.01 and for [0.1, 0.9] to be with dx = 0.1. Try this: mesh = Grid1D(dx=[0.01]*10 + [0.1]*8 + [0.01]*10)
