1 /// @file 2 /// Test Simultaneous Diagonalization 3 /// \test Simultaneous Diagonalization 4 #include <ceed.h> 5 6 int main(int argc, char **argv) { 7 Ceed ceed; 8 CeedScalar M[16] = {0.19996678, 0.0745459, -0.07448852, 0.0332866, 9 0.0745459, 1., 0.16666509, -0.07448852, 10 -0.07448852, 0.16666509, 1., 0.0745459, 11 0.0332866, -0.07448852, 0.0745459, 0.19996678}; 12 CeedScalar K[16] = {3.03344425, -3.41501767, 0.49824435, -0.11667092, 13 -3.41501767, 5.83354662, -2.9167733, 0.49824435, 14 0.49824435, -2.9167733, 5.83354662, -3.41501767, 15 -0.11667092, 0.49824435, -3.41501767, 3.03344425}; 16 CeedScalar x[16], lambda[4]; 17 18 CeedInit(argv[1], &ceed); 19 20 CeedSimultaneousDiagonalization(ceed, K, M, x, lambda, 4); 21 fprintf(stdout, "x:\n"); 22 for (int i=0; i<4; i++) { 23 for (int j=0; j<4; j++) { 24 if (x[j+4*i] <= 1E-14 && x[j+4*i] >= -1E-14) x[j+4*i] = 0; 25 fprintf(stdout, "%12.8f\t", x[j+4*i]); 26 } 27 fprintf(stdout, "\n"); 28 } 29 fprintf(stdout, "lambda:\n"); 30 for (int i=0; i<4; i++) { 31 if (lambda[i] <= 1E-14 && lambda[i] >= -1E-14) lambda[i] = 0; 32 fprintf(stdout, "%12.8f\n", lambda[i]); 33 } 34 CeedDestroy(&ceed); 35 return 0; 36 } 37