2 /// Test div with a 2D Quad non-tensor H(div) basis
3 /// \test Test div with a 2D Quad non-tensor H(div) basis
8 #include "t330-basis.h"
12   CeedVector    u, v;  in main()  local
29 …if (fabs(div[i] - div_in_basis[i]) > 100. * CEED_EPSILON) printf("%f != %f\n", div[i], div_in_basi…  in main()
35   CeedVectorCreate(ceed, num_qpts, &v);  in main()
36   CeedVectorSetValue(v, 0.0);  in main()
38   CeedBasisApply(basis, 1, CEED_NOTRANSPOSE, CEED_EVAL_DIV, u, v);  in main()
43     CeedVectorGetArrayRead(v, CEED_MEM_HOST, &v_array);  in main()
45       if (fabs(p * 0.25 - v_array[i]) > 100. * CEED_EPSILON) printf("%f != %f\n", 2.0, v_array[i]);  in main()
47     CeedVectorRestoreArrayRead(v, &v_array);  in main()
50   CeedVectorSetValue(v, 1.0);  in main()
53   CeedBasisApply(basis, 1, CEED_TRANSPOSE, CEED_EVAL_DIV, v, u);  in main()
60 …if (fabs(num_qpts * 0.25 - u_array[i]) > 100. * CEED_EPSILON) printf("%f != %f\n", 2.0, u_array[i]…  in main()
67   CeedVectorDestroy(&v);  in main()