44     integer fe_dim  local
51     num_elem = prob_size/(degree**fe_dim)
62   r = mod(s, fe_dim)
64   do d = 1, fe_dim
65     sd = s/fe_dim
81     integer fe_dim  local
105     num_nodes = p**fe_dim
106     elem_qpts = num_qpts**fe_dim
110     do i = 1, fe_dim
128       do j = 1, fe_dim
138         do k = 1, fe_dim
160     integer fe_dim  local
172     scalar_size = mesh_size/fe_dim
173     select case (fe_dim)
197     integer fe_dim  local
217     do i = 1, fe_dim
222     coords_size = scalar_size * fe_dim
237       do j = 1, fe_dim
245     call transformmeshcoords(fe_dim, coords_size, coords, exact_volume, err)
258     integer fe_dim, num_comp_x, mesh_degree, sol_degree, num_qpts  variable
283     fe_dim      = 3
310           read(arg_value, '(I10)') fe_dim
311           num_comp_x = fe_dim
355     write (*, *) '  Mesh dimension         [-d] : ', fe_dim
379 …call ceedbasiscreatetensorh1lagrange(ceed, fe_dim, num_comp_x, mesh_degree + 1, num_qpts, ceed_gau…
381 …call ceedbasiscreatetensorh1lagrange(ceed, fe_dim, 1, sol_degree + 1, num_qpts, ceed_gauss, sol_ba…
384     call getcartesianmeshsize(fe_dim, sol_degree, prob_size, num_xyz)
399 …call buildcartesianrestriction(ceed, fe_dim, num_xyz, mesh_degree, num_comp_x, mesh_size, num_qpts…
401 …call buildcartesianrestriction(ceed, fe_dim, num_xyz, sol_degree, 1, sol_size, num_qpts, sol_restr…
406       write (*, *) 'Number of mesh nodes     : ', mesh_size/fe_dim
414     call setcartesianmeshcoords(fe_dim, num_xyz, mesh_degree, mesh_coords, exact_volume, err)
417     build_ctx_data(1) = fe_dim
427       select case (fe_dim)
441       call ceedqfunctionaddinput(qf_build, 'dx', num_comp_x * fe_dim, ceed_eval_grad, err)
455     elem_qpts = num_qpts**fe_dim
456     do i = 1, fe_dim
526       if (fe_dim == 1) then