make bempp optional
the mesher takes number of cells not nodes for creation of cuboids
enable domain creation in demag-field mesh method