is.submodel function

Identify correctness-preserving submodel relations