14411cf47Sjeremylt /// @file 2*52bfb9bbSJeremy L Thompson /// Test Simultaneous Diagonalization 3*52bfb9bbSJeremy L Thompson /// \test Simultaneous Diagonalization 457c64913Sjeremylt #include <ceed.h> 557c64913Sjeremylt 657c64913Sjeremylt int main(int argc, char **argv) { 757c64913Sjeremylt Ceed ceed; 8*52bfb9bbSJeremy L Thompson CeedScalar M[16] = {0.19996678, 0.0745459, -0.07448852, 0.0332866, 9*52bfb9bbSJeremy L Thompson 0.0745459, 1., 0.16666509, -0.07448852, 10*52bfb9bbSJeremy L Thompson -0.07448852, 0.16666509, 1., 0.0745459, 11*52bfb9bbSJeremy L Thompson 0.0332866, -0.07448852, 0.0745459, 0.19996678}; 12*52bfb9bbSJeremy L Thompson CeedScalar K[16] = {3.03344425, -3.41501767, 0.49824435, -0.11667092, 13*52bfb9bbSJeremy L Thompson -3.41501767, 5.83354662, -2.9167733, 0.49824435, 14*52bfb9bbSJeremy L Thompson 0.49824435, -2.9167733, 5.83354662, -3.41501767, 15*52bfb9bbSJeremy L Thompson -0.11667092, 0.49824435, -3.41501767, 3.03344425}; 16*52bfb9bbSJeremy L Thompson CeedScalar x[16], lambda[4]; 1757c64913Sjeremylt 1857c64913Sjeremylt CeedInit(argv[1], &ceed); 19aedaa0e5Sjeremylt 20*52bfb9bbSJeremy L Thompson CeedSimultaneousDiagonalization(ceed, K, M, x, lambda, 4); 21*52bfb9bbSJeremy L Thompson fprintf(stdout, "x:\n"); 22*52bfb9bbSJeremy L Thompson for (int i=0; i<4; i++) { 23*52bfb9bbSJeremy L Thompson for (int j=0; j<4; j++) { 24*52bfb9bbSJeremy L Thompson if (x[j+4*i] <= 1E-14 && x[j+4*i] >= -1E-14) x[j+4*i] = 0; 25*52bfb9bbSJeremy L Thompson fprintf(stdout, "%12.8f\t", x[j+4*i]); 26*52bfb9bbSJeremy L Thompson } 27*52bfb9bbSJeremy L Thompson fprintf(stdout, "\n"); 28*52bfb9bbSJeremy L Thompson } 29*52bfb9bbSJeremy L Thompson fprintf(stdout, "lambda:\n"); 30*52bfb9bbSJeremy L Thompson for (int i=0; i<4; i++) { 31*52bfb9bbSJeremy L Thompson if (lambda[i] <= 1E-14 && lambda[i] >= -1E-14) lambda[i] = 0; 32*52bfb9bbSJeremy L Thompson fprintf(stdout, "%12.8f\n", lambda[i]); 33*52bfb9bbSJeremy L Thompson } 3457c64913Sjeremylt CeedDestroy(&ceed); 3557c64913Sjeremylt return 0; 3657c64913Sjeremylt } 37