docstring update for materials
don't use assemble_system() for subdomains with Neumann boundaries