better documentation of load_mesh method, consider moving some part of this method directly into BDDC if something is problem independent.