  We extend the construction of so-called encapsulated global
summation-by-parts operators to the general case of a mesh which is not
boundary conforming. Owing to this development, energy stable discretizations
of nonlinear and variable coefficient initial boundary value problems can be
formulated in simple and straightforward ways using high-order accurate
operators of generalized summation-by-parts type. Encapsulated features on a
single computational block or element may include polynomial bases, tensor
products as well as curvilinear coordinate transformations. Moreover, through
the use of inner product preserving interpolation or projection, the global
summation-by-parts property in extended to arbitrary multi-block or
multi-element meshes with non-conforming nodal interfaces.
