disable pariter for now (conflicts with libmatrix)
rectilinear nsubdomains default to prop.nprocs
reenable pariter for Topology.elem_eval