#include <math.h>#include <string.h>#include <gandalf/linalg/3x3matrix.h>#include <gandalf/common/misc_error.h>#include <gandalf/linalg/mat_symmetric.h>#include <gandalf/linalg/mat_diagonal.h>#include <gandalf/linalg/mat_scaledI.h>#include <gandalf/linalg/mat_triangular.h>#include <gandalf/common/endian_io.h>