Add sumpy qbmax support for volume expansion (self eval), add tests