14411cf47Sjeremylt /// @file 2*2e95cc35Svaleriabarra /// Test that length of BasisApply input/output vectors is incompatible with basis dimensions 352bfb9bbSJeremy L Thompson /// \test Test that topological and geometric dimensions of basis match 457c64913Sjeremylt #include <ceed.h> 557c64913Sjeremylt #include <math.h> 657c64913Sjeremylt 757c64913Sjeremylt int main(int argc, char **argv) { 857c64913Sjeremylt Ceed ceed; 952bfb9bbSJeremy L Thompson CeedBasis b; 1052bfb9bbSJeremy L Thompson CeedVector U, V; 1152bfb9bbSJeremy L Thompson CeedInt Q = 8, P = 2, ncomp = 1, dim = 3, 1252bfb9bbSJeremy L Thompson len = pow((double)(Q), dim); 1357c64913Sjeremylt 1457c64913Sjeremylt CeedInit(argv[1], &ceed); 1557c64913Sjeremylt 1652bfb9bbSJeremy L Thompson CeedVectorCreate(ceed, len, &U); 1752bfb9bbSJeremy L Thompson CeedVectorCreate(ceed, len+1, &V); 18aedaa0e5Sjeremylt 1952bfb9bbSJeremy L Thompson CeedBasisCreateTensorH1Lagrange(ceed, dim, ncomp, P, Q, CEED_GAUSS, &b); 20aedaa0e5Sjeremylt 2152bfb9bbSJeremy L Thompson // Basis apply will error because dimensions don't agree 2252bfb9bbSJeremy L Thompson CeedBasisApply(b, 1, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, U, V); 23aedaa0e5Sjeremylt 24a2546046Sjeremylt // LCOV_EXCL_START 2552bfb9bbSJeremy L Thompson CeedBasisDestroy(&b); 26aedaa0e5Sjeremylt CeedVectorDestroy(&U); 2752bfb9bbSJeremy L Thompson CeedVectorDestroy(&V); 2857c64913Sjeremylt CeedDestroy(&ceed); 2957c64913Sjeremylt return 0; 3052bfb9bbSJeremy L Thompson // LCOV_EXCL_STOP 3157c64913Sjeremylt } 32