20 constexpr static std::array<std::array<int, 3>, 22> PQ_VALUES = {
28 constexpr static std::array<std::pair<int, int>, 7> N_VALUES = {
32 using Clock    = std::chrono::steady_clock;
33 using Duration = std::chrono::duration<double>;
38   std::random_device               rand_device;  in main()
39   std::default_random_engine       rand_engine(rand_device());  in main()
40   std::uniform_real_distribution<> rand_dist(0.0, 1.0);  in main()
54     std::vector<double> q_ref(dim * Q, 0.0), q_weight(Q, 0.0), interp(P * Q), grad(P * Q * dim);  in main()
55     std::generate(interp.begin(), interp.end(), generate_random);  in main()
56     std::generate(grad.begin(), grad.end(), generate_random);  in main()
79           data_interp_n = std::chrono::duration_cast<Duration>(Clock::now() - start).count();  in main()
91           data_interp_t = std::chrono::duration_cast<Duration>(Clock::now() - start).count();  in main()
112           data_grad_n = std::chrono::duration_cast<Duration>(Clock::now() - start).count();  in main()
124           data_grad_t = std::chrono::duration_cast<Duration>(Clock::now() - start).count();  in main()
136       std::printf("%-*d%-*d%-*d%-*d%-*d%*.*f\n",  in main()
139       std::printf("%-*d%-*d%-*d%-*d%-*d%*.*f\n",  in main()
142       std::printf("%-*d%-*d%-*d%-*d%-*d%*.*f\n",  in main()
145       std::printf("%-*d%-*d%-*d%-*d%-*d%*.*f\n",  in main()