Hi!
I want to generate meshes of checkerboard`s structure.
I tried mshr:
from mshr import *
domain = Rectangle(Point(0., 0.), Point(1., 1.))
n = 10
for k in range(n):
for m in range(n):
domain.set_subdomain(n*k + m + 1, Rectangle(Point(float(k)/n, float(m)/n), Point(float(k+1)/n, float(m+1)/n)))
It works when n is small enough:
But there is some troubles when n>8
Is it bug, or I chose not the right way?
Thanks a lot!