Efficient discrete implementations for a dynamic problem of linear elasticity

For a dynamic three-dimensional linear elasticity problem in velocities-stresses, we construct efficient difference schemes on the basis of various additive decompositions of the original spatial operator. They include a difference scheme whose efficient implementation at the “predictor” stage has the property of complete conservativeness. Another class of efficient difference schemes is related to the representation of the operator as a product of triangular operators, that is, an operator analog of the LU-decomposition. The parallelism degree of these difference schemes is the same as of explicit schemes.