Hi all,
I have the following problem when I tried to refine mesh. Here is the code:
mesh1=UnitSquareMesh(2,2)
f=CellFunction("bool",mesh)
f.set_all(False)
f.array()[0]=True
mesh2=refine(mesh,f)
I expected the mesh2 is the one obtained after the only standard refine the first cell in mesh1, which means the first cell is split up into 4 cells, and other does not change. However, it is not. why is this and is there any way to get the mesh just refining the first one?(I guess adapt would be one choose, but it does not work now for me.) Any help would be appreciated. Thanks very much.
