Bulk waves

SlownessBuddy computes slowness curves for bulk waves with a high degree of accuracy and for any angle (provided the material constants define tensors). The algorithm used is a reformulation of the Christoffel equation as an eigenvalue problem (click here for more... but in French). Due to the properties of the linear problem, eigenvalues are real and a fast and accurate algorithm can be used.

Additionally, the eigenvalues are sorted so as to identify to which slowness curve they actually belong (there are 3 piezoelectrically stiffened acoustic slowness curves). The algorithm used tries to ensure the continuity of the curves; it works well in most practical cases.

The 3 acoustic slowness curves are saved under the names volume1.dat, volume2.dat and volume3.dat.