It is obvious from the previous two paragraphs that the equation structure is identical for different type of media. One can either treat the acoustic case as a completely self contained problem or as an elastic problem where the material properties happen to be acoustic. Differences in spatial dimensionality or component dimensionality are reflected just in the bounds of the subscripts. From that we can conclude that if we have a given modeling structure we can use it for all our modeling purposes. It is just a matter of keeping track of the indices and to make sure those indices are efficiently mapped in each case. Thus all the building blocks can reuse the same pieces of code if they are written in a general form (i.e. using the full anisotropic formulation). Clearly for small spatial variations in density we could use a more compact Laplacian in case of the cascaded first order derivative operators.