Finite Element Analysis, Python, Sfepy, Motivation.
def get_circle_outer(coors, domain=None): r = np.sqrt(coors[:,0]**2.0 + coors[:,1]**2.0) return np.where(r > 3.9)[0] def get_circle_inner(coors, domain=None): r = np.sqrt(coors[:,0]**2.0 + (coors[:,1]-1.0)**2.0) return np.where(r < 1.1)[0] functions = { 'get_out_cir' : (get_circle_outer,),